Replication timing of the human genome

KMT2C/KMT2D-dependent H3K4me1 mediates changes in DNA replication timing and origin activity during a cell fate transition

Mammalian genomes replicate in a cell-type-specific order during the S phase, correlated to transcriptional activity, histone modifications, and chromatin structure. The causal relationships between these features and DNA replication timing (RT), especially during cell fate changes, are largely unknown. Using machine learning, we quantify 21 chromatin features predicting local RT and RT changes during differentiation in embryonic stem cells (ESCs). About one-third of the genome shows RT changes during differentiation. Chromatin features accurately predict both steady-state RT and RT changes. Histone H3 lysine 4 monomethylation (H3K4me1), catalyzed by KMT2C and KMT2D (KMT2C/D), emerges as a top predictor. Loss of KMT2C/D or their enzymatic activities impairs RT changes during differentiation. This correlates with local H3K4me1 loss and reduced replication origin firing, while transcription remains largely unaffected. Our findings reveal KMT2C/D-dependent H3K4me1 as a key regulator of RT and replication initiation, a role that likely impacts diseases associated with KMT2C/D mutations.

Genome organization and stability in mammalian pre-implantation development

A largely stable genome is required for normal development, even as genetic change is an integral aspect of reproduction, genetic adaptation and evolution. Recent studies highlight a critical window of mammalian development with intrinsic DNA replication stress and genome instability in the first cell divisions after fertilization. Patterns of DNA replication and genome stability are established very early in mammals, alongside patterns of nuclear organization, and before the emergence of gene expression patterns, and prior to cell specification and germline formation. The study of DNA replication and genome stability in the mammalian embryo provides a unique cellular system due to the resetting of the epigenome to a totipotent state, and the de novo establishment of the patterns of nuclear organization, gene expression, DNA methylation, histone modifications and DNA replication. Studies on DNA replication and genome stability in the early mammalian embryo is relevant for understanding both normal and disease-causing genetic variation, and to uncover basic principles of genome regulation.

MLL3/MLL4 enzymatic activity shapes DNA replication timing

Mammalian genomes are replicated in a precise order during S phase, which is cell-type-specific1-3 and correlates with local transcriptional activity2,4-8, chromatin modifications9 and chromatin architecture1,10,11,12. However, the causal relationships between these features and the key regulators of DNA replication timing (RT) are largely unknown. Here, machine learning was applied to quantify chromatin features, including epigenetic marks, histone variants and chromatin architectural factors, best predicting local RT under steady-state and RT changes during early embryonic stem (ES) cell differentiation. About one-third of genome exhibited RT changes during the differentiation. Combined, chromatin features predicted steady-state RT and RT changes with high accuracy. Of these features, histone H3 lysine 4 monomethylation (H3K4me1) catalyzed by MLL3/4 (also known as KMT2C/D) emerged as a top predictor. Loss of Mll3/4 (but not Mll3 alone) or their enzymatic activity resulted in erasure of genome-wide RT dynamics during ES cell differentiation. Sites that normally gain H3K4me1 in a MLL3/4-dependent fashion during the transition failed to transition towards earlier RT, often with transcriptional activation unaffected. Further analysis revealed a requirement for MLL3/4 in promoting DNA replication initiation zones through MCM2 recruitment, providing a direct link for its role in regulating RT. Our results uncover MLL3/4-dependent H3K4me1 as a functional regulator of RT and highlight a causal relationship between the epigenome and RT that is largely uncoupled from transcription. These findings uncover a previously unknown role for MLL3/4-dependent chromatin functions which is likely relevant to the numerous diseases associated with MLL3/4 mutations.

Where and when to start: Regulating DNA replication origin activity in eukaryotic genomes

In eukaryotic genomes, hundreds to thousands of potential start sites of DNA replication named origins are dispersed across each of the linear chromosomes. During S-phase, only a subset of origins is selected in a stochastic manner to assemble bidirectional replication forks and initiate DNA synthesis. Despite substantial progress in our understanding of this complex process, a comprehensive 'identity code' that defines origins based on specific nucleotide sequences, DNA structural features, the local chromatin environment, or 3D genome architecture is still missing. In this article, we review the genetic and epigenetic features of replication origins in yeast and metazoan chromosomes and highlight recent insights into how this flexibility in origin usage contributes to nuclear organization, cell growth, differentiation, and genome stability.

Genome-Wide Extrachromosomal Circular DNA Profiling of Paired Hepatocellular Carcinoma and Adjacent Liver Tissues

Hepatocellular carcinoma (HCC) develops through multiple mechanisms. While recent studies have shown the presence of extrachromosomal circular DNA (eccDNA) in most cancer types, the eccDNA expression pattern and its association with HCC remain obscure. We aimed to investigate this problem. The genome-wide eccDNA profiles of eight paired HCC and adjacent non-tumor tissue samples were comprehensively elucidated based on Circle-seq, and they were further cross-analyzed with the RNA sequencing data to determine the association between eccDNA expression and transcriptome dysregulation. A total of 60,423 unique eccDNA types were identified. Most of the detected eccDNAs were smaller than 1 kb, with a length up to 182,363 bp and a mean sizes of 674 bp (non-tumor) and 813 bp (tumor), showing a greater association with gene-rich rather than with gene-poor regions. Although there was no statistical difference in length and chromosome distribution, the eccDNA patterns between HCC and adjacent non-tumor tissues showed significant differences at both the chromosomal and single gene levels. Five of the eight HCC tissues showed significantly higher amounts of chromosome 22-derived eccDNA expression compared to the non-tumor tissue. Furthermore, two genes, SLC16A3 and BAIAP2L2, with a higher transcription level in tumor tissues, were related to eccDNAs exclusively detected in three HCC samples and were negatively associated with survival rates in HCC cohorts from public databases. These results indicate the existence and massive heterogeneity of eccDNAs in HCC and adjacent liver tissues, and suggest their potential association with dysregulated gene expression.

Emergence and influence of sequence bias in evolutionarily malleable, mammalian tandem arrays

BACKGROUND

The radiation of mammals at the extinction of the dinosaurs produced a plethora of new forms-as diverse as bats, dolphins, and elephants-in only 10-20 million years. Behind the scenes, adaptation to new niches is accompanied by extensive innovation in large families of genes that allow animals to contact the environment, including chemosensors, xenobiotic enzymes, and immune and barrier proteins. Genes in these "outward-looking" families are allelically diverse among humans and exhibit tissue-specific and sometimes stochastic expression.

RESULTS

Here, we show that these tandem arrays of outward-looking genes occupy AT-biased isochores and comprise the "tissue-specific" gene class that lack CpG islands in their promoters. Models of mammalian genome evolution have not incorporated the sharply different functions and transcriptional patterns of genes in AT- versus GC-biased regions. To examine the relationship between gene family expansion, sequence content, and allelic diversity, we use population genetic data and comparative analysis. First, we find that AT bias can emerge during evolutionary expansion of gene families in cis. Second, human genes in AT-biased isochores or with GC-poor promoters experience relatively low rates of de novo point mutation today but are enriched for non-synonymous variants. Finally, we find that isochores containing gene clusters exhibit low rates of recombination.

CONCLUSIONS

Our analyses suggest that tolerance of non-synonymous variation and low recombination are two forces that have produced the depletion of GC bases in outward-facing gene arrays. In turn, high AT content exerts a profound effect on their chromatin organization and transcriptional regulation.

Genomic mutation landscape of skin cancers from DNA repair-deficient xeroderma pigmentosum patients

Xeroderma pigmentosum (XP) is a genetic disorder caused by mutations in genes of the Nucleotide Excision Repair (NER) pathway (groups A-G) or in Translesion Synthesis DNA polymerase η (V). XP is associated with an increased skin cancer risk, reaching, for some groups, several thousand-fold compared to the general population. Here, we analyze 38 skin cancer genomes from five XP groups. We find that the activity of NER determines heterogeneity of the mutation rates across skin cancer genomes and that transcription-coupled NER extends beyond the gene boundaries reducing the intergenic mutation rate. Mutational profile in XP-V tumors and experiments with POLH knockout cell line reveal the role of polymerase η in the error-free bypass of (i) rare TpG and TpA DNA lesions, (ii) 3' nucleotides in pyrimidine dimers, and (iii) TpT photodimers. Our study unravels the genetic basis of skin cancer risk in XP and provides insights into the mechanisms reducing UV-induced mutagenesis in the general population.

Strand asymmetries across genomic processes

Across biological systems, a number of genomic processes, including transcription, replication, DNA repair, and transcription factor binding, display intrinsic directionalities. These directionalities are reflected in the asymmetric distribution of nucleotides, motifs, genes, transposon integration sites, and other functional elements across the two complementary strands. Strand asymmetries, including GC skews and mutational biases, have shaped the nucleotide composition of diverse organisms. The investigation of strand asymmetries often serves as a method to understand underlying biological mechanisms, including protein binding preferences, transcription factor interactions, retrotransposition, DNA damage and repair preferences, transcription-replication collisions, and mutagenesis mechanisms. Research into this subject also enables the identification of functional genomic sites, such as replication origins and transcription start sites. Improvements in our ability to detect and quantify DNA strand asymmetries will provide insights into diverse functionalities of the genome, the contribution of different mutational mechanisms in germline and somatic mutagenesis, and our knowledge of genome instability and evolution, which all have significant clinical implications in human disease, including cancer. In this review, we describe key developments that have been made across the field of genomic strand asymmetries, as well as the discovery of associated mechanisms.

Optimized Repli-seq: improved DNA replication timing analysis by next-generation sequencing

The human genome is divided into functional units that replicate at specific times during S-phase. This temporal program is known as replication timing (RT) and is coordinated with the spatial organization of the genome and transcriptional activity. RT is also cell type-specific, dynamically regulated during development, and alterations in RT are observed in multiple diseases. Thus, the precise measure of RT is critical to understand the role of RT in gene function regulation. Distinct methods for assaying the RT program exist; however, conventional methods require thousands of cells as input, prohibiting its applicability to samples with limited cell numbers such as those from disease patients or from early developing embryos. Although single-cell RT analyses have been developed, these methods are low throughput, require generation of numerous libraries, increased sequencing costs, and produce low resolution data. Here, we developed an improved method to measure RT genome-wide that enables high-resolution analysis of low input samples. This method incorporates direct cell sorting into lysis buffer, as well as DNA fragmentation and library preparation in a single tube, resulting in higher yields, increased quality, and reproducibility with decreased costs. We also performed a systematic data processing analysis to provide standardized parameters for RT measurement. This optimized method facilitates RT analysis and will enable its application to a broad range of studies investigating the role of RT in gene expression, nuclear architecture, and disease.

Replication stress generates distinctive landscapes of DNA copy number alterations and chromosome scale losses

BACKGROUND

A major driver of cancer chromosomal instability is replication stress, the slowing or stalling of DNA replication. How replication stress and genomic instability are connected is not known. Aphidicolin-induced replication stress induces breakages at common fragile sites, but the exact causes of fragility are debated, and acute genomic consequences of replication stress are not fully explored.

RESULTS

We characterize DNA copy number alterations (CNAs) in single, diploid non-transformed cells, caused by one cell cycle in the presence of either aphidicolin or hydroxyurea. Multiple types of CNAs are generated, associated with different genomic regions and features, and observed copy number landscapes are distinct between aphidicolin and hydroxyurea-induced replication stress. Coupling cell type-specific analysis of CNAs to gene expression and single-cell replication timing analyses pinpointed the causative large genes of the most recurrent chromosome-scale CNAs in aphidicolin. These are clustered on chromosome 7 in RPE1 epithelial cells but chromosome 1 in BJ fibroblasts. Chromosome arm level CNAs also generate acentric lagging chromatin and micronuclei containing these chromosomes.

CONCLUSIONS

Chromosomal instability driven by replication stress occurs via focal CNAs and chromosome arm scale changes, with the latter confined to a very small subset of chromosome regions, potentially heavily skewing cancer genome evolution. Different inducers of replication stress lead to distinctive CNA landscapes providing the opportunity to derive copy number signatures of specific replication stress mechanisms. Single-cell CNA analysis thus reveals the impact of replication stress on the genome, providing insights into the molecular mechanisms which fuel chromosomal instability in cancer.

MTBP and MYC: A Dynamic Duo in Proliferation, Cancer, and Aging

The oncogenic transcription factor c-MYC (MYC) is highly conserved across species and is frequently overexpressed or dysregulated in human cancers. MYC regulates a wide range of critical cellular and oncogenic activities including proliferation, metabolism, metastasis, apoptosis, and differentiation by transcriptionally activating or repressing the expression of a large number of genes. This activity of MYC is not carried out in isolation, instead relying on its association with a myriad of protein cofactors. We determined that MDM Two Binding Protein (MTBP) indirectly binds MYC and is a novel MYC transcriptional cofactor. MTBP promotes MYC-mediated transcriptional activity, proliferation, and cellular transformation by binding in a protein complex with MYC at MYC-bound promoters. This discovery provided critical context for data linking MTBP to aging as well as a rapidly expanding body of evidence demonstrating MTBP is overexpressed in many human malignancies, is often linked to poor patient outcomes, and is necessary for cancer cell survival. As such, MTBP represents a novel and potentially broad reaching oncologic drug target, particularly when MYC is dysregulated. Here we have reviewed the discovery of MTBP and the initial controversy with its function as well as its associations with proliferation, MYC, DNA replication, aging, and human cancer.

Transmission ratio distortion of mutations in the master regulator of centriole biogenesis PLK4

The evolutionary conserved Polo-like kinase 4 (PLK4) is essential for centriole duplication, spindle assembly, and de novo centriole formation. In man, homozygous mutations in PLK4 lead to primary microcephaly, altered PLK4 expression is associated with aneuploidy in human embryos. Here, we report on a consanguineous four-generation family with 8 affected individuals compound heterozygous for a novel missense variant, c.881 T > G, and a deletion of the PLK4 gene. The clinical phenotype of the adult patients is mild compared to individuals with previously described PLK4 mutations. One individual was homozygous for the variant c.881G and phenotypically unaffected. The deletion was inherited by 14 of 16 offspring and thus exhibits transmission ratio distortion (TRD). Moreover, based on the already published families with PLK4 mutations, it could be shown that due to the preferential transmission of the mutant alleles, the number of affected offspring is significantly increased. It is assumed that reduced expression of PLK4 decreases the intrinsically high error rate of the first cell divisions after fertilization, increases the number of viable embryos and thus leads to preferential transmission of the deleted/mutated alleles.

Chemically Induced Chromosomal Interaction (CICI) method to study chromosome dynamics and its biological roles

Numerous intra- and inter-chromosomal contacts have been mapped in eukaryotic genomes, but it remains challenging to link these 3D structures to their regulatory functions. To establish the causal relationships between chromosome conformation and genome functions, we  develop a method, Chemically Induced Chromosomal Interaction (CICI), to selectively perturb the chromosome conformation at targeted loci. Using this method, long-distance chromosomal interactions can be induced dynamically between intra- or inter-chromosomal loci pairs, including the ones with very low Hi-C contact frequencies. Measurement of CICI formation time allows us to probe chromosome encounter dynamics between different loci pairs across the cell cycle. We also conduct two functional tests of CICI. We perturb the chromosome conformation near a DNA double-strand break and observe altered donor preference in homologous recombination; we force interactions between early and late-firing DNA replication origins and find no significant changes in replication timing. These results suggest that chromosome conformation plays a deterministic role in homology-directed DNA repair, but not in the establishment of replication timing. Overall, our study demonstrates that CICI is a powerful tool to study chromosome dynamics and 3D genome function.

Methods to selectively manipulate specific long-distance chromosomal interactions are limited. Here the authors develop a method called Chemically Induced Chromosomal Interaction (CICI) to engineer interactions and demonstrate that 3D conformation plays a causal role in establishing donor DNA preference during DNA repair.

Transcription-Replication Coordination

Transcription and replication are the two most essential processes that a cell does with its DNA: they allow cells to express the genomic content that is required for their functions and to create a perfect copy of this genomic information to pass on to the daughter cells. Nevertheless, these two processes are in a constant ambivalent relationship. When transcription and replication occupy the same regions, there is the possibility of conflicts between transcription and replication as transcription can impair DNA replication progression leading to increased DNA damage. Nevertheless, DNA replication origins are preferentially located in open chromatin next to actively transcribed regions, meaning that the possibility of conflicts is potentially an accepted incident for cells. Data in the literature point both towards the existence or not of coordination between these two processes to avoid the danger of collisions. Several reviews have been published on transcription-replication conflicts, but we focus here on the most recent findings that relate to how these two processes are coordinated in eukaryotes, considering advantages and disadvantages from coordination, how likely conflicts are at any given time, and which are their potential hotspots in the genome.

Chromatin loading of MCM hexamers is associated with di-/tri-methylation of histone H4K20 toward S phase entry

DNA replication is a key step in initiating cell proliferation. Loading hexameric complexes of minichromosome maintenance (MCM) helicase onto DNA replication origins during the G1 phase is essential for initiating DNA replication. Here, we examined MCM hexamer states during the cell cycle in human hTERT-RPE1 cells using multicolor immunofluorescence-based, single-cell plot analysis, and biochemical size fractionation. Experiments involving cell-cycle arrest at the G1 phase and release from the arrest revealed that a double MCM hexamer was formed via a single hexamer during G1 progression. A single MCM hexamer was recruited to chromatin in the early G1 phase. Another single hexamer was recruited to form a double hexamer in the late G1 phase. We further examined relationship between the MCM hexamer states and the methylation levels at lysine 20 of histone H4 (H4K20) and found that the double MCM hexamer state was correlated with di/trimethyl-H4K20 (H4K20me2/3). Inhibiting the conversion from monomethyl-H4K20 (H4K20me1) to H4K20me2/3 retained the cells in the single MCM hexamer state. Non-proliferative cells, including confluent cells or Cdk4/6 inhibitor-treated cells, also remained halted in the single MCM hexamer state. We propose that the single MCM hexamer state is a halting step in the determination of cell cycle progression.

Chromosomal Rearrangements and Altered Nuclear Organization: Recent Mechanistic Models in Cancer

Simple Summary

New methodologies and technologies developed in the last few decades have highlighted the precise spatial organization of the genome into the cell nucleus, with chromatin architecture playing a central role in controlling several genome functions. Genes are expressed in a well-defined way and at a well-defined time during cell differentiation, and alterations in genome organization can lead to genetic diseases, such as cancers. Here we review how the genome is organized in the cell nucleus and the evidence of genome misorganization leading to cancer diseases.

Abstract

The last decade has seen significant progress in understanding how the genome is organized spatially within interphase nuclei. Recent analyses have confirmed earlier molecular cytogenetic studies on chromosome positioning within interphase nuclei and provided new information about the topologically associated domains (TADs). Examining the nuances of how genomes are organized within interphase nuclei will provide information fundamental to understanding gene regulation and expression in health and disease. Indeed, the radial spatial positioning of individual gene loci within nuclei has been associated with up- and down-regulation of specific genes, and disruption of normal genome organization within nuclei will result in compromised cellular health. In cancer cells, where reorganization of the nuclear architecture may occur in the presence of chromosomal rearrangements such as translocations, inversions, or deletions, gene repositioning can change their expression. To date, very few studies have focused on radial gene positioning and the correlation to gene expression in cancers. Further investigations would improve our understanding of the biological mechanisms at the basis of cancer and, in particular, in leukemia initiation and progression, especially in those cases where the molecular consequences of chromosomal rearrangements are still unclear. In this review, we summarize the main milestones in the field of genome organization in the nucleus and the alterations to this organization that can lead to cancer diseases.

Deposition Bias of Chromatin Proteins Inverts under DNA Replication Stress Conditions

Following DNA replication, equal amounts of chromatin proteins are distributed over sister chromatids by re-deposition of parental chromatin proteins and deposition of newly synthesized chromatin proteins. Molecular mechanisms balancing the allocation of new and old chromatin proteins remain largely unknown. Here, we studied the genome-wide distribution of new chromatin proteins relative to parental DNA template strands and replication initiation zones using the double-click-seq. Under control conditions, new chromatin proteins were preferentially found on DNA replicated by the lagging strand machinery. Strikingly, replication stress induced by hydroxyurea or curaxin treatment and inhibition of ataxia telangiectasia and Rad3-related protein (ATR) or p53 inactivation inverted the observed chromatin protein deposition bias to the strand replicated by the leading strand polymerase in line with previously reported effects on replication protein A occupancy. We propose that asymmetric deposition of newly synthesized chromatin proteins onto sister chromatids reflects differences in the processivity of leading and lagging strand synthesis.

The 3D Genome: From Structure to Function

The genome is the most functional part of a cell, and genomic contents are organized in a compact three-dimensional (3D) structure. The genome contains millions of nucleotide bases organized in its proper frame. Rapid development in genome sequencing and advanced microscopy techniques have enabled us to understand the 3D spatial organization of the genome. Chromosome capture methods using a ligation approach and the visualization tool of a 3D genome browser have facilitated detailed exploration of the genome. Topologically associated domains (TADs), lamin-associated domains, CCCTC-binding factor domains, cohesin, and chromatin structures are the prominent identified components that encode the 3D structure of the genome. Although TADs are the major contributors to 3D genome organization, they are absent in Arabidopsis. However, a few research groups have reported the presence of TAD-like structures in the plant kingdom.

Trisomy 21-associated increases in chromosomal instability are unmasked by comparing isogenic trisomic/disomic leukocytes from people with mosaic Down syndrome

Down syndrome, which results from a trisomic imbalance for chromosome 21, has been associated with 80+ phenotypic traits. However, the cellular changes that arise in somatic cells due to this aneuploid condition are not fully understood. The primary aim of this study was to determine if germline trisomy 21 is associated with an increase in spontaneous somatic cell chromosomal instability frequencies (SCINF). To achieve this aim, we quantified SCINF in people with mosaic Down syndrome using a cytokinesis-blocked micronucleus assay. By comparing values in their isogenic trisomic/disomic cells, we obtained a measure of differences in SCINF that are directly attributable to a trisomy 21 imbalance, since differential effects attributable to "background" genetic factors and environmental exposures could be eliminated. A cross-sectional assessment of 69 people with mosaic Down syndrome (ages 1 to 44; mean age of 12.84 years) showed a significantly higher frequency of micronuclei in their trisomic (0.37 ± 0.35 [mean ± standard deviation]) compared to disomic cells (0.18 ± 0.11)(P <0.0001). The daughter binucleates also showed significantly higher levels of abnormal patterns in the trisomic (1.68 ± 1.21) compared to disomic (0.35 ± 0.45) cells (P <0.0001). Moreover, a significant Age x Cell Type interaction was noted (P = 0.0113), indicating the relationship between age and SCINF differed between the trisomic and disomic cells. Similarly, a longitudinal assessment (mean time interval of 3.9 years; range of 2 to 6 years) of 18 participants showed a mean 1.63-fold increase in SCINF within individuals over time for their trisomic cells (P = 0.0186), compared to a 1.13-fold change in their disomic cells (P = 0.0464). In summary, these results showed a trisomy 21-associated, age-related increase in SCINF. They also underscore the strength of the isogenic mosaic Down syndrome model system for "unmasking" cellular changes arising from a trisomy 21 imbalance.

Analysis of gene expression and mutation data points on contribution of transcription to the mutagenesis by APOBEC enzymes

Since the discovery of the role of the APOBEC enzymes in human cancers, the mechanisms of this type of mutagenesis remain little understood. Theoretically, targeting of single-stranded DNA by the APOBEC enzymes could occur during cellular processes leading to the unwinding of DNA double-stranded structure. Some evidence points to the importance of replication in the APOBEC mutagenesis, while the role of transcription is still underexplored. Here, we analyzed gene expression and whole genome sequencing data from five types of human cancers with substantial APOBEC activity to estimate the involvement of transcription in the APOBEC mutagenesis and compare its impact with that of replication. Using the TCN motif as the mutation signature of the APOBEC enzymes, we observed a correlation of active APOBEC mutagenesis with gene expression, confirmed the increase of APOBEC-induced mutations in early-replicating regions and estimated the relative impact of transcription and replication on the APOBEC mutagenesis. We also found that the known effect of higher density of APOBEC-induced mutations on the lagging strand was highest in middle-replicating regions and observed higher APOBEC mutation density on the sense strand, the latter bias positively correlated with the gene expression level.

Low Replicative Stress Triggers Cell-Type Specific Inheritable Advanced Replication Timing

DNA replication timing (RT), reflecting the temporal order of origin activation, is known as a robust and conserved cell-type specific process. Upon low replication stress, the slowing of replication forks induces well-documented RT delays associated to genetic instability, but it can also generate RT advances that are still uncharacterized. In order to characterize these advanced initiation events, we monitored the whole genome RT from six independent human cell lines treated with low doses of aphidicolin. We report that RT advances are cell-type-specific and involve large heterochromatin domains. Importantly, we found that some major late to early RT advances can be inherited by the unstressed next-cellular generation, which is a unique process that correlates with enhanced chromatin accessibility, as well as modified replication origin landscape and gene expression in daughter cells. Collectively, this work highlights how low replication stress may impact cellular identity by RT advances events at a subset of chromosomal domains.

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.

Replication Stress, Genomic Instability, and Replication Timing: A Complex Relationship

The replication-timing program constitutes a key element of the organization and coordination of numerous nuclear processes in eukaryotes. This program is established at a crucial moment in the cell cycle and occurs simultaneously with the organization of the genome, thus indicating the vital significance of this process. With recent technological achievements of high-throughput approaches, a very strong link has been confirmed between replication timing, transcriptional activity, the epigenetic and mutational landscape, and the 3D organization of the genome. There is also a clear relationship between replication stress, replication timing, and genomic instability, but the extent to which they are mutually linked to each other is unclear. Recent evidence has shown that replication timing is affected in cancer cells, although the cause and consequence of this effect remain unknown. However, in-depth studies remain to be performed to characterize the molecular mechanisms of replication-timing regulation and clearly identify different cis- and trans-acting factors. The results of these studies will potentially facilitate the discovery of new therapeutic pathways, particularly for personalized medicine, or new biomarkers. This review focuses on the complex relationship between replication timing, replication stress, and genomic instability.

Developmental differences in genome replication program and origin activation

To ensure error-free duplication of all (epi)genetic information once per cell cycle, DNA replication follows a cell type and developmental stage specific spatio-temporal program. Here, we analyze the spatio-temporal DNA replication progression in (un)differentiated mouse embryonic stem (mES) cells. Whereas telomeres replicate throughout S-phase, we observe mid S-phase replication of (peri)centromeric heterochromatin in mES cells, which switches to late S-phase replication upon differentiation. This replication timing reversal correlates with and depends on an increase in condensation and a decrease in acetylation of chromatin. We further find synchronous duplication of the Y chromosome, marking the end of S-phase, irrespectively of the pluripotency state. Using a combination of single-molecule and super-resolution microscopy, we measure molecular properties of the mES cell replicon, the number of replication foci active in parallel and their spatial clustering. We conclude that each replication nanofocus in mES cells corresponds to an individual replicon, with up to one quarter representing unidirectional forks. Furthermore, with molecular combing and genome-wide origin mapping analyses, we find that mES cells activate twice as many origins spaced at half the distance than somatic cells. Altogether, our results highlight fundamental developmental differences on progression of genome replication and origin activation in pluripotent cells.

The mutational signature profile of known and suspected human carcinogens in mice

Epidemiological studies have identified many environmental agents that appear to significantly increase cancer risk in human populations. By analyzing tumor genomes from mice chronically exposed to 1 of 20 known or suspected human carcinogens, we reveal that most agents do not generate distinct mutational signatures or increase mutation burden, with most mutations, including driver mutations, resulting from tissue-specific endogenous processes. We identify signatures resulting from exposure to cobalt and vinylidene chloride and link distinct human signatures (SBS19 and SBS42) with 1,2,3-trichloropropane, a haloalkane and pollutant of drinking water, and find these and other signatures in human tumor genomes. We define the cross-species genomic landscape of tumors induced by an important compendium of agents with relevance to human health.

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.

DNA copy-number measurement of genome replication dynamics by high-throughput sequencing: the sort-seq, sync-seq and MFA-seq family

Genome replication follows a defined temporal programme that can change during cellular differentiation and disease onset. DNA replication results in an increase in DNA copy number that can be measured by high-throughput sequencing. Here we present a protocol to determine genome replication dynamics using DNA copy-number measurements. Cell populations can be obtained in three variants of the method. First, sort-seq reveals the average replication dynamics across S phase in an unperturbed cell population; FACS is used to isolate replicating and non-replicating subpopulations from asynchronous cells. Second, sync-seq measures absolute replication time at specific points during S phase using a synchronized cell population. Third, marker frequency analysis can be used to reveal the average replication dynamics using copy-number analysis in any proliferating asynchronous cell culture. These approaches have been used to reveal genome replication dynamics in prokaryotes, archaea and a wide range of eukaryotes, including yeasts and mammalian cells. We have found this approach straightforward to apply to other organisms and highlight example studies from across the three domains of life. Here we present a Saccharomyces cerevisiae version of the protocol that can be performed in 7-10 d. It requires basic molecular and cellular biology skills, as well as a basic understanding of Unix and R.

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.

Clonal Decomposition and DNA Replication States Defined by Scaled Single-Cell Genome Sequencing

Accurate measurement of clonal genotypes, mutational processes, and replication states from individual tumor-cell genomes will facilitate improved understanding of tumor evolution. We have developed DLP+, a scalable single-cell whole-genome sequencing platform implemented using commodity instruments, image-based object recognition, and open source computational methods. Using DLP+, we have generated a resource of 51,926 single-cell genomes and matched cell images from diverse cell types including cell lines, xenografts, and diagnostic samples with limited material. From this resource we have defined variation in mitotic mis-segregation rates across tissue types and genotypes. Analysis of matched genomic and image measurements revealed correlations between cellular morphology and genome ploidy states. Aggregation of cells sharing copy number profiles allowed for calculation of single-nucleotide resolution clonal genotypes and inference of clonal phylogenies and avoided the limitations of bulk deconvolution. Finally, joint analysis over the above features defined clone-specific chromosomal aneuploidy in polyclonal populations.

Synchronized replication of genes encoding the same protein complex in fast-proliferating cells

DNA replication perturbs the dosage balance among genes; at mid-S phase, early-replicating genes have doubled their copies while late-replicating ones have not. Dosage imbalance among genes, especially within members of a protein complex, is toxic to cells. However, the molecular mechanisms that cells use to deal with such imbalance remain not fully understood. Here, we validate at the genomic scale that the dosage between early- and late-replicating genes is imbalanced in HeLa cells. We propose the synchronized replication hypothesis that genes sensitive to stoichiometric relationships will be replicated simultaneously to maintain stoichiometry. In support of this hypothesis, we observe that genes encoding the same protein complex have similar replication timing but mainly in fast-proliferating cells such as embryonic stem cells and cancer cells. We find that the synchronized replication observed in cancer cells, but not in slow-proliferating differentiated cells, is due to convergent evolution during tumorigenesis that restores synchronized replication timing within protein complexes. Taken together, our study reveals that the demand for dosage balance during S phase plays an important role in the optimization of the replication-timing program; this selection is relaxed during differentiation as the cell cycle prolongs and is restored during tumorigenesis as the cell cycle shortens.

Germline DNA replication timing shapes mammalian genome composition

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

Developmental and cancer-associated plasticity of DNA replication preferentially targets GC-poor, lowly expressed and late-replicating regions

The spatiotemporal program of metazoan DNA replication is regulated during development and altered in cancers. We have generated novel OK-seq, Repli-seq and RNA-seq data to compare the DNA replication and gene expression programs of twelve cancer and non-cancer human cell types. Changes in replication fork directionality (RFD) determined by OK-seq are widespread but more frequent within GC-poor isochores and largely disconnected from transcription changes. Cancer cell RFD profiles cluster with non-cancer cells of similar developmental origin but not with different cancer types. Importantly, recurrent RFD changes are detected in specific tumour progression pathways. Using a model for establishment and early progression of chronic myeloid leukemia (CML), we identify 1027 replication initiation zones (IZs) that progressively change efficiency during long-term expression of the BCR-ABL1 oncogene, being twice more often downregulated than upregulated. Prolonged expression of BCR-ABL1 results in targeting of new IZs and accentuation of previous efficiency changes. Targeted IZs are predominantly located in GC-poor, late replicating gene deserts and frequently silenced in late CML. Prolonged expression of BCR-ABL1 results in massive deletion of GC-poor, late replicating DNA sequences enriched in origin silencing events. We conclude that BCR-ABL1 expression progressively affects replication and stability of GC-poor, late-replicating regions during CML progression.

Rapid Irreversible Transcriptional Reprogramming in Human Stem Cells Accompanied by Discordance between Replication Timing and Chromatin Compartment

The temporal order of DNA replication is regulated during development and is highly correlated with gene expression, histone modifications and 3D genome architecture. We tracked changes in replication timing, gene expression, and chromatin conformation capture (Hi-C) A/B compartments over the first two cell cycles during differentiation of human embryonic stem cells to definitive endoderm. Remarkably, transcriptional programs were irreversibly reprogrammed within the first cell cycle and were largely but not universally coordinated with replication timing changes. Moreover, changes in A/B compartment and several histone modifications that normally correlate strongly with replication timing showed weak correlation during the early cell cycles of differentiation but showed increased alignment in later differentiation stages and in terminally differentiated cell lines. Thus, epigenetic cell fate transitions during early differentiation can occur despite dynamic and discordant changes in otherwise highly correlated genomic properties.

Next-Generation Sequencing Enables Spatiotemporal Resolution of Human Centromere Replication Timing

Centromeres serve a critical function in preserving genome integrity across sequential cell divisions, by mediating symmetric chromosome segregation. The repetitive, heterochromatic nature of centromeres is thought to be inhibitory to DNA replication, but has also led to their underrepresentation in human reference genome assemblies. Consequently, centromeres have been excluded from genomic replication timing analyses, leaving their time of replication unresolved. However, the most recent human reference genome, hg38, included models of centromere sequences. To establish the experimental requirements for achieving replication timing profiles for centromeres, we sequenced G₁- and S-phase cells from five human cell lines, and aligned the sequence reads to hg38. We were able to infer DNA replication timing profiles for the centromeres in each of the five cell lines, which showed that centromere replication occurs in mid-to-late S phase. Furthermore, we found that replication timing was more variable between cell lines in the centromere regions than expected, given the distribution of variation in replication timing genome-wide. These results suggest the potential of these, and future, sequence models to enable high-resolution studies of replication in centromeres and other heterochromatic regions.

DNA Replication Timing Enters the Single-Cell Era

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

DNA replication dynamics of vole genome and its epigenetic regulation

Background

The genome of some vole rodents exhibit large blocks of heterochromatin coupled to their sex chromosomes. The DNA composition and transcriptional activity of these heterochromatin blocks have been studied, but little is known about their DNA replication dynamics and epigenetic composition.

Results

Here, we show prominent epigenetic marks of the heterochromatic blocks in the giant sex chromosomes of female Microtus cabrerae cells. While the X chromosomes are hypoacetylated and cytosine hypomethylated, they are either enriched for macroH2A and H3K27me3 typical for facultative heterochromatin or for H3K9me3 and HP1 beta typical for constitutive heterochromatin. Using pulse-chase replication labeling and time-lapse microscopy, we found that the heterochromatic block enriched for macroH2A/H3K27me3 of the X chromosome is replicated during mid-S-phase, prior to the heterochromatic block enriched for H3K9me3/HP1 beta, which is replicated during late S-phase. To test whether histone acetylation level regulates its replication dynamics, we induced either global hyperacetylation by pharmacological inhibition or by targeting a histone acetyltransferase to the heterochromatic region of the X chromosomes. Our data reveal that histone acetylation level affects DNA replication dynamics of the sex chromosomes’ heterochromatin and leads to a global reduction in replication fork rate genome wide.

Conclusions

In conclusion, we mapped major epigenetic modifications controlling the structure of the sex chromosome-associated heterochromatin and demonstrated the occurrence of differences in the molecular mechanisms controlling the replication timing of the heterochromatic blocks at the sex chromosomes in female Microtus cabrerae cells. Furthermore, we highlighted a conserved role of histone acetylation level on replication dynamics across mammalian species.

Electronic supplementary material

The online version of this article (10.1186/s13072-019-0262-0) contains supplementary material, which is available to authorized users.

Genome-wide stability of the DNA replication program in single mammalian cells

Here, we report a single-cell DNA replication sequencing method, scRepli-seq, a genome-wide methodology that measures copy number differences between replicated and unreplicated DNA. Using scRepli-seq, we demonstrate that replication-domain organization is conserved among individual mouse embryonic stem cells (mESCs). Differentiated mESCs exhibited distinct profiles, which were also conserved among cells. Haplotype-resolved scRepli-seq revealed similar replication profiles of homologous autosomes, while the inactive X chromosome was clearly replicated later than its active counterpart. However, a small degree of cell-to-cell replication-timing heterogeneity was present, which was smallest at the beginning and the end of S phase. In addition, developmentally regulated domains were found to deviate from others and showed a higher degree of heterogeneity, thus suggesting a link to developmental plasticity. Moreover, allelic expression imbalance was found to strongly associate with replication-timing asynchrony. Our results form a foundation for single-cell-level understanding of DNA replication regulation and provide insights into three-dimensional genome organization.

DNA replication and repair kinetics of Alu, LINE-1 and satellite III genomic repetitive elements

Background

Preservation of genome integrity by complete, error-free DNA duplication prior to cell division and by correct DNA damage repair is paramount for the development and maintenance of an organism. This holds true not only for protein-encoding genes, but also it applies to repetitive DNA elements, which make up more than half of the human genome. Here, we focused on the replication and repair kinetics of interspersed and tandem repetitive DNA elements.

Results

We integrated genomic population level data with a single cell immunofluorescence in situ hybridization approach to simultaneously label replication/repair and repetitive DNA elements. We found that: (1) the euchromatic Alu element was replicated during early S-phase; (2) LINE-1, which is associated with AT-rich genomic regions, was replicated throughout S-phase, with the majority being replicated according to their particular histone marks; (3) satellite III, which constitutes pericentromeric heterochromatin, was replicated exclusively during the mid-to-late S-phase. As for the DNA double-strand break repair process, we observed that Alu elements followed the global genome repair kinetics, while LINE-1 elements repaired at a slower rate. Finally, satellite III repeats were repaired at later time points.

Conclusions

We conclude that the histone modifications in the specific repeat element predominantly determine its replication and repair timing. Thus, Alu elements, which are characterized by euchromatic chromatin features, are repaired and replicated the earliest, followed by LINE-1 elements, including more variegated eu/heterochromatic features and, lastly, satellite tandem repeats, which are homogeneously characterized by heterochromatic features and extend over megabase-long genomic regions. Altogether, this work reemphasizes the need for complementary approaches to achieve an integrated and comprehensive investigation of genomic processes.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0226-9) contains supplementary material, which is available to authorized users.

Functional dissection of Drosophila melanogaster SUUR protein influence on H3K27me3 profile

BACKGROUND

In eukaryotes, heterochromatin replicates late in S phase of the cell cycle and contains specific covalent modifications of histones. SuUR mutation found in Drosophila makes heterochromatin replicate earlier than in wild type and reduces the level of repressive histone modifications. SUUR protein was shown to be associated with moving replication forks, apparently through the interaction with PCNA. The biological process underlying the effects of SUUR on replication and composition of heterochromatin remains unknown.

RESULTS

Here we performed a functional dissection of SUUR protein effects on H3K27me3 level. Using hidden Markow model-based algorithm we revealed SuUR-sensitive chromosomal regions that demonstrated unusual characteristics: They do not contain Polycomb and require SUUR function to sustain H3K27me3 level. We tested the role of SUUR protein in the mechanisms that could affect H3K27me3 histone levels in these regions. We found that SUUR does not affect the initial H3K27me3 pattern formation in embryogenesis or Polycomb distribution in the chromosomes. We also ruled out the possible effect of SUUR on histone genes expression and its involvement in DSB repair.

CONCLUSIONS

Obtained results support the idea that SUUR protein contributes to the heterochromatin maintenance during the chromosome replication. A model that explains major SUUR-associated phenotypes is proposed.

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

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

Recent Advancements in DNA Damage-Transcription Crosstalk and High-Resolution Mapping of DNA Breaks

Until recently, DNA damage arising from physiological DNA metabolism was considered a detrimental by-product for cells. However, an increasing amount of evidence has shown that DNA damage could have a positive role in transcription activation. In particular, DNA damage has been detected in transcriptional elements following different stimuli. These physiological DNA breaks are thought to be instrumental for the correct expression of genomic loci through different mechanisms. In this regard, although a plethora of methods are available to precisely map transcribed regions and transcription start sites, commonly used techniques for mapping DNA breaks lack sufficient resolution and sensitivity to draw a robust correlation between DNA damage generation and transcription. Recently, however, several methods have been developed to map DNA damage at single-nucleotide resolution, thus providing a new set of tools to correlate DNA damage and transcription. Here, we review how DNA damage can positively regulate transcription initiation, the current techniques for mapping DNA breaks at high resolution, and how these techniques can benefit future studies of DNA damage and transcription.

Replicating Large Genomes: Divide and Conquer

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

Repliscan: a tool for classifying replication timing regions

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

The three-dimensional genome organization of Drosophila melanogaster through data integration

BACKGROUND

Genome structures are dynamic and non-randomly organized in the nucleus of higher eukaryotes. To maximize the accuracy and coverage of three-dimensional genome structural models, it is important to integrate all available sources of experimental information about a genome's organization. It remains a major challenge to integrate such data from various complementary experimental methods. Here, we present an approach for data integration to determine a population of complete three-dimensional genome structures that are statistically consistent with data from both genome-wide chromosome conformation capture (Hi-C) and lamina-DamID experiments.

RESULTS

Our structures resolve the genome at the resolution of topological domains, and reproduce simultaneously both sets of experimental data. Importantly, this data deconvolution framework allows for structural heterogeneity between cells, and hence accounts for the expected plasticity of genome structures. As a case study we choose Drosophila melanogaster embryonic cells, for which both data types are available. Our three-dimensional genome structures have strong predictive power for structural features not directly visible in the initial data sets, and reproduce experimental hallmarks of the D. melanogaster genome organization from independent and our own imaging experiments. Also they reveal a number of new insights about genome organization and its functional relevance, including the preferred locations of heterochromatic satellites of different chromosomes, and observations about homologous pairing that cannot be directly observed in the original Hi-C or lamina-DamID data.

CONCLUSIONS

Our approach allows systematic integration of Hi-C and lamina-DamID data for complete three-dimensional genome structure calculation, while also explicitly considering genome structural variability.

Genomic properties of chromosomal bands are linked to evolutionary rearrangements and new centromere formation in primates

Chromosomal rearrangements in humans are largely related to pathological conditions, and phenotypic effects are also linked to alterations in the expression profile following nuclear relocation of genes between functionally different compartments, generally occupying the periphery or the inner part of the cell nuclei. On the other hand, during evolution, chromosomal rearrangements may occur apparently without damaging phenotypic effects and are visible in currently phylogenetically related species. To increase our insight into chromosomal reorganisation in the cell nucleus, we analysed 18 chromosomal regions endowed with different genomic properties in cell lines derived from eight primate species covering the entire evolutionary tree. We show that homologous loci, in spite of their evolutionary relocation along the chromosomes, generally remain localised to the same functional compartment of the cell nuclei. We conclude that evolutionarily successful chromosomal rearrangements are those that leave the nuclear position of the regions involved unchanged. On the contrary, in pathological situations, the effect typically observed is on gene structure alteration or gene nuclear reposition. Moreover, our data indicate that new centromere formation could potentially occur everywhere in the chromosomes, but only those emerging in very GC-poor/gene-poor regions, generally located in the nuclear periphery, have a high probability of being retained through evolution. This suggests that, in the cell nucleus of related species, evolutionary chromosomal reshufflings or new centromere formation does not alter the functionality of the regions involved or the interactions between different loci, thus preserving the expression pattern of orthologous genes.

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.

Rapid evolution of the human mutation spectrum

DNA is a remarkably precise medium for copying and storing biological information. This high fidelity results from the action of hundreds of genes involved in replication, proofreading, and damage repair. Evolutionary theory suggests that in such a system, selection has limited ability to remove genetic variants that change mutation rates by small amounts or in specific sequence contexts. Consistent with this, using SNV variation as a proxy for mutational input, we report here that mutational spectra differ substantially among species, human continental groups and even some closely related populations. Close examination of one signal, an increased TCC→TTC mutation rate in Europeans, indicates a burst of mutations from about 15,000 to 2000 years ago, perhaps due to the appearance, drift, and ultimate elimination of a genetic modifier of mutation rate. Our results suggest that mutation rates can evolve markedly over short evolutionary timescales and suggest the possibility of mapping mutational modifiers.

DOI:http://dx.doi.org/10.7554/eLife.24284.001

DNA is a molecule that contains the information needed to build an organism. This information is stored as a code made up of four chemicals: adenine (A), guanine (G), cytosine (C), and thymine (T). Every time a cell divides and copies its DNA, it accidentally introduces ‘typos’ into the code, known as mutations. Most mutations are harmless, but some can cause damage. All cells have ways to proofread DNA, and the more resources are devoted to proofreading, the less mutations occur. Simple organisms such as bacteria use less energy to reduce mutations, because their genomes may tolerate more damage. More complex organisms, from yeast to humans, instead need to proofread their genomes more thoroughly.

Recent research has shown that humans have a lower mutation rate than chimpanzees and gorillas, their closest living relatives. Humans and other apes copy and proofread their DNA with basically the same biological machinery as yeast, which is about a billion years old. Yet, humans and apes have only existed for a small fraction of this time, a few million years. Why then do humans need to replicate and proofread their DNA differently from apes, and could it be that the way mutations arise is still evolving?

Previous research revealed that European people experience more mutations within certain DNA motifs (specifically, the DNA sequences ‘TCC’, ‘TCT’, ‘CCC’ and ‘ACC’) than Africans or East Asians do.

Now, Harris (who conducted the previous research) and Pritchard have compared how various human ethnic groups accumulate mutations and how these processes differ in different groups.

Statistical analysis of the genomes of thousands of people from all over the world did indeed show that the mutation rates of many different three-letter DNA motifs have changed during the past 20,000 years of human evolution. Harris and Pritchard report that when groups of humans left Africa and settled in isolated populations across different continents, each population quickly became better at avoiding mutations in some genomic contexts, but worse in others. This suggests that the risk of passing on harmful mutations to future generations is changing and evolving at an even faster rate than was originally suspected.

The results suggest that every human ethnic group carries specific variants of the genes which ensure that DNA replication and repair are accurate. These differences appear to influence which types of mutations are frequently passed down to future generations. An important next step will be to identify the genetic variants that could be controlling mutational patterns and how they affect human health.

DOI:http://dx.doi.org/10.7554/eLife.24284.002

How does chromatin package DNA within nucleus and regulate gene expression?

The human body is made up of 60 trillion cells, each cell containing 2 millions of genomic DNA in its nucleus. How is this genomic deoxyribonucleic acid [DNA] organised into nuclei? Around 1880, W. Flemming discovered a nuclear substance that was clearly visible on staining under primitive light microscopes and named it 'chromatin'; this is now thought to be the basic unit of genomic DNA organization. Since long before DNA was known to carry genetic information, chromatin has fascinated biologists. DNA has a negatively charged phosphate backbone that produces electrostatic repulsion between adjacent DNA regions, making it difficult for DNA to fold upon itself. In this article, we will try to shed light on how does chromatin package DNA within nucleus and regulate gene expression?

Large-scale heterochromatin remodeling linked to overreplication-associated DNA damage

Previously, we have shown that loss of the histone 3 lysine 27 (H3K27) monomethyltransferases ARABIDOPSIS TRITHORAX-RELATED 5 (ATXR5) and ATXR6 (ATXR6) results in the overreplication of heterochromatin. Here we show that the overreplication results in DNA damage and extensive chromocenter remodeling into unique structures we have named "overreplication-associated centers" (RACs). RACs have a highly ordered structure with an outer layer of condensed heterochromatin, an inner layer enriched in the histone variant H2AX, and a low-density core containing foci of phosphorylated H2AX (a marker of double-strand breaks) and the DNA-repair enzyme RAD51. atxr5,6 mutants are strongly affected by mutations in DNA repair, such as ATM and ATR. Because of its dense packaging and repetitive DNA sequence, heterochromatin is a challenging environment in which to repair DNA damage. Previous work in animals has shown that heterochromatic breaks are translocated out of the heterochromatic domain for repair. Our results show that atxr5,6 mutants use a variation on this strategy for repairing heterochromatic DNA damage. Rather than being moved to adjacent euchromatic regions, as in animals, heterochromatin undergoes large-scale remodeling to create a compartment with low chromatin density.