Intragenomic variability and extended sequence patterns in the mutational signature of ultraviolet light

Case report: Is severe toxicity the price to pay for high sensitivity to checkpoint inhibitors immunotherapy in desmoplastic melanoma?

Background

Desmoplastic melanoma (DM) is a rare subtype of melanoma characterized by high immunogenicity which makes it particularly suitable for immune checkpoint inhibitors (ICIs) treatment.

Case presentation

We report the case of a 53-year-old man with metastatic DM successfully treated with the combination of anti-CTLA-4 and anti-PD-1 antibodies, who developed serious immune-related adverse events (irAEs). The primary tumor was characterized by absent PD-L1 expression and no-brisk lymphocytes infiltration. NGS showed absence of BRAF mutation, a high tumor mutational burden, and an UV-induced DNA damage signature. Metastatic lesions regressed rapidly after few cycles of ICIs until complete response, however the patient developed serious irAEs including hypothyroidism, adrenal deficiency, and acute interstitial nephritis which led to the definitive suspension of treatment. Currently, the patient has normal renal functionality and no disease relapse after 26 months from starting immunotherapy, and after 9 months from its definitive suspension.

Conclusion

Efficacy and toxicity are two sides of the same coin of high sensitivity to ICIs in DM. For this reason, these patients should be closely monitored during ICIs therapy to promptly identify serious side effects and to correctly manage them.

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.

Quantitative and qualitative mutational impact of ionizing radiation on normal cells

The comprehensive genomic impact of ionizing radiation (IR), a carcinogen, on healthy somatic cells remains unclear. Using large-scale whole-genome sequencing (WGS) of clones expanded from irradiated murine and human single cells, we revealed that IR induces a characteristic spectrum of short insertions or deletions (indels) and structural variations (SVs), including balanced inversions, translocations, composite SVs (deletion-insertion, deletion-inversion, and deletion-translocation composites), and complex genomic rearrangements (CGRs), including chromoplexy, chromothripsis, and SV by breakage-fusion-bridge cycles. Our findings suggest that 1 Gy IR exposure causes an average of 2.33 mutational events per Gb genome, comprising 2.15 indels, 0.17 SVs, and 0.01 CGRs, despite a high level of inter-cellular stochasticity. The mutational burden was dependent on total irradiation dose, regardless of dose rate or cell type. The findings were further validated in IR-induced secondary cancers and single cells without clonalization. Overall, our study highlights a comprehensive and clear picture of IR effects on normal mammalian genomes.

Flexible model-based non-negative matrix factorization with application to mutational signatures

Somatic mutations in cancer can be viewed as a mixture distribution of several mutational signatures, which can be inferred using non-negative matrix factorization (NMF). Mutational signatures have previously been parametrized using either simple mono-nucleotide interaction models or general tri-nucleotide interaction models. We describe a flexible and novel framework for identifying biologically plausible parametrizations of mutational signatures, and in particular for estimating di-nucleotide interaction models. Our novel estimation procedure is based on the expectation-maximization (EM) algorithm and regression in the log-linear quasi-Poisson model. We show that di-nucleotide interaction signatures are statistically stable and sufficiently complex to fit the mutational patterns. Di-nucleotide interaction signatures often strike the right balance between appropriately fitting the data and avoiding over-fitting. They provide a better fit to data and are biologically more plausible than mono-nucleotide interaction signatures, and the parametrization is more stable than the parameter-rich tri-nucleotide interaction signatures. We illustrate our framework in a large simulation study where we compare to state of the art methods, and show results for three data sets of somatic mutation counts from patients with cancer in the breast, Liver and urinary tract.

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.

Is Cutaneous T-Cell Lymphoma Caused by Ultraviolet Radiation? A Comparison of UV Mutational Signatures in Malignant Melanoma and Mycosis Fungoides

Ultraviolet (UV) radiation is a strong environmental carcinogen responsible for the pathogenesis of most skin cancers, including malignant melanoma (MM) and non-melanoma (keratinocyte) skin cancers. The carcinogenic role of UV was firmly established based on epidemiological evidence and molecular findings of the characteristic mutation signatures which occur during the excision repair of cyclobutane pyrimidine dimers and 6,4-photoproducts. The role of UV in the pathogenesis of mycosis fungoides (MF), the most common type of primary cutaneous T-cell lymphoma, remains controversial. Here, we performed whole-exome sequencing of 61 samples of MF cells microdissected from cutaneous lesions, and compared their mutational signatures to 340 MMs. The vast majority of MM mutations had a typical UV mutational signature (SBS 7, SBS 38, or DSB 1), underscoring the key role of ultraviolet as a mutagen. In contrast, the SBS 7 signature in MF comprised < 5% of all mutations. SBS 7 was higher in the intraepidermal MF cells (when compared to the dermal cells) and in the cells from tumors as compared to that in early-stage plaques. In conclusion, our data do not support the pathogenic role of UV in the pathogenesis of MF and suggest that the UV mutations are the result of the cumulative environmental ultraviolet exposure of cutaneous lesions rather than an early mutagenic event.

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.

Somatic mutation distribution across tumour cohorts provides a signal for positive selection in cancer

Cancer gene discovery is reliant on distinguishing driver mutations from a multitude of passenger mutations in tumour genomes. While driver genes may be revealed based on excess mutation recurrence or clustering, there is a need for orthogonal principles. Here, we take advantage of the fact that non-cancer genes, containing only passenger mutations under neutral selection, exhibit a likelihood of mutagenesis in a given tumour determined by the tumour's mutational signature and burden. This relationship can be disrupted by positive selection, leading to a difference in the distribution of mutated cases across a cohort for driver and passenger genes. We apply this principle to detect cancer drivers independently of recurrence in large pan-cancer cohorts, and show that our method (SEISMIC) performs comparably to traditional approaches and can provide resistance to known confounding mutational phenomena. Being based on a different principle, the approach provides a much-needed complement to existing methods for detecting signals of selection.

Somatic mutations during rapid clonal domestication of Populus alba var. pyramidalis

For many clonally propagated species, the accumulation of somatic mutations is the principal driver of declines in yield and quality. However, somatic mutations may also promote genetic diversification. Thus, elucidating somatic mutation rates and patterns is important to understand the genetic basis undergirding the emergence of commercially valuable traits and developmental processes. In this study, we studied the effect of short-time clonal domestication of Populus alba var. pyramidalis, a species that has been propagated by cutting for the last 67 years. We found that: (1) the somatic mutation rate for P. alba var. pyramidalis is 9.24 × 10-9, which is higher than rates observed in related species; (2) there were more mutations near heterozygous regions, and a larger proportion of CpG and CHG sites were associated with somatic mutations, which may be related to the blocking of DNA repair by methylation; and (3) deleterious mutations were not shared by multiple individuals, and all occurred in heterozygous states, demonstrating the strong selective pressures that act against deleterious mutations. Taken together, the results of our study provide a global view of somatic mutation that will aid efforts to understand the genetic basis of commercially valuable traits and to improve clonally breeding species.

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.

The major mechanism of melanoma mutations is based on deamination of cytosine in pyrimidine dimers as determined by circle damage sequencing

Sunlight-associated melanomas carry a unique C-to-T mutation signature. UVB radiation induces cyclobutane pyrimidine dimers (CPDs) as the major form of DNA damage, but the mechanism of how CPDs cause mutations is unclear. To map CPDs at single-base resolution genome wide, we developed the circle damage sequencing (circle-damage-seq) method. In human cells, CPDs form preferentially in a tetranucleotide sequence context (5'-Py-T<>Py-T/A), but this alone does not explain the tumor mutation patterns. To test whether mutations arise at CPDs by cytosine deamination, we specifically mapped UVB-induced cytosine-deaminated CPDs. Transcription start sites (TSSs) were protected from CPDs and deaminated CPDs, but both lesions were enriched immediately upstream of the TSS, suggesting a mutation-promoting role of bound transcription factors. Most importantly, the genomic dinucleotide and trinucleotide sequence specificity of deaminated CPDs matched the prominent mutation signature of melanomas. Our data identify the cytosine-deaminated CPD as the leading premutagenic lesion responsible for mutations in melanomas.

Undifferentiated large cell/rhabdoid carcinoma presenting in the intestines of patients with concurrent or recent non-small cell lung cancer (NSCLC): clinicopathologic and molecular analysis of 14 cases indicates an unusual pattern of dedifferentiated metastases

Undifferentiated carcinoma metastatic to the bowel is uncommon in surgical pathology practice and might be confused with primary gastrointestinal carcinoma, melanoma, lymphoma, and others. We present 14 cases of uni- (n = 9) or multifocal (n = 5) undifferentiated large cell/rhabdoid carcinoma presenting in the bowel of patients with concurrent (n = 9) or recent (diagnosed 1 to 25 months earlier; median, 4) non-small cell lung cancer (NSCLC). Patients were 6 females and 8 males, aged 52 to 85 years. Primary NSCLC was verified histologically in 10 cases and by imaging in 4. The undifferentiated histology was present in the lung biopsy in 4/10 patients (as sole pattern in 3 and combined with adenocarcinoma in 1) and was limited to the intestinal metastases in the remainder. PDL1 was strongly expressed in 7/9 cases (CPS: 41 to 100). Loss of at least one SWI/SNF subunit was detected in 7/13 cases (54%). SMARCA2 loss (n = 6) was most frequent and was combined with SMARCA4 loss in one case. PBRM1 loss was observed in one tumor. Successful molecular testing of 11 cases revealed BRAF mutations in 4 (3 were non-V600E variants), KRAS mutations in 3, and wildtype in 4. None had EGFR mutations. Analysis of 4 paired samples revealed concordant KRAS (2) and BRAF (1) mutations or wildtype (1). Our study indicates that undifferentiated carcinoma within the intestines of patients with concurrent/recent NSCLC represents dedifferentiated metastasis from the NSCLC. Recognition of this unusual presentation is cardinal to avoid misdiagnosis with inappropriate therapeutic and prognostic implications.

Human MettL3-MettL14 RNA adenine methyltransferase complex is active on double-stranded DNA containing lesions

MettL3-MettL14 methyltransferase complex has been studied widely for its role in RNA adenine methylation. This complex is also recruited to UV- and X-ray exposed DNA damaged sites, and its methyltransfer activity is required for subsequent DNA repair, though in theory this could result from RNA methylation of short transcripts made at the site of damage. We report here that MettL3-MettL14 is active in vitro on double-stranded DNA containing a cyclopyrimidine dimer – a major lesion of UV radiation-induced products – or an abasic site or mismatches. Furthermore, N6-methyladenine (N6mA) decreases misincorporation of 8-oxo-guanine (8-oxoG) opposite to N6mA by repair DNA polymerases. When 8-oxoG is nevertheless incorporated opposite N6mA, the methylation inhibits N6mA excision from the template (correct) strand by the adenine DNA glycosylase (MYH), implying that the methylation decreases inappropriate misrepair. Finally, we observed that the N6mA reader domain of YTHDC1, which is also recruited to sites of DNA damage, binds N6mA that is located across from a single-base gap between two canonical DNA helices. This YTHDC1 complex with a gapped duplex is structurally similar to DNA complexes with FEN1 and GEN1 – two members of the nuclease family that act in nucleotide excision repair, mismatch repair and homologous recombination, and which incise distinct non-B DNA structures. Together, the parts of our study provide a plausible mechanism for N6mA writer and reader proteins acting directly on lesion-containing DNA, and suggest in vivo experiments to test the mechanisms involving methylation of adenine.

A streamlined solution for processing, elucidating and quality control of cyclobutane pyrimidine dimer sequencing data

UV radiation may lead to melanoma and nonmelanoma skin cancers by causing helix-distorting DNA damage such as cyclobutane pyrimidine dimers (CPDs). These DNA lesions, if located in important genes and not repaired promptly, are mutagenic and may eventually result in carcinogenesis. Examining CPD formation and repair processes across the genome can shed light on the mutagenesis mechanisms associated with UV damage in relevant cancers. We recently developed CPD-Seq, a high-throughput and single-nucleotide resolution sequencing technique that can specifically capture UV-induced CPD lesions across the genome. This novel technique has been increasingly used in studies of UV damage and can be adapted to sequence other clinically relevant DNA lesions. Although the library preparation protocol has been established, a systematic protocol to analyze CPD-Seq data has not been described yet. To streamline the various general or specific analysis steps, we developed a protocol named CPDSeqer to assist researchers with CPD-Seq data processing. CPDSeqer can accommodate both a single- and multiple-sample experimental design, and it allows both genome-wide analyses and regional scrutiny (such as of suspected UV damage hotspots). The runtime of CPDSeqer scales with raw data size and takes roughly 4 h per sample with the possibility of acceleration by parallel computing. Various guiding graphics are generated to help diagnose the performance of the experiment and inform regional enrichment of CPD formation. UV damage comparison analyses are set forth in three analysis scenarios, and the resulting HTML pages report damage directional trends and statistical significance. CPDSeqer can be accessed at https://github.com/shengqh/cpdseqer .

Dedifferentiated and Undifferentiated Melanomas: Report of 35 New Cases With Literature Review and Proposal of Diagnostic Criteria

Dedifferentiated melanoma (DM) and undifferentiated melanoma (UM) is defined as a primary or metastatic melanoma showing transition between conventional and undifferentiated components (DM) or lacking histologic and immunophenotypic features of melanoma altogether (UM). The latter is impossible to verify as melanoma by conventional diagnostic tools alone. We herein describe our experience with 35 unpublished cases to expand on their morphologic, phenotypic, and genotypic spectrum, along with a review of 50 previously reported cases (total: 85) to establish the diagnostic criteria. By definition, the dedifferentiated/undifferentiated component lacked expression of 5 routinely used melanoma markers (S100, SOX10, Melan-A, HMB45, Pan-melanoma). Initial diagnoses (known in 66 cases) were undifferentiated/unclassified pleomorphic sarcoma (n=30), unclassified epithelioid malignancy (n=7), pleomorphic rhabdomyosarcoma (n=5), other specific sarcoma types (n=6), poorly differentiated carcinoma (n=2), collision tumor (n=2), atypical fibroxanthoma (n=2), and reactive osteochondromatous lesion (n=1). In only 11 cases (16.6%) was a diagnosis of melanoma considered. Three main categories were identified: The largest group (n=56) comprised patients with a history of verified previous melanoma who presented with metastatic DM or UM. Axillary or inguinal lymph nodes, soft tissue, bone, and lung were mainly affected. A melanoma-compatible mutation was detected in 35 of 48 (73%) evaluable cases: BRAF (n=20; 40.8%), and NRAS (n=15; 30.6%). The second group (n=15) had clinicopathologic features similar to group 1, but a melanoma history was lacking. Axillary lymph nodes (n=6) was the major site in this group followed by the lung, soft tissue, and multiple site involvement. For this group, NRAS mutation was much more frequent (n=9; 60%) than BRAF (n=3; 20%) and NF1 (n=1; 6.6%). The third category (n=14) comprised primary DM (12) or UM (2). A melanoma-compatible mutation was detected in only 7 cases: BRAF (n=2), NF1 (n=2), NRAS (n=2), and KIT exon 11 (n=1). This extended follow-up study highlights the high phenotypic plasticity of DM/UM and indicates significant underrecognition of this aggressive disease among general surgical pathologists. The major clues to the diagnosis of DM and UM are: (1) presence of minimal differentiated clone in DM, (2) earlier history of melanoma, (3) undifferentiated histology that does not fit any defined entity, (4) locations at sites that are unusual for undifferentiated/unclassified pleomorphic sarcoma (axilla, inguinal, neck, digestive system, etc.), (5) unusual multifocal disease typical of melanoma spread, (6) detection of a melanoma-compatible gene mutation, and (7) absence of another genuine primary (eg, anaplastic carcinoma) in other organs.

Human genes differ by their UV sensitivity estimated through analysis of UV-induced silent mutations in melanoma

We hypothesized that human genes differ by their sensitivity to ultraviolet (UV) exposure. We used somatic mutations detected by genome-wide screens in melanoma and reported in the Catalog Of Somatic Mutations In Cancer. As a measure of UV sensitivity, we used the number of silent mutations generated by C>T transitions in pyrimidine dimers of a given transcript divided by the number of potential sites for this type of mutations in the transcript. We found that human genes varied by UV sensitivity by two orders of magnitude. We noted that the melanoma-associated tumor suppressor gene CDKN2A was among the top five most UV-sensitive genes in the human genome. Melanoma driver genes have a higher UV-sensitivity compared with other genes in the human genome. The difference was more prominent for tumor suppressors compared with oncogene. The results of this study suggest that differential sensitivity of human transcripts to UV light may explain melanoma specificity of some driver genes. Practical significance of the study relates to the fact that differences in UV sensitivity among human genes need to be taken into consideration whereas predicting melanoma-associated genes by the number of somatic mutations detected in a given gene.

New insights into the functions of Cox-2 in skin and esophageal malignancies

Understanding the cellular and molecular mechanisms of tumor initiation and progression for each cancer type is central to making improvements in both prevention and therapy. Identifying the cancer cells of origin and the necessary and sufficient mechanisms of transformation and progression provide opportunities for improved specific clinical interventions. In the last few decades, advanced genetic manipulation techniques have facilitated rapid progress in defining the etiologies of cancers and their cells of origin. Recent studies driven by various groups have provided experimental evidence indicating the cellular origins for each type of skin and esophageal cancer and have identified underlying mechanisms that stem/progenitor cells use to initiate tumor development. Specifically, cyclooxygenase-2 (Cox-2) is associated with tumor initiation and progression in many cancer types. Recent studies provide data demonstrating the roles of Cox-2 in skin and esophageal malignancies, especially in squamous cell carcinomas (SCCs) occurring in both sites. Here, we review experimental evidence aiming to define the origins of skin and esophageal cancers and discuss how Cox-2 contributes to tumorigenesis and differentiation.

Intragenomic variability and extended sequence patterns in the mutational signature of ultraviolet light

Mutational signatures have emerged as essential tools in cancer genomics, providing clinically relevant insights as well as accurate background models needed when assessing signals of selection in cancer. Here, we observe that the mutational signature of ultraviolet (UV) light varies across chromatin states, highlighting a previously unappreciated aspect of mutational signatures. Our results imply that locally derived, rather than genome-wide or exome-wide, signatures are more accurate, which is of relevance in situations such as cancer driver gene detection, where correct modelling of signatures and expected mutation rates is critical. We also show that incorporation of longer contextual patterns into the signature further improves modeling of UV mutations.

Mutational signatures can reveal properties of underlying mutational processes and are important when assessing signals of selection in cancer. Here, we describe the sequence characteristics of mutations induced by ultraviolet (UV) light, a major mutagen in several human cancers, in terms of extended (longer than trinucleotide) patterns as well as variability of the signature across chromatin states. Promoter regions display a distinct UV signature with reduced TCG > TTG transitions, and genome-wide mapping of UVB-induced DNA photoproducts (pyrimidine dimers) showed that this may be explained by decreased damage formation at hypomethylated promoter CpG sites. Further, an extended signature model encompassing additional information from longer contextual patterns improves modeling of UV mutations, which may enhance discrimination between drivers and passenger events. Our study presents a refined picture of the UV signature and underscores that the characteristics of a single mutational process may vary across the genome.