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

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.

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.

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.

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.

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.

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.

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.

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.

Spatially resolved transcriptomics reveals the architecture of the tumor-microenvironment interface

During tumor progression, cancer cells come into contact with various non-tumor cell types, but it is unclear how tumors adapt to these new environments. Here, we integrate spatially resolved transcriptomics, single-cell RNA-seq, and single-nucleus RNA-seq to characterize tumor-microenvironment interactions at the tumor boundary. Using a zebrafish model of melanoma, we identify a distinct "interface" cell state where the tumor contacts neighboring tissues. This interface is composed of specialized tumor and microenvironment cells that upregulate a common set of cilia genes, and cilia proteins are enriched only where the tumor contacts the microenvironment. Cilia gene expression is regulated by ETS-family transcription factors, which normally act to suppress cilia genes outside of the interface. A cilia-enriched interface is conserved in human patient samples, suggesting it is a conserved feature of human melanoma. Our results demonstrate the power of spatially resolved transcriptomics in uncovering mechanisms that allow tumors to adapt to new environments.

CTCF binding modulates UV damage formation to promote mutation hot spots in melanoma

Somatic mutations in DNA-binding sites for CCCTC-binding factor (CTCF) are significantly elevated in many cancers. Prior analysis has suggested that elevated mutation rates at CTCF-binding sites in skin cancers are a consequence of the CTCF-cohesin complex inhibiting repair of UV damage. Here, we show that CTCF binding modulates the formation of UV damage to induce mutation hot spots. Analysis of genome-wide CPD-seq data in UV-irradiated human cells indicates that formation of UV-induced cyclobutane pyrimidine dimers (CPDs) is primarily suppressed by CTCF binding but elevated at specific locations within the CTCF motif. Locations of CPD hot spots in the CTCF-binding motif coincide with mutation hot spots in melanoma. A similar pattern of damage formation is observed at CTCF-binding sites in vitro, indicating that UV damage modulation is a direct consequence of CTCF binding. We show that CTCF interacts with binding sites containing UV damage and inhibits repair by a model repair enzyme in vitro. Structural analysis and molecular dynamic simulations reveal the molecular mechanism for how CTCF binding modulates CPD formation.

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.

Genome-wide mapping of genomic DNA damage: methods and implications

Exposures from the external and internal environments lead to the modification of genomic DNA, which is implicated in the cause of numerous diseases, including cancer, cardiovascular, pulmonary and neurodegenerative diseases, together with ageing. However, the precise mechanism(s) linking the presence of damage, to impact upon cellular function and pathogenesis, is far from clear. Genomic location of specific forms of damage is likely to be highly informative in understanding this process, as the impact of downstream events (e.g. mutation, microsatellite instability, altered methylation and gene expression) on cellular function will be positional—events at key locations will have the greatest impact. However, until recently, methods for assessing DNA damage determined the totality of damage in the genomic location, with no positional information. The technique of “mapping DNA adductomics” describes the molecular approaches that map a variety of forms of DNA damage, to specific locations across the nuclear and mitochondrial genomes. We propose that integrated comparison of this information with other genome-wide data, such as mutational hotspots for specific genotoxins, tumour-specific mutation patterns and chromatin organisation and transcriptional activity in non-cancerous lesions (such as nevi), pre-cancerous conditions (such as polyps) and tumours, will improve our understanding of how environmental toxins lead to cancer. Adopting an analogous approach for non-cancer diseases, including the development of genome-wide assays for other cellular outcomes of DNA damage, will improve our understanding of the role of DNA damage in pathogenesis more generally.

Comparative analyses of two primate species diverged by more than 60 million years show different rates but similar distribution of genome-wide UV repair events

Background

Nucleotide excision repair is the primary DNA repair mechanism that removes bulky DNA adducts such as UV-induced pyrimidine dimers. Correspondingly, genome-wide mapping of nucleotide excision repair with eXcision Repair sequencing (XR-seq), provides comprehensive profiling of DNA damage repair. A number of XR-seq experiments at a variety of conditions for different damage types revealed heterogenous repair in the human genome. Although human repair profiles were extensively studied, how repair maps vary between primates is yet to be investigated. Here, we characterized the genome-wide UV-induced damage repair in gray mouse lemur, Microcebus murinus, in comparison to human.

Results

We derived fibroblast cell lines from mouse lemur, exposed them to UV irradiation, and analyzed the repair events genome-wide using the XR-seq protocol. Mouse lemur repair profiles were analyzed in comparison to the equivalent human fibroblast datasets. We found that overall UV sensitivity, repair efficiency, and transcription-coupled repair levels differ between the two primates. Despite this, comparative analysis of human and mouse lemur fibroblasts revealed that genome-wide repair profiles of the homologous regions are highly correlated, and this correlation is stronger for highly expressed genes. With the inclusion of an additional XR-seq sample derived from another human cell line in the analysis, we found that fibroblasts of the two primates repair UV-induced DNA lesions in a more similar pattern than two distinct human cell lines do.

Conclusion

Our results suggest that mouse lemurs and humans, and possibly primates in general, share a homologous repair mechanism as well as genomic variance distribution, albeit with their variable repair efficiency. This result also emphasizes the deep homologies of individual tissue types across the eukaryotic phylogeny.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07898-3.

Non-coding driver mutations in human cancer

Tumour formation involves random mutagenic events and positive evolutionary selection acting on a subset of such events, referred to as driver mutations. A decade of careful surveying of tumour DNA using exome-based analyses has revealed a multitude of protein-coding somatic driver mutations, some of which are clinically actionable. Today, a transition towards whole-genome analysis is well under way, technically enabling the discovery of potential driver mutations occurring outside protein-coding sequences. Mutations are abundant in this vast non-coding space, which is more than 50 times larger than the coding exome, but reliable identification of selection signals in non-coding DNA remains a challenge. In this Review, we discuss recent findings in the field, where the emerging landscape is one in which non-coding driver mutations appear to be relatively infrequent. Nevertheless, we highlight several notable discoveries. We consider possible reasons for the relative absence of non-coding driver events, as well as the difficulties associated with detecting signals of positive selection in non-coding DNA.

The landscape of driver mutations in cutaneous squamous cell carcinoma

Cutaneous squamous cell carcinoma is a form of skin cancer originating from keratinocytes in the skin. It is the second most common type of cancer and is responsible for an estimated 8000 deaths per year in the United States. Compared to other cancer subtypes with similar incidences and death tolls, our understanding of the somatic mutations driving cutaneous squamous cell carcinoma is limited. The main challenge is that these tumors have high mutation burdens, primarily a consequence of UV-radiation-induced DNA damage from sunlight, making it difficult to distinguish driver mutations from passenger mutations. We overcame this challenge by performing a meta-analysis of publicly available sequencing data covering 105 tumors from 10 different studies. Moreover, we eliminated tumors with issues, such as low neoplastic cell content, and from the tumors that passed quality control, we utilized multiple strategies to reveal genes under selection. In total, we nominated 30 cancer genes. Among the more novel genes, mutations frequently affected EP300, PBRM1, USP28, and CHUK. Collectively, mutations in the NOTCH and p53 pathways were ubiquitous, and to a lesser extent, mutations affected genes in the Hippo pathway, genes in the Ras/MAPK/PI3K pathway, genes critical for cell-cycle checkpoint control, and genes encoding chromatin remodeling factors. Taken together, our study provides a catalog of driver genes in cutaneous squamous cell carcinoma, offering points of therapeutic intervention and insights into the biology of cutaneous squamous cell carcinoma.