Allele-specific analysis of DNA replication origins in mammalian cells

Joint optimization of replicative rate and information storage set the letter size of primordial genetic alphabet

The simplest possible informational heteropolymer requires only a two-letter alphabet to be able to store information. The evolutionary choice of four monomers in the informational biomolecules RNA/DNA or their progenitors is intriguing, given the inherent difficulties in the simultaneous and localized prebiotic synthesis of all four monomers of progenitors of RNA/DNA from common precursors on early Earth. Excluding the scenario where a two-letter alphabet genome eventually expanded to include two more letters to code for more amino acids on teleological grounds, we show here that a replicatively superior heteropolymer sequence in an RNA-world-like scenario would have to be composed of at least four letters in order to predictably fold into a specific secondary structure, and hence must have out-competed the two-letter alphabet genomes. As a consequence of our earlier demonstration of the replicative rate advantage of maximal-nucleotide-skew sequences, in this follow-up article, we show that the competing constraints of maximum replicative rate and predictable secondary structure formation can be simultaneously satisfied only by maximally-skewed palindromic heteropolymer sequences composed of a minimum of four letters.

Statistical framework for calling allelic imbalance in high-throughput sequencing data

High-throughput sequencing facilitates large-scale studies of gene regulation and allows tracing the associations of individual genomic variants with changes in gene regulation and expression. Compared to classic association studies, the assessment of an allelic imbalance at heterozygous variants captures functional variant effects with smaller sample sizes, higher sensitivity, and better resolution. Yet, identification of allele-specific variants from allelic read counts remains challenging due to data-dependent biases and overdispersion arising from technical and biological variability. We present MIXALIME, a novel computational framework for calling allele-specific variants in diverse omics data with a repertoire of statistical models accounting for read mapping bias and copy number variation. We benchmark MIXALIME with DNase-Seq, ATAC-Seq, and CAGE-Seq data, and we demonstrate that the allelic imbalance highlights causal variants in GWAS results. Finally, as a showcase of the large-scale practical application of MIXALIME, we present an atlas of variants exhibiting allele-specific chromatin accessibility, built from thousands of available datasets obtained from diverse cell types.

High Nucleotide Skew Palindromic DNA Sequences Function as Potential Replication Origins due to their Unzipping Propensity

Locations of DNA replication initiation in prokaryotes, called "origins of replication", are well-characterized. However, a mechanistic understanding of the sequence dependence of the local unzipping of double-stranded DNA, the first step towards replication initiation, is lacking. Here, utilizing a Markov chain model that was created to address the directional nature of DNA unzipping and replication, we model the sequence dependence of local melting of double-stranded linear DNA segments. We show that generalized palindromic sequences with high nucleotide skews have a low kinetic barrier for local melting near melting temperatures. This allows for such sequences to function as potential replication origins. We support our claim with evidence for high-skew palindromic sequences within the replication origins of mitochondrial DNA, bacteria, archaea and plasmids.

Parent-of-origin-specific DNA replication timing is confined to large imprinted regions

Genomic imprinting involves differential DNA methylation and gene expression between homologous paternal and maternal loci. It remains unclear, however, whether DNA replication also shows parent-of-origin-specific patterns at imprinted or other genomic regions. Here, we investigate genome-wide asynchronous DNA replication utilizing uniparental human embryonic stem cells containing either maternal-only (parthenogenetic) or paternal-only (androgenetic) DNA. Four clusters of imprinted genes exhibited differential replication timing based on parent of origin, while the remainder of the genome, 99.82%, showed no significant replication asynchrony between parental origins. Active alleles in imprinted gene clusters replicated earlier than their inactive counterparts. At the Prader-Willi syndrome locus, replication asynchrony spanned virtually the entirety of S phase. Replication asynchrony was carried through differentiation to neuronal precursor cells in a manner consistent with gene expression. This study establishes asynchronous DNA replication as a hallmark of large imprinted gene clusters.

Replication initiation sites and zones in the mammalian genome: Where are they located and how are they defined?

Eukaryotic DNA replication is a tightly controlled process that occurs in two main steps, i.e., licensing and firing, which take place in the G1 and S phases of the cell cycle, respectively. In Saccharomyces cerevisiae, the budding yeast, replication origins contain consensus sequences that are recognized and bound by the licensing factor Orc1-6, which then recruits the replicative Mcm2-7 helicase. By contrast, mammalian initiation sites lack such consensus sequences, and the mammalian ORC does not exhibit sequence specificity. Studies performed over the past decades have identified replication initiation sites in the mammalian genome using sequencing-based assays, raising the question of whether replication initiation occurs at confined sites or in broad zones across the genome. Although recent reports have shown that the licensed MCMs in mammalian cells are broadly distributed, suggesting that ORC-dependent licensing may not determine the initiation sites/zones, they are predominantly located upstream of actively transcribed genes. This review compares the mechanism of replication initiation in yeast and mammalian cells, summarizes the sequencing-based technologies used for the identification of initiation sites/zones, and proposes a possible mechanism of initiation-site/zone selection in mammalian cells. Future directions and challenges in this field are also discussed.

Dynamic alternative DNA structures in biology and disease

Repetitive elements in the human genome, once considered 'junk DNA', are now known to adopt more than a dozen alternative (that is, non-B) DNA structures, such as self-annealed hairpins, left-handed Z-DNA, three-stranded triplexes (H-DNA) or four-stranded guanine quadruplex structures (G4 DNA). These dynamic conformations can act as functional genomic elements involved in DNA replication and transcription, chromatin organization and genome stability. In addition, recent studies have revealed a role for these alternative structures in triggering error-generating DNA repair processes, thereby actively enabling genome plasticity. As a driving force for genetic variation, non-B DNA structures thus contribute to both disease aetiology and evolution.

Replication initiation: Implications in genome integrity

In every cell cycle, billions of nucleotides need to be duplicated within hours, with extraordinary precision and accuracy. The molecular mechanism by which cells regulate the replication event is very complicated, and the entire process begins way before the onset of S phase. During the G1 phase of the cell cycle, cells prepare by assembling essential replication factors to establish the pre-replicative complex at origins, sites that dictate where replication would initiate during S phase. During S phase, the replication process is tightly coupled with the DNA repair system to ensure the fidelity of replication. Defects in replication and any error must be recognized by DNA damage response and checkpoint signaling pathways in order to halt the cell cycle before cells are allowed to divide. The coordination of these processes throughout the cell cycle is therefore critical to achieve genomic integrity and prevent diseases. In this review, we focus on the current understanding of how the replication initiation events are regulated to achieve genome stability.

Chromosomal coordination and differential structure of asynchronous replicating regions

Stochastic asynchronous replication timing (AS-RT) is a phenomenon in which the time of replication of each allele is different, and the identity of the early allele varies between cells. By taking advantage of stable clonal pre-B cell populations derived from C57BL6/Castaneous mice, we have mapped the genome-wide AS-RT loci, independently of genetic differences. These regions are characterized by differential chromatin accessibility, mono-allelic expression and include new gene families involved in specifying cell identity. By combining population level mapping with single cell FISH, our data reveal the existence of a novel regulatory program that coordinates a fixed relationship between AS-RT regions on any given chromosome, with some loci set to replicate in a parallel and others set in the anti-parallel orientation. Our results show that AS-RT is a highly regulated epigenetic mark established during early embryogenesis that may be used for facilitating the programming of mono-allelic choice throughout development.

Most regions of the mammalian genome replicate both alleles in a synchronous manner, but some loci have been found to replicate asynchronously and the time of replication of each allele is different. Here the authors, by employing clonal mouse cells from a hybrid strain chart replication timing over the entire genome, using polymorphisms to distinguish between the paternal and maternal alleles.

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.

Evolutionary advantage of anti-parallel strand orientation of duplex DNA

DNA in all living systems shares common properties that are remarkably well suited to its function, suggesting refinement by evolution. However, DNA also shares some counter-intuitive properties which confer no obvious benefit, such as strand directionality and anti-parallel strand orientation, which together result in the complicated lagging strand replication. The evolutionary dynamics that led to these properties of DNA remain unknown but their universality suggests that they confer as yet unknown selective advantage to DNA. In this article, we identify an evolutionary advantage of anti-parallel strand orientation of duplex DNA, within a given set of plausible premises. The advantage stems from the increased rate of replication, achieved by dividing the DNA into predictable, independently and simultaneously replicating segments, as opposed to sequentially replicating the entire DNA, thereby parallelizing the replication process. We show that anti-parallel strand orientation is essential for such a replicative organization of DNA, given our premises, the most important of which is the assumption of the presence of sequence-dependent asymmetric cooperativity in DNA.

Reciprocal monoallelic expression of ASAR lncRNA genes controls replication timing of human chromosome 6

DNA replication occurs on mammalian chromosomes in a cell-type distinctive temporal order known as the replication timing program. We previously found that disruption of the noncanonical lncRNA genes ASAR6 and ASAR15 results in delayed replication timing and delayed mitotic chromosome condensation of human chromosomes 6 and 15, respectively. ASAR6 and ASAR15 display random monoallelic expression and display asynchronous replication between alleles that is coordinated with other random monoallelic genes on their respective chromosomes. Disruption of the expressed allele, but not the silent allele, of ASAR6 leads to delayed replication, activation of the previously silent alleles of linked monoallelic genes, and structural instability of human chromosome 6. In this report, we describe a second lncRNA gene (ASAR6-141) on human chromosome 6 that when disrupted results in delayed replication timing in cisASAR6-141 is subject to random monoallelic expression and asynchronous replication and is expressed from the opposite chromosome 6 homolog as ASAR6 ASAR6-141 RNA, like ASAR6 and ASAR15 RNAs, contains a high L1 content and remains associated with the chromosome territory where it is transcribed. Three classes of cis-acting elements control proper chromosome function in mammals: origins of replication, centromeres, and telomeres, which are responsible for replication, segregation, and stability of all chromosomes. Our work supports a fourth type of essential chromosomal element, the "Inactivation/Stability Center," which expresses ASAR lncRNAs responsible for proper replication timing, monoallelic expression, and structural stability of each chromosome.

Genomic methods for measuring DNA replication dynamics

Genomic DNA replicates according to a defined temporal program in which early-replicating loci are associated with open chromatin, higher gene density, and increased gene expression levels, while late-replicating loci tend to be heterochromatic and show higher rates of genomic instability. The ability to measure DNA replication dynamics at genome scale has proven crucial for understanding the mechanisms and cellular consequences of DNA replication timing. Several methods, such as quantification of nucleotide analog incorporation and DNA copy number analyses, can accurately reconstruct the genomic replication timing profiles of various species and cell types. More recent developments have expanded the DNA replication genomic toolkit to assays that directly measure the activity of replication origins, while single-cell replication timing assays are beginning to reveal a new level of replication timing regulation. The combination of these methods, applied on a genomic scale and in multiple biological systems, promises to resolve many open questions and lead to a holistic understanding of how eukaryotic cells replicate their genomes accurately and efficiently.

H2A.Z facilitates licensing and activation of early replication origins

DNA replication is a tightly regulated process that ensures the precise duplication of the genome during the cell cycle¹. In eukaryotes, the licensing and activation of replication origins are regulated by both DNA sequence and chromatin features². However, the chromatin-based regulatory mechanisms remain largely uncharacterized. Here we show that, in HeLa cells, nucleosomes containing the histone variant H2A.Z are enriched with histone H4 that is dimethylated on its lysine 20 residue (H4K20me2) and with bound origin-recognition complex (ORC). In vitro studies show that H2A.Z-containing nucleosomes bind directly to the histone lysine methyltransferase enzyme SUV420H1, promoting H4K20me2 deposition, which is in turn required for ORC1 binding. Genome-wide studies show that signals from H4K20me2, ORC1 and nascent DNA strands co-localize with H2A.Z, and that depletion of H2A.Z results in decreased H4K20me2, ORC1 and nascent-strand signals throughout the genome. H2A.Z-regulated replication origins have a higher firing efficiency and early replication timing compared with other origins. Our results suggest that the histone variant H2A.Z epigenetically regulates the licensing and activation of early replication origins and maintains replication timing through the SUV420H1-H4K20me2-ORC1 axis.

Control of DNA replication timing in the 3D genome

The 3D organization of mammalian chromatin was described more than 30 years ago by visualizing sites of DNA synthesis at different times during the S phase of the cell cycle. These early cytogenetic studies revealed structurally stable chromosome domains organized into subnuclear compartments. Active-gene-rich domains in the nuclear interior replicate early, whereas more condensed chromatin domains that are largely at the nuclear and nucleolar periphery replicate later. During the past decade, this spatiotemporal DNA replication programme has been mapped along the genome and found to correlate with epigenetic marks, transcriptional activity and features of 3D genome architecture such as chromosome compartments and topologically associated domains. But the causal relationship between these features and DNA replication timing and the regulatory mechanisms involved have remained an enigma. The recent identification of cis-acting elements regulating the replication time and 3D architecture of individual replication domains and of long non-coding RNAs that coordinate whole chromosome replication provide insights into such mechanisms.

Mechanisms of DNA Replication and Repair: Insights from the Study of G-Quadruplexes

G-quadruplexes are four-stranded guanine-rich structures that have been demonstrated to occur across the genome in humans and other organisms. They provide regulatory functions during transcription, translation and immunoglobulin gene rearrangement, but there is also a large amount of evidence that they can present a potent barrier to the DNA replication machinery. This mini-review will summarize recent advances in understanding the many strategies nature has evolved to overcome G-quadruplex-mediated replication blockage, including removal of the structure by helicases or nucleases, or circumventing the deleterious effects on the genome through homologous recombination, alternative end-joining or synthesis re-priming. Paradoxically, G-quadruplexes have also recently been demonstrated to provide a positive role in stimulating the initiation of DNA replication. These recent studies have not only illuminated the many roles and consequences of G-quadruplexes, but have also provided fundamental insights into the general mechanisms of DNA replication and its links with genetic and epigenetic stability.

Involvement of G-quadruplex regions in mammalian replication origin activity

Genome-wide studies of DNA replication origins revealed that origins preferentially associate with an Origin G-rich Repeated Element (OGRE), potentially forming G-quadruplexes (G4). Here, we functionally address their requirements for DNA replication initiation in a series of independent approaches. Deletion of the OGRE/G4 sequence strongly decreased the corresponding origin activity. Conversely, the insertion of an OGRE/G4 element created a new replication origin. This element also promoted replication of episomal EBV vectors lacking the viral origin, but not if the OGRE/G4 sequence was deleted. A potent G4 ligand, PhenDC3, stabilized G4s but did not alter the global origin activity. However, a set of new, G4-associated origins was created, whereas suppressed origins were largely G4-free. In vitro Xenopus laevis replication systems showed that OGRE/G4 sequences are involved in the activation of DNA replication, but not in the pre-replication complex formation. Altogether, these results converge to the functional importance of OGRE/G4 elements in DNA replication initiation.

Metazoan DNA replication origins

DNA replication starts with the opening of DNA at sites called DNA replication origins. From the single sequence-specific DNA replication origin of the small Escherichia coli genome, up to thousands of origins that are necessary to replicate the large human genome, strict sequence specificity has been lost. Nevertheless, genome-wide analyses performed in the recent years, using different mapping methods, demonstrated that there are precise locations along the metazoan genome from which replication initiates. These sites contain relaxed sequence consensus and epigenetic features. There is flexibility in the choice of origins to be used during a given cell cycle, probably imposed by evolution and developmental constraints. Here, we will briefly describe their main features.

Mechanisms of establishment and functional significance of DNA demethylation during erythroid differentiation

Erythroid differentiation is associated with global DNA demethylation, but a complete methylome was lacking in the erythroid lineage. We have generated allele-specific base resolution methylomes of primary basophilic erythroblasts (BasoEs) and compared these with 8 other cell types. We found that DNA demethylation during differentiation from hematopoietic stem/progenitor cells (HSPCs) to BasoEs occurred predominantly in intergenic sequences and in inactive gene bodies causing the formation of partially methylated domains (PMDs) in 74% of the BasoE methylome. Moreover, differentially methylated regions (DMRs) between HSPCs and BasoEs occurred mostly in putative enhancer regions and were most often associated with GATA, EKLF, and AP1 binding motifs. Surprisingly, promoters silent in both HSPCs and BasoEs exhibited much more dramatic chromatin changes during differentiation than activated promoters. Unmethylated silent promoters were often associated with active chromatin states in highly methylated domains (HMDs) but with polycomb-repression in PMDs, indicating that silent promoters are generally regulated differently in HMDs and PMDs. We show that long PMDs replicate late, but that short PMDs replicate early and therefore that the partial methylation of DNA after replication during erythroid expansion occurs throughout S phase of the cell cycle. We propose that baseline maintenance methylation following replication decreases during erythroid differentiation resulting in PMD formation and that the presence of HMDs in the BasoE methylome results from transcription-associated DNA methylation of gene bodies. We detected ∼700 large allele-specific DMRs that were enriched in single-nucleotide polymorphisms, suggesting that primary DNA sequence might be a determinant of DNA methylation levels within PMDs.

The Protective Role of Dormant Origins in Response to Replicative Stress

Genome stability requires tight regulation of DNA replication to ensure that the entire genome of the cell is duplicated once and only once per cell cycle. In mammalian cells, origin activation is controlled in space and time by a cell-specific and robust program called replication timing. About 100,000 potential replication origins form on the chromatin in the gap 1 (G1) phase but only 20⁻30% of them are active during the DNA replication of a given cell in the synthesis (S) phase. When the progress of replication forks is slowed by exogenous or endogenous impediments, the cell must activate some of the inactive or "dormant" origins to complete replication on time. Thus, the many origins that may be activated are probably key to protect the genome against replication stress. This review aims to discuss the role of these dormant origins as safeguards of the human genome during replicative stress.

Replication timing and nuclear structure

DNA replication proceeds along spatially and temporally coordinated patterns within the nucleus, thus protecting the genome during the synthesis of new genetic material. While we have been able to visualize replication patterns on DNA fibers for 50 years, recent developments and discoveries have provided a greater insight into how DNA replication is controlled. In this review, we highlight many of these discoveries. Of great interest are the physiological role of the replication timing program, cis and trans-acting factors that modulate replication timing and the effects of chromatin structure on the replication timing program. We also discuss future directions in the study of replication timing.

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

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

Genome-wide analysis of replication timing by next-generation sequencing with E/L Repli-seq

This protocol is an extension to: Nat. Protoc. 6, 870-895 (2014); doi:10.1038/nprot.2011.328; published online 02 June 2011Cycling cells duplicate their DNA content during S phase, following a defined program called replication timing (RT). Early- and late-replicating regions differ in terms of mutation rates, transcriptional activity, chromatin marks and subnuclear position. Moreover, RT is regulated during development and is altered in diseases. Here, we describe E/L Repli-seq, an extension of our Repli-chip protocol. E/L Repli-seq is a rapid, robust and relatively inexpensive protocol for analyzing RT by next-generation sequencing (NGS), allowing genome-wide assessment of how cellular processes are linked to RT. Briefly, cells are pulse-labeled with BrdU, and early and late S-phase fractions are sorted by flow cytometry. Labeled nascent DNA is immunoprecipitated from both fractions and sequenced. Data processing leads to a single bedGraph file containing the ratio of nascent DNA from early versus late S-phase fractions. The results are comparable to those of Repli-chip, with the additional benefits of genome-wide sequence information and an increased dynamic range. We also provide computational pipelines for downstream analyses, for parsing phased genomes using single-nucleotide polymorphisms (SNPs) to analyze RT allelic asynchrony, and for direct comparison to Repli-chip data. This protocol can be performed in up to 3 d before sequencing, and requires basic cellular and molecular biology skills, as well as a basic understanding of Unix and R.

L1 retrotransposon antisense RNA within ASAR lncRNAs controls chromosome-wide replication timing

Proper chromosome duplication is critical for genome integrity and normal cellular function. Platt et al. show that the lncRNA genes ASAR6 and ASAR15 control chromosome-wide replication timing via the antisense strand of L1 retrotransposons located within ASAR6 and ASAR15 RNAs.

Mammalian cells replicate their chromosomes via a temporal replication program. The ASAR6 and ASAR15 genes were identified as loci that when disrupted result in delayed replication and condensation of entire human chromosomes. ASAR6 and ASAR15 are monoallelically expressed long noncoding RNAs that remain associated with the chromosome from which they are transcribed. The chromosome-wide effects of ASAR6 map to the antisense strand of an L1 retrotransposon within ASAR6 RNA, deletion or inversion of which delayed replication of human chromosome 6. Furthermore, ectopic integration of ASAR6 or ASAR15 transgenes into mouse chromosomes resulted in delayed replication and condensation, an increase in H3K27me3, coating of the mouse chromosome with ASAR RNA, and a loss of mouse Cot-1 RNA expression in cis. Targeting the antisense strand of the L1 within ectopically expressed ASAR6 RNA restored normal replication timing. Our results provide direct evidence that L1 antisense RNA plays a functional role in chromosome-wide replication timing of mammalian chromosomes.

Genomic Analysis of the DNA Replication Timing Program during Mitotic S Phase in Maize (Zea mays) Root Tips

All plants and animals must replicate their DNA, using a regulated process to ensure that their genomes are completely and accurately replicated. DNA replication timing programs have been extensively studied in yeast and animal systems, but much less is known about the replication programs of plants. We report a novel adaptation of the "Repli-seq" assay for use in intact root tips of maize (Zea mays) that includes several different cell lineages and present whole-genome replication timing profiles from cells in early, mid, and late S phase of the mitotic cell cycle. Maize root tips have a complex replication timing program, including regions of distinct early, mid, and late S replication that each constitute between 20 and 24% of the genome, as well as other loci corresponding to ∼32% of the genome that exhibit replication activity in two different time windows. Analyses of genomic, transcriptional, and chromatin features of the euchromatic portion of the maize genome provide evidence for a gradient of early replicating, open chromatin that transitions gradually to less open and less transcriptionally active chromatin replicating in mid S phase. Our genomic level analysis also demonstrated that the centromere core replicates in mid S, before heavily compacted classical heterochromatin, including pericentromeres and knobs, which replicate during late S phase.

Biological action in Read-Write genome evolution

Many of the most important evolutionary variations that generated phenotypic adaptations and originated novel taxa resulted from complex cellular activities affecting genome content and expression. These activities included (i) the symbiogenetic cell merger that produced the mitochondrion-bearing ancestor of all extant eukaryotes, (ii) symbiogenetic cell mergers that produced chloroplast-bearing ancestors of photosynthetic eukaryotes, and (iii) interspecific hybridizations and genome doublings that generated new species and adaptive radiations of higher plants and animals. Adaptive variations also involved horizontal DNA transfers and natural genetic engineering by mobile DNA elements to rewire regulatory networks, such as those essential to viviparous reproduction in mammals. In the most highly evolved multicellular organisms, biological complexity scales with 'non-coding' DNA content rather than with protein-coding capacity in the genome. Coincidentally, 'non-coding' RNAs rich in repetitive mobile DNA sequences function as key regulators of complex adaptive phenotypes, such as stem cell pluripotency. The intersections of cell fusion activities, horizontal DNA transfers and natural genetic engineering of Read-Write genomes provide a rich molecular and biological foundation for understanding how ecological disruptions can stimulate productive, often abrupt, evolutionary transformations.

Replicating Large Genomes: Divide and Conquer

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

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.

The gastrula transition reorganizes replication-origin selection in Caenorhabditis elegans

Although some features underlying replication-origin activation in metazoan cells have been determined, little is known about their regulation during metazoan development. Using the nascent-strand purification method, here we identified replication origins throughout Caenorhabditis elegans embryonic development and found that the origin repertoire is thoroughly reorganized after gastrulation onset. During the pluripotent embryonic stages (pregastrula), potential cruciform structures and open chromatin are determining factors that establish replication origins. The observed enrichment of replication origins in transcription factor-binding sites and their presence in promoters of highly transcribed genes, particularly operons, suggest that transcriptional activity contributes to replication initiation before gastrulation. After the gastrula transition, when embryonic differentiation programs are set, new origins are selected at enhancers, close to CpG-island-like sequences, and at noncoding genes. Our findings suggest that origin selection coordinates replication initiation with transcriptional programs during metazoan development.

Regulation of Replication Origins

In eukaryotes, genome duplication starts concomitantly at many replication initiation sites termed replication origins. The replication initiation program is spatially and temporally coordinated to ensure accurate, efficient DNA synthesis that duplicates the entire genome while maintaining other chromatin-dependent functions. Unlike in prokaryotes, not all potential replication origins in eukaryotes are needed for complete genome duplication during each cell cycle. Instead, eukaryotic cells vary the use of initiation sites so that only a fraction of potential replication origins initiate replication each cell cycle. Flexibility in origin choice allows each eukaryotic cell type to utilize different initiation sites, corresponding to unique nuclear DNA packaging patterns. These patterns coordinate replication with gene expression and chromatin condensation. Budding yeast replication origins share a consensus sequence that marks potential initiation sites. Metazoan origins, on the other hand, lack a consensus sequence. Rather, they are associated with a collection of structural features, chromatin packaging features, histone modifications, transcription, and DNA-DNA/DNA-protein interactions. These features confer cell type-specific replication and expression and play an essential role in maintaining genomic stability.

Genome-wide identification and characterisation of human DNA replication origins by initiation site sequencing (ini-seq)

Next-generation sequencing has enabled the genome-wide identification of human DNA replication origins. However, different approaches to mapping replication origins, namely (i) sequencing isolated small nascent DNA strands (SNS-seq); (ii) sequencing replication bubbles (bubble-seq) and (iii) sequencing Okazaki fragments (OK-seq), show only limited concordance. To address this controversy, we describe here an independent high-resolution origin mapping technique that we call initiation site sequencing (ini-seq). In this approach, newly replicated DNA is directly labelled with digoxigenin-dUTP near the sites of its initiation in a cell-free system. The labelled DNA is then immunoprecipitated and genomic locations are determined by DNA sequencing. Using this technique we identify >25,000 discrete origin sites at sub-kilobase resolution on the human genome, with high concordance between biological replicates. Most activated origins identified by ini-seq are found at transcriptional start sites and contain G-quadruplex (G4) motifs. They tend to cluster in early-replicating domains, providing a correlation between early replication timing and local density of activated origins. Origins identified by ini-seq show highest concordance with sites identified by SNS-seq, followed by OK-seq and bubble-seq. Furthermore, germline origins identified by positive nucleotide distribution skew jumps overlap with origins identified by ini-seq and OK-seq more frequently and more specifically than do sites identified by either SNS-seq or bubble-seq.

Distinct epigenetic features of differentiation-regulated replication origins

BACKGROUND

Eukaryotic genome duplication starts at discrete sequences (replication origins) that coordinate cell cycle progression, ensure genomic stability and modulate gene expression. Origins share some sequence features, but their activity also responds to changes in transcription and cellular differentiation status.

RESULTS

To identify chromatin states and histone modifications that locally mark replication origins, we profiled origin distributions in eight human cell lines representing embryonic and differentiated cell types. Consistent with a role of chromatin structure in determining origin activity, we found that cancer and non-cancer cells of similar lineages exhibited highly similar replication origin distributions. Surprisingly, our study revealed that DNase hypersensitivity, which often correlates with early replication at large-scale chromatin domains, did not emerge as a strong local determinant of origin activity. Instead, we found that two distinct sets of chromatin modifications exhibited strong local associations with two discrete groups of replication origins. The first origin group consisted of about 40,000 regions that actively initiated replication in all cell types and preferentially colocalized with unmethylated CpGs and with the euchromatin markers, H3K4me3 and H3K9Ac. The second group included origins that were consistently active in cells of a single type or lineage and preferentially colocalized with the heterochromatin marker, H3K9me3. Shared origins replicated throughout the S-phase of the cell cycle, whereas cell-type-specific origins preferentially replicated during late S-phase.

CONCLUSIONS

These observations are in line with the hypothesis that differentiation-associated changes in chromatin and gene expression affect the activation of specific replication origins.

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

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

Replication origins: determinants or consequences of nuclear organization?

Chromosome replication, gene expression and chromatin assembly all occur on the same template, necessitating a tight spatial and temporal coordination to maintain genomic stability. The distribution of replication initiation events is responsive to local and global changes in chromatin structure and is affected by transcriptional activity. Concomitantly, replication origin sequences, which determine the locations of replication initiation events, can affect chromatin structure and modulate transcriptional efficiency. The flexibility observed in the replication initiation landscape might help achieve complete and accurate genome duplication while coordinating the DNA replication program with transcription and other nuclear processes in a cell-type specific manner. This review discusses the relationships among replication origin distribution, local and global chromatin structures and concomitant nuclear metabolic processes.

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

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

Quadruplex formation by both G-rich and C-rich DNA strands of the C9orf72 (GGGGCC)8•(GGCCCC)8 repeat: effect of CpG methylation

Unusual DNA/RNA structures of the C9orf72 repeat may participate in repeat expansions or pathogenesis of amyotrophic lateral sclerosis and frontotemporal dementia. Expanded repeats are CpG methylated with unknown consequences. Typically, quadruplex structures form by G-rich but not complementary C-rich strands. Using CD, UV and electrophoresis, we characterized the structures formed by (GGGGCC)8 and (GGCCCC)8 strands with and without 5-methylcytosine (5mCpG) or 5-hydroxymethylcytosine (5hmCpG) methylation. All strands formed heterogenous mixtures of structures, with features of quadruplexes (at pH 7.5, in K(+), Na(+) or Li(+)), but no feature typical of i-motifs. C-rich strands formed quadruplexes, likely stabilized by G•C•G•C-tetrads and C•C•C•C-tetrads. Unlike G•G•G•G-tetrads, some G•C•G•C-tetrad conformations do not require the N7-Guanine position, hence C9orf72 quadruplexes still formed when N7-deazaGuanine replace all Guanines. 5mCpG and 5hmCpG increased and decreased the thermal stability of these structures. hnRNPK, through band-shift analysis, bound C-rich but not G-rich strands, with a binding preference of unmethylated > 5hmCpG > 5mCpG, where methylated DNA-protein complexes were retained in the wells, distinct from unmethylated complexes. Our findings suggest that for C-rich sequences interspersed with G-residues, one must consider quadruplex formation and that methylation of quadruplexes may affect epigenetic processes.