Mutational signatures: experimental design and analytical framework

Funder: Dr. Josef Steiner Krebsstiftung; doi: http://dx.doi.org/10.13039/100008669Abstract: Mutational signatures provide a powerful alternative for understanding the pathophysiology of cancer. Currently, experimental efforts aimed at validating and understanding the etiologies of cancer-derived mutational signatures are underway. In this review, we highlight key aspects of mutational signature experimental design and describe the analytical framework. We suggest guidelines and quality control measures for handling whole-genome sequencing data for mutational signature analyses and discuss pitfalls in interpretation. We envision that improved next-generation sequencing technologies and molecular cell biology approaches will usher in the next generation of studies into the etiologies and mechanisms of mutational patterns uncovered in cancers

Short inverted repeats contribute to localized mutability in human somatic cells.

Selected repetitive sequences termed short inverted repeats (SIRs) have the propensity to form secondary DNA structures called hairpins. SIRs comprise palindromic arm sequences separated by short spacer sequences that form the hairpin stem and loop respectively. Here, we show that SIRs confer an increase in localized mutability in breast cancer, which is domain-dependent with the greatest mutability observed within spacer sequences (∼1.35-fold above background). Mutability is influenced by factors that increase the likelihood of formation of hairpins such as loop lengths (of 4-5 bp) and stem lengths (of 7-15 bp). Increased mutability is an intrinsic property of SIRs as evidenced by how almost all mutational processes demonstrate a higher rate of mutagenesis of spacer sequences. We further identified 88 spacer sequences showing enrichment from 1.8- to 90-fold of local mutability distributed across 283 sites in the genome that intriguingly, can be used to inform the biological status of a tumor

Epigenetic modifiers DNMT3A and BCOR are recurrently mutated in CYLD cutaneous syndrome

Abstract: Patients with CYLD cutaneous syndrome (CCS; syn. Brooke-Spiegler syndrome) carry germline mutations in the tumor suppressor CYLD and develop multiple skin tumors with diverse histophenotypes. Here, we comprehensively profile the genomic landscape of 42 benign and malignant tumors across 13 individuals from four multigenerational families and discover recurrent mutations in epigenetic modifiers DNMT3A and BCOR in 29% of benign tumors. Multi-level and microdissected sampling strikingly reveal that many clones with different DNMT3A mutations exist in these benign tumors, suggesting that intra-tumor heterogeneity is common. Integrated genomic, methylation and transcriptomic profiling in selected tumors suggest that isoform-specific DNMT3A2 mutations are associated with dysregulated methylation. Phylogenetic and mutational signature analyses confirm cylindroma pulmonary metastases from primary skin tumors. These findings contribute to existing paradigms of cutaneous tumorigenesis and metastasis

Whole-Genome Sequencing of Retinoblastoma Reveals the Diversity of Rearrangements Disrupting RB1 and Uncovers a Treatment-Related Mutational Signature.

The development of retinoblastoma is thought to require pathological genetic changes in both alleles of the RB1 gene. However, cases exist where RB1 mutations are undetectable, suggesting alternative pathways to malignancy. We used whole-genome sequencing (WGS) and transcriptomics to investigate the landscape of sporadic retinoblastomas derived from twenty patients, sought RB1 and other driver mutations and investigated mutational signatures. At least one RB1 mutation was identified in all retinoblastomas, including new mutations in addition to those previously identified by clinical screening. Ten tumours carried structural rearrangements involving RB1 ranging from relatively simple to extremely complex rearrangement patterns, including a chromothripsis-like pattern in one tumour. Bilateral tumours obtained from one patient harboured conserved germline but divergent somatic RB1 mutations, indicating independent evolution. Mutational signature analysis showed predominance of signatures associated with cell division, an absence of ultraviolet-related DNA damage and a profound platinum-related mutational signature in a chemotherapy-exposed tumour. Most RB1 mutations are identifiable by clinical screening. However, the increased resolution and ability to detect otherwise elusive rearrangements by WGS have important repercussions on clinical management and advice on recurrence risks

Detection and quantification of methylation in DNA using solid-state nanopores.

Epigenetic modifications in eukaryotic genomes occur primarily in the form of 5-methylcytosine (5 mC). These modifications are heavily involved in transcriptional repression, gene regulation, development and the progression of diseases including cancer. We report a new single-molecule assay for the detection of DNA methylation using solid-state nanopores. Methylation is detected by selectively labeling methylation sites with MBD1 (MBD-1x) proteins, the complex inducing a 3 fold increase in ionic blockage current relative to unmethylated DNA. Furthermore, the discrimination of methylated and unmethylated DNA is demonstrated in the presence of only a single bound protein, thereby giving a resolution of a single methylated CpG dinucleotide. The extent of methylation of a target molecule could also be coarsely quantified using this novel approach. This nanopore-based methylation sensitive assay circumvents the need for bisulfite conversion, fluorescent labeling, and PCR and could therefore prove very useful in studying the role of epigenetics in human disease

A systematic CRISPR screen defines mutational mechanisms underpinning signatures caused by replication errors and endogenous DNA damage.

Mutational signatures are imprints of pathophysiological processes arising through tumorigenesis. We generated isogenic CRISPR-Cas9 knockouts (Δ) of 43 genes in human induced pluripotent stem cells, cultured them in the absence of added DNA damage, and performed whole-genome sequencing of 173 subclones. ΔOGG1, ΔUNG, ΔEXO1, ΔRNF168, ΔMLH1, ΔMSH2, ΔMSH6, ΔPMS1, and ΔPMS2 produced marked mutational signatures indicative of being critical mitigators of endogenous DNA modifications. Detailed analyses revealed mutational mechanistic insights, including how 8-oxo-dG elimination is sequence-context-specific while uracil clearance is sequence-context-independent. Mismatch repair (MMR) deficiency signatures are engendered by oxidative damage (C>A transversions), differential misincorporation by replicative polymerases (T>C and C>T transitions), and we propose a 'reverse template slippage' model for T>A transversions. ΔMLH1, ΔMSH6, and ΔMSH2 signatures were similar to each other but distinct from ΔPMS2. Finally, we developed a classifier, MMRDetect, where application to 7,695 WGS cancers showed enhanced detection of MMR-deficient tumors, with implications for responsiveness to immunotherapies

The topography of mutational processes in breast cancer genomes.

Somatic mutations in human cancers show unevenness in genomic distribution that correlate with aspects of genome structure and function. These mutations are, however, generated by multiple mutational processes operating through the cellular lineage between the fertilized egg and the cancer cell, each composed of specific DNA damage and repair components and leaving its own characteristic mutational signature on the genome. Using somatic mutation catalogues from 560 breast cancer whole-genome sequences, here we show that each of 12 base substitution, 2 insertion/deletion (indel) and 6 rearrangement mutational signatures present in breast tissue, exhibit distinct relationships with genomic features relating to transcription, DNA replication and chromatin organization. This signature-based approach permits visualization of the genomic distribution of mutational processes associated with APOBEC enzymes, mismatch repair deficiency and homologous recombinational repair deficiency, as well as mutational processes of unknown aetiology. Furthermore, it highlights mechanistic insights including a putative replication-dependent mechanism of APOBEC-related mutagenesis

Recognition of methylated DNA through methyl-CpG binding domain proteins

DNA methylation is a key regulatory control route in epigenetics, involving gene silencing and chromosome inactivation. It has been recognized that methyl-CpG binding domain (MBD) proteins play an important role in interpreting the genetic information encoded by methylated DNA (mDNA). Although the function of MBD proteins has attracted considerable attention and is well characterized, the mechanism underlying mDNA recognition by MBD proteins is still poorly understood. In this article, we demonstrate that the methyl-CpG dinucleotides are recognized at the MBD–mDNA interface by two MBD arginines through an interplay of hydrogen bonding and cation-π interaction. Through molecular dynamics and quantum-chemistry calculations we investigate the methyl-cytosine recognition process and demonstrate that methylation enhances MBD–mDNA binding by increasing the hydrophobic interfacial area and by strengthening the interaction between mDNA and MBD proteins. Free-energy perturbation calculations also show that methylation yields favorable contribution to the binding free energy for MBD–mDNA complex

A somatic-mutational process recurrently duplicates germline susceptibility loci and tissue-specific super-enhancers in breast cancers

Somatic rearrangements contribute to the mutagenized landscape of cancer genomes. Here, we systematically interrogated rearrangements in 560 breast cancers by using a piecewise constant fitting approach. We identified 33 hotspots of large (>100 kb) tandem duplications, a mutational signature associated with homologous-recombination-repair deficiency. Notably, these tandem-duplication hotspots were enriched in breast cancer germline susceptibility loci (odds ratio (OR) = 4.28) and breast-specific 'super-enhancer' regulatory elements (OR = 3.54). These hotspots may b