ETS transcription factors induce a unique UV damage signature that drives recurrent mutagenesis in melanoma

YBX1 Modulates 8-Oxoguanine Recognition and Repair in DNA

8-Oxoguanine (8-oxoG) is not only a biomarker of oxidative DNA damage but also an epigenetic-like regulator in mammalian cells. The identification and characterization of 8-oxoG-binding proteins would be crucial for further understanding the biological consequences of 8-oxoG. Here, we identified human Y-box-binding protein 1 (YBX1) as a novel binding protein for 8-oxoG modification in DNA by using a quantitative proteomic approach. Moreover, we found that the deficiency of YBX1 can substantially decrease the cellular sensitivity to oxidative stress and facilitate the repair of 8-oxoG embedded in DNA. These findings provided new insight into the biological significance of the functional interplay between YBX1 and 8-oxoG modification in DNA.

The epigenetic landscape shapes smoking-induced mutagenesis by modulating DNA damage susceptibility and repair efficiency

Lung cancer sequencing efforts have uncovered mutational signatures that are attributed to exposure to the cigarette smoke carcinogen benzo[a]pyrene. Benzo[a]pyrene metabolizes in cells to benzo[a]pyrene diol epoxide (BPDE) and reacts with guanine nucleotides to form bulky BPDE adducts. These DNA adducts block transcription and replication, compromising cell function and survival, and are repaired in human cells by the nucleotide excision repair pathway. Here, we applied high-resolution genomic assays to measure BPDE-induced damage formation and mutagenesis in human cells. We integrated the new damage and mutagenesis data with previous repair, DNA methylation, RNA expression, DNA replication, and chromatin component measurements in the same cell lines, along with lung cancer mutagenesis data. BPDE damage formation is significantly enhanced by DNA methylation and in accessible chromatin regions, including transcribed and early-replicating regions. Binding of transcription factors is associated primarily with reduced, but also enhanced damage formation, depending on the factor. While DNA methylation does not appear to influence repair efficiency, this repair was significantly elevated in accessible chromatin regions, which accumulated fewer mutations. Thus, when damage and repair drive mutagenesis in opposing directions, the final mutational patterns appear to be dictated by the efficiency of repair rather than the frequency of underlying damages.

Transcription factors, nucleotide excision repair, and cancer: A review of molecular interplay

Bulky DNA adducts are mostly formed by external factors such as UV irradiation, smoking or treatment with DNA crosslinking agents. If such DNA adducts are not removed by nucleotide excision repair, they can lead to formation of driver mutations that contribute to cancer formation. Transcription factors (TFs) may critically affect both DNA adduct formation and repair efficiency at the binding site to DNA. For example, "hotspot" mutations in melanoma coincide with UV-induced accumulated cyclobutane pyrimidine dimer (CPD) adducts and/or inhibited repair at the binding sites of some TFs. Similarly, anticancer treatment with DNA cross-linkers may additionally generate DNA adducts leading to secondary mutations and the formation of malignant subclones. In addition, some TFs are overexpressed in response to UV irradiation or chemotherapeutic treatment, activating oncogenic and anti-oncogenic pathways independently of nucleotide excision repair itself. This review focuses on the interplay between TFs and nucleotide excision repair during cancer development and progression.

Contributing factors to the oxidation-induced mutational landscape in human cells

8-oxoguanine (8-oxoG) is a common oxidative DNA lesion that causes G > T substitutions. Determinants of local and regional differences in 8-oxoG-induced mutability across genomes are currently unknown. Here, we show DNA oxidation induces G > T substitutions and insertion/deletion (INDEL) mutations in human cells and cancers. Potassium bromate (KBrO3)-induced 8-oxoGs occur with similar sequence preferences as their derived substitutions, indicating that the reactivity of specific oxidants dictates mutation sequence specificity. While 8-oxoG occurs uniformly across chromatin, 8-oxoG-induced mutations are elevated in compact genomic regions, within nucleosomes, and at inward facing guanines within strongly positioned nucleosomes. Cryo-electron microscopy structures of OGG1-nucleosome complexes indicate that these effects originate from OGG1's ability to flip outward positioned 8-oxoG lesions into the catalytic pocket while inward facing lesions are occluded by the histone octamer. Mutation spectra from human cells with DNA repair deficiencies reveals contributions of a DNA repair network limiting 8-oxoG mutagenesis, where OGG1- and MUTYH-mediated base excision repair is supplemented by the replication-associated factors Pol η and HMCES. Transcriptional asymmetry of KBrO3-induced mutations in OGG1- and Pol η-deficient cells also demonstrates transcription-coupled repair can prevent 8-oxoG-induced mutation. Thus, oxidant chemistry, chromatin structures, and DNA repair processes combine to dictate the oxidative mutational landscape in human genomes.

Methods and applications of genome-wide profiling of DNA damage and rare mutations

DNA damage is a threat to genome integrity and can be a cause of many human diseases, owing to either changes in the chemical structure of DNA or conversion of the damage into a mutation, that is, a permanent change in DNA sequence. Determining the exact positions of DNA damage and ensuing mutations in the genome are important for identifying mechanisms of disease aetiology when characteristic mutations are prevalent and probably causative in a particular disease. However, this approach is challenging particularly when levels of DNA damage are low, for example, as a result of chronic exposure to environmental agents or certain endogenous processes, such as the generation of reactive oxygen species. Over the past few years, a comprehensive toolbox of genome-wide methods has been developed for the detection of DNA damage and rare mutations at single-nucleotide resolution in mammalian cells. Here, we review and compare these methods, describe their current applications and discuss future research questions that can now be addressed.

Mechanistic basis of atypical TERT promoter mutations

Non-coding mutations in the TERT promoter (TERTp), typically at one of two bases -124 and -146 bp upstream of the start codon, are among the most prevalent driver mutations in human cancer. Several additional recurrent TERTp mutations have been reported but their functions and origins remain largely unexplained. Here, we show that atypical TERTp mutations arise secondary to canonical TERTp mutations in a two-step process. Canonical TERTp mutations create de novo binding sites for ETS family transcription factors that induce favourable conditions for DNA damage formation by UV light, thus creating a hotspot effect but only after a first mutational hit. In agreement, atypical TERTp mutations co-occur with canonical driver mutations in large cancer cohorts and arise subclonally specifically on the TERTp driver mutant chromosome homolog of melanoma cells treated with UV light in vitro. Our study gives an in-depth view of TERTp mutations in cancer and provides a mechanistic explanation for atypical TERTp mutations.

Nucleotide excision repair of aflatoxin-induced DNA damage within the 3D human genome organization

Aflatoxin B1 (AFB1), a potent mycotoxin, is one of the environmental risk factors that cause liver cancer. In the liver, the bioactivated AFB1 intercalates into the DNA double helix to form a bulky DNA adduct which will lead to mutation if left unrepaired. Here, we adapted the tXR-seq method to measure the nucleotide excision repair of AFB1-induced DNA adducts at single-nucleotide resolution on a genome-wide scale, and compared it with repair data obtained from conventional UV-damage XR-seq. Our results showed that transcription-coupled repair plays a major role in the damage removal process. We further analyzed the distribution of nucleotide excision repair sites for AFB1-induced DNA adducts within the 3D human genome organization. Our analysis revealed a heterogeneous AFB1-dG repair across four different organization levels, including chromosome territories, A/B compartments, TADs, and chromatin loops. We found that chromosomes positioned closer to the nuclear center and regions within A compartments have higher levels of nucleotide excision repair. Notably, we observed high repair activity around both TAD boundaries and loop anchors. These findings provide insights into the complex interplay between AFB1-induced DNA damage repair, transcription, and 3D genome organization, shedding light on the mechanisms underlying AFB1-induced mutagenesis.

Pervasive Induction of Regulatory Mutation Microclones in Sun-exposed Skin

Carcinogen-induced mutations are thought near-random, with rare cancer-driver mutations underlying clonal expansion. Using high-fidelity Duplex Sequencing to reach a mutation frequency sensitivity of 4×10 -9 per nt, we report that sun exposure creates pervasive mutations at sites with ∼100-fold UV-sensitivity in RNA-processing gene promoters - cyclobutane pyrimidine dimer (CPD) hyperhotspots - and these mutations have a mini-driver clonal expansion phenotype. Numerically, human skin harbored 10-fold more genuine mutations than previously reported, with neonatal skin containing 90,000 per cell; UV signature mutations increased 8,000-fold in sun-exposed skin, averaging 3×10 -5 per nt. Clonal expansion by neutral drift or passenger formation was nil. Tumor suppressor gene hotspots reached variant allele frequency 0.1-10% via 30-3,000 fold clonal expansion, in occasional biopsies. CPD hyperhotspots reached those frequencies in every biopsy, with modest clonal expansion. In vitro, tumor hotspot mutations arose occasionally over weeks of chronic low-dose exposure, whereas CPD hyperhotspot mutations arose in days at 1000-fold higher frequencies, growing exponentially. UV targeted mini-drivers in every skin cell.

ASH1L guards cis-regulatory elements against cyclobutane pyrimidine dimer induction

The histone methyltransferase ASH1L, first discovered for its role in transcription, has been shown to accelerate the removal of ultraviolet (UV) light-induced cyclobutane pyrimidine dimers (CPDs) by nucleotide excision repair. Previous reports demonstrated that CPD excision is most efficient at transcriptional regulatory elements, including enhancers, relative to other genomic sites. Therefore, we analyzed DNA damage maps in ASH1L-proficient and ASH1L-deficient cells to understand how ASH1L controls enhancer stability. This comparison showed that ASH1L protects enhancer sequences against the induction of CPDs besides stimulating repair activity. ASH1L reduces CPD formation at C-containing but not at TT dinucleotides, and no protection occurs against pyrimidine-(6,4)-pyrimidone photoproducts or cisplatin crosslinks. The diminished CPD induction extends to gene promoters but excludes retrotransposons. This guardian role against CPDs in regulatory elements is associated with the presence of H3K4me3 and H3K27ac histone marks, which are known to interact with the PHD and BRD motifs of ASH1L, respectively. Molecular dynamics simulations identified a DNA-binding AT hook of ASH1L that alters the distance and dihedral angle between neighboring C nucleotides to disfavor dimerization. The loss of this protection results in a higher frequency of C->T transitions at enhancers of skin cancers carrying ASH1L mutations compared to ASH1L-intact counterparts.

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.

Deciphering the dynamic code: DNA recognition by transcription factors in the ever-changing genome

Transcription factors (TFs) intricately navigate the vast genomic landscape to locate and bind specific DNA sequences for the regulation of gene expression programs. These interactions occur within a dynamic cellular environment, where both DNA and TF proteins experience continual chemical and structural perturbations, including epigenetic modifications, DNA damage, mechanical stress, and post-translational modifications (PTMs). While many of these factors impact TF-DNA binding interactions, understanding their effects remains challenging and incomplete. This review explores the existing literature on these dynamic changes and their potential impact on TF-DNA interactions.

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.

Genomic analysis of human brain metastases treated with stereotactic radiosurgery reveals unique signature based on treatment failure

Stereotactic radiosurgery (SRS) has been shown to be efficacious for the treatment of limited brain metastasis (BM); however, the effects of SRS on human brain metastases have yet to be studied. We performed genomic analysis on resected brain metastases from patients whose resected lesion was previously treated with SRS. Our analyses demonstrated for the first time that patients possess a distinct genomic signature based on type of treatment failure including local failure, leptomeningeal spread, and radio-necrosis. Examination of the center and peripheral edge of the tumors treated with SRS indicated differential DNA damage distribution and an enrichment for tumor suppressor mutations and DNA damage repair pathways along the peripheral edge. Furthermore, the two clinical modalities used to deliver SRS, LINAC and GK, demonstrated differential effects on the tumor landscape even between controlled primary sites. Our study provides, in human, biological evidence of differential effects of SRS across BM's.

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.

Gene architecture is a determinant of the transcriptional response to bulky DNA damages

Bulky DNA damages block transcription and compromise genome integrity and function. The cellular response to these damages includes global transcription shutdown. Still, active transcription is necessary for transcription-coupled repair and for induction of damage-response genes. To uncover common features of a general bulky DNA damage response, and to identify response-related transcripts that are expressed despite damage, we performed a systematic RNA-seq study comparing the transcriptional response to three independent damage-inducing agents: UV, the chemotherapy cisplatin, and benzo[a]pyrene, a component of cigarette smoke. Reduction in gene expression after damage was associated with higher damage rates, longer gene length, and low GC content. We identified genes with relatively higher expression after all three damage treatments, including NR4A2, a potential novel damage-response transcription factor. Up-regulated genes exhibit higher exon content that is associated with preferential repair, which could enable rapid damage removal and transcription restoration. The attenuated response to BPDE highlights that not all bulky damages elicit the same response. These findings frame gene architecture as a major determinant of the transcriptional response that is hardwired into the human genome.

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.

High UV damage and low repair, but not cytosine deamination, stimulate mutation hotspots at ETS binding sites in melanoma

Noncoding mutation hotspots have been identified in melanoma and many of them occur at the binding sites of E26 transformation-specific (ETS) proteins; however, their formation mechanism and functional impacts are not fully understood. Here, we used UV (Ultraviolet) damage sequencing data and analyzed cyclobutane pyrimidine dimer (CPD) formation, DNA repair, and CPD deamination in human cells at single-nucleotide resolution. Our data show prominent CPD hotspots immediately after UV irradiation at ETS binding sites, particularly at sites with a conserved TTCCGG motif, which correlate with mutation hotspots identified in cutaneous melanoma. Additionally, CPDs are repaired slower at ETS binding sites than in flanking DNA. Cytosine deamination in CPDs to uracil is suggested as an important step for UV mutagenesis. However, we found that CPD deamination is significantly suppressed at ETS binding sites, particularly for the CPD hotspot on the 5' side of the ETS motif, arguing against a role for CPD deamination in promoting ETS-associated UV mutations. Finally, we analyzed a subset of frequently mutated promoters, including the ribosomal protein genes RPL13A and RPS20, and found that mutations in the ETS motif can significantly reduce the promoter activity. Thus, our data identify high UV damage and low repair, but not CPD deamination, as the main mechanism for ETS-associated mutations in melanoma and uncover important roles of often-overlooked mutation hotspots in perturbing gene transcription.

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.

A mutation rate model at the basepair resolution identifies the mutagenic effect of polymerase III transcription

De novo mutations occur at substantially different rates depending on genomic location, sequence context and DNA strand. The success of methods to estimate selection intensity, infer demographic history and map rare disease genes, depends strongly on assumptions about the local mutation rate. Here we present Roulette, a genome-wide mutation rate model at basepair resolution that incorporates known determinants of local mutation rate. Roulette is shown to be more accurate than existing models. We use Roulette to refine the estimates of population growth within Europe by incorporating the full range of human mutation rates. The analysis of significant deviations from the model predictions revealed a tenfold increase in mutation rate in nearly all genes transcribed by polymerase III (Pol III), suggesting a new mutagenic mechanism. We also detected an elevated mutation rate within transcription factor binding sites restricted to sites actively used in testis and residing in promoters.

Functional analysis of recurrent CDC20 promoter variants in human melanoma

Small nucleotide variants in non-coding regions of the genome can alter transcriptional regulation, leading to changes in gene expression which can activate oncogenic gene regulatory networks. Melanoma is heavily burdened by non-coding variants, representing over 99% of total genetic variation, including the well-characterized TERT promoter mutation. However, the compendium of regulatory non-coding variants is likely still functionally under-characterized. We developed a pipeline to identify hotspots, i.e. recurrently mutated regions, in melanoma containing putatively functional non-coding somatic variants that are located within predicted melanoma-specific regulatory regions. We identified hundreds of statistically significant hotspots, including the hotspot containing the TERT promoter variants, and focused on a hotspot in the promoter of CDC20. We found that variants in the promoter of CDC20, which putatively disrupt an ETS motif, lead to lower transcriptional activity in reporter assays. Using CRISPR/Cas9, we generated an indel in the CDC20 promoter in human A375 melanoma cell lines and observed decreased expression of CDC20, changes in migration capabilities, increased growth of xenografts, and an altered transcriptional state previously associated with a more proliferative and less migratory state. Overall, our analysis prioritized several recurrent functional non-coding variants that, through downregulation of CDC20, led to perturbation of key melanoma phenotypes.

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.

Coupling of conformation and CPD damage in nucleosomal DNA

UV-light can cause photodimerization and hence damages in DNA. Most frequent are cyclobutane pyrimidine dimer (CPD) damages, which predominantly form at TpT (thymine-thymine) steps. It is well known that CPD damage probability is different for single-stranded or double stranded DNA and depends on the sequence context. However, DNA deformation due to packing in nucleosomes can also influence CPD formation. Quantum mechanical calculations and Molecular Dynamics simulations indicate little CPD damage probability for DNA's equilibrium structure. We find that DNA needs to be deformed in a specific way to allow the HOMO → LUMO transition required for CPD damage formation. The simulation studies further show that the periodic CPD damage patterns measured in chromosomes and nucleosomes can be directly explained by the periodic deformation pattern of the DNA in the nucleosome complex. It supports previous findings on characteristic deformation patterns found in experimental nucleosome structures that relate to CPD damage formation. The result may have important implications for our understanding of UV-induced DNA mutations in human cancers.

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.

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.

Base-resolution UV footprinting by sequencing reveals distinctive damage signatures for DNA-binding proteins

Decades ago, it was shown that proteins binding to DNA can quantitatively alter the formation of DNA damage by UV light. This established the principle of UV footprinting for non-intrusive study of protein-DNA contacts in living cells, albeit at limited scale and precision. Here, we perform deep base-resolution quantification of the principal UV damage lesion, the cyclobutane pyrimidine dimer (CPD), at select human promoter regions using targeted CPD sequencing. Several transcription factors exhibited distinctive and repeatable damage signatures indicative of site occupancy, involving strong (up to 17-fold) position-specific elevations and reductions in CPD formation frequency relative to naked DNA. Positive damage modulation at some ETS transcription factor binding sites coincided at base level with melanoma somatic mutation hotspots. Our work provides proof of concept for the study of protein-DNA interactions at individual loci using light and sequencing, and reveals widespread and potent modulation of UV damage in regulatory regions.

Detecting recurrent passenger mutations in melanoma by targeted UV damage sequencing

Sequencing of melanomas has identified hundreds of recurrent mutations in both coding and non-coding DNA. These include a number of well-characterized oncogenic driver mutations, such as coding mutations in the BRAF and NRAS oncogenes, and non-coding mutations in the promoter of telomerase reverse transcriptase (TERT). However, the molecular etiology and significance of most of these mutations is unknown. Here, we use a new method known as CPD-capture-seq to map UV-induced cyclobutane pyrimidine dimers (CPDs) with high sequencing depth and single nucleotide resolution at sites of recurrent mutations in melanoma. Our data reveal that many previously identified drivers and other recurrent mutations in melanoma occur at CPD hotspots in UV-irradiated melanocytes, often associated with an overlapping binding site of an E26 transformation-specific (ETS) transcription factor. In contrast, recurrent mutations in the promoters of a number of known or suspected cancer genes are not associated with elevated CPD levels. Our data indicate that a subset of recurrent protein-coding mutations are also likely caused by ETS-induced CPD hotspots. This analysis indicates that ETS proteins profoundly shape the mutation landscape of melanoma and reveals a method for distinguishing potential driver mutations from passenger mutations whose recurrence is due to elevated UV damage.

Genomic Analysis of Human Brain Metastases Treated with Stereotactic Radiosurgery Under the Phase-II Clinical Trial (NCT03398694) Reveals DNA Damage Repair at the Peripheral Tumor Edge

Stereotactic Radiosurgery (SRS) is one of the leading treatment modalities for oligo brain metastasis (BM), however no comprehensive genomic data assessing the effect of radiation on BM in humans exist. Leveraging a unique opportunity, as part of the clinical trial (NCT03398694), we collected post-SRS, delivered via Gamma-knife or LINAC, tumor samples from core and peripheral-edges of the resected tumor to characterize the genomic effects of overall SRS as well as the SRS delivery modality. Using these rare patient samples, we show that SRS results in significant genomic changes at DNA and RNA levels throughout the tumor. Mutations and expression profiles of peripheral tumor samples indicated interaction with surrounding brain tissue as well as elevated DNA damage repair. Central samples show GSEA enrichment for cellular apoptosis while peripheral samples carried an increase in tumor suppressor mutations. There are significant differences in the transcriptomic profile at the periphery between Gamma-knife vs LINAC.

UV irradiation remodels the specificity landscape of transcription factors

Somatic mutations are highly enriched at transcription factor (TF) binding sites, with the strongest trend being observed for ultraviolet light (UV)-induced mutations in melanomas. One of the main mechanisms proposed for this hypermutation pattern is the inefficient repair of UV lesions within TF-binding sites, caused by competition between TFs bound to these lesions and the DNA repair proteins that must recognize the lesions to initiate repair. However, TF binding to UV-irradiated DNA is poorly characterized, and it is unclear whether TFs maintain specificity for their DNA sites after UV exposure. We developed UV-Bind, a high-throughput approach to investigate the impact of UV irradiation on protein-DNA binding specificity. We applied UV-Bind to ten TFs from eight structural families, and found that UV lesions significantly altered the DNA-binding preferences of all the TFs tested. The main effect was a decrease in binding specificity, but the precise effects and their magnitude differ across factors. Importantly, we found that despite the overall reduction in DNA-binding specificity in the presence of UV lesions, TFs can still compete with repair proteins for lesion recognition, in a manner consistent with their specificity for UV-irradiated DNA. In addition, for a subset of TFs, we identified a surprising but reproducible effect at certain nonconsensus DNA sequences, where UV irradiation leads to a high increase in the level of TF binding. These changes in DNA-binding specificity after UV irradiation, at both consensus and nonconsensus sites, have important implications for the regulatory and mutagenic roles of TFs in the cell.

Genome-wide maps of rare and atypical UV photoproducts reveal distinct patterns of damage formation and mutagenesis in yeast chromatin

Ultraviolet (UV) light induces different classes of mutagenic photoproducts in DNA, namely cyclobutane pyrimidine dimers (CPDs), 6-4 photoproducts (6-4PPs), and atypical thymine-adenine photoproducts (TA-PPs). CPD formation is modulated by nucleosomes and transcription factors (TFs), which has important ramifications for Ultraviolet (UV) mutagenesis. How chromatin affects the formation of 6-4PPs and TA-PPs is unclear. Here, we use UV damage endonuclease-sequencing (UVDE-seq) to map these UV photoproducts across the yeast genome. Our results indicate that nucleosomes, the fundamental building block of chromatin, have opposing effects on photoproduct formation. Nucleosomes induce CPDs and 6-4PPs at outward rotational settings in nucleosomal DNA but suppress TA-PPs at these settings. Our data also indicate that DNA binding by different classes of yeast TFs causes lesion-specific hotspots of 6-4PPs or TA-PPs. For example, DNA binding by the TF Rap1 generally suppresses CPD and 6-4PP formation but induces a TA-PP hotspot. Finally, we show that 6-4PP formation is strongly induced at the binding sites of TATA-binding protein (TBP), which is correlated with higher mutation rates in UV-exposed yeast. These results indicate that the formation of 6-4PPs and TA-PPs is modulated by chromatin differently than CPDs and that this may have important implications for UV mutagenesis.

ETS transcription factors: Multifaceted players from cancer progression to tumor immunity

The E26 transformation specific (ETS) family comprises 28 transcription factors, the majority of which are involved in tumor initiation and development. Serving as a group of functionally heterogeneous gene regulators, ETS factors possess a structurally conserved DNA-binding domain. As one of the most prominent families of transcription factors that control diverse cellular functions, ETS activation is modulated by multiple intracellular signaling pathways and post-translational modifications. Disturbances in ETS activity often lead to abnormal changes in oncogenicity, including cancer cell survival, growth, proliferation, metastasis, genetic instability, cell metabolism, and tumor immunity. This review systematically addresses the basics and advances in studying ETS factors, from their tumor relevance to clinical translational utility, with a particular focus on elucidating the role of ETS family in tumor immunity, aiming to decipher the vital role and clinical potential of regulation of ETS factors in the cancer field.

Widespread perturbation of ETS factor binding sites in cancer

Although >90% of somatic mutations reside in non-coding regions, few have been reported as cancer drivers. To predict driver non-coding variants (NCVs), we present a transcription factor (TF)-aware burden test based on a model of coherent TF function in promoters. We apply this test to NCVs from the Pan-Cancer Analysis of Whole Genomes cohort and predict 2555 driver NCVs in the promoters of 813 genes across 20 cancer types. These genes are enriched in cancer-related gene ontologies, essential genes, and genes associated with cancer prognosis. We find that 765 candidate driver NCVs alter transcriptional activity, 510 lead to differential binding of TF-cofactor regulatory complexes, and that they primarily impact the binding of ETS factors. Finally, we show that different NCVs within a promoter often affect transcriptional activity through shared mechanisms. Our integrated computational and experimental approach shows that cancer NCVs are widespread and that ETS factors are commonly disrupted.

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.

Circle Damage Sequencing for Whole-Genome Analysis of DNA Damage

There are many sources of endogenous and exogenous DNA damage. Damaged bases represent a threat to genome integrity and may interfere with normal cellular processes such as replication and transcription. To understand the specificity and biological consequences of DNA damage, it is essential to employ methods that are sensitive enough to detect damaged DNA bases at the level of single nucleotide resolution and genome-wide. Here we describe in detail a method we developed for this purpose, circle damage sequencing (CD-seq). This method is based on the circularization of genomic DNA that contains damaged bases and conversion of the damaged sites into double-strand breaks using specific DNA repair enzymes. Library sequencing of the opened circles yields the precise positions of the DNA lesions that are present. CD-seq can be adopted to various types of DNA damage as long as a specific cleavage scheme can be designed.

Mutually exclusive genetic interactions and gene essentiality shape the genomic landscape of primary melanoma

Melanoma is a heterogenous malignancy with an unpredictable clinical course. Most patients who present in the clinic are diagnosed with primary melanoma, yet large-scale sequencing efforts have focused primarily on metastatic disease. In this study we sequence-profiled 524 American Joint Committee on Cancer Stage I-III primary tumours. Our analysis of these data reveals recurrent driver mutations, mutually exclusive genetic interactions, where two genes were never or rarely co-mutated, and an absence of co-occurring genetic events. Further, we intersected copy number calls from our primary melanoma data with whole-genome CRISPR screening data to identify the transcription factor interferon regulatory factor 4 (IRF4) as a melanoma-associated dependency. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

UVC radiation intensity dependence of pathogen decontamination rate: semiclassical theory and experiment

A semiclassical (light classical and molecule quantum) model describing the dependence of DNA/RNA dimerization rate as function of the ultraviolet C (UVC) radiation's intensity is proposed. Particularly, a nonlinear model is developed based on the Raman-like processes in quantum optics. The main result of the theory shows that the process of dimerization in the DNA/RNA depends strongly on the UVC light's intensity, thus proving a possible quantum microscopical mechanism of the interaction of UV light with the DNA. To corroborate the theoretical findings, we realize some experiments, by which want to investigate how the inactivation rate of the yeast colonies depends on the intensity of the UVC irradiation. The experimental results evidence a nonlinear decreasing of the residual yeast colonies as a function of the intensity in the irradiation process. The possibilities to optimize the intensity of UVC radiation in the considered decontamination equipment by using metamaterials are studied. The application of such equipment in disinfection of fluids (air, water, droplets, etc.), as well for the SARS-CoV-2-infected aerosols, is discussed.

Cyclobutane Pyrimidine Dimer Hyperhotspots as Sensitive Indicators of Keratinocyte UV Exposure†

The dominant DNA damage generated by UV exposure is the cyclobutane pyrimidine dimer (CPD), which alters skin cell physiology and induces cell death and mutation. Genome-wide nucleotide-resolution analysis of CPDs in melanocytes and fibroblasts has identified "CPD hyperhotspots", pyrimidine-pyrimidine sites hundreds of fold more susceptible to the generation of CPDs than the genomic average. Identifying hyperhotspots in keratinocytes could enable measuring individual past UV exposure in small skin samples and predicting future skin cancer risk. We therefore exposed neonatal human epidermal keratinocytes to narrowband UVB and quantified CPDs using the adductSeq high-throughput DNA sequencing method. Keratinocytes contained thousands of CPD hyperhotspots, with a UVB-sensitivity up to 550 fold greater than the genomic average. As with melanocytes, the most sensitive sites were located in promoter regions at ETS-family transcription factor binding sequence motifs, near RNA processing genes. Moreover, they lay at sequence motifs bound to ETS1 in CpG islands. These genes were specifically upregulated in skin and the CPD hyperhotspots were mutated in a fraction of keratinocyte cancers. Crucially for their biological importance and practical application, CPD hyperhotspot locations and UV-sensitivity ranking demonstrated high reproducibility across experiments and across skin donors. CPD hyperhotspots are therefore sensitive indicators of UV exposure.

MTBP enhances the activation of transcription factor ETS-1 and promotes the proliferation of hepatocellular carcinoma cells

Increasing evidence indicates that the oncoprotein murine double minute (MDM2) binding protein (MTBP) can be considered a pro-oncogene of human malignancies; however, its function and mechanisms in hepatocellular carcinoma (HCC) are still not clear. In the present work, our results demonstrate that MTBP could function as a co-activator of transcription factor E26 transformation-specific sequence (ETS-1), which plays an important role in HCC cell proliferation and/or metastasis and promotes proliferation of HCC cells. Using luciferase and real-time polymerase chain reaction (qPCR) assays, MTBP was found to enhance the transcription factor activation of ETS-1. The results from chromatin co-immunoprecipitation showed that MTBP enhanced the recruitment of ETS-1 to its downstream gene’s (mmp1’s) promoter region with ETS-1 binding sites. In cellular and nude mice models, overexpression of MTBP was shown to promote the proliferation of MHCC97-L cells with low endogenous MTBP levels, whereas the knockdown of MTBP led to inhibition of the proliferation of MHCC97-H cells that possessed high endogenous levels of MTBP. The effect of MTBP on ETS-1 was confirmed in the clinical specimens; the expression of MTBP was positively correlated with the downstream genes of ETS-1, mmp3, mmp9, and uPA. Therefore, by establishing the role of MTBP as a novel co-activator of ETS-1, this work expands our knowledge of MTBP or ETS-1 and helps to provide new ideas concerning HCC-related research.

Whole genome analysis reveals the genomic complexity in metastatic cutaneous squamous cell carcinoma

Metastatic cutaneous squamous cell carcinoma (CSCC) is a highly morbid disease requiring radical surgery and adjuvant therapy, which is associated with a poor prognosis. Yet, compared to other advanced malignancies, relatively little is known of the genomic landscape of metastatic CSCC. We have previously reported the mutational signatures and mutational patterns of CCCTC-binding factor (CTCF) regions in metastatic CSCC. However, many other genomic components (indel signatures, non-coding drivers, and structural variants) of metastatic CSCC have not been reported. To this end, we performed whole genome sequencing on lymph node metastases and blood DNA from 25 CSCC patients with regional metastases of the head and neck. We designed a multifaceted computational analysis at the whole genome level to provide a more comprehensive perspective of the genomic landscape of metastatic CSCC. In the non-coding genome, 3′ untranslated region (3′UTR) regions of EVC (48% of specimens), PPP1R1A (48% of specimens), and ABCA4 (20% of specimens) along with the tumor-suppressing long non-coding RNA (lncRNA) LINC01003 (64% of specimens) were significantly functionally altered (Q-value < 0.05) and represent potential non-coding biomarkers of CSCC. Recurrent copy number loss in the tumor suppressor gene PTPRD was observed. Gene amplification was much less frequent, and few genes were recurrently amplified. Single nucleotide variants driver analyses from three tools confirmed TP53 and CDKN2A as recurrently mutated genes but also identified C9 as a potential novel driver in this disease. Furthermore, indel signature analysis highlighted the dominance of ID signature 13 (ID13) followed by ID8 and ID9. ID9 has previously been shown to have no association with skin melanoma, unlike ID13 and ID8, suggesting a novel pattern of indel variation in metastatic CSCC. The enrichment analysis of various genetically altered candidates shows enrichment of “TGF-beta regulation of extracellular matrix” and “cell cycle G1 to S check points.” These enriched terms are associated with genetic instability, cell proliferation, and migration as mechanisms of genomic drivers of metastatic CSCC.

Mediating and maintaining methylation while minimizing mutation: Recent advances on mammalian DNA methyltransferases

Mammalian genomes are methylated on carbon-5 of many cytosines, mostly in CpG dinucleotides. Methylation patterns are maintained during mitosis via DNMT1, and regulatory factors involved in processes that include histone modifications. Methylation in a sequence longer than CpG can influence the binding of sequence-specific transcription factors, thus affecting gene expression. 5-Methylcytosine deamination results in C-to-T transition. While some mutations are beneficial, most are not; so boosting C-to-T transitions can be dangerous. Given the role of DNMT3A in establishing de novo DNA methylation during development, it is this CpG methylation and deamination that provide the major mutagenic impetus in the DNMT3A gene itself, including the R882H dominant-negative substitution associated with two diseases: germline mutations in DNMT3A overgrowth syndrome, and somatic mutations in clonal hematopoiesis that can initiate acute myeloid leukemia. We discuss recent developments in therapeutics targeting DNMT1, the role of noncatalytic isoform DNMT3B3 in regulating de novo methylation by DNMT3A, and structural characterization of DNMT3A in various configurations.

Differential ETS1 binding to T:G mismatches within a CpG dinucleotide contributes to C-to-T somatic mutation rate of the IDH2 hotspot at codon Arg140

Cytosine to thymine (C>T) somatic mutation is highly enriched in certain types of cancer, and most commonly occurs via deamination of a 5-methylcytosine (5mC) to thymine, in the context of a CpG dinucleotide. In theory, deamination should occur at equal rates to both 5mC nucleotides on opposite strands. In most cases, the resulting T:G or G:T mismatch can be repaired by thymine DNA glycosylase activities. However, while some hotspot-associated CpG mutations have approximately equal numbers of mutations that resulted either from C>T or G>A in a CpG dinucleotide, many showed strand bias, being skewed toward C>T of the first base pair or G>A of the second base pair. Using the IDH2 Arg140 codon as a case study, we show that the two possible T:G mismatches at the codon-specific CpG site have differing effects on transcription factor ETS1 binding affinity, differentially affecting access of a repair enzyme (MBD4) to the deamination-caused T:G mismatch. Our study thus provides a plausible mechanism for exclusion of repair enzymes by the differential binding of transcription factors affecting the rate at which the antecedent opposite-strand mutations occur.

A new technique for genome-wide mapping of nucleotide excision repair without immunopurification of damaged DNA

Nucleotide excision repair functions to protect genome integrity, and ongoing studies using excision repair sequencing (XR-seq) have contributed to our understanding of how cells prioritize repair across the genome. In this method, the products of excision repair bearing damaged DNA are captured, sequenced, and then mapped genome-wide at single-nucleotide resolution. However, reagent requirements and complex procedures have limited widespread usage of this technique. In addition to the expense of these reagents, it has been hypothesized that the immunoprecipitation step using antibodies directed against damaged DNA may introduce bias in different sequence contexts. Here, we describe a newly developed adaptation called dA-tailing and adaptor ligation (ATL)–XR-seq, a relatively simple XR-seq method that avoids the use of immunoprecipitation targeting damaged DNA. ATL-XR-seq captures repair products by 3′-dA-tailing and 5′-adapter ligation instead of the original 5′- and 3′-dual adapter ligation. This new approach avoids adapter dimer formation during subsequent PCR, omits inefficient and time-consuming purification steps, and is very sensitive. In addition, poly(dA) tail length heterogeneity can serve as a molecular identifier, allowing more repair hotspots to be mapped. Importantly, a comparison of both repair mapping methods showed that no major bias is introduced by the anti-UV damage antibodies used in the original XR-seq procedure. Finally, we also coupled the described dA-tailing approach with quantitative PCR in a new method to quantify repair products. These new methods provide powerful and user-friendly tools to qualitatively and quantitatively measure excision repair.

Mapping atypical UV photoproducts in vitro and across the S. cerevisiae genome

Exposure to ultraviolet (UV) light induces DNA damage, predominantly cyclobutane pyrimidine dimers (CPD) and 6,4-photoproducts (6,4-PP), as well as rare, atypical photoproducts at thymidine-adenine (TA) sequences. We have recently shown 'TA' photoproducts are induced in UV-irradiated oligonucleotides and across the budding yeast genome. Here, we describe a protocol for mapping atypical 'TA' photoproducts in vitro and in vivo. This protocol overcomes the technical challenges involved in accurately mapping such rare photoproducts by using ultraviolet damage endonuclease (UVDE) enzymes. For complete details on the use and execution of this protocol, please refer to Laughery et al. (2020).

High-resolution mapping demonstrates inhibition of DNA excision repair by transcription factors

DNA base damage arises frequently in living cells and needs to be removed by base excision repair (BER) to prevent mutagenesis and genome instability. Both the formation and repair of base damage occur in chromatin and are conceivably affected by DNA-binding proteins such as transcription factors (TFs). However, to what extent TF binding affects base damage distribution and BER in cells is unclear. Here, we used a genome-wide damage mapping method, N-methylpurine-sequencing (NMP-seq), and characterized alkylation damage distribution and BER at TF binding sites in yeast cells treated with the alkylating agent methyl methanesulfonate (MMS). Our data show that alkylation damage formation was mainly suppressed at the binding sites of yeast TFs ARS binding factor 1 (Abf1) and rDNA enhancer binding protein 1 (Reb1), but individual hotspots with elevated damage levels were also found. Additionally, Abf1 and Reb1 binding strongly inhibits BER in vivo and in vitro, causing slow repair both within the core motif and its adjacent DNA. Repair of ultraviolet (UV) damage by nucleotide excision repair (NER) was also inhibited by TF binding. Interestingly, TF binding inhibits a larger DNA region for NER relative to BER. The observed effects are caused by the TF–DNA interaction, because damage formation and BER can be restored by depletion of Abf1 or Reb1 protein from the nucleus. Thus, our data reveal that TF binding significantly modulates alkylation base damage formation and inhibits repair by the BER pathway. The interplay between base damage formation and BER may play an important role in affecting mutation frequency in gene regulatory regions.

Biallelic mutations in cancer genomes reveal local mutational determinants

The infinite sites model of molecular evolution posits that every position in the genome is mutated at most once1. By restricting the number of possible mutation histories, haplotypes and alleles, it forms a cornerstone of tumor phylogenetic analysis2 and is often implied when calling, phasing and interpreting variants3,4 or studying the mutational landscape as a whole5. Here we identify 18,295 biallelic mutations, where the same base is mutated independently on both parental copies, in 559 (21%) bulk sequencing samples from the Pan-Cancer Analysis of Whole Genomes study. Biallelic mutations reveal ultraviolet light damage hotspots at E26 transformation-specific (ETS) and nuclear factor of activated T cells (NFAT) binding sites, and hypermutable motifs in POLE-mutant and other cancers. We formulate recommendations for variant calling and provide frameworks to model and detect biallelic mutations. These results highlight the need for accurate models of mutation rates and tumor evolution, as well as their inference from sequencing data.

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.

CTCF puts a new twist on UV damage and repair in skin cancer

Somatic mutations in skin cancers are highly enriched at binding sites for CCCTC-binding factor (CTCF). We have discovered that CTCF binding alters the DNA structure to render it more susceptible to UV damage. Elevated UV damage formation at CTCF binding sites, in conjunction with subsequent repair inhibition, promotes UV mutagenesis.

Hypermethylation of the Promoter of miR-338-5p Mediates Aberrant Expression of ETS-1 and Is Correlated With Disease Severity Of Astrocytoma Patients

The pro-oncogene ETS-1 (E26 transformation-specific sequence 1) is a key regulator of the proliferation and invasion of cancer cells. The present work examined the correlation of the aberrant expression of ETS-1 with histological or clinical classification of astrocytoma: grade I (pilocytic astrocytoma), grade II (diffuse astrocytoma), grade III (anaplastic astrocytoma), and grade IV (glioblastoma multiforme). MicroRNA, miR-338-5p, was predicted by an online tool (miRDB) to potentially target the 3' untranslated region of ETS-1; this was confirmed by multi-assays, including western blot experiments or the point mutation of the targeting sites of miR-338-5p in ETS-1's 3'untralation region (3'UTR). The expression of miR-338-5p was negatively associated with that of ETS-1 in astrocytoma, and deficiency of miR-338-5p would mediate aberrant expression of ETS-1 in astrocytoma. Mechanistically, hypermethylation of miR-338-5p by DNA methyltransferase 1 (DNMT1) resulted in repression of miR-338-5p expression and the aberrant expression of ETS-1. Knockdown or deactivation of DNMT1 decreased the methylation rate of the miR-338-5p promoter, increased the expression of miR-338-5p, and repressed the expression of ETS-1 in astrocytoma cell lines U251 and U87. These results indicate that hypermethylation of the miR-338-5p promoter by DNMT1 mediates the aberrant expression of ETS-1 related to disease severity of patients with astrocytoma.