Somatic and Germline Mutation Periodicity Follow the Orientation of the DNA Minor Groove around Nucleosomes

Strand-resolved mutagenicity of DNA damage and repair

DNA base damage is a major source of oncogenic mutations1. Such damage can produce strand-phased mutation patterns and multiallelic variation through the process of lesion segregation2. Here we exploited these properties to reveal how strand-asymmetric processes, such as replication and transcription, shape DNA damage and repair. Despite distinct mechanisms of leading and lagging strand replication3,4, we observe identical fidelity and damage tolerance for both strands. For small alkylation adducts of DNA, our results support a model in which the same translesion polymerase is recruited on-the-fly to both replication strands, starkly contrasting the strand asymmetric tolerance of bulky UV-induced adducts5. The accumulation of multiple distinct mutations at the site of persistent lesions provides the means to quantify the relative efficiency of repair processes genome wide and at single-base resolution. At multiple scales, we show DNA damage-induced mutations are largely shaped by the influence of DNA accessibility on repair efficiency, rather than gradients of DNA damage. Finally, we reveal specific genomic conditions that can actively drive oncogenic mutagenesis by corrupting the fidelity of nucleotide excision repair. These results provide insight into how strand-asymmetric mechanisms underlie the formation, tolerance and repair of DNA damage, thereby shaping cancer genome evolution.

Future in the past: paternal reprogramming of offspring phenotype and the epigenetic mechanisms

That certain preconceptual paternal exposures reprogram the developmental phenotypic plasticity in future generation(s) has conceptualized the "paternal programming of offspring health" hypothesis. This transgenerational effect is transmitted primarily through sperm epigenetic mechanisms-DNA methylation, non-coding RNAs (ncRNAs) and associated RNA modifications, and histone modifications-and potentially through non-sperm-specific mechanisms-seminal plasma and circulating factors-that create 'imprinted' memory of ancestral information. The epigenetic landscape in sperm is highly responsive to environmental cues, due to, in part, the soma-to-germline communication mediated by epididymosomes. While human epidemiological studies and experimental animal studies have provided solid evidences in support of transgenerational epigenetic inheritance, how ancestral information is memorized as epigenetic codes for germline transmission is poorly understood. Particular elusive is what the downstream effector pathways that decode those epigenetic codes into persistent phenotypes. In this review, we discuss the paternal reprogramming of offspring phenotype and the possible underlying epigenetic mechanisms. Cracking these epigenetic mechanisms will lead to a better appreciation of "Paternal Origins of Health and Disease" and guide innovation of intervention algorithms to achieve 'healthier' outcomes in future generations. All this will revolutionize our understanding of human disease etiology.

Mutation rate heterogeneity at the sub-gene scale due to local DNA hypomethylation

Local mutation rates in human are highly heterogeneous, with known variability at the scale of megabase-sized chromosomal domains, and, on the other extreme, at the scale of oligonucleotides. The intermediate, kilobase-scale heterogeneity in mutation risk is less well characterized. Here, by analyzing thousands of somatic genomes, we studied mutation risk gradients along gene bodies, representing a genomic scale spanning roughly 1-10 kb, hypothesizing that different mutational mechanisms are differently distributed across gene segments. The main heterogeneity concerns several kilobases at the transcription start site and further downstream into 5' ends of gene bodies; these are commonly hypomutated with several mutational signatures, most prominently the ubiquitous C > T changes at CpG dinucleotides. The width and shape of this mutational coldspot at 5' gene ends is variable across genes, and corresponds to variable interval of lowered DNA methylation depending on gene activity level and regulation. Such hypomutated loci, at 5' gene ends or elsewhere, correspond to DNA hypomethylation that can associate with various landmarks, including intragenic enhancers, Polycomb-marked regions, or chromatin loop anchor points. Tissue-specific DNA hypomethylation begets tissue-specific local hypomutation. Of note, direction of mutation risk is inverted for AID/APOBEC3 cytosine deaminase activity, whose signatures are enriched in hypomethylated regions.

DNA Repair in Nucleosomes: Insights from Histone Modifications and Mutants

DNA repair pathways play a critical role in genome stability, but in eukaryotic cells, they must operate to repair DNA lesions in the compact and tangled environment of chromatin. Previous studies have shown that the packaging of DNA into nucleosomes, which form the basic building block of chromatin, has a profound impact on DNA repair. In this review, we discuss the principles and mechanisms governing DNA repair in chromatin. We focus on the role of histone post-translational modifications (PTMs) in repair, as well as the molecular mechanisms by which histone mutants affect cellular sensitivity to DNA damage agents and repair activity in chromatin. Importantly, these mechanisms are thought to significantly impact somatic mutation rates in human cancers and potentially contribute to carcinogenesis and other human diseases. For example, a number of the histone mutants studied primarily in yeast have been identified as candidate oncohistone mutations in different cancers. This review highlights these connections and discusses the potential importance of DNA repair in chromatin to human health.

Genome-wide impact of cytosine methylation and DNA sequence context on UV-induced CPD formation

Exposure to ultraviolet (UV) light is the primary etiological agent for skin cancers because UV damages cellular DNA. The most frequent form of UV damage is the cyclobutane pyrimidine dimer (CPD), which consists of covalent linkages between neighboring pyrimidine bases in DNA. In human cells, the 5' position of cytosine bases in CG dinucleotides is frequently methylated, and methylated cytosines in the TP53 tumor suppressor are often sites of mutation hotspots in skin cancers. It has been argued that this is because cytosine methylation promotes UV-induced CPD formation; however, the effects of cytosine methylation on CPD formation are controversial, with conflicting results from previous studies. Here, we use a genome-wide method known as CPD-seq to map UVB- and UVC-induced CPDs across the yeast genome in the presence or absence in vitro methylation by the CpG methyltransferase M.SssI. Our data indicate that cytosine methylation increases UVB-induced CPD formation nearly 2-fold relative to unmethylated DNA, but the magnitude of induction depends on the flanking sequence context. Sequence contexts with a 5' guanine base (e.g., GCCG and GTCG) show the strongest induction due to cytosine methylation, potentially because these sequence contexts are less efficient at forming CPD lesions in the absence of methylation. We show that cytosine methylation also modulates UVC-induced CPD formation, albeit to a lesser extent than UVB. These findings can potentially reconcile previous studies, and define the impact of cytosine methylation on UV damage across a eukaryotic genome.

Synthetic reversed sequences reveal default genomic states

Pervasive transcriptional activity is observed across diverse species. The genomes of extant organisms have undergone billions of years of evolution, making it unclear whether these genomic activities represent effects of selection or 'noise'1-4. Characterizing default genome states could help understand whether pervasive transcriptional activity has biological meaning. Here we addressed this question by introducing a synthetic 101-kb locus into the genomes of Saccharomyces cerevisiae and Mus musculus and characterizing genomic activity. The locus was designed by reversing but not complementing human HPRT1, including its flanking regions, thus retaining basic features of the natural sequence but ablating evolved coding or regulatory information. We observed widespread activity of both reversed and native HPRT1 loci in yeast, despite the lack of evolved yeast promoters. By contrast, the reversed locus displayed no activity at all in mouse embryonic stem cells, and instead exhibited repressive chromatin signatures. The repressive signature was alleviated in a locus variant lacking CpG dinucleotides; nevertheless, this variant was also transcriptionally inactive. These results show that synthetic genomic sequences that lack coding information are active in yeast, but inactive in mouse embryonic stem cells, consistent with a major difference in 'default genomic states' between these two divergent eukaryotic cell types, with implications for understanding pervasive transcription, horizontal transfer of genetic information and the birth of new genes.

Analysis of cytosine deamination events in excision repair sequencing reads reveals mechanisms of incision site selection in NER

Nucleotide excision repair (NER) removes helix-distorting DNA lesions and is therefore critical for genome stability. During NER, DNA is unwound on either side of the lesion and excised, but the rules governing incision site selection, particularly in eukaryotic cells, are unclear. Excision repair-sequencing (XR-seq) sequences excised NER fragments, but analysis has been limited because the lesion location is unknown. Here, we exploit accelerated cytosine deamination rates in UV-induced CPD (cyclobutane pyrimidine dimer) lesions to precisely map their locations at C to T mismatches in XR-seq reads, revealing general and species-specific patterns of incision site selection during NER. Our data indicate that the 5' incision site occurs preferentially in HYV (i.e. not G; C/T; not T) sequence motifs, a pattern that can be explained by sequence preferences of the XPF-ERCC1 endonuclease. In contrast, the 3' incision site does not show strong sequence preferences, once truncated reads arising from mispriming events are excluded. Instead, the 3' incision is partially determined by the 5' incision site distance, indicating that the two incision events are coupled. Finally, our data reveal unique and coupled NER incision patterns at nucleosome boundaries. These findings reveal key principles governing NER incision site selection in eukaryotic cells.

Evaluating topography of mutational signatures with SigProfilerTopography

The mutations found in a cancer genome are shaped by diverse processes, each displaying a characteristic mutational signature that may be influenced by the genome's architecture. While prior analyses have evaluated the effect of topographical genomic features on mutational signatures, there has been no computational tool that can comprehensively examine this interplay. Here, we present SigProfilerTopography, a Python package that allows evaluating the effect of chromatin organization, histone modifications, transcription factor binding, DNA replication, and DNA transcription on the activities of different mutational processes. SigProfilerTopography elucidates the unique topographical characteristics of mutational signatures, unveiling their underlying biological and molecular mechanisms.

Hotspot propensity across mutational processes

The sparsity of mutations observed across tumours hinders our ability to study mutation rate variability at nucleotide resolution. To circumvent this, here we investigated the propensity of mutational processes to form mutational hotspots as a readout of their mutation rate variability at single base resolution. Mutational signatures 1 and 17 have the highest hotspot propensity (5-78 times higher than other processes). After accounting for trinucleotide mutational probabilities, sequence composition and mutational heterogeneity at 10 Kbp, most (94-95%) signature 17 hotspots remain unexplained, suggesting a significant role of local genomic features. For signature 1, the inclusion of genome-wide distribution of methylated CpG sites into models can explain most (80-100%) of the hotspot propensity. There is an increased hotspot propensity of signature 1 in normal tissues and de novo germline mutations. We demonstrate that hotspot propensity is a useful readout to assess the accuracy of mutation rate models at nucleotide resolution. This new approach and the findings derived from it open up new avenues for a range of somatic and germline studies investigating and modelling mutagenesis.

Detection of new pioneer transcription factors as cell-type specific nucleosome binders

Wrapping of DNA into nucleosomes restricts accessibility to the DNA and may affect the recognition of binding motifs by transcription factors. A certain class of transcription factors, the pioneer transcription factors, can specifically recognize their DNA binding sites on nucleosomes, may initiate local chromatin opening and facilitate the binding of co-factors in a cell-type-specific manner. For the majority of human pioneer transcription factors, the locations of their binding sites, mechanisms of binding and regulation remain unknown. We have developed a computational method to predict the cell-type-specific ability of transcription factors to bind nucleosomes by integrating ChIP-seq, MNase-seq and DNase-seq data with details of nucleosome structure. We have demonstrated the ability of our approach in discriminating pioneer from canonical transcription factors and predicted new potential pioneer transcription factors in H1, K562, HepG2 and HeLa cell lines. Lastly, we systemically analyzed the interaction modes between various pioneer transcription factors and detected several clusters of distinctive binding sites on nucleosomal DNA.

Novel insights into bulky DNA damage formation and nucleotide excision repair from high-resolution genomics

DNA damages compromise cell function and fate. Cells of all organisms activate a global DNA damage response that includes a signaling stress response, activation of checkpoints, and recruitment of repair enzymes. Especially deleterious are bulky, helix-distorting damages that block transcription and replication. Due to their miscoding nature, these damages lead to mutations and cancer. In human cells, bulky DNA damages are repaired by nucleotide excision repair (NER). To date, the basic mechanism of NER in naked DNA is well defined. Still, there is a fundamental gap in our understanding of how repair is orchestrated despite the packaging of DNA in chromatin, and how it is coordinated with active transcription and replication. The last decade has brought forth huge advances in our ability to detect and assay bulky DNA damages and their repair at single nucleotide resolution across the human genome. Here we review recent findings on the effect of chromatin and DNA-binding proteins on the formation of bulky DNA damages, and novel insights on NER, provided by the recent application of genomic methods.

Predicting regional somatic mutation rates using DNA motifs

How the locus-specificity of epigenetic modifications is regulated remains an unanswered question. A contributing mechanism is that epigenetic enzymes are recruited to specific loci by DNA binding factors recognizing particular sequence motifs (referred to as epi-motifs). Using these motifs to predict biological outputs depending on local epigenetic state such as somatic mutation rates would confirm their functionality. Here, we used DNA motifs including known TF motifs and epi-motifs as a surrogate of epigenetic signals to predict somatic mutation rates in 13 cancers at an average 23kbp resolution. We implemented an interpretable neural network model, called contextual regression, to successfully learn the universal relationship between mutations and DNA motifs, and uncovered motifs that are most impactful on the regional mutation rates such as TP53 and epi-motifs associated with H3K9me3. Furthermore, we identified genomic regions with significantly higher mutation rates than the expected values in each individual tumor and demonstrated that such cancer-related regions can accurately predict cancer types. Interestingly, we found that the same mutation signatures often have different contributions to cancer-related and cancer-independent regions, and we also identified the motifs with the most contribution to each mutation signature.

A half century of exploring DNA excision repair in chromatin

DNA in eukaryotic cells is packaged into the compact and dynamic structure of chromatin. This packaging is a double-edged sword for DNA repair and genomic stability. Chromatin restricts the access of repair proteins to DNA lesions embedded in nucleosomes and higher order chromatin structures. However, chromatin also serves as a signaling platform in which post-translational modifications of histones and other chromatin-bound proteins promote lesion recognition and repair. Similarly, chromatin modulates the formation of DNA damage, promoting or suppressing lesion formation depending on the chromatin context. Therefore, the modulation of DNA damage and its repair in chromatin is crucial to our understanding of the fate of potentially mutagenic and carcinogenic lesions in DNA. Here, we survey many of the landmark findings on DNA damage and repair in chromatin over the last 50 years (i.e., since the beginning of this field), focusing on excision repair, the first repair mechanism studied in the chromatin landscape. For example, we highlight how the impact of chromatin on these processes explains the distinct patterns of somatic mutations observed in cancer genomes.

Topography of mutational signatures in human cancer

The somatic mutations found in a cancer genome are imprinted by different mutational processes. Each process exhibits a characteristic mutational signature, which can be affected by the genome architecture. However, the interplay between mutational signatures and topographical genomic features has not been extensively explored. Here, we integrate mutations from 5,120 whole-genome-sequenced tumors from 40 cancer types with 516 topographical features from ENCODE to evaluate the effect of nucleosome occupancy, histone modifications, CTCF binding, replication timing, and transcription/replication strand asymmetries on the cancer-specific accumulation of mutations from distinct mutagenic processes. Most mutational signatures are affected by topographical features, with signatures of related etiologies being similarly affected. Certain signatures exhibit periodic behaviors or cancer-type-specific enrichments/depletions near topographical features, revealing further information about the processes that imprinted them. Our findings, disseminated via the COSMIC (Catalog of Somatic Mutations in Cancer) signatures database, provide a comprehensive online resource for exploring the interactions between mutational signatures and topographical features across human cancer.

Sequence dependencies and mutation rates of localized mutational processes in cancer

BACKGROUND

Cancer mutations accumulate through replication errors and DNA damage coupled with incomplete repair. Individual mutational processes often show nucleotide sequence and functional region preferences. As a result, some sequence contexts mutate at much higher rates than others, with additional variation found between functional regions. Mutational hotspots, with recurrent mutations across cancer samples, represent genomic positions with elevated mutation rates, often caused by highly localized mutational processes.

METHODS

We count the 11-mer genomic sequences across the genome, and using the PCAWG set of 2583 pan-cancer whole genomes, we associate 11-mers with mutational signatures, hotspots of single nucleotide variants, and specific genomic regions. We evaluate the mutation rates of individual and combined sets of 11-mers and derive mutational sequence motifs.

RESULTS

We show that hotspots generally identify highly mutable sequence contexts. Using these, we show that some mutational signatures are enriched in hotspot sequence contexts, corresponding to well-defined sequence preferences for the underlying localized mutational processes. This includes signature 17b (of unknown etiology) and signatures 62 (POLE deficiency), 7a (UV), and 72 (linked to lymphomas). In some cases, the mutation rate and sequence preference increase further when focusing on certain genomic regions, such as signature 62 in transcribed regions, where the mutation rate is increased up to 9-folds over cancer type and mutational signature average.

CONCLUSIONS

We summarize our findings in a catalog of localized mutational processes, their sequence preferences, and their estimated mutation rates.

Structural underpinnings of mutation rate variations in the human genome

Single nucleotide mutation rates have critical implications for human evolution and genetic diseases. Importantly, the rates vary substantially across the genome and the principles underlying such variations remain poorly understood. A recent model explained much of this variation by considering higher-order nucleotide interactions in the 7-mer sequence context around mutated nucleotides. This model's success implicates a connection between DNA shape and mutation rates. DNA shape, i.e. structural properties like helical twist and tilt, is known to capture interactions between nucleotides within a local context. Thus, we hypothesized that changes in DNA shape features at and around mutated positions can explain mutation rate variations in the human genome. Indeed, DNA shape-based models of mutation rates showed similar or improved performance over current nucleotide sequence-based models. These models accurately characterized mutation hotspots in the human genome and revealed the shape features whose interactions underlie mutation rate variations. DNA shape also impacts mutation rates within putative functional regions like transcription factor binding sites where we find a strong association between DNA shape and position-specific mutation rates. This work demonstrates the structural underpinnings of nucleotide mutations in the human genome and lays the groundwork for future models of genetic variations to incorporate DNA shape.

Mutational signature dynamics shaping the evolution of oesophageal adenocarcinoma

A variety of mutational processes drive cancer development, but their dynamics across the entire disease spectrum from pre-cancerous to advanced neoplasia are poorly understood. We explore the mutagenic processes shaping oesophageal adenocarcinoma tumorigenesis in 997 instances comprising distinct stages of this malignancy, from Barrett Oesophagus to primary tumours and advanced metastatic disease. The mutational landscape is dominated by the C[T > C/G]T substitution enriched signatures SBS17a/b, which are linked with TP53 mutations, increased proliferation, genomic instability and disease progression. The APOBEC mutagenesis signature is a weak but persistent signal amplified in primary tumours. We also identify prevalent alterations in DNA damage repair pathways, with homologous recombination, base and nucleotide excision repair and translesion synthesis mutated in up to 50% of the cohort, and surprisingly uncoupled from transcriptional activity. Among these, the presence of base excision repair deficiencies show remarkably poor prognosis in the cohort. In this work, we provide insights on the mutational aetiology and changes enabling the transition from pre-neoplastic to advanced oesophageal adenocarcinoma.

Structural and mechanistic insights into the DNA glycosylase AAG-mediated base excision in nucleosome

The engagement of a DNA glycosylase with a damaged DNA base marks the initiation of base excision repair. Nucleosome-based packaging of eukaryotic genome obstructs DNA accessibility, and how DNA glycosylases locate the substrate site on nucleosomes is currently unclear. Here, we report cryo-electron microscopy structures of nucleosomes bearing a deoxyinosine (DI) in various geometric positions and structures of them in complex with the DNA glycosylase AAG. The apo nucleosome structures show that the presence of a DI alone perturbs nucleosomal DNA globally, leading to a general weakening of the interface between DNA and the histone core and greater flexibility for the exit/entry of the nucleosomal DNA. AAG makes use of this nucleosomal plasticity and imposes further local deformation of the DNA through formation of the stable enzyme-substrate complex. Mechanistically, local distortion augmentation, translation/rotational register shift and partial opening of the nucleosome are employed by AAG to cope with substrate sites in fully exposed, occluded and completely buried positions, respectively. Our findings reveal the molecular basis for the DI-induced modification on the structural dynamics of the nucleosome and elucidate how the DNA glycosylase AAG accesses damaged sites on the nucleosome with different solution accessibility.

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.

Damage mapping techniques and the light they have shed on canonical and atypical UV photoproducts

Ultraviolet (UV) light is a pervasive threat to the DNA of terrestrial organisms. UV light induces helix-distorting DNA lesions, primarily cyclobutane pyrimidine dimers (CPDs) that form between neighboring pyrimidine bases. Unrepaired CPD lesions cause cytosine-to-thymine (C>T) substitutions in dipyrimidine sequences, which is the predominant mutation class in skin cancer genomes. However, many driver mutations in melanoma (e.g., in the BRAF and NRAS oncogenes) do not fit this UV mutation signature. Recent studies have brought to light the intriguing hypothesis that these driver mutations may be induced by infrequent or atypical UV photoproducts, including pyrimidine 6-4 pyrimidone photoproducts (6-4PP) and thymine-adenine (TA) photoproducts. Here, we review innovative methods for mapping both canonical and atypical UV-induced photoproducts across the genome.

DNA replication initiation shapes the mutational landscape and expression of the human genome

The interplay between active biological processes and DNA repair is central to mutagenesis. Here, we show that the ubiquitous process of replication initiation is mutagenic, leaving a specific mutational footprint at thousands of early and efficient replication origins. The observed mutational pattern is consistent with two distinct mechanisms, reflecting the two-step process of origin activation, triggering the formation of DNA breaks at the center of origins and local error-prone DNA synthesis in their immediate vicinity. We demonstrate that these replication initiation-dependent mutational processes exert an influence on phenotypic diversity in humans that is disproportionate to the origins' genomic size: By increasing mutational loads at gene promoters and splice junctions, the presence of an origin significantly influences both gene expression and mRNA isoform usage. Last, we show that mutagenesis at origins not only drives the evolution of origin sequences but also contributes to sculpting regulatory domains of the human genome.

The co-evolution of the genome and epigenome in colorectal cancer

Colorectal malignancies are a leading cause of cancer-related death1 and have undergone extensive genomic study2,3. However, DNA mutations alone do not fully explain malignant transformation4-7. Here we investigate the co-evolution of the genome and epigenome of colorectal tumours at single-clone resolution using spatial multi-omic profiling of individual glands. We collected 1,370 samples from 30 primary cancers and 8 concomitant adenomas and generated 1,207 chromatin accessibility profiles, 527 whole genomes and 297 whole transcriptomes. We found positive selection for DNA mutations in chromatin modifier genes and recurrent somatic chromatin accessibility alterations, including in regulatory regions of cancer driver genes that were otherwise devoid of genetic mutations. Genome-wide alterations in accessibility for transcription factor binding involved CTCF, downregulation of interferon and increased accessibility for SOX and HOX transcription factor families, suggesting the involvement of developmental genes during tumourigenesis. Somatic chromatin accessibility alterations were heritable and distinguished adenomas from cancers. Mutational signature analysis showed that the epigenome in turn influences the accumulation of DNA mutations. This study provides a map of genetic and epigenetic tumour heterogeneity, with fundamental implications for understanding colorectal cancer biology.

Evaluation of Mutagenicity and Anti-Mutagenicity of Various Bean Milks Using Drosophila with High Bioactivation

The consumption of a nutritious diet including phytochemicals can minimize mutations as the primary cause of carcinogenesis. Bean consumption supplies calories, minerals and phytochemicals but their anti-mutagenic properties in vivo remain little understood. Hence, the present study aimed to study the mutagenicity and anti-mutagenic properties of five bean milks using the somatic mutation and recombination test (SMART) involving Drosophila with high bioactivation. Milk derived from five bean varieties, namely black bean (Phaseolus vulgaris), red kidney bean (Phaseolus vulgaris), mung bean (Phaseolus aureus), peanut (Arachis hypogaea) and soybean (Glycine max) did not induce DNA mutations in Drosophila with high bioactivation, indicating their genome-safe properties. All bean milks showed anti-mutagenicity against the food-derived mutagen, urethane, in vivo with different degrees of inhibition. In the co-administration study, larvae were treated with each bean milk together with urethane. Soybean milk showed the highest anti-mutagenicity at 27.75%; peanut milk exhibited the lowest at 7.51%. In the pre-feeding study, the larvae received each bean milk followed by urethane. Soybean milk exhibited the highest anti-mutagenic potential, followed by red kidney bean and black bean milks. Total phenolic and antioxidant data revealed that the anti-mutagenicity of both red kidney bean milk and black bean milk might be derived from their phenolic or antioxidant properties; other phytochemicals may contribute to the high anti-mutagenicity observed in soybean milk. Further investigations on the anti-mutagenicity of bean milks against other dietary mutagens are required to develop bean-based products with potent anti-mutagenic properties.

Effects of replication domains on genome-wide UV-induced DNA damage and repair

Nucleotide excision repair is the primary repair mechanism that removes UV-induced DNA lesions in placentals. Unrepaired UV-induced lesions could result in mutations during DNA replication. Although the mutagenesis of pyrimidine dimers is reasonably well understood, the direct effects of replication fork progression on nucleotide excision repair are yet to be clarified. Here, we applied Damage-seq and XR-seq techniques and generated replication maps in synchronized UV-treated HeLa cells. The results suggest that ongoing replication stimulates local repair in both early and late replication domains. Additionally, it was revealed that lesions on lagging strand templates are repaired slower in late replication domains, which is probably due to the imbalanced sequence context. Asymmetric relative repair is in line with the strand bias of melanoma mutations, suggesting a role of exogenous damage, repair, and replication in mutational strand asymmetry.

Structural basis for APE1 processing DNA damage in the nucleosome

Genomic DNA is continually exposed to endogenous and exogenous factors that promote DNA damage. Eukaryotic genomic DNA is packaged into nucleosomes, which present a barrier to accessing and effectively repairing DNA damage. The mechanisms by which DNA repair proteins overcome this barrier to repair DNA damage in the nucleosome and protect genomic stability is unknown. Here, we determine how the base excision repair (BER) endonuclease AP-endonuclease 1 (APE1) recognizes and cleaves DNA damage in the nucleosome. Kinetic assays determine that APE1 cleaves solvent-exposed AP sites in the nucleosome with 3 − 6 orders of magnitude higher efficiency than occluded AP sites. A cryo-electron microscopy structure of APE1 bound to a nucleosome containing a solvent-exposed AP site reveal that APE1 uses a DNA sculpting mechanism for AP site recognition, where APE1 bends the nucleosomal DNA to access the AP site. Notably, additional biochemical and structural characterization of occluded AP sites identify contacts between the nucleosomal DNA and histone octamer that prevent efficient processing of the AP site by APE1. These findings provide a rationale for the position-dependent activity of BER proteins in the nucleosome and suggests the ability of BER proteins to sculpt nucleosomal DNA drives efficient BER in chromatin.

AP endonuclease 1 (APE1) processes genomic AP sites during base excision repair. Here, the authors determine the structural mechanism used by APE1 to process nucleosomal AP sites, providing new insight into DNA repair in chromatin.

The shaping of cancer genomes with the regional impact of mutation processes

Mutation signature analysis has been used to infer the contributions of various DNA mutagenic-repair events in individual cancer genomes. Here, we build a statistical framework using a multinomial distribution to assign individual mutations to their cognate mutation signatures. We applied it to 47 million somatic mutations in 1925 publicly available cancer genomes to obtain a mutation signature map at the resolution of individual somatic mutations. Based on mutation signature-level genetic-epigenetic correlative analyses, mutations with transcriptional and replicative strand asymmetries show different enrichment patterns across genomes, and “transcribed” chromatin states and gene boundaries are particularly vulnerable to transcription-coupled repair activities. While causative processes of cancer-driving mutations can be diverse, as shown for converging effects of multiple mutational processes on TP53 mutations, the substantial fraction of recurrently mutated amino acids points to specific mutational processes, e.g., age-related C-to-T transition for KRAS p.G12 mutations. Our investigation of evolutionary trajectories with respect to mutation signatures further revealed that candidate pairs of early- vs. late-operative mutation processes in cancer genomes represent evolutionary dynamics of multiple mutational processes in the shaping of cancer genomes. We also observed that the local mutation clusters of kataegis often include mutations arising from multiple mutational processes, suggestive of a locally synchronous impact of multiple mutational processes on cancer genomes. Taken together, our examination of the genome-wide landscape of mutation signatures at the resolution of individual somatic mutations shows the spatially and temporally distinct mutagenesis-repair-replication histories of various mutational processes and their effects on shaping cancer genomes.

A statistical model that assigns non-hereditary DNA alterations known as somatic mutations to mutation “signatures” (groups of mutations arising from a specific biological process) on cancer genomes provides novel insights into disease evolution. Somatic mutations result from exposure to factors often linked to cancer development, such as tobacco or ultraviolet radiation. However, assigning a somatic mutation to a particular mutation “signature” remains challenging. The model created by Ruibin Xi (Peking University, China) and Tae-Min Kim (Catholic University of Korea, Seoul, South Korea) and co-workers grouped 47 million somatic mutations in 1925 cancer genomes into localized clusters before connecting them with mutation signatures. This strategy highlights the spatial and temporal patterns related to the origins of mutations, how the DNA strands are repaired and replicated, and how this influences the emerging cancer genome.

DNA methylation cues in nucleosome geometry, stability and unwrapping

Cytosine methylation at the 5-carbon position is an essential DNA epigenetic mark in many eukaryotic organisms. Although countless structural and functional studies of cytosine methylation have been reported, our understanding of how it influences the nucleosome assembly, structure, and dynamics remains obscure. Here, we investigate the effects of cytosine methylation at CpG sites on nucleosome dynamics and stability. By applying long molecular dynamics simulations on several microsecond time scale, we generate extensive atomistic conformational ensembles of full nucleosomes. Our results reveal that methylation induces pronounced changes in geometry for both linker and nucleosomal DNA, leading to a more curved, under-twisted DNA, narrowing the adjacent minor grooves, and shifting the population equilibrium of sugar-phosphate backbone geometry. These DNA conformational changes are associated with a considerable enhancement of interactions between methylated DNA and the histone octamer, doubling the number of contacts at some key arginines. H2A and H3 tails play important roles in these interactions, especially for DNA methylated nucleosomes. This, in turn, prevents a spontaneous DNA unwrapping of 3-4 helical turns for the methylated nucleosome with truncated histone tails, otherwise observed in the unmethylated system on several microseconds time scale.

NucPosDB: a database of nucleosome positioning in vivo and nucleosomics of cell-free DNA

Nucleosome positioning is involved in many gene regulatory processes happening in the cell, and it may change as cells differentiate or respond to the changing microenvironment in a healthy or diseased organism. One important implication of nucleosome positioning in clinical epigenetics is its use in the “nucleosomics” analysis of cell-free DNA (cfDNA) for the purpose of patient diagnostics in liquid biopsies. The rationale for this is that the apoptotic nucleases that digest chromatin of the dying cells mostly cut DNA between nucleosomes. Thus, the short pieces of DNA in body fluids reflect the positions of nucleosomes in the cells of origin. Here, we report a systematic nucleosomics database — NucPosDB — curating published nucleosome positioning datasets in vivo as well as datasets of sequenced cell-free DNA (cfDNA) that reflect nucleosome positioning in situ in the cells of origin. Users can select subsets of the database by a number of criteria and then obtain raw or processed data. NucPosDB also reports the originally determined regions with stable nucleosome occupancy across several individuals with a given condition. An additional section provides a catalogue of computational tools for the analysis of nucleosome positioning or cfDNA experiments and theoretical algorithms for the prediction of nucleosome positioning preferences from DNA sequence. We provide an overview of the field, describe the structure of the database in this context, and demonstrate data variability using examples of different medical conditions. NucPosDB is useful both for the analysis of fundamental gene regulation processes and the training of computational models for patient diagnostics based on cfDNA. The database currently curates ~ 400 publications on nucleosome positioning in cell lines and in situ as well as cfDNA from > 10,000 patients and healthy volunteers. For open-access cfDNA datasets as well as key MNase-seq datasets in human cells, NucPosDB allows downloading processed mapped data in addition to the regions with stable nucleosome occupancy. NucPosDB is available at https://generegulation.org/nucposdb/.

Molecular mechanism of UV damage modulation in nucleosomes

Exposure to ultraviolet (UV) light causes the formation of mutagenic cyclobutane pyrimidine dimers (CPDs) in cellular DNA. Previous studies have revealed that CPD formation in nucleosomes, the building blocks of chromatin, shows a striking ∼10 base pair (bp) periodic pattern. CPD formation is suppressed at positions where the DNA minor groove faces toward the histone octamer (minor-in) and elevated CPD formation at positions where the minor groove faces away from the histone octamer (minor-out). However, the molecular mechanism underlying this nucleosome photofootprint is unclear. Here, we analyzed ∼180 high-resolution nucleosome structures to characterize whether differences in DNA mobility or conformation are responsible for the CPD modulation in nucleosomes. Our results indicate that differences in DNA mobility cannot explain CPD modulation in nucleosome. Instead, we find that the sharp DNA bending around the histone octamer results in DNA conformations with structural parameters more susceptible to UV damage formation at minor-out positions and more resistant to CPD formation at minor-in positions. This analysis reveals the molecular mechanism responsible for periodic modulation of CPD formation and UV mutagenesis in nucleosomal DNA.

Evolution and progression of Barrett's oesophagus to oesophageal cancer

Cancer cells are shaped through an evolutionary process of DNA mutation, cell selection and population expansion. Early steps in this process are driven by a set of mutated driver genes and structural alterations to the genome through copy number gains or losses. Oesophageal adenocarcinoma (EAC) and the pre-invasive tissue, Barrett's oesophagus (BE), provide an ideal example in which to observe and study this evolution. BE displays early genomic instability, specifically in copy number changes that may later be observed in EAC. Furthermore, these early changes result in patterns of progression (that is, 'born bad', gradual or catastrophic) that may help to describe the evolution of EAC. As only a small proportion of patients with BE will go on to develop cancer, a better understanding of these patterns and the resulting genomic changes should improve early detection in EAC and may provide clues for the evolution of cancer more broadly.

Chaperoning histones at the DNA repair dance

Unlike all other biological molecules that are degraded and replaced if damaged, DNA must be repaired as chromosomes cannot be replaced. Indeed, DNA endures a wide variety of structural damage that need to be repaired accurately to maintain genomic stability and proper functioning of cells and to prevent mutation leading to disease. Given that the genome is packaged into chromatin within eukaryotic cells, it has become increasingly evident that the chromatin context of DNA both facilitates and regulates DNA repair processes. In this review, we discuss mechanisms involved in removal of histones (chromatin disassembly) from around DNA lesions, by histone chaperones and chromatin remodelers, that promotes accessibility of the DNA repair machinery. We also elaborate on how the deposition of core histones and specific histone variants onto DNA (chromatin assembly) during DNA repair promotes repair processes, the role of histone post translational modifications in these processes and how chromatin structure is reestablished after DNA repair is complete.

The origin of human mutation in light of genomic data

Despite years of active research into the role of DNA repair and replication in mutagenesis, surprisingly little is known about the origin of spontaneous human mutation in the germ line. With the advent of high-throughput sequencing, genome-scale data have revealed statistical properties of mutagenesis in humans. These properties include variation of the mutation rate and spectrum along the genome at different scales in relation to epigenomic features and dependency on parental age. Moreover, mutations originated in mothers are less frequent than mutations originated in fathers and have a distinct genomic distribution. Statistical analyses that interpret these patterns in the context of known biochemistry can provide mechanistic models of mutagenesis in humans.

Mutperiod: Analysis of periodic mutation rates in nucleosomes

Nucleosomes modulate DNA damage and repair, resulting in periodic mutation rates in nucleosomal DNA. Previous research has characterized these patterns in many sequenced tumor genomes; however, computational tools to identify and quantify these periodicities have not been developed for the broader scientific community. Here, we describe mutperiod, a Python and R based toolset that quantifies nucleosome mutational periodicities and compares them across different genetic and cellular backgrounds. We use mutperiod to demonstrate that DNA mismatch repair contributes to the nucleosome mutational periodicity observed in esophageal adenocarcinomas, and that the strength of this mutational periodicity varies in different chromatin states.

Binding of regulatory proteins to nucleosomes is modulated by dynamic histone tails

Little is known about the roles of histone tails in modulating nucleosomal DNA accessibility and its recognition by other macromolecules. Here we generate extensive atomic level conformational ensembles of histone tails in the context of the full nucleosome, totaling 65 microseconds of molecular dynamics simulations. We observe rapid conformational transitions between tail bound and unbound states, and characterize kinetic and thermodynamic properties of histone tail-DNA interactions. Different histone types exhibit distinct binding modes to specific DNA regions. Using a comprehensive set of experimental nucleosome complexes, we find that the majority of them target mutually exclusive regions with histone tails on nucleosomal/linker DNA around the super-helical locations ± 1, ± 2, and ± 7, and histone tails H3 and H4 contribute most to this process. These findings are explained within competitive binding and tail displacement models. Finally, we demonstrate the crosstalk between different histone tail post-translational modifications and mutations; those which change charge, suppress tail-DNA interactions and enhance histone tail dynamics and DNA accessibility.

The intrinsic disorder of histone tails poses challenges in their characterization. Here the authors apply extensive molecular dynamics simulations of the full nucleosome to show reversible binding to DNA with specific binding modes of different types of histone tails, where charge-altering modifications suppress tail-DNA interactions and may boost interactions between nucleosomes and nucleosome-binding proteins.

In silico saturation mutagenesis of cancer genes

Despite the existence of good catalogues of cancer genes^1,2, identifying the specific mutations of those genes that drive tumorigenesis across tumour types is still a largely unsolved problem. As a result, most mutations identified in cancer genes across tumours are of unknown significance to tumorigenesis³. We propose that the mutations observed in thousands of tumours-natural experiments testing their oncogenic potential replicated across individuals and tissues-can be exploited to solve this problem. From these mutations, features that describe the mechanism of tumorigenesis of each cancer gene and tissue may be computed and used to build machine learning models that encapsulate these mechanisms. Here we demonstrate the feasibility of this solution by building and validating 185 gene-tissue-specific machine learning models that outperform experimental saturation mutagenesis in the identification of  driver and passenger mutations. The models and their assessment of each mutation are designed to be interpretable, thus avoiding a black-box prediction device. Using these models, we outline the blueprints of potential driver mutations in cancer genes, and demonstrate the role of mutation probability in shaping the landscape of observed driver mutations. These blueprints will support the interpretation of newly sequenced tumours in patients and the study of the mechanisms of tumorigenesis of cancer genes across tissues.

Mutational Signatures: From Methods to Mechanisms

Mutations are the driving force of evolution, yet they underlie many diseases, in particular, cancer. They are thought to arise from a combination of stochastic errors in DNA processing, naturally occurring DNA damage (e.g., the spontaneous deamination of methylated CpG sites), replication errors, and dysregulation of DNA repair mechanisms. High-throughput sequencing has made it possible to generate large datasets to study mutational processes in health and disease. Since the emergence of the first mutational process studies in 2012, this field is gaining increasing attention and has already accumulated a host of computational approaches and biomedical applications.

Learning mutational signatures and their multidimensional genomic properties with TensorSignatures

We present TensorSignatures, an algorithm to learn mutational signatures jointly across different variant categories and their genomic localisation and properties. The analysis of 2778 primary and 3824 metastatic cancer genomes of the PCAWG consortium and the HMF cohort shows that all signatures operate dynamically in response to genomic states. The analysis pins differential spectra of UV mutagenesis found in active and inactive chromatin to global genome nucleotide excision repair. TensorSignatures accurately characterises transcription-associated mutagenesis in 7 different cancer types. The algorithm also extracts distinct signatures of replication- and double strand break repair-driven mutagenesis by APOBEC3A and 3B with differential numbers and length of mutation clusters. Finally, TensorSignatures reproduces a signature of somatic hypermutation generating highly clustered variants at transcription start sites of active genes in lymphoid leukaemia, distinct from a general and less clustered signature of Polη-driven translesion synthesis found in a broad range of cancer types. In summary, TensorSignatures elucidates complex mutational footprints by characterising their underlying processes with respect to a multitude of genomic variables.

Coupling between Sequence-Mediated Nucleosome Organization and Genome Evolution

The nucleosome is a major modulator of DNA accessibility to other cellular factors. Nucleosome positioning has a critical importance in regulating cell processes such as transcription, replication, recombination or DNA repair. The DNA sequence has an influence on the position of nucleosomes on genomes, although other factors are also implicated, such as ATP-dependent remodelers or competition of the nucleosome with DNA binding proteins. Different sequence motifs can promote or inhibit the nucleosome formation, thus influencing the accessibility to the DNA. Sequence-encoded nucleosome positioning having functional consequences on cell processes can then be selected or counter-selected during evolution. We review the interplay between sequence evolution and nucleosome positioning evolution. We first focus on the different ways to encode nucleosome positions in the DNA sequence, and to which extent these mechanisms are responsible of genome-wide nucleosome positioning in vivo. Then, we discuss the findings about selection of sequences for their nucleosomal properties. Finally, we illustrate how the nucleosome can directly influence sequence evolution through its interactions with DNA damage and repair mechanisms. This review aims to provide an overview of the mutual influence of sequence evolution and nucleosome positioning evolution, possibly leading to complex evolutionary dynamics.

The impact of epigenetic information on genome evolution

Epigenetic information affects gene function by interacting with chromatin, while not changing the DNA sequence itself. However, it has become apparent that the interactions between epigenetic information and chromatin can, in fact, indirectly lead to DNA mutations and ultimately influence genome evolution. This review evaluates the ways in which epigenetic information affects genome sequence and evolution. We discuss how DNA methylation has strong and pervasive effects on DNA sequence evolution in eukaryotic organisms. We also review how the physical interactions arising from the connections between histone proteins and DNA affect DNA mutation and repair. We then discuss how a variety of epigenetic mechanisms exert substantial effects on genome evolution by suppressing the movement of transposable elements. Finally, we examine how genome expansion through gene duplication is also partially controlled by epigenetic information. Overall, we conclude that epigenetic information has widespread indirect effects on DNA sequences in eukaryotes and represents a potent cause and constraint of genome evolution. This article is part of the theme issue 'How does epigenetics influence the course of evolution?'

Molecular Mechanisms of Oncogenesis through the Lens of Nucleosomes and Histones

At the cellular level, cancer is the disease of both the genome and the epigenome, and the interplay between genetic mutations and epigenetic states may occur at the level of elementary chromatin units, the nucleosomes. They are formed by a segment of DNA wrapped around an octamer of histone proteins. In this review, we survey various mechanisms of cancer etiology and progression mediated by histones and nucleosomes. In particular, we discuss the effects of mutations in histones, changes in their expression and slicing on epigenetic dysregulation and carcinogenesis. The links between cancer phenotypes and differential expression of histone variants and isoforms are summarized. Finally, we discourse the geometric and steric effects of DNA compaction in nucleosomes on DNA mutation rate, interactions with transcription factors, including pioneer transcription factors, and prospects of cancer cells' genome and epigenome editing.

Epigenetic Biomarkers for the Diagnosis and Treatment of Liver Disease

Research in the last decades has demonstrated the relevance of epigenetics in controlling gene expression to maintain cell homeostasis, and the important role played by epigenome alterations in disease development. Moreover, the reversibility of epigenetic marks can be harnessed as a therapeutic strategy, and epigenetic marks can be used as diagnosis biomarkers. Epigenetic alterations in DNA methylation, histone post-translational modifications (PTMs), and non-coding RNA (ncRNA) expression have been associated with the process of hepatocarcinogenesis. Here, we summarize epigenetic alterations involved in the pathogenesis of chronic liver disease (CLD), particularly focusing on DNA methylation. We also discuss their utility as epigenetic biomarkers in liquid biopsy for the diagnosis and prognosis of hepatocellular carcinoma (HCC). Finally, we discuss the potential of epigenetic therapeutic strategies for HCC treatment.

Exons and introns exhibit transcriptional strand asymmetry of dinucleotide distribution, damage formation and DNA repair

Recent cancer sequencing efforts have uncovered asymmetry in DNA damage induced mutagenesis between the transcribed and non-transcribed strands of genes. Here, we investigate the major type of damage induced by ultraviolet (UV) radiation, the cyclobutane pyrimidine dimers (CPDs), which are formed primarily in TT dinucleotides. We reveal that a transcriptional asymmetry already exists at the level of TT dinucleotide frequency and therefore also in CPD damage formation. This asymmetry is conserved in vertebrates and invertebrates and is completely reversed between introns and exons. We show the asymmetry in introns is linked to the transcription process itself, and is also found in enhancer elements. In contrast, the asymmetry in exons is not correlated to transcription, and is associated with codon usage preferences. Reanalysis of nucleotide excision repair, normalizing repair to the underlying TT frequencies, we show repair of CPDs is more efficient in exons compared to introns, contributing to the maintenance and integrity of coding regions. Our results highlight the importance of considering the primary sequence of the DNA in determining DNA damage sensitivity and mutagenic potential.

Variable interplay of UV-induced DNA damage and repair at transcription factor binding sites

An abnormally high rate of UV-light related mutations appears at transcription factor binding sites (TFBS) across melanomas. The binding of transcription factors (TFs) to the DNA impairs the repair of UV-induced lesions and certain TFs have been shown to increase the rate of generation of these lesions at their binding sites. However, the precise contribution of these two elements to the increase in mutation rate at TFBS in these malignant cells is not understood. Here, exploiting nucleotide-resolution data, we computed the rate of formation and repair of UV-lesions within the binding sites of TFs of different families. We observed, at certain dipyrimidine positions within the binding site of TFs in the Tryptophan Cluster family, an increased rate of formation of UV-induced lesions, corroborating previous studies. Nevertheless, across most families of TFs, the observed increased mutation rate within the entire DNA region covered by the protein results from the decreased repair efficiency. While the rate of mutations across all TFBS does not agree with the amount of UV-induced lesions observed immediately after UV exposure, it strongly agrees with that observed after 48 h. This corroborates the determinant role of the impaired repair in the observed increase of mutation rate.

Chromatin and other obstacles to base excision repair: potential roles in carcinogenesis

DNA is comprised of chemically reactive nucleobases that exist under a constant barrage from damaging agents. Failure to repair chemical modifications to these nucleobases can result in mutations that can cause various diseases, including cancer. Fortunately, the base excision repair (BER) pathway can repair modified nucleobases and prevent these deleterious mutations. However, this pathway can be hindered through several mechanisms. For instance, mutations to the enzymes in the BER pathway have been identified in cancers. Biochemical characterisation of these mutants has elucidated various mechanisms that inhibit their activity. Furthermore, the packaging of DNA into chromatin poses another obstacle to the ability of BER enzymes to function properly. Investigations of BER in the base unit of chromatin, the nucleosome core particle (NCP), have revealed that the NCP acts as a complex substrate for BER enzymes. The constituent proteins of the NCP, the histones, also have variants that can further impact the structure of the NCP and may modulate access of enzymes to the packaged DNA. These histone variants have also displayed significant clinical effects both in carcinogenesis and patient prognosis. This review focuses on the underlying molecular mechanisms that present obstacles to BER and the relationship of these obstacles to cancer. In addition, several chemotherapeutics induce DNA damage that can be repaired by the BER pathway and understanding obstacles to BER can inform how resistance and/or sensitivity to these therapies may occur. With the understanding of these molecular mechanisms, current chemotherapeutic treatment regiments may be improved, and future therapies developed.

Nucleosome Core Particles Lacking H2B or H3 Tails Are Altered Structurally and Have Differential Base Excision Repair Fingerprints

A recently discovered post-translational modification of histone proteins is the irreversible proteolytic clipping of the histone N-terminal tail domains. This modification is involved in the regulation of various biological processes, including the DNA damage response. In this work, we used chemical footprinting to characterize the structural alterations to nucleosome core particles (NCPs) that result from a lack of a histone H2B or H3 tail. We also examine the influence of these histone tails on excision of the mutagenic lesion 1,N⁶-ethenoadenine (εA) by the repair enzyme alkyladenine DNA glycosylase. We found that the absence of the H2B or H3 tail results in altered DNA periodicity relative to that of native NCPs. We correlated these structural alterations to εA excision by utilizing a global analysis of 21 εA sites in NCPs and unincorporated duplex DNA. In comparison to native NCPs, there is enhanced excision of εA in tailless H2B NCPs in regions that undergo DNA unwrapping. This enhanced excision is not observed for tailless H3 NCPs; rather, excision is inhibited in more static areas of the NCP not prone to unwrapping. Our results support in vivo observations of alkylation damage profiles and the potential role of tail clipping as a mechanism for overcoming physical obstructions caused by packaging in NCPs but also reveal the potential inhibition of repair by tail clipping in some locations. Taken together, these results further our understanding of how base excision repair can be facilitated or diminished by histone tail removal and contribute to our understanding of the underlying mechanism that leads to mutational hot spots.

Cancer-associated mutations in endometriosis: shedding light on the pathogenesis and pathophysiology

BACKGROUND

Endometriosis is a benign gynaecological disease. Thus, it came as a complete surprise when it was reported recently that the majority of deep endometriosis lesions harbour somatic mutations and a sizeable portion of them contain known cancer-associated mutations (CAMs). Four more studies have since been published, all demonstrating the existence of CAMs in different subtypes of endometriosis. While the field is still evolving, the confirmation of CAMs has raised many questions that were previously overlooked.

OBJECTIVE AND RATIONALE

A comprehensive overview of CAMs in endometriosis has been produced. In addition, with the recently emerged understanding of the natural history of endometriotic lesions as well as CAMs in normal and apparently healthy tissues, this review attempts to address the following questions: Why has there been such a wild discrepancy in reported mutation frequencies? Why does ectopic endometrium have a higher mutation rate than that of eutopic endometrium? Would the presence of CAMs in endometriotic lesions increase the risk of cancer to the bearers? Why do endometriotic epithelial cells have much higher mutation frequencies than their stromal counterpart? What clinical implications, if any, do the CAMs have for the bearers? Do these CAMs tell us anything about the pathogenesis and/or pathophysiology of endometriosis?

SEARCH METHODS

The PubMed database was searched, from its inception to September 2019, for all papers in English using the term 'endometriosis and CAM', 'endometriosis and cancer-driver mutation', 'somatic mutations', 'fibrosis', 'fibrosis and epigenetic', 'CAMs and tumorigenesis', 'somatic mutation and normal tissues', 'oestrogen receptor and fibrosis', 'oxidative stress and fibrosis', 'ARID1A mutation', and 'Kirsten rat sarcoma mutation and therapeutics'. All retrieved papers were read and, when relevant, incorporated into the review results.

OUTCOMES

Seven papers that identified CAMs in endometriosis using various sequencing methods were retrieved, and their results were somewhat different. Yet, it is apparent that those using microdissection techniques and more accurate sequencing methods found more CAMs, echoing recent discoveries that apparently healthy tissues also harbour CAMs as a result of the replicative aging process. Hence endometriotic lesions, irrespective of subtype, if left intact, would generate CAMs as part of replicative aging, oxidative stress and perhaps other factors yet to be identified and, in some rare cases, develop cancer. The published data still are unable to paint a clear picture on pathogenesis of endometriosis. However, since endometriotic epithelial cells have a higher turnover than their stromal counterpart due to cyclic bleeding, and since the endometriotic stromal component can be formed by refresh influx of mesenchymal cells through epithelial-mesenchymal transition, endothelial-mesenchymal transition, mesothelial-mesenchymal transition and other processes as well as recruitment of bone-marrow-derived stem cells and outflow due to smooth muscle metaplasia, endometriotic epithelial cells have much higher mutation frequencies than their stromal counterpart. The epithelial and stromal cellular components develop in a dependent and co-evolving manner. Genes involved in CAMs are likely to be active players in lesional fibrogenesis, and hyperestrogenism and oxidative stress are likely drivers of both CAMs and fibrogenesis. Finally, endometriotic lesions harbouring CAMs would conceivably be more refractory to medical treatment, due, in no small part, to their high fibrotic content and reduced vascularity and cellularity.

WIDER IMPLICATIONS

The accumulating data on CAMs in endometriosis have shed new light on the pathogenesis and pathophysiology of endometriosis. They also suggest new challenges in management. The distinct yet co-evolving developmental trajectories of endometriotic stroma and epithelium underscore the importance of the lesional microenvironment and ever-changing cellular identity. Mutational profiling of normal endometrium from women of different ages and reproductive history is needed in order to gain a deeper understanding of the pathogenesis. Moreover, one area that has conspicuously received scant attention is the epigenetic landscape of ectopic, eutopic and normal endometrium.

DNA mismatches reveal conformational penalties in protein-DNA recognition

Transcription factors recognize specific genomic sequences to regulate complex gene-expression programs. Although it is well-established that transcription factors bind to specific DNA sequences using a combination of base readout and shape recognition, some fundamental aspects of protein-DNA binding remain poorly understood1,2. Many DNA-binding proteins induce changes in the structure of the DNA outside the intrinsic B-DNA envelope. However, how the energetic cost that is associated with distorting the DNA contributes to recognition has proven difficult to study, because the distorted DNA exists in low abundance in the unbound ensemble3-9. Here we use a high-throughput assay that we term SaMBA (saturation mismatch-binding assay) to investigate the role of DNA conformational penalties in transcription factor-DNA recognition. In SaMBA, mismatched base pairs are introduced to pre-induce structural distortions in the DNA that are much larger than those induced by changes in the Watson-Crick sequence. Notably, approximately 10% of mismatches increased transcription factor binding, and for each of the 22 transcription factors that were examined, at least one mismatch was found that increased the binding affinity. Mismatches also converted non-specific sites into high-affinity sites, and high-affinity sites into 'super sites' that exhibit stronger affinity than any known canonical binding site. Determination of high-resolution X-ray structures, combined with nuclear magnetic resonance measurements and structural analyses, showed that many of the DNA mismatches that increase binding induce distortions that are similar to those induced by protein binding-thus prepaying some of the energetic cost incurred from deforming the DNA. Our work indicates that conformational penalties are a major determinant of protein-DNA recognition, and reveals mechanisms by which mismatches can recruit transcription factors and thus modulate replication and repair activities in the cell10,11.

Distinctive regulatory architectures of germline-active and somatic genes in C. elegans

RNA profiling has provided increasingly detailed knowledge of gene expression patterns, yet the different regulatory architectures that drive them are not well understood. To address this, we profiled and compared transcriptional and regulatory element activities across five tissues of Caenorhabditis elegans, covering ∼90% of cells. We find that the majority of promoters and enhancers have tissue-specific accessibility, and we discover regulatory grammars associated with ubiquitous, germline, and somatic tissue–specific gene expression patterns. In addition, we find that germline-active and soma-specific promoters have distinct features. Germline-active promoters have well-positioned +1 and −1 nucleosomes associated with a periodic 10-bp WW signal (W = A/T). Somatic tissue–specific promoters lack positioned nucleosomes and this signal, have wide nucleosome-depleted regions, and are more enriched for core promoter elements, which largely differ between tissues. We observe the 10-bp periodic WW signal at ubiquitous promoters in other animals, suggesting it is an ancient conserved signal. Our results show fundamental differences in regulatory architectures of germline and somatic tissue–specific genes, uncover regulatory rules for generating diverse gene expression patterns, and provide a tissue-specific resource for future studies.

A compendium of mutational cancer driver genes

A fundamental goal in cancer research is to understand the mechanisms of cell transformation. This is key to developing more efficient cancer detection methods and therapeutic approaches. One milestone towards this objective is the identification of all the genes with mutations capable of driving tumours. Since the 1970s, the list of cancer genes has been growing steadily. Because cancer driver genes are under positive selection in tumorigenesis, their observed patterns of somatic mutations across tumours in a cohort deviate from those expected from neutral mutagenesis. These deviations, which constitute signals of positive selection, may be detected by carefully designed bioinformatics methods, which have become the state of the art in the identification of driver genes. A systematic approach combining several of these signals could lead to a compendium of mutational cancer genes. In this Review, we present the Integrative OncoGenomics (IntOGen) pipeline, an implementation of such an approach to obtain the compendium of mutational cancer drivers. Its application to somatic mutations of more than 28,000 tumours of 66 cancer types reveals 568 cancer genes and points towards their mechanisms of tumorigenesis. The application of this approach to the ever-growing datasets of somatic tumour mutations will support the continuous refinement of our knowledge of the genetic basis of cancer.

MASS SPECTROMETRY-BASED PERSONALIZED DRUG THERAPY

Personalized drug therapy aims to provide tailored treatment for individual patient. Mass spectrometry (MS) is revolutionarily involved in this area because MS is a rapid, customizable, cost-effective, and easy to be used high-throughput method with high sensitivity, specificity, and accuracy. It is driving the formation of a new field, MS-based personalized drug therapy, which currently mainly includes five subfields: therapeutic drug monitoring (TDM), pharmacogenomics (PGx), pharmacomicrobiomics, pharmacoepigenomics, and immunopeptidomics. Gas chromatography-MS (GC-MS) and liquid chromatography-MS (LC-MS) are considered as the gold standard for TDM, which can be used to optimize drug dosage. Matrix-assisted laser desorption ionization-time of flight-MS (MALDI-TOF-MS) significantly improves the capability of detecting biomacromolecule, and largely promotes the application of MS in PGx. It is becoming an indispensable tool for genotyping, which is used to discover and validate genetic biomarkers. In addition, MALDI-TOF-MS also plays important roles in identity of human microbiome whose diversity can explain interindividual differences of drug response. Pharmacoepigenetics is to study the role of epigenetic factors in individualized drug treatment. MS can be used to discover and validate pharmacoepigenetic markers (DNA methylation, histone modification, and noncoding RNA). For the emerging cancer immunotherapy, personalized cancer vaccine has effective immunotherapeutic activity in the clinic. MS-based immunopeptidomics can effectively discover and screen neoantigens. This article systematically reviewed MS-based personalized drug therapy in the above mentioned five subfields. © 2020 John Wiley & Sons Ltd. Mass Spec Rev.