Data from Cancer Mutational Processes Vary in Their Association with Replication Timing and Chromatin Accessibility

Article Figures and data Abstract Editor's evaluation eLife digest Introduction Results Discussion Materials and methods Data availability References Decision letter Author response Article and author information Metrics Abstract Recent studies have suggested that the human germline mutation rate and spectrum evolve rapidly. Variation in generation time has been linked to these changes, though its contribution remains unclear. We develop a framework to characterize temporal changes in polymorphisms within and between populations, while controlling for the effects of natural selection and biased gene conversion. Application to the 1000 Genomes Project dataset reveals multiple independent changes that arose after the split of continental groups, including a previously reported, transient elevation in TCC>TTC mutations in Europeans and novel signals of divergence in C>Gand T>A mutation rates among population samples. We also find a significant difference between groups sampled in and outside of Africa in old T>C polymorphisms that predate the out-of-Africa migration. This surprising signal is driven by TpG>CpG mutations and stems in part from mis-polarized CpG transitions, which are more likely to undergo recurrent mutations. Finally, by relating the mutation spectrum of polymorphisms to parental age effects on de novo mutations, we show that plausible changes in the generation time cannot explain the patterns observed for different mutation types jointly. Thus, other factors – genetic modifiers or environmental exposures – must have had a non-negligible impact on the human mutation landscape. Editor's evaluation This important study investigates temporal variation in patterns of germline mutation during the evolution of human populations. Using a compelling approach that controls for the effects of selection and biased gene conversion the authors show that changes in generation time alone cannot explain the joint patterns observed for different mutation types, suggesting that other factors such as genetic modifiers or environmental exposures must have played a role as well. This work will be of broad interest to population geneticists and evolutionary biologists. https://doi.org/10.7554/eLife.81188.sa0 Decision letter Reviews on Sciety eLife's review process eLife digest Each human has 23 pairs of chromosomes, one set inherited from each parent. But the child's chromosomes are not an exact copy of their parents' chromosomes. Spontaneous changes or mutations in the DNA during the formation of the egg or sperm cells, or early development of the embryo, can change a small fraction of the nucleotides or 'letters' that make up the DNA. These modifications are an important source of genetic diversity in human populations and contribute to the evolution of new traits. Each genetic variant in present-day human populations represents a mutation in one of their ancestors. The types and frequencies of variants vary across human populations and have changed over time, suggesting that mutation patterns have evolved in the past. But the processes driving these population-level differences remain elusive. One possible factor may be changes in the average age of reproduction or the generation time in a population . For example, older parents contribute more – and also different types of – new mutations to their children than younger parents do. Populations, where it is customary to have children at older ages, may therefore have a different mutation landscape. To find out if this is indeed the case, Gao et al. used computer algorithms to analyze the genomes of hundreds of people living on three continents who participated in 'the 1,000 Genomes Project'. The analysis identified differences in mutation patterns across continental groups and estimated when these changes occurred. Further, they showed that although the age of reproduction had an impact on the mutation landscape, differences in generation time alone could not explain the observed changes in the human mutation spectrum. Factors other than generation time, such as environmental exposures, may have played a role in shifting these patterns. The study provides new insights into the changes in the mutation landscape over the course of human evolution. Mapping these patterns in humans worldwide may help scientists understand the causes underlying these changes. The techniques used by Gao et al. may also help analyze changes in mutation patterns in other organisms. Introduction Recent advances in high-throughput sequencing have enabled large-scale surveys of genetic variation in thousands of humans, providing a rich resource for understanding the source and mechanisms shaping the mutation landscape over time. Comparisons of polymorphism patterns across geographic population samples have uncovered numerous differences in the mutation rates and spectra (i.e., relative proportions of different types of mutations) (DeWitt et al., 2021; Goldberg and Harris, 2022; Harris, 2015; Harris and Pritchard, 2017; Hwang and Green, 2004; Mathieson and Reich, 2017; Moorjani et al., 2016a; Narasimhan et al., 2017; Speidel et al., 2019). A notable signal in humans is the enrichment of TCC>TTC variants in polymorphism data from Europeans relative to Africans and Asians (Harris and Pritchard, 2017). This signal is also observed in South Asians to a lesser degree and has been suggested to originate in ancient Neolithic farmers (Harris and Pritchard, 2017; Speidel et al., 2021). Many other subtle but statistically significant signals have also been detected; given the recent common ancestry of human populations, this finding indicates that the mutational spectrum in humans has been evolving rapidly. Several genetic and nongenetic factors have been implicated as affecting mutation rates and acting as potential drivers of observed interpopulation differences in the mutation spectrum of polymorphisms. First, some environmental exposures can increase mutation rates, especially of particular types. As humans in different geographic locations and environments may have experienced differential exposures over the past 50,000–100,000 years since the out-of-Africa (OOA) migration, rates of specific mutation types could have diverged between populations (Harris, 2015; Mathieson and Reich, 2017). Second, genetic modifiers of mutation rates, such as variants in genes that copy or repair DNA, could segregate at different frequencies across populations. Despite the deleterious effects of alleles that modify mutation rates, in recombining species, they could be nearly neutral and maintained for a long time, leading to genome-wide differences across populations (Milligan et al., 2022; Seoighe and Scally, 2017). In addition, direct sequencing of human pedigrees has revealed the effects of the parental ages at reproduction on the relative fractions of mutation types (Goldmann et al., 2018; Jónsson et al., 2017). For example, as parents age, fathers pass on disproportionally more T>C mutations, and mothers contribute a higher fraction of C>G mutations (Jónsson et al., 2017). Thus, differences in the average reproductive ages, or equivalently 'generation times,' alone can lead to differences in mutation spectrum across populations; indeed, such differences have been invoked to explain a large fraction of observed variation in types of polymorphisms among population samples (Coll Macià et al., 2021). The joint distribution of mutation type and frequency of polymorphisms, however, depends not only on the mutational input, but also on other evolutionary forces such as natural selection, biased gene conversion, and demography. In particular, natural selection distorts the allele frequency distribution and fixation probability of non-neutral variants, and the average effect of natural selection can differ across mutation types (Wakeley, 2010). As an example, genic regions tend to be more GC-rich, so mutations at G:C base pairs may be subject to stronger purifying or background selection compared to mutations at A:T base pairs (Lander et al., 2001; McVicker et al., 2009). GC-biased gene conversion (gBGC) is another process that exerts differential effects across mutation types by effectively acting like positive selection favoring mutations from weak alleles (A or T) to strong alleles (C or G) and negative selection against mutations from strong to weak alleles (Duret and Galtier, 2009). The strengths of selection and gBGC depend on the effective population size and thus on the demographic history of a population. Demographic history also influences allele frequencies for a given allele age (Kimura, 1969). This poses a challenge in interpreting previous studies (Harris and Pritchard, 2017; Mathieson and Reich, 2017) aimed at learning about when changes in mutational processes may have occurred by using allele frequencies, as mutations of the same frequency can have drastically different distributions of ages in distinct populations (e.g., doubletons in Africans are substantially older than doubletons in Europeans or Asians; Mathieson and McVean, 2014). Beyond the biological processes that shape polymorphism data, the characterization of the mutational spectrum can be biased by many technical issues. For instance, a recent study showed that some interpopulation differences discovered in low-coverage 1000 Genomes data may be driven by cell line artifacts or errors in PCR amplification (Anderson-Trocmé et al., 2020). Further, comparisons of mutation patterns across datasets are sensitive to differences in the accessible genomic regions across studies. Because there is large variation in mutation rates and base pair composition across genomic regions, differences in the regions sequenced across studies can have a non-negligible impact on comparisons of mutation spectrum across datasets (Monroe et al., 2022; Seplyarskiy et al., 2021). In addition, the number of genomes surveyed, in combination with the specific population demographic history, influences the chance of observing repeated mutations at the same site, and thus the observed polymorphism patterns (Lek et al., 2016). Given these challenges, it remains unclear whether the numerous observed differences across human populations stem from rapid evolution of the mutation process itself, other evolutionary processes, or technical factors. Motivated by these considerations, we propose a new framework to compare the mutation spectrum over time and across human populations. First, we infer the age of each derived allele observed in a population using a newly developed approach, Relate, which reconstructs local genealogies and estimates allele ages (Speidel et al., 2019). This approach allows us to perform more reliable comparisons across populations as well as to investigate changes in mutation processes across time. Next, we minimize confounding effects of selection by removing constrained regions and known targets of selection in the genome. We also control the effects of biased gene conversion by focusing on comparison of pairs of mutations (e.g., T>C and T>G) that are subject to similar effects of gBGC. This pairwise comparison further mitigates the issue of interdependencies in comparing mutation fractions (i.e., an increased contribution of one mutation type necessarily lowers the contribution of other mutation types). Based on this new framework, we re-evaluate the evidence for evolution of the mutation spectrum in human populations and investigate when, how, and in which populations significant changes have occurred over the course of human evolution. Finally, by relating parental age effects on the mutation spectrum estimated in contemporary pedigrees to the observed patterns of polymorphisms of varying ages, we evaluate the role of changes in generation times in shaping the human mutation landscape. Results Variation in the spectrum of human polymorphisms over time We analyzed single-nucleotide polymorphisms (SNP) identified in high-coverage whole-genome sequencing data from the 1000 Genomes Project, including 178 individuals of West African ancestry living in Ibadan, Nigeria (YRI), 179 individuals of Northern European ancestry living in the United States (CEU), and 103 individuals of East Asian ancestry living in Beijing, China (CHB) (Byrska-Bishop et al., 2022). To focus on putatively neutral mutations, we removed exons and phylogenetically conserved regions as previous studies (Harris and Pritchard, 2017; Moorjani et al., 2016a). To perform reliable comparison between datasets in downstream analysis and ensure the results are not driven by local genomic differences in mutation rate, we focused on regions that were accessible in both population and pedigree datasets ( hereafter, referred to as 'commonly accessible regions') ('Materials and Methods'). We inferred the age of each derived variant (with the ancestral allele determined based on the six primate EPO (Enredo, Pecan, Ortheus) alignment) in YRI, CEU, and CHB using Relate, a method to reconstruct local genealogies based on phased haplotype sequences (Speidel et al., 2019). We then divided all SNPs into 15 bins based on the ages of the derived allele inferred by Relate, accounting for uncertainty in the inferred mutation age by assuming a uniform distribution of ages between the inferred lower and upper bounds for each variant ('Materials and methods'). We classified each SNP into six disjoint classes based on the type of base pair substitution: T>A, T>C, T>G, C>A, C>G, and C>T (each including the corresponding substitution on the reverse complement strand, e.g., T>C includes both T>C and A>G substitutions). Given the well-characterized hypermutability of methylated CpG sites (Duncan and Miller, 1980; Kong et al., 2012), we further divided C>T SNPs into subtypes occurring in CpG and non-CpG contexts by considering the flanking base pair on either side of the variant. We find marked differences in the relative proportions of different mutation types (i.e., the mutation spectrum) across varying allele age bins within CEU (Figure 1) as well as in YRI and in CHB (Figure 1—figure supplement 1), as seen earlier in the low-coverage 1000G data (Speidel et al., 2019). We obtain qualitatively similar results when considering other 1000G populations of TSI, LWK, and JPT (Figure 1—figure supplement 1). This observation echoes previous findings about the evolution of the mutation spectrum comparing polymorphisms across allele frequencies (Carlson et al., 2018; Harris and Pritchard, 2017; Mathieson and Reich, 2017). As noted previously, however, differences in mutation spectrum across frequencies alone are weak evidence for the evolution of the mutation process itself because patterns of standing polymorphisms can be affected by repeat mutations and other evolutionary forces, including selection and gene conversion. Figure 1 with 8 supplements see all Download asset Open asset Changes in the mutation spectrum of polymorphisms in CEU over evolutionary time. Figure 1—source data 1 Bedfile for the commonly accessible region excluding exons and phylogenetically conserved elements. https://cdn.elifesciences.org/articles/81188/elife-81188-fig1-data1-v3.zip Download elife-81188-fig1-data1-v3.zip Figure 1—source data 2 Text files with (pseudo-)counts of different types of mutations in YRI, LWK, CEU, TSI, CHB, JPT in each time window. https://cdn.elifesciences.org/articles/81188/elife-81188-fig1-data2-v3.zip Download elife-81188-fig1-data2-v3.zip Notably, the infinite sites model is a reasonable assumption for small sample sizes (Kimura, 1969), but recurrent mutations become highly likely in large datasets, especially at sites with higher mutation rates (Harpak et al., 2016; Lek et al., 2016). Recurrent, multi-allelic, and back mutations violate the model assumptions of Relate and are often excluded from its output. For instance, given the higher mutation rate of transitions at CpG sites, such SNPs are more likely to be subject to recurrent mutations in a large sample and thus may map to multiple branches in the tree, leading to their exclusion from Relate's output (Speidel et al., 2019). As expected from these considerations, the fraction of CpG C>T SNPs in young mutations (i.e., those estimated to have occurred in the past ~50 generations) is lower than proportions in de novo mutations (DNMs) in present-day pedigree studies (Figure 1—figure supplement 2). Differences in mutation spectrum across age bins in modern humans persist even after excluding CpG C>T mutations (Figure 1—figure supplement 3), however, indicating that other mutation types are also changing in relative frequency over time and the observed patterns are not driven solely by recurrent mutation at CpG sites. Next, we examined the effect of linked selection on different mutation types. While we excluded direct targets of selection from analysis (i.e., exons and conserved regions), much of the genome is linked to non-neutral variants and subject, notably, to background selection (Charlesworth et al., 1993; McVicker et al., 2009; Murphy et al., 2022). A common measure of the effects of background selection is the B-statistic or B-score that estimates the reduction in nucleotide diversity levels compared to the neutral expectation (McVicker et al., 2009). To characterize the impact of linked selection, we calculated the average genome-wide B-score of each mutation type. We find nearly identical average B-scores and similar distributions for all mutation types (Figure 1—figure supplement 4). Further, comparing the mutation spectrum over time in CEU, YRI, and CHB, we obtain qualitatively similar results when restricting to regions with weak background selection (B-score > 800, where the genetic diversity is reduced by <20% compared to the neutral expectation; Figure 1—figure supplement 5). These analyses suggest that although linked selection has pervasive effects, its average impact is relatively uniform across the seven mutation types in commonly accessible regions ('Materials and methods). Gene conversion is another evolutionary process that can have a profound impact on the mutation spectrum of polymorphisms. gBGC acts like selection for certain mutation types by causing the preferential transmission of strong (S) alleles (C or G) over weak (W) alleles (A or T) in heterozygotes (Duret and Galtier, 2009). Accordingly, we observe enrichments of W>S mutations (T>C and T>G) in common variants and of S>W mutations (C>A and C>T) in rare variants (Figure 1—figure supplement 6A). Moreover, gBGC violates model assumptions of Relate (for both neutrality and infinite-sites mutation model) and could lead to subtle biases in estimated allele ages (Speidel et al., 2019). Due to the effect of gBGC, W>S mutations are expected to be enriched in older variants compared to S>W variants, and this enrichment is expected to be stronger in regions with high recombination rates (Glémin et al., 2015). Indeed, we observe such enrichment and the expected correlation with recombination rate (Figure 1—figure supplement 7A), supporting the effect of gBGC on the mutation spectrum of variants across mutation ages. Furthermore, the effect of gBGC is expected to vary across populations as its strength depends on the effective population size. Accordingly, we observe that the trends of the ratio of W>S to S>W over time differ across human populations (Figure 1—figure supplement 6B, Figure 1—figure supplement 7B). These results highlight the need to account for gBGC in order to reliably interpret the source of observed differences within and between populations (whether using allele frequency bins or allele age estimates). Pairwise comparisons of mutation types accounting for gBGC In light of the impact of gBGC on the mutation spectrum, we focused on comparisons of pairs of mutation types subject to similar effects of gBGC (i.e., in which both are favored, disfavored, or unaffected by gBGC). Specifically, we focused on four pairwise comparisons including (1) C>T at non-CpGs vs. C>A at non-CpGs; (2) C>T CpGs vs. C>A CpGs; (3) C>G vs. T>A; and (4) T>C vs. T>G. In principle, it is possible that the strength of gBGC is distinct for different types of variantsinvolving S and W alleles (Tsai-Wu et al., 1992). However, in mice, roughly similar conversion rates are observed for C>A and C>T non-crossover gene conversion events as well as for T>C and T>G events (Li et al., 2019), lending support to using pairwise mutation ratios for controlling the effects of gBGC at least to a first approximation. Three of the four pairwise comparisons involve mutation types with the same mutational opportunity (e.g., both T>C and T>G mutations involve changes at ancestral T bases in the genome), which further minimizes the confounding effects of regional variation on the chance of recurrent mutation or strength of background selection. Moreover, the pairwise ratios impose no co-dependency among mutation types as the four comparisons are mathematically independent of each other (although they may be biologically dependent if multiple ratios are affected simultaneously by some change in the mutational process). Investigating the mutation spectrum using these four pairwise comparisons, we observe marked differences in the ratios both over evolutionary time and across populations. Specifically, we find multiple independent signals of mutation rate evolution, reflected by both temporal variation within a population and differences between YRI, CHB, and CEU (p<0.01 by chi-square test after correcting for multiple hypothesis testing; 'Materials and methods'; Figure 2A). These differences may represent geographic or population differences as we these findings in other population samples from the same continents – LWK, TSI, and JPT – from the 1000 Genomes Project (Figure supplement 1). Figure 2 with supplements see all Download asset Open asset of pairwise mutation ratios for polymorphisms in different time pairwise mutation ratios are each of which mutation types that are for mutational opportunity and effects of GC-biased gene conversion The indicates the with the out-of-Africa (OOA) migration. The represent the observed polymorphism while the the assuming a distribution of polymorphism in are three ratios that show significant interpopulation with into each in lower to the from a chi-square test after a for and of levels were used in Figure in ratio in CEU at non-CpG sites, after excluding the four contexts and previously identified to be with the in Europeans by Harris and Pritchard, as well as contexts affected by of in mutational of and (Harris, 2015; Mathieson and Reich, 2017). divergence in ratio among three population ratio in YRI than CEU and CHB samples among old variants, driven by Figure data 1 Text files with (pseudo-)counts of mutations classified into types in genomic regions and within the C>G mutation in YRI, LWK, CEU, TSI, CHB, and JPT in each time window. Download We multiple to out technical artifacts or to the mutation process in to the observed interpopulation We obtain qualitatively similar results when restricting the analysis to putatively neutral regions with B-score > (Figure supplement or comparing regions with high and recombination rates (Figure supplement 3), that the of pairwise comparisons effectively controls for the effects of selection and gBGC. In to account for potential in mutation ages estimated by Relate, we variants by allele frequencies of inferred mutation ages and the signals in mutation age analysis (Figure supplement 4). We observe similar results with different for allele age, based on inferred mutation ages in YRI or CHB (Figure supplement 5). these results strong evidence that the human germline mutation spectrum has evolved over time and across populations. we the and population of the mutation rate changes to each of the signals we in of non-CpG ratio in Europeans The signal that we observed is the transient elevation in the ratio of mutations at non-CpG sites in CEU compared to the ancestral the in the non-CpG ratios of CHB and YRI not a similar at recent This signal the previously enrichment of C>T polymorphisms in a in as well as other contexts (Harris and Pritchard, 2017; Mathieson and Reich, 2017; Speidel et al., 2019). Investigating the temporal patterns in CEU, we find that the increase in the ratio of mutations at non-CpG sites from the time the years or and and in the recent age of (Figure 2A). Because there is large uncertainty in inferred allele ages and approach often the contribution of each variant over or more age bins ('Materials and the and of variation be the transient change in non-CpG C>T mutations likely and of higher than results However, the temporal and geographic enrichment patterns from analysis are with previous based on low-coverage or other datasets (Harris and Pritchard, 2017; Mathieson and Reich, 2017; Speidel et al., 2019). all non-CpG the interpopulation differences are in the four previously contexts and Harris and Pritchard, but are in other non-CpG contexts as well (Figure analysis that these mutational contexts are enriched in of the mutational from mutations in the of in and with exposures to light and et al., Harris, 2015; Mathieson and Reich, 2017). To test whether one of these mutational may be for the observed differences in polymorphism data, we the mutation ratio at non-CpG sites after excluding the contexts affected by or ('Materials and methods'). While we observe some reduction in the of non-CpG ratio in the interpopulation differences remain significant (Figure These results suggest the transient change in non-CpG ratio is not driven by the mutational mechanisms corresponding to either or Thus, the of this signal in Europeans remains unclear. of ratio among populations The interpopulation difference is in the ratio (Figure the among all three populations. In the past both YRI and CEU samples show an increase in the ratio of different while in the CHB, the ratio and then relatively for roughly (Figure 2A). the previous interpopulation differences in remain highly significant for the recent variants as well to factors the relative rates of C>G and T>A mutations at The fraction of C>G in de novo germline mutations is sensitive to parental ages, with the age at (Jónsson et al., 2017). This the that the interpopulation differences in ratio are driven by different average reproductive ages among populations (Coll Macià et al., 2021). To test this we the regional enrichment of C>G mutations – enriched – as of the genome with the C>G SNP that to of the age effect (i.e., the increase in with Jónsson et al., 2017). The ratio within the C>G enriched regions not show significant interpopulation differences (Figure reduced to the much lower SNP in these regions of all but see Figure supplement for of the C>G enriched regions, the three populations differ as much as they genome-wide (Figure indicating that the differential of C>G mutations with ages is not the of the differences

There are two frameworks for characterizing mutational signatures which are commonly used to describe the nucleotide patterns that arise from mutational processes. Estimated mutational signatures from fitting these two methods in human cancer can be found online, in the Catalogue Of Somatic Mutations In Cancer (COSMIC) website or a GitHub repository. The two frameworks make differing assumptions regarding independence of base pairs and for that reason may produce different results. Consequently, there is a need to compare and contrast the results of the two methods, but no such tool currently exists. In this paper, we provide a simple and intuitive interface that allows comparisons of pairs of mutational signatures to be easily performed. Cosine similarity measures the extent of signature similarity. To compare mutational signatures of different formats, one signature type (COSMIC or pmsignature) is converted to the format of the other before the signatures are compared. iMutSig provides a simple and user-friendly web application allowing researchers to download published mutational signatures of either type and to compare signatures from COSMIC to those from pmsignature, and vice versa. Furthermore, iMutSig allows users to input a self-defined mutational signature and examine its similarity to published signatures from both data sources. iMutSig is accessible online and source code is available for download from GitHub.

Introduction: Identifying the causes of early-onset colorectal cancer (EOCRC) is integral to reducing its rising global burden. Approximately 25% of EOCRC cases have a first degree relative (FDR) who has also developed CRC. Determining if familial EOCRC is the result of shared hereditary or environmental factors or chance familial aggregation has important implications for clinical management and cancer prevention. We investigated the similarity of tumor mutational and experimental signature profiles from EOCRC cases and their EOCRC-affected FDRs to identify evidence of common hereditary or environmental etiologies. Methods: We identified 13 people with a mismatch repair proficient (MMRp) EOCRC who had at least one FDR with a MMRp EOCRC with tumor tissue available (familial EOCRC pairs). To calibrate mutational signature profile “similarity”, a reference group of CRCs from FDR pairs with known hereditary cancer syndromes were included (hereditary CRCs: Lynch syndrome (n=4), biallelic MUTYH (n=1), biallelic NTHL1 (n=1), BRCA1 (n=1)). In addition, a group of 254 MMRp EOCRC cases were analyzed to assess the prevalence of any familial signatures. All participants were recruited through the ANGELS and CCFR studies. COSMIC v3.4 tumour mutational signatures (TMS) and COSMIC v1.0 experimental mutational signatures (EMS) were calculated from whole exome sequencing. Spearman’s rank correlation coefficient was calculated across tumours against which familial EOCRCs and hereditary CRCs were assessed for similarity. Hierarchical clustering with bootstrap resampling (pvclust, R) was applied to correlation-based distances, supporting clusters with an approximately unbiased (AU) value ≥95% as significant. Results: Hereditary CRCs had a Spearman’s correlation range of 0.34-0.84, derived from TMS and EMS. Of the 13 familial EOCRC pairs, five (38%) shared a TMS/EMS profile (Spearman’s correlation range 0.37-0.94). One familial EOCRC pair exhibited a dominant EMS related to n-nitrosopyrrolidine (NPYR) exposure (correlation=0.94). This familial EOCRC pair, a mother and daughter, were diagnosed within a year of each other at ages 49 and 18 years, respectively. Hierarchical bootstrap clustering identified 11 EOCRCs (11/254 or 4%) with a dominant EMS profile related to NPYR exposure, two of which had a CRC-affected FDR. Conclusions: TMS and EMS profiling of unexplained familial EOCRC-affected FDR pairs showed heterogeneity with regards to shared mutational processes among the FDR pairs. Only 38% of familial EOCRC pairs demonstrated evidence of similar mutational processes in their EOCRCs that were suggestive of a shared etiology. One experimental signature related to NPYR exposure, a carcinogen found in food, tobacco smoke and occupational settings, was present as a dominant mutational process in 4% of the EOCRCs tested, highlighting the potential utility of mutational signature analysis for elucidating the etiology of EOCRC. Citation Format: Alysha Prisc, Peter Georgeson, Jihoon E. Joo, Romy Walker, Mark Clendenning, Julia Como, Natalie Diepenhorst, Julie McDonald, Steven Gallinger, Robert C. Grant, Dylan E. O'Sullivan, Darren R. Brenner, Finlay A. Macrae, Christophe Rosty, Mark A. Jenkins, Ingrid M. Winship, Daniel D. Buchanan. Unexplained familial early-onset colorectal cancer: Determining a shared etiology by tumor mutational signature profiling [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: The Rise in Early-Onset Cancers—Knowledge Gaps and Research Opportunities; 2025 Dec 10-13; Montreal, QC, Canada. Philadelphia (PA): AACR; Clin Cancer Res 2025;31(23_Suppl):Abstract nr C014.

I come from a family in which there have been no scientists or doctors. I was interested, however, in biology at school and started my scientific career by training in medicine at Oxford University and Guys Hospital, London. Practising as a doctor reinforced my curiosity about the biological processes underlying human disease. As a consequence, I pursued a clinical vocation in histopathology, a discipline that couples exposure to the sights and smells of the autopsy room with a daily journey into the often beautiful, sometimes ugly world of healthy and diseased human tissues under the microscope. After an introduction to general histopathology in Nick Wright's department at the Hammersmith Hospital, London, I completed my postgraduate medical training in neuropathology with Peter Lantos at the Maudsley Hospital, London.Peering at the nuclei of cancer cells under the microscope, for me it was a matter of fascination that hidden within them were the key events converting normal cells into cancer cells, and frustration because they were out of reach. Many of the tissue samples examined by pathologists are from cancers. The clonal theory of cancer development and the general role of DNA mutations in generating cancer cell clones had been established by 1986 when I was working as a junior pathologist. Indeed, the first mutated cancer gene, HRAS , had recently been identified through application of the, then new, technologies of recombinant DNA technology. Peering at the nuclei of cancer cells under the microscope, for me it was a matter of fascination that hidden within them were the key events converting normal cells into cancer cells, and frustration because they were out of reach. So, I took 3 years break from medicine to study for a PhD, learning the methods and thinking of molecular oncology in Colin Cooper's laboratory at the Institute …

Complex correlations: replication timing and mutational landscapes during cancer and genome evolution

DNA replication is a dynamic process that occurs in a temporal order along each of the chromosomes. A consequence of the temporally coordinated activation of replication origins is the establishment of broad domains (>100 kb) that replicate either early or late in S phase. This partitioning of the genome into early and late replication domains is important for maintaining genome stability, gene dosage, and epigenetic inheritance; however, the molecular mechanisms that define and establish these domains are poorly understood. The modENCODE Project provided an opportunity to investigate the chromatin features that define the Drosophila replication timing program in multiple cell lines. The majority of early and late replicating domains in the Drosophila genome were static across all cell lines; however, a small subset of domains was dynamic and exhibited differences in replication timing between the cell lines. Both origin selection and activation contribute to defining the DNA replication program. Our results suggest that static early and late replicating domains were defined at the level of origin selection (ORC binding) and likely mediated by chromatin accessibility. In contrast, dynamic domains exhibited low ORC densities in both cell types, suggesting that origin activation and not origin selection governs the plasticity of the DNA replication program. Finally, we show that the male-specific early replication of the X chromosome is dependent on the dosage compensation complex (DCC), suggesting that the transcription and replication programs respond to the same chromatin cues. Specifically, MOF-mediated hyperacetylation of H4K16 on the X chromosome promotes both the up-regulation of male-specific transcription and origin activation.

BackgroundThere has been a growing appreciation recently that mutagenic processes can be studied through the lenses of mutational signatures, which represent characteristic mutation patterns attributed to individual mutagens. However, the causal links between mutagens and observed mutation patterns as well as other types of interactions between mutagenic processes and molecular pathways are not fully understood, limiting the utility of mutational signatures.MethodsTo gain insights into these relationships, we developed a network-based method, named GeneSigNet that constructs an influence network among genes and mutational signatures. The approach leverages sparse partial correlation among other statistical techniques to uncover dominant influence relations between the activities of network nodes.ResultsApplying GeneSigNet to cancer data sets, we uncovered important relations between mutational signatures and several cellular processes that can shed light on cancer-related processes. Our results are consistent with previous findings, such as the impact of homologous recombination deficiency on clustered APOBEC mutations in breast cancer. The network identified by GeneSigNet also suggest an interaction between APOBEC hypermutation and activation of regulatory T Cells (Tregs), as well as a relation between APOBEC mutations and changes in DNA conformation. GeneSigNet also exposed a possible link between the SBS8 signature of unknown etiology and the Nucleotide Excision Repair (NER) pathway.ConclusionsGeneSigNet provides a new and powerful method to reveal the relation between mutational signatures and gene expression. The GeneSigNet method was implemented in python, and installable package, source codes and the data sets used for and generated during this study are available at the Github site https://github.com/ncbi/GeneSigNet.

Recent evidence suggests that the timing of DNA replication is coordinated across megabase-scale domains in metazoan genomes, yet the importance of this aspect of genome organization is unclear. Here we show that replication timing is remarkably conserved between human and mouse, uncovering large regions that may have been governed by similar replication dynamics since these species have diverged. This conservation is both tissue-specific and independent of the genomic G+C content conservation. Moreover, we show that time of replication is globally conserved despite numerous large-scale genome rearrangements. We systematically identify rearrangement fusion points and demonstrate that replication time can be locally diverged at these loci. Conversely, rearrangements are shown to be correlated with early replication and physical chromosomal proximity. These results suggest that large chromosomal domains of coordinated replication are shuffled by evolution while conserving the large-scale nuclear architecture of the genome.

Analysis of Mutational Signatures Suggest That Aid Has an Early and Driver Role in Multiple Myeloma

The human beta globin locus replicates late in most cell types, but becomes early replicating in erythroid cells. Using FISH to map DNA replication timing around the endogenous beta globin locus and by applying a genetic approach in transgenic mice, we have demonstrated that both the late and early replication states are controlled by regulatory elements within the locus control region. These results also show that the pattern of replication timing is set up by mechanisms that work independently of gene transcription.

Breast cancer, a complex and heterogeneous disease driven by genetic mutations in breast tissue cells, remains a leading cause of cancer-related mortality among women globally. A mutational signature can reveal the genomic landscape and history of breast cancer as it reflects the cumulative effect of various mutational processes that operate in cancer cells. This review provides an overview of the concept and classification of mutational signatures and discusses their clinical implications for breast cancer. We highlight how mutational signatures can provide insights into the therapeutic strategies, prognostic indicators, resistance mechanisms, and evolution of mutational signatures during treatment. Besides, we explore the potential applications of mutational signatures in personalized medicine for breast cancer, such as their integration with genomic profiling, prediction of treatment response, monitoring of treatment progression, and tailoring of therapeutic regimens based on signature analysis. We also address the challenges and limitations that need to be overcome before mutational signatures can be fully exploited for clinical benefit, such as the technical issues of data interpretation and standardization, the clinical translation of signature-based biomarkers, the exploration of emerging mutational signatures, and the longitudinal study of signature evolution. Future directions in mutational signature research encompass the exploration of emerging signatures, longitudinal studies to capture signature evolution, and the application of artificial intelligence to enhance signature detection and interpretation. While challenges remain, mutational signatures in breast cancer stand as a powerful tool that can revolutionize diagnosis and treatment, ultimately advancing our understanding and management of this complex disease.

BackgroundAllelic counterparts of biallelically expressed genes display an epigenetic symmetry normally manifested by synchronous replication, different from genes subjected to monoallelic expression, which normally are characterized by an asynchronous mode of replication (well exemplified by the SNRPN imprinted locus). Malignancy was documented to be associated with gross modifications in the inherent replication-timing coordination between allelic counterparts of imprinted genes as well as of biallelically expressed loci. The cancer-related allelic replication timing aberrations are non-disease specific and appear in peripheral blood cells of cancer patients, including those with solid tumors. As such they offer potential blood markers for non-invasive cancer test. The present study was aimed to gain some insight into the mechanism leading to the replication timing alterations of genes in blood lymphocytes of cancer patients.MethodsPeripheral blood samples derived from patients with prostate cancer were chosen to represent the cancerous status, and samples taken from patients with no cancer but with benign prostate hyperplasia were used to portray the normal status. Fluorescence In Situ Hybridization (FISH) replication assay, applied to phytohemagglutinin (PHA)-stimulated blood lymphocytes, was used to evaluate the temporal order (either synchronous or asynchronous) of genes in the patients' cells.ResultsWe demonstrated that: (i) the aberrant epigenetic profile, as delineated by the cancer status, is a reversible modification, evidenced by our ability to restore the normal patterns of replication in three unrelated loci (CEN15, SNRPN and RB1) by introducing an archetypical demethylating agent, 5-azacytidine; (ii) following the rehabilitating effect of demethylation, an imprinted gene (SNRPN) retains its original parental imprint; and (iii) the choice of an allele between early or late replication in the aberrant asynchronous replication, delineated by the cancer status, is not random but is independent of the parental origin.ConclusionThe non-disease specific aberrant epigenetic profile displayed in peripheral blood cells of patients with a solid tumour (unlike genetic aberrations) can be reversed, by an epigenetic drug applied in vitro, to the normal. It appears that the cancerous status differentiates between two allelic counterparts in a non-random manner, but independent of the parental origin

The timing of DNA replication generally correlates with transcription, gene density and sequence composition. How is the timing affected if a genomic region has a combination of features that individually correlate with either early or late replication? The Major Histocompatibility Complex (MHC) class II region is an AT-rich isochore that would be expected to replicate late, but it also contains coordinately regulated genes that are highly expressed in antigen-presenting cells and are strongly inducible in other cell types. Using cytological and biochemical assays, we find that the entire MHC replicates within the first half of S-phase, and that the class II region replicates slightly later than the adjacent regions irrespective of gene expression. These data suggest that despite AT-richness, an early-to-middle replication time in the class II region is defined by an open chromatin conformation that allows rapid transcriptional activation as a defence against pathogens.

Genomes in eukaryotes normally undergo DNA replication in a choreographed temporal order, resulting in early and late replicating chromosome compartments. Leishmania, a human protozoan parasite, displays an unconventional DNA replication program in which the timing of DNA replication completion is chromosome size-dependent: larger chromosomes complete replication later then smaller ones. Here we show that both R-loops and RNase H1, a ribonuclease that resolves RNA-DNA hybrids, accumulate in Leishmania major chromosomes in a pattern that reflects their replication timing. Furthermore, we demonstrate that such differential organisation of R-loops, RNase H1 and DNA replication timing across the parasite’s chromosomes correlates with size-dependent differences in chromatin accessibility, G quadruplex distribution and sequence content. Using conditional gene excision, we show that loss of RNase H1 leads to transient growth perturbation and permanently abrogates the differences in DNA replication timing across chromosomes, as well as altering levels of aneuploidy and increasing chromosome instability in a size-dependent manner. This work provides a link between R-loop homeostasis and DNA replication timing in a eukaryotic parasite and demonstrates that orchestration of DNA replication dictates levels of genome plasticity in Leishmania.

Driver mutations are the genetic variants responsible for oncogenesis, but how specific somatic mutational events arise in cells remains poorly understood. Mutational signatures derive from the frequency of mutated trinucleotides in a given cancer sample, and they provide an avenue for investigating the underlying mutational processes that operate in cancer. Here we analyse somatic mutations from 7,815 cancer exomes from The Cancer Genome Atlas (TCGA) across 26 cancer types. We curate a list of 50 known cancer driver mutations by analysing recurrence in our cohort and annotations of known cancer-associated genes from the Cancer Gene Census, IntOGen database and Cancer Genome Interpreter. We then use these datasets to perform binary univariate logistic regression and establish the statistical relationship between individual driver mutations and known mutational signatures across different cancer types. Our analysis led to the identification of 39 significant associations between driver mutations and mutational signatures (P < 0.004, with a false discovery rate of < 5%). We first validate our methodology by establishing statistical links for known and novel associations between driver mutations and the mutational signature arising from Polymerase Epsilon proofreading deficiency. We then examine associations between driver mutations and mutational signatures for AID/APOBEC enzyme activity and deficient mismatch repair. We also identify negative associations (odds ratio < 1) between mutational signatures and driver mutations, and here we examine the role of aging and cigarette smoke mutagenesis in the generation of driver mutations in IDH1 and KRAS in brain cancers and lung adenocarcinomas respectively. Our study provides statistical foundations for hypothesised links between otherwise independent biological processes and we uncover previously unexplored relationships between driver mutations and mutagenic processes during cancer development. These associations give insights into how cancers acquire advantageous mutations and can provide direction to guide further mechanistic studies into cancer pathogenesis.