Search

BackgroundCancer somatic mutations are the product of multiple mutational and repair processes, which are tightly associated with DNA replication. Distinctive patterns of somatic mutations accumulation in tumors, termed mutational signatures, are indicative of processes the tumors underwent. While tumor mutational load is correlated with late replicating regions and spatial genome organization, much is unknown about the association of many different mutational processes and replication timing, and the interplay with chromatin structure remains an open question.MethodsWe systematically analyzed the mutational landscape of 2,787 WGS tumors from 32 different tumor types separately for early and late replicating regions. We used sequence context normalization and chromatin data to account for sequence and chromatin accessibility differences between early and late replicating regions. Moreover, we expanded the signature analyses to doublet base substitutions and small insertions and deletions by developing an artificial genomes-based approach to account for sequence differences between various genomic regions.ResultsWe revealed the replication timing (RT) association of single base, doublet base and small insertions and deletions mutational signatures. The association is signature specific: some are associated with early or late replication (such as UV-exposure signatures SBS7b and SBS7a, respectively) and others have no association. Most associations exist even after normalizing for genome accessibility. We further developed a focused mutational signature identification approach, which uses RT information to improve signature identification, and found that SBS16, which is biased towards early replication, is strongly associated with better survival rates in liver cancer.ConclusionsOur comprehensive analyses enabled a more robust classification of RT association of single base, doublet base and indels signatures. By doing so, we demonstrated a variation in the association with RT, as many mutational processes biased towards either early or late replication timing, and others have an equal RT distribution. These associations were independent from chromatin accessibility in most cases. This work highlights that restricting signatures analyses to concise genomic regions improves identification of signatures, such as SBS16, and demonstrates its clinically relevance as a predictor of improved survival of liver cancer patients.

Funder: National Institute of Health (NIH) DP2-GM123495, DNA replication follows a strict spatiotemporal program that intersects with chromatin structure but has a poorly understood genetic basis. To systematically identify genetic regulators of replication timing, we exploited inter-individual variation in human pluripotent stem cells from 349 individuals. We show that the human genome���s replication program is broadly encoded in DNA and identify 1,617 cis-acting replication timing quantitative trait loci (rtQTLs) ��� sequence determinants of replication initiation. rtQTLs function individually, or in combinations of proximal and distal regulators, and are enriched at sites of histone H3 trimethylation of lysines 4, 9, and 36 together with histone hyperacetylation. H3 trimethylation marks are individually repressive yet synergistically associate with early replication. We identify pluripotency-related transcription factors and boundary elements as positive and negative regulators of replication timing, respectively. Taken together, human replication timing is controlled by a multi-layered mechanism with dozens of effectors working combinatorially and following principles analogous to transcription regulation.

Cellular proliferation depends on the accurate and timely replication of the genome. Several genetic diseases are caused by mutations in key DNA replication genes; however, it remains unclear whether these genes influence the normal program of DNA replication timing. Similarly, the factors that regulate DNA replication dynamics are poorly understood. To systematically identify trans-acting modulators of replication timing, we profiled replication in 184 cell lines from three cell types, encompassing 60 different gene knockouts or genetic diseases. Through a rigorous approach that considers the background variability of replication timing, we concluded that most samples displayed normal replication timing. However, mutations in two genes showed consistently abnormal replication timing. The first gene was RIF1, a known modulator of replication timing. The second was MCM10, a highly conserved member of the pre-replication complex. MCM10 mutant cells demonstrated replication timing variability comprising 46% of the genome and at different locations than RIF1 knockouts. Replication timing alterations in MCM10-mutant cells was predominantly comprised of replication initiation defects. Taken together, this study demonstrates the remarkable robustness of the human replication timing program and reveals MCM10 as a novel modulator of DNA replication timing.

Motivation Genomic DNA replicates according to a reproducible spatiotemporal program, with some loci replicating early in S phase while others replicate late. Despite being a central cellular process, DNA replication timing studies have been limited in scale due to technical challenges. Results We present TIGER (Timing Inferred from Genome Replication), a computational approach for extracting DNA replication timing information from whole genome sequence data obtained from proliferating cell samples. The presence of replicating cells in a biological specimen leads to non-uniform representation of genomic DNA that depends on the timing of replication of different genomic loci. Replication dynamics can hence be observed in genome sequence data by analyzing DNA copy number along chromosomes while accounting for other sources of sequence coverage variation. TIGER is applicable to any species with a contiguous genome assembly and rivals the quality of experimental measurements of DNA replication timing. It provides a straightforward approach for measuring replication timing and can readily be applied at scale. Availability and implementation TIGER is available at https://github.com/TheKorenLab/TIGER. Supplementary information Supplementary data are available at Bioinformatics online.

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.

Sensitive detection of somatic copy number alterations (SCNA) in cancer genomes is confounded by “waviness” in read depth data. We present dryclean, a signal processing algorithm to optimize SCNA detection in whole genome (WGS) and targeted sequencing platforms through foreground detection and background subtraction of read depth data. Application of dryclean to WGS demonstrates that WGS waviness is driven by replication timing. Re-analysis of thousands of tumor profiles reveals that dryclean provides superior detection of biologically relevant SCNAs relative to state-of-the-art algorithms. Applied to in silico tumor dilutions, dryclean improves the sensitivity of relapse detection 10-fold relative to current standards. dryclean is available as an R package in the GitHub repository https://github.com/mskilab/dryclean

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 analogue 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.

The DNA replication timing program is modulated throughout development and is also one of the main factors influencing the distribution of mutation rates across the genome. However, the relationship between the mutagenic influence of replication timing and its developmental plasticity remains unexplored. Here, we studied the distribution of copy number variations (CNVs) and single nucleotide polymorphisms across the zebrafish genome in relation to changes in DNA replication timing during embryonic development in this model vertebrate species. We show that CNV sites exhibit strong replication timing plasticity during development, replicating significantly early during early development but significantly late during more advanced developmental stages. Reciprocally, genomic regions that changed their replication timing during development contained a higher proportion of CNVs than developmentally constant regions. Developmentally plastic CNV sites, in particular those that become delayed in their replication timing, were enriched for the clustered protocadherins, a set of genes important for neuronal development that have undergone extensive genetic and epigenetic diversification during zebrafish evolution. In contrast, single nucleotide polymorphism sites replicated consistently early throughout embryonic development, highlighting a unique aspect of the zebrafish genome. Our results uncover a hitherto unrecognized interface between development and evolution.

Human B-lineage precursor acute lymphoid leukemias (BCP-ALLs) comprise a group of genetically and clinically distinct disease entities with features of differentiation arrest at known stages of normal B-lineage differentiation. We previously showed BCP-ALL cells display unique and clonally heritable DNA-replication timing (RT) programs; i.e., programs describing the variable order of replication of megabase-scale chromosomal units of DNA in different cell types. To determine the extent to which BCP-ALL RT programs mirror or deviate from specific stages of normal human B-cell differentiation, we transplanted immunodeficient mice with quiescent normal human CD34+ cord blood cells and obtained RT signatures of the regenerating B-lineage populations. We then compared these with RT signatures for leukemic cells from a large cohort of BCP-ALL patients. The results identify BCP-ALL subtype-specific features that resemble specific stages of B-cell differentiation and features that appear associated with relapse. These results suggest the genesis of BCP-ALL involves alterations in RT that reflect clinically relevant leukemia-specific genetic and/or epigenetic changes.SUMMARYGenome-wide DNA replication timing profiles of >100 pediatric leukemic samples and normally differentiating human B-lineage cells isolated from xenografted immunodeficient mice were generated. Comparison of these identified potentially clinically relevant features that both match and deviate from the normal profiles.

Chromatin structure has a major influence on the cell-specific density of somatic mutations along the cancer genome. Here, we present a pan-cancer study in which we searched for the putative cancer cell-of-origin of 2,550 whole genomes, representing 32 cancer types by matching their mutational landscape to the regional patterns of chromatin modifications ascertained in 104 normal tissue types. We found that, in almost all cancer types, the cell-of-origin can be predicted solely from their DNA sequences. Our analysis validated the hypothesis that high-grade serous ovarian cancer originates in the fallopian tube and identified distinct origins of breast cancer subtypes. We also demonstrated that the technique is equally capable of identifying the cell-of-origin for a series of 2,044 metastatic samples from 22 of the tumor types available as primaries. Moreover, cancer drivers, whether inherited or acquired, reside in active chromatin regions in the respective cell-of-origin. Taken together, our findings highlight that many somatic mutations accumulate while the chromatin structure of the cell-of-origin is maintained and that this historical record, captured in the DNA, can be used to identify the often elusive cancer cell-of-origin.

In dividing cells, DNA replication occurs in a precise order, but many questions remain regarding the mechanisms of replication timing establishment and regulation. We now have generated genome-wide, high-resolution replication timing maps throughout zebrafish development. Unexpectedly, in the rapid cell cycles preceding the midblastula transition, a defined timing program was present that predicted the initial wave of zygotic transcription. Replication timing was thereafter progressively and continuously remodeled across the majority of the genome, and epigenetic changes involved in enhancer activation frequently paralleled developmental changes in replication timing. The long arm of Chromosome 4 underwent a dramatic developmentally regulated switch to late replication during gastrulation, reminiscent of mammalian X Chromosome inactivation. This study reveals that replication timing is dynamic and tightly linked to epigenetic and transcriptional changes throughout early zebrafish development. These data provide insight into the regulation and functions of replication timing and will enable further mechanistic studies.

SummaryGenomic DNA replicates in a choreographed temporal order that impacts the distribution of mutations along the genome. We show here that DNA replication timing is shaped by genetic polymorphisms that act in cis upon megabase-scale DNA segments. In genome sequences from proliferating cells, read depth along chromosomes reflected DNA replication activity in those cells. We used this relationship to analyze variation in replication timing among 161 individuals sequenced by the 1000 Genomes Project. Genome-wide association of replication timing with genetic variation identified 16 loci at which inherited alleles associate with replication timing. We call these “replication timing quantitative trait loci” (rtQTLs). rtQTLs involved the differential use of replication origins, exhibited allele-specific effects on replication timing, and associated with gene expression variation at megabase scales. Our results show replication timing to be shaped by genetic polymorphism and identify a means by which inherited polymorphism regulates the mutability of nearby sequences.

Major international projects are underway that are aimed at creating a comprehensive catalogue of all the genes responsible for the initiation and progression of cancer. These studies involve the sequencing of matched tumour-normal samples followed by mathematical analysis to identify those genes in which mutations occur more frequently than expected by random chance. Here we describe a fundamental problem with cancer genome studies: as the sample size increases, the list of putatively significant genes produced by current analytical methods burgeons into the hundreds. The list includes many implausible genes (such as those encoding olfactory receptors and the muscle protein titin), suggesting extensive false-positive findings that overshadow true driver events. We show that this problem stems largely from mutational heterogeneity and provide a novel analytical methodology, MutSigCV, for resolving the problem. We apply MutSigCV to exome sequences from 3,083 tumour-normal pairs and discover extraordinary variation in mutation frequency and spectrum within cancer types, which sheds light on mutational processes and disease aetiology, and in mutation frequency across the genome, which is strongly correlated with DNA replication timing and also with transcriptional activity. By incorporating mutational heterogeneity into the analyses, MutSigCV is able to eliminate most of the apparent artefactual findings and enable the identification of genes truly associated with cancer.

A new analysis of cancer genomes identifies a decrease in the mutation burden of exons, but not introns, as compared to expectation. This difference can be explained by preferential recruitment of the DNA mismatch repair machinery to a protein modification that marks exons.

SummaryDe novo mutation plays an important role in autism spectrum disorders (ASDs). Notably, pathogenic copy number variants (CNVs) are characterized by high mutation rates. We hypothesize that hypermutability is a property of ASD genes and may also include nucleotide-substitution hot spots. We investigated global patterns of germline mutation by whole-genome sequencing of monozygotic twins concordant for ASD and their parents. Mutation rates varied widely throughout the genome (by 100-fold) and could be explained by intrinsic characteristics of DNA sequence and chromatin structure. Dense clusters of mutations within individual genomes were attributable to compound mutation or gene conversion. Hypermutability was a characteristic of genes involved in ASD and other diseases. In addition, genes impacted by mutations in this study were associated with ASD in independent exome-sequencing data sets. Our findings suggest that regional hypermutation is a significant factor shaping patterns of genetic variation and disease risk in humans.PaperFlick

Many overlapping surveillance and repair mechanisms operate in eukaryotic cells to ensure the stability of the genome. We have screened to isolate yeast mutants exhibiting increased levels of recombination between repeated sequences. Here we characterize one of these mutants, elg1 . Strains lacking Elg1p exhibit elevated levels of recombination between homologous and nonhomologous chromosomes, as well as between sister chromatids and direct repeats. These strains also exhibit increased levels of chromosome loss. The Elg1 protein shares sequence homology with the large subunit of the clamp loader replication factor C (RFC) and with the product of two additional genes involved in checkpoint functions and genome maintenance: RAD24 and CTF18 . Elg1p forms a complex with the Rfc2–5 subunits of RFC that is distinct from the previously described RFC-like complexes containing Rad24 and Ctf18. Genetic data indicate that the Elg1, Ctf18, and Rad24 RFC-like complexes work in three separate pathways important for maintaining the integrity of the genome and for coping with various genomic stresses.