Initiation of DNA replication at CpG islands in mammalian chromosomes

The genetic landscape of origins of replication in P. falciparum

Various origin mapping approaches have enabled genome-wide identification of origins of replication (ORI) in model organisms, but only a few studies have focused on divergent organisms. By employing three complementary approaches we provide a high-resolution map of ORIs in Plasmodium falciparum, the deadliest human malaria parasite. We profiled the distribution of origin of recognition complex (ORC) binding sites by ChIP-seq of two PfORC subunits and mapped active ORIs using NFS and SNS-seq. We show that ORIs lack sequence specificity but are not randomly distributed, and group in clusters. Licensing is biased towards regions of higher GC content and associated with G-quadruplex forming sequences (G4FS). While strong transcription likely enhances firing, active origins are depleted from transcription start sites. Instead, most accumulate in transcriptionally active gene bodies. Single molecule analysis of nanopore reads containing multiple initiation events, which could have only come from individual nuclei, showed a relationship between the replication fork pace and the distance to the nearest origin. While some similarities were drawn with the canonic eukaryote model, the distribution of ORIs in P. falciparum is likely shaped by unique genomic features such as extreme AT-richness-a product of evolutionary pressure imposed by the parasitic lifestyle.

CHEX-seq detects single-cell genomic single-stranded DNA with catalytical potential

Genomic DNA (gDNA) undergoes structural interconversion between single- and double-stranded states during transcription, DNA repair and replication, which is critical for cellular homeostasis. We describe "CHEX-seq" which identifies the single-stranded DNA (ssDNA) in situ in individual cells. CHEX-seq uses 3'-terminal blocked, light-activatable probes to prime the copying of ssDNA into complementary DNA that is sequenced, thereby reporting the genome-wide single-stranded chromatin landscape. CHEX-seq is benchmarked in human K562 cells, and its utilities are demonstrated in cultures of mouse and human brain cells as well as immunostained spatially localized neurons in brain sections. The amount of ssDNA is dynamically regulated in response to perturbation. CHEX-seq also identifies single-stranded regions of mitochondrial DNA in single cells. Surprisingly, CHEX-seq identifies single-stranded loci in mouse and human gDNA that catalyze porphyrin metalation in vitro, suggesting a catalytic activity for genomic ssDNA. We posit that endogenous DNA enzymatic activity is a function of genomic ssDNA.

AAV genome modification for efficient AAV production

The adeno-associated virus (AAV) is one of the most potent vectors in gene therapy. The experimental profile of this vector shows its efficiency and accepted safety, which explains its increased usage by scientists for the research and treatment of a wide range of diseases. These studies require using functional, pure, and high titers of vector particles. In fact, the current knowledge of AAV structure and genome helps improve the scalable production of AAV vectors. In this review, we summarize the latest studies on the optimization of scalable AAV production through modifying the AAV genome or biological processes inside the cell.

Detection and characterization of constitutive replication origins defined by DNA polymerase epsilon

BACKGROUND

Despite the process of DNA replication being mechanistically highly conserved, the location of origins of replication (ORI) may vary from one tissue to the next, or between rounds of replication in eukaryotes, suggesting flexibility in the choice of locations to initiate replication. Lists of human ORI therefore vary widely in number and location, and there are currently no methods available to compare them. Here, we propose a method of detection of ORI based on somatic mutation patterns generated by the mutator phenotype of damaged DNA polymerase epsilon (POLE).

RESULTS

We report the genome-wide localization of constitutive ORI in POLE-mutated human tumors using whole genome sequencing data. Mutations accumulated after many rounds of replication of unsynchronized dividing cell populations in tumors allow to identify constitutive origins, which we show are shared with high fidelity between individuals and tumor types. Using a Smith-Waterman-like dynamic programming approach, we compared replication origin positions obtained from multiple different methods. The comparison allowed us to define a consensus set of replication origins, identified consistently by multiple ORI detection methods. Many DNA features co-localized with the consensus set of ORI, including chromatin loop anchors, G-quadruplexes, S/MARs, and CpGs. Among all features, the H2A.Z histone exhibited the most significant association.

CONCLUSIONS

Our results show that mutation-based detection of replication origins is a viable approach to determining their location and associated sequence features.

MTA2 sensitizes gastric cancer cells to PARP inhibition by induction of DNA replication stress

Poly (ADP-ribose) polymerase (PARP) inhibitor olaparib selectively kills cancer cells with BRCA-deficiency and is approved for BRCA-mutated breast, ovarian and pancreatic cancers by FDA. However, phase III study of olaparib failed to show a significant improvement in overall survival in patients with gastric cancer (GC). To discover an effective biomarker for GC patient-selection in olaparib treatment, we analyzed proteomic profiling of 12 GC cell lines. MTA2 was identified to confer sensitivity to olaparib by aggravating olaparib-induced replication stress in cancer cells. Mechanistically, we applied Cleavage Under Targets and Tagmentation assay to find that MTA2 proteins preferentially bind regions of replication origin-associated DNA sequences, which could be enhanced by olaparib treatment. Furthermore, MTA2 was validated here to render cancer cells susceptible to combination of olaparib with ATR inhibitor AZD6738. In general, our study identified MTA2 as a potential biomarker for olaparib sensitivity by aggravating olaparib-induced replication stress.

Adeno-Associated Virus Serotype-Specific Inverted Terminal Repeat Sequence Role in Vector Transgene Expression

Adeno-associated viral vectors have been successfully used in laboratory and clinical settings for efficient gene delivery. In these vectors, 96% of the adeno-associated virus (AAV) genome is replaced with a gene cassette of interest, leaving only the 145 bp inverted terminal repeat (ITR) sequences. These cis-elements, primarily from AAV serotype 2, are required for genome rescue, replication, packaging, and vector persistence. Previous work from our lab and others have demonstrated that the AAV ITR2 sequence has inherent transcriptional activity, which may confound intended transgene expression in therapeutic applications. Currently, AAV capsids are extensively study for vector contribution; however, a comprehensive analysis of ITR promoter activity of various AAV serotypes has not been described to date. Here, the transcriptional activity of AAV ITRs from different serotypes (1-4, 6, and 7) was compared in numerous cell lines and a mouse model. Under the conditions used here, all ITRs tested were capable of promoting transgene expression both in vitro and in vivo. However, we observed three classes of AAV ITR expression in vitro. Class I ITRs (AAV2 and 3) generated the highest level, whereas class II (AAV 4) had intermediate levels, and class III (AAV1 and 6) had the lowest levels. These expression levels were consistent across multiple cell lines. Only ITR7 demonstrated cell-type dependent transcriptional activity. In vivo, all classes had promoter activity. Next-generation sequencing revealed multiple transcriptional start sites that originated from the ITR sequence, with most arising from within the Rep binding element. The collective results demonstrate that the serotype ITR sequence may have multiple levels of influence on transgene expression cassettes independent of promoter selection.

Binding of the Treslin-MTBP Complex to Specific Regions of the Human Genome Promotes the Initiation of DNA Replication

The processes that control where higher eukaryotic cells initiate DNA replication throughout the genome are not understood clearly. In metazoans, the Treslin-MTBP complex mediates critical final steps in formation of the activated replicative helicase prior to initiation of replication. Here, we map the genome-wide distribution of the MTBP subunit of this complex in human cells. Our results indicate that MTBP binds to at least 30,000 sites in the genome. A majority of these sites reside in regions of open chromatin that contain transcriptional-regulatory elements (e.g., promoters, enhancers, and super-enhancers), which are known to be preferred areas for initiation of replication. Furthermore, many binding sites encompass two genomic features: a nucleosome-free DNA sequence (e.g., G-quadruplex DNA or AP-1 motif) and a nucleosome bearing histone marks characteristic of open chromatin, such as H3K4me2. Taken together, these findings indicate that Treslin-MTBP associates coordinately with multiple genomic signals to promote initiation of replication.

Kumagai and Dunphy show that Treslin-MTBP, activator of the replicative helicase, binds to at least 30,000 sites in the human genome. Many sites contain a nucleosome with active chromatin marks and nucleosome-free DNA (G-quadruplex or AP-1 site). Thus, Treslin-MTBP associates with multiple genomic elements to promote initiation of DNA replication.

A predictable conserved DNA base composition signature defines human core DNA replication origins

DNA replication initiates from multiple genomic locations called replication origins. In metazoa, DNA sequence elements involved in origin specification remain elusive. Here, we examine pluripotent, primary, differentiating, and immortalized human cells, and demonstrate that a class of origins, termed core origins, is shared by different cell types and host ~80% of all DNA replication initiation events in any cell population. We detect a shared G-rich DNA sequence signature that coincides with most core origins in both human and mouse genomes. Transcription and G-rich elements can independently associate with replication origin activity. Computational algorithms show that core origins can be predicted, based solely on DNA sequence patterns but not on consensus motifs. Our results demonstrate that, despite an attributed stochasticity, core origins are chosen from a limited pool of genomic regions. Immortalization through oncogenic gene expression, but not normal cellular differentiation, results in increased stochastic firing from heterochromatin and decreased origin density at TAD borders.

In metazoan the DNA sequence elements characterizing origin specification are unknown. By generating and analysing 19 SNS-seq datasets from different human cell types, the authors reveal a class and features of Core origins of replication which can be predicted by an algorithm.

Delving into Eukaryotic Origins of Replication Using DNA Structural Features

DNA replication in eukaryotes is an intricate process, which is precisely synchronized by a set of regulatory proteins, and the replication fork emanates from discrete sites on chromatin called origins of replication (Oris). These spots are considered as the gateway to chromosomal replication and are stereotyped by sequence motifs. The cognate sequences are noticeable in a small group of entire origin regions or totally absent across different metazoans. Alternatively, the use of DNA secondary structural features can provide additional information compared to the primary sequence. In this article, we report the trends in DNA sequence-based structural properties of origin sequences in nine eukaryotic systems representing different families of life. Biologically relevant DNA secondary structural properties, namely, stability, propeller twist, flexibility, and minor groove shape were studied in the sequences flanking replication start sites. Results indicate that Oris in yeasts show lower stability, more rigidity, and narrow minor groove preferences compared to genomic sequences surrounding them. Yeast Oris also show preference for A-tracts and the promoter element TATA box in the vicinity of replication start sites. On the contrary, Drosophila melanogaster, humans, and Arabidopsis thaliana do not have such features in their Oris, and instead, they show high preponderance of G-rich sequence motifs such as putative G-quadruplexes or i-motifs and CpG islands. Our extensive study applies the DNA structural feature computation to delve into origins of replication across organisms ranging from yeasts to mammals and including a plant. Insights from this study would be significant in understanding origin architecture and help in designing new algorithms for predicting DNA trans-acting factor recognition events.

Origins of DNA replication

In all kingdoms of life, DNA is used to encode hereditary information. Propagation of the genetic material between generations requires timely and accurate duplication of DNA by semiconservative replication prior to cell division to ensure each daughter cell receives the full complement of chromosomes. DNA synthesis of daughter strands starts at discrete sites, termed replication origins, and proceeds in a bidirectional manner until all genomic DNA is replicated. Despite the fundamental nature of these events, organisms have evolved surprisingly divergent strategies that control replication onset. Here, we discuss commonalities and differences in replication origin organization and recognition in the three domains of life.

5-hydroxymethylcytosine Marks Mammalian Origins Acting as a Barrier to Replication

In most mammalian cells, DNA replication occurs once, and only once between cell divisions. Replication initiation is a highly regulated process with redundant mechanisms that prevent errant initiation events. In lower eukaryotes, replication is initiated from a defined consensus sequence, whereas a consensus sequence delineating mammalian origin of replication has not been identified. Here we show that 5-hydroxymethylcytosine (5hmC) is present at mammalian replication origins. Our data support the hypothesis that 5hmC has a role in cell cycle regulation. We show that 5hmC level is inversely proportional to proliferation; indeed, 5hmC negatively influences cell division by increasing the time a cell resides in G1. Our data suggest that 5hmC recruits replication-licensing factors, then is removed prior to or during origin firing. Later we propose that TET2, the enzyme catalyzing 5mC to 5hmC conversion, acts as barrier to rereplication. In a broader context, our results significantly advance the understating of 5hmC involvement in cell proliferation and disease states.

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.

Genome-wide mapping and profiling of γH2AX binding hotspots in response to different replication stress inducers

BACKGROUND

Replication stress (RS) gives rise to DNA damage that threatens genome stability. RS can originate from different sources that stall replication by diverse mechanisms. However, the mechanism underlying how different types of RS contribute to genome instability is unclear, in part due to the poor understanding of the distribution and characteristics of damage sites induced by different RS mechanisms.

RESULTS

We use ChIP-seq to map γH2AX binding sites genome-wide caused by aphidicolin (APH), hydroxyurea (HU), and methyl methanesulfonate (MMS) treatments in human lymphocyte cells. Mapping of γH2AX ChIP-seq reveals that APH, HU, and MMS treatments induce non-random γH2AX chromatin binding at discrete regions, suggesting that there are γH2AX binding hotspots in the genome. Characterization of the distribution and sequence/epigenetic features of γH2AX binding sites reveals that the three treatments induce γH2AX binding at largely non-overlapping regions, suggesting that RS may cause damage at specific genomic loci in a manner dependent on the fork stalling mechanism. Nonetheless, γH2AX binding sites induced by the three treatments share common features including compact chromatin, coinciding with larger-than-average genes, and depletion of CpG islands and transcription start sites. Moreover, we observe significant enrichment of SINEs in γH2AX sites in all treatments, indicating that SINEs may be a common barrier for replication polymerases.

CONCLUSIONS

Our results identify the location and common features of genome instability hotspots induced by different types of RS, and help in deciphering the mechanisms underlying RS-induced genetic diseases and carcinogenesis.

Epigenetic Regulations in Neuropsychiatric Disorders

Precise genetic and epigenetic spatiotemporal regulation of gene expression is critical for proper brain development, function and circuitry formation in the mammalian central nervous system. Neuronal differentiation processes are tightly regulated by epigenetic mechanisms including DNA methylation, histone modifications, chromatin remodelers and non-coding RNAs. Dysregulation of any of these pathways is detrimental to normal neuronal development and functions, which can result in devastating neuropsychiatric disorders, such as depression, schizophrenia and autism spectrum disorders. In this review, we focus on the current understanding of epigenetic regulations in brain development and functions, as well as their implications in neuropsychiatric disorders.

The Promoter Toolbox for Recombinant Gene Expression in Trichoderma reesei

The ascomycete Trichoderma reesei is one of the main fungal producers of cellulases and xylanases based on its high production capacity. Its enzymes are applied in food, feed, and textile industry or in lignocellulose hydrolysis in biofuel and biorefinery industry. Over the last years, the demand to expand the molecular toolbox for T. reesei to facilitate genetic engineering and improve the production of heterologous proteins grew. An important instrument to modify the expression of key genes are promoters to initiate and control their transcription. To date, the most commonly used promoter for T. reesei is the strong inducible promoter of the main cellobiohydrolase cel7a. Beside this one, there is a number of alternative inducible promoters derived from other cellulase- and xylanase encoding genes and a few constitutive promoters. With the advances in genomics and transcriptomics the identification of new constitutive and tunable promoters with different expression strength was simplified. In this review, we will discuss new developments in the field of promoters and compare their advantages and disadvantages. Synthetic expression systems constitute a new option to control gene expression and build up complex gene circuits. Therefore, we will address common structural features of promoters and describe options for promoter engineering and synthetic design of promoters. The availability of well-characterized gene expression control tools is essential for the analysis of gene function, detection of bottlenecks in gene networks and yield increase for biotechnology applications.

DNA Topoisomerase I differentially modulates R-loops across the human genome

BACKGROUND

Co-transcriptional R-loops are abundant non-B DNA structures in mammalian genomes. DNA Topoisomerase I (Top1) is often thought to regulate R-loop formation owing to its ability to resolve both positive and negative supercoils. How Top1 regulates R-loop structures at a global level is unknown.

RESULTS

Here, we perform high-resolution strand-specific R-loop mapping in human cells depleted for Top1 and find that Top1 depletion results in both R-loop gains and losses at thousands of transcribed loci, delineating two distinct gene classes. R-loop gains are characteristic for long, highly transcribed, genes located in gene-poor regions anchored to Lamin B1 domains and in proximity to H3K9me3-marked heterochromatic patches. R-loop losses, by contrast, occur in gene-rich regions overlapping H3K27me3-marked active replication initiation regions. Interestingly, Top1 depletion coincides with a block of the cell cycle in G0/G1 phase and a trend towards replication delay.

CONCLUSIONS

Our findings reveal new properties of Top1 in regulating R-loop homeostasis in a context-dependent manner and suggest a potential role for Top1 in modulating the replication process via R-loop formation.

Methylation-Demethylation Dynamics: Implications of Changes in Acute Kidney Injury

Over the years, the epigenetic landscape has grown increasingly complex. Until recently, methylation of DNA and histones was considered one of the most important epigenetic modifications. However, with the discovery of enzymes involved in the demethylation process, several exciting prospects have emerged that focus on the dynamic regulation of methylation and its crucial role in development and disease. An interplay of the methylation-demethylation machinery controls the process of gene expression. Since acute kidney injury (AKI), a major risk factor for chronic kidney disease and death, is characterised by aberrant expression of genes, understanding the dynamics of methylation and demethylation will provide new insights into the intricacies of the disease. Research on epigenetics in AKI has only made its mark in the recent years but has provided compelling evidence that implicates the involvement of methylation and demethylation changes in its pathophysiology. In this review, we explore the role of methylation and demethylation machinery in cellular epigenetic control and further discuss the contribution of methylomic changes and histone modifications to the pathophysiology of AKI.

Targeting GLI by GANT61 involves mechanisms dependent on inhibition of both transcription and DNA licensing

The GLI genes are transcription factors and in cancers are oncogenes, aberrantly and constitutively activated. GANT61, a specific GLI inhibitor, has induced extensive cytotoxicity in human models of colon cancer. The FOXM1 promoter was determined to be a transcriptional target of GLI1. In HT29 cells, inhibition of GLI1 binding at the GLI consensus sequence by GANT61 led to inhibited binding of Pol II, the pause-release factors DSIF, NELF and p-TEFb. The formation of R-loops (RNA:DNA hybrids, ssDNA), were reduced by GANT61 at the FOXM1 promoter. Pretreatment of HT29 cells with α-amanitin reduced GANT61-induced γH2AX foci. Co-localization of GLI1 and BrdU foci, inhibited by GANT61, indicated GLI1 and DNA replication to be linked. By co-immunoprecipitation and confocal microscopy, GLI1 co-localized with the DNA licensing factors ORC4, CDT1, and MCM2. Significant co-localization of GLI1 and ORC4 was inhibited by GANT61, and enrichment of ORC4 occurred at the GLI binding site in the FOXM1 promoter. CDT1 was found to be a transcription target of GLI1. Overexpression of CDT1 in HT29 and SW480 cells reduced GANT61-induced cell death, gH2AX foci, and cleavage of caspase-3. Data demonstrate involvement of transcription and of DNA replication licensing factors by non-transcriptional and transcriptional mechanisms in the GLI-dependent mechanism of action of GANT61.

Conservation and Variation in Strategies for DNA Replication of Kinetoplastid Nuclear Genomes

Introduction:

Understanding how the nuclear genome of kinetoplastid parasites is replicated received experimental stimulus from sequencing of the Leishmania major, Trypanosoma brucei and Trypanosoma cruzi genomes around 10 years ago. Gene annotations suggested key players in DNA replication initiation could not be found in these organisms, despite considerable conservation amongst characterised eukaryotes. Initial studies that indicated trypanosomatids might possess an archaeal-like Origin Recognition Complex (ORC), composed of only a single factor termed ORC1/CDC6, have been supplanted by the more recent identification of an ORC in T. brucei. However, the constituent subunits of T. brucei ORC are highly diverged relative to other eukaryotic ORCs and the activity of the complex appears subject to novel, positive regulation. The availability of whole genome sequences has also allowed the deployment of genome-wide strategies to map DNA replication dynamics, to date in T. brucei and Leishmania. ORC1/CDC6 binding and function in T. brucei displays pronounced overlap with the unconventional organisation of gene expression in the genome. Moreover, mapping of sites of replication initiation suggests pronounced differences in replication dynamics in Leishmania relative to T. brucei.

Conclusion:

Here we discuss what implications these emerging data may have for parasite and eukaryotic biology of DNA replication.

Methylation-Demethylation Dynamics: Implications of Changes in Acute Kidney Injury

Over the years, the epigenetic landscape has grown increasingly complex. Until recently, methylation of DNA and histones was considered one of the most important epigenetic modifications. However, with the discovery of enzymes involved in the demethylation process, several exciting prospects have emerged that focus on the dynamic regulation of methylation and its crucial role in development and disease. An interplay of the methylation-demethylation machinery controls the process of gene expression. Since acute kidney injury (AKI), a major risk factor for chronic kidney disease and death, is characterised by aberrant expression of genes, understanding the dynamics of methylation and demethylation will provide new insights into the intricacies of the disease. Research on epigenetics in AKI has only made its mark in the recent years but has provided compelling evidence that implicates the involvement of methylation and demethylation changes in its pathophysiology. In this review, we explore the role of methylation and demethylation machinery in cellular epigenetic control and further discuss the contribution of methylomic changes and histone modifications to the pathophysiology of AKI.

CpG Islands: A Historical Perspective

The discovery of CpG islands (CGIs) and the study of their structure and properties run parallel to the development of molecular biology in the last two decades of the twentieth century and to the development of high-throughput genomic technologies at the turn of the millennium. First identified as discrete G + C-rich regions of unmethylated DNA in several vertebrates, CGIs were soon found to display additional distinctive chromatin features from the rest of the genome in terms of accessibility and of the epigenetic modifications of their histones. These features, together with their colocalization with promoters and with origins of DNA replication in mammals, highlighted their relevance in the regulation of genomic processes. Recent approaches have shown with unprecedented detail the dynamics and diversity of the epigenetic landscape of CGIs during normal development and under pathological conditions. Also, comparative analyses across species have started revealing how CGIs evolve and contribute to the evolution of the vertebrate genome.

DNA replication origins-where do we begin?

In this review, Prioleau and MacAlpine summarize recent advances in our understanding of how primary sequence, chromatin environment, and nuclear architecture contribute to the dynamic selection and activation of replication origins across diverse cell types and developmental stages.

For more than three decades, investigators have sought to identify the precise locations where DNA replication initiates in mammalian genomes. The development of molecular and biochemical approaches to identify start sites of DNA replication (origins) based on the presence of defining and characteristic replication intermediates at specific loci led to the identification of only a handful of mammalian replication origins. The limited number of identified origins prevented a comprehensive and exhaustive search for conserved genomic features that were capable of specifying origins of DNA replication. More recently, the adaptation of origin-mapping assays to genome-wide approaches has led to the identification of tens of thousands of replication origins throughout mammalian genomes, providing an unprecedented opportunity to identify both genetic and epigenetic features that define and regulate their distribution and utilization. Here we summarize recent advances in our understanding of how primary sequence, chromatin environment, and nuclear architecture contribute to the dynamic selection and activation of replication origins across diverse cell types and developmental stages.

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.

Mechanisms and regulation of DNA replication initiation in eukaryotes

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

Transcriptionally Driven DNA Replication Program of the Human Parasite Leishmania major

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

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.

The chromatin environment shapes DNA replication origin organization and defines origin classes

To unveil the still-elusive nature of metazoan replication origins, we identified them genome-wide and at unprecedented high-resolution in mouse ES cells. This allowed initiation sites (IS) and initiation zones (IZ) to be differentiated. We then characterized their genetic signatures and organization and integrated these data with 43 chromatin marks and factors. Our results reveal that replication origins can be grouped into three main classes with distinct organization, chromatin environment, and sequence motifs. Class 1 contains relatively isolated, low-efficiency origins that are poor in epigenetic marks and are enriched in an asymmetric AC repeat at the initiation site. Late origins are mainly found in this class. Class 2 origins are particularly rich in enhancer elements. Class 3 origins are the most efficient and are associated with open chromatin and polycomb protein-enriched regions. The presence of Origin G-rich Repeated elements (OGRE) potentially forming G-quadruplexes (G4) was confirmed at most origins. These coincide with nucleosome-depleted regions located upstream of the initiation sites, which are associated with a labile nucleosome containing H3K64ac. These data demonstrate that specific chromatin landscapes and combinations of specific signatures regulate origin localization. They explain the frequently observed links between DNA replication and transcription. They also emphasize the plasticity of metazoan replication origins and suggest that in multicellular eukaryotes, the combination of distinct genetic features and chromatin configurations act in synergy to define and adapt the origin profile.

Capturing snapshots of APE1 processing DNA damage

DNA apurinic-apyrimidinic (AP) sites are prevalent noncoding threats to genomic stability and are processed by AP endonuclease 1 (APE1). APE1 incises the AP-site phosphodiester backbone, generating a DNA-repair intermediate that is potentially cytotoxic. The molecular events of the incision reaction remain elusive, owing in part to limited structural information. We report multiple high-resolution human APE1-DNA structures that divulge new features of the APE1 reaction, including the metal-binding site, the nucleophile and the arginine clamps that mediate product release. We also report APE1-DNA structures with a T-G mismatch 5' to the AP site, representing a clustered lesion occurring in methylated CpG dinucleotides. These structures reveal that APE1 molds the T-G mismatch into a unique Watson-Crick-like geometry that distorts the active site, thus reducing incision. These snapshots provide mechanistic clarity for APE1 while affording a rational framework to manipulate biological responses to DNA damage.

DNA replication origin activation in space and time

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

HsOrc4-Dependent Dna Remodeling of the ori-β Dhfr Replicator

Replication of DNA in multicellular organisms initiates from origin of replication (ori) sequences, which significantly differ in length and complexity. One of the best characterized is hamster dihydrofolate reductase (DHFR), which contains the ori-β sequence with several functionally relevant domains, such as an AT-rich region, dinucleotide repeat element (DNR), sequence-induced bend DNA (BEND) and a RIP60 protein-binding site (RIP60). Prior to initiation, ori sequences are recognized by origin recognition complex (ORC), which is a hetero hexamer complex that serves as the landing pad for proteins of the pre-replication complex. The function of each ORC subunit is still unclear. In this study, we analyze the function of subunit 4 of the human ORC complex (HsOrc4) in interaction with a plasmid bearing the ori-β DHFR sequence. We show that the topologically closed DHFR ori-β replicator contains a bubble-like structure within its AT-rich region and that it is reversibly modified in the interaction with HsOrc4. The non-canonical structure of the AT-rich region in the topologically closed ori sequence is recognized and changed by HsOrc4 using the energy of supercoiled DNA. These findings could help to further elucidate DNA replication and its possible association with human genetic diseases.

Allele-specific analysis of DNA replication origins in mammalian cells

The mechanisms that control the location and timing of firing of replication origins are poorly understood. Using a novel functional genomic approach based on the analysis of SNPs and indels in phased human genomes, we observe that replication asynchrony is associated with small cumulative variations in the initiation efficiency of multiple origins between the chromosome homologues, rather than with the activation of dormant origins. Allele-specific measurements demonstrate that the presence of G-quadruplex-forming sequences does not correlate with the efficiency of initiation. Sequence analysis reveals that the origins are highly enriched in sequences with profoundly asymmetric G/C and A/T nucleotide distributions and are almost completely depleted of antiparallel triplex-forming sequences. We therefore propose that although G4-forming sequences are abundant in replication origins, an asymmetry in nucleotide distribution, which increases the propensity of origins to unwind and adopt non-B DNA structure, rather than the ability to form G4, is directly associated with origin activity.

The hunt for origins of DNA replication in multicellular eukaryotes

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

Characterizing and controlling intrinsic biases of lambda exonuclease in nascent strand sequencing reveals phasing between nucleosomes and G-quadruplex motifs around a subset of human replication origins

Nascent strand sequencing (NS-seq) is used to discover DNA replication origins genome-wide, allowing identification of features for their specification. NS-seq depends on the ability of lambda exonuclease (λ-exo) to efficiently digest parental DNA while leaving RNA-primer protected nascent strands intact. We used genomics and biochemical approaches to determine if λ-exo digests all parental DNA sequences equally. We report that λ-exo does not efficiently digest G-quadruplex (G4) structures in a plasmid. Moreover, λ-exo digestion of nonreplicating genomic DNA (LexoG0) enriches GC-rich DNA and G4 motifs genome-wide. We used LexoG0 data to control for nascent strand–independent λ-exo biases in NS-seq and validated this approach at the rDNA locus. The λ-exo–controlled NS-seq peaks are not GC-rich, and only 35.5% overlap with 6.8% of all G4s, suggesting that G4s are not general determinants for origin specification but may play a role for a subset. Interestingly, we observed a periodic spacing of G4 motifs and nucleosomes around the peak summits, suggesting that G4s may position nucleosomes at this subset of origins. Finally, we demonstrate that use of Na⁺ instead of K⁺ in the λ-exo digestion buffer reduced the effect of G4s on λ-exo digestion and discuss ways to increase both the sensitivity and specificity of NS-seq.

Possible role of H1 histone in replication timing

AT-rich repetitive DNA sequences become late replicating during cell differentiation. Replication timing is not correlated with LINE density in human cells (Ryba et al. 2010). However, short and properly spaced runs of oligo dA or dT present in nuclear matrix attachment regions (MARs) of the genome are good candidates for elements of AT-rich repetitive late replicating DNA. MAR attachment to the nuclear matrix is negatively regulated by chromatin binding of H1 histone, but this is counteracted by H1 phosphorylation, high mobility group proteins or, indirectly, core histone acetylation. Fewer MAR attachments correlates positively with longer average DNA loop size, longer replicons and an increase of late replicating DNA.

Evidence for evolutionary and nonevolutionary forces shaping the distribution of human genetic variants near transcription start sites

The regions surrounding transcription start sites (TSSs) of genes play a critical role in the regulation of gene expression. At the same time, current evidence indicates that these regions are particularly stressed by transcription-related mutagenic phenomena. In this work we performed a genome-wide analysis of the distribution of single nucleotide polymorphisms (SNPs) inside the 10 kb region flanking human TSSs by dividing SNPs into four classes according to their frequency (rare, two intermediate classes, and common). We found that, in this 10 kb region, the distribution of variants depends on their frequency and on their localization relative to the TSS. We found that the distribution of variants is generally different for TSSs located inside or outside of CpG islands. We found a significant relationship between the distribution of rare variants and nucleosome occupancy scores. Furthermore, our analysis suggests that evolutionary (purifying selection) and nonevolutionary (biased gene conversion) forces both play a role in determining the relative SNP frequency around TSSs. Finally, we analyzed the potential pathogenicity of each class of variant using the Combined Annotation Dependent Depletion score. In conclusion, this study provides a novel and detailed view of the distribution of genomic variants around TSSs, providing insight into the forces that instigate and maintain variability in such critical regions.

Allele-specific genome-wide profiling in human primary erythroblasts reveal replication program organization

We have developed a new approach to characterize allele-specific timing of DNA replication genome-wide in human primary basophilic erythroblasts. We show that the two chromosome homologs replicate at the same time in about 88% of the genome and that large structural variants are preferentially associated with asynchronous replication. We identified about 600 megabase-sized asynchronously replicated domains in two tested individuals. The longest asynchronously replicated domains are enriched in imprinted genes suggesting that structural variants and parental imprinting are two causes of replication asynchrony in the human genome. Biased chromosome X inactivation in one of the two individuals tested was another source of detectable replication asynchrony. Analysis of high-resolution TimEX profiles revealed small variations termed timing ripples, which were undetected in previous, lower resolution analyses. Timing ripples reflect highly reproducible, variations of the timing of replication in the 100 kb-range that exist within the well-characterized megabase-sized replication timing domains. These ripples correspond to clusters of origins of replication that we detected using novel nascent strands DNA profiling methods. Analysis of the distribution of replication origins revealed dramatic differences in initiation of replication frequencies during S phase and a strong association, in both synchronous and asynchronous regions, between origins of replication and three genomic features: G-quadruplexes, CpG Islands and transcription start sites. The frequency of initiation in asynchronous regions was similar in the two homologs. Asynchronous regions were richer in origins of replication than synchronous regions.

DNA replication in mammalian cells proceeds according to a distinct order. Genes that are expressed tend to replicate before genes that are not expressed. We report here that we have developed a method to measure the timing of replication of the maternal and paternal chromosomes separately. We found that the paternal and maternal chromosomes replicate at exactly the same time in the large majority of the genome and that the 12% of the genome that replicated asynchronously was enriched in imprinted genes and in structural variants. Previous experiments have shown that chromosomes could be divided into replication timing domains that are a few hundred thousand to a few megabases in size. We show here that these domains can be divided into sub-domains defined by ripples in the timing profile. These ripples corresponded to clusters of origins of replication. Finally, we show that the frequency of initiation in asynchronous regions was similar in the two homologs.

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

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

Epigenetic landscape for initiation of DNA replication

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

DNA replication origins

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

GC skew at the 5' and 3' ends of human genes links R-loop formation to epigenetic regulation and transcription termination

Strand asymmetry in the distribution of guanines and cytosines, measured by GC skew, predisposes DNA sequences toward R-loop formation upon transcription. Previous work revealed that GC skew and R-loop formation associate with a core set of unmethylated CpG island (CGI) promoters in the human genome. Here, we show that GC skew can distinguish four classes of promoters, including three types of CGI promoters, each associated with unique epigenetic and gene ontology signatures. In particular, we identify a strong and a weak class of CGI promoters and show that these loci are enriched in distinct chromosomal territories reflecting the intrinsic strength of their protection against DNA methylation. Interestingly, we show that strong CGI promoters are depleted from the X chromosome while weak CGIs are enriched, a property consistent with the acquisition of DNA methylation during dosage compensation. Furthermore, we identify a third class of CGI promoters based on its unique GC skew profile and show that this gene set is enriched for Polycomb group targets. Lastly, we show that nearly 2000 genes harbor GC skew at their 3′ ends and that these genes are preferentially located in gene-dense regions and tend to be closely arranged. Genomic profiling of R-loops accordingly showed that a large proportion of genes with terminal GC skew form R-loops at their 3′ ends, consistent with a role for these structures in permitting efficient transcription termination. Altogether, we show that GC skew and R-loop formation offer significant insights into the epigenetic regulation, genomic organization, and function of human genes.

[Regulation of DNA replication timing]

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

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.