Examining clustered somatic mutations with SigProfilerClusters

APOBEC affects tumor evolution and age at onset of lung cancer in smokers

Most solid tumors harbor somatic mutations attributed to off-target activities of APOBEC3A (A3A) and/or APOBEC3B (A3B). However, how APOBEC3A/B enzymes affect tumor evolution in the presence of exogenous mutagenic processes is largely unknown. Here, multi-omics profiling of 309 lung cancers from smokers identifies two subtypes defined by low (LAS) and high (HAS) APOBEC mutagenesis. LAS are enriched for A3B-like mutagenesis and KRAS mutations; HAS for A3A-like mutagenesis and TP53 mutations. Compared to LAS, HAS have older age at onset and high proportions of newly generated progenitor-like cells likely due to the combined tobacco smoking- and APOBEC3A-associated DNA damage and apoptosis. Consistently, HAS exhibit high expression of pulmonary healing signaling pathway, stemness markers, distal cell-of-origin, more neoantigens, slower clonal expansion, but no smoking-associated genomic/epigenomic changes. With validation in 184 lung tumor samples, these findings show how heterogeneity in mutational burden across co-occurring mutational processes and cell types contributes to tumor development.

Tobacco smoke carcinogens exacerbate APOBEC mutagenesis and carcinogenesis

Mutations in somatic cells are inflicted by both extrinsic and intrinsic sources and contribute over time to cancer. Tobacco smoke contains chemical carcinogens that have been causatively implicated with cancers of the lung and head & neck1,2. APOBEC family DNA cytosine deaminases have emerged as endogenous sources of mutation in cancer, with hallmark mutational signatures (SBS2/SBS13) that often co-occur in tumors of tobacco smokers with an equally diagnostic mutational signature (SBS4)3,4. Here we challenge the dogma that mutational processes are thought to occur independently and with additive impact by showing that 4-nitroquinoline 1-oxide (NQO), a model carcinogen for tobacco exposure, sensitizes cells to APOBEC3B (A3B) mutagenesis and leads to synergistic increases in both SBS2 mutation loads and oral carcinomas in vivo. NQO-exposed/A3B-expressing animals exhibit twice as many head & neck lesions as carcinogen-exposed wildtype animals. This increase in carcinogenesis is accompanied by a synergistic increase in mutations from APOBEC signature SBS2, but not from NQO signature SBS4. Interestingly, a large proportion of A3B-catalyzed SBS2 mutations occurs as strand-coordinated pairs within 32 nucleotides of each other in transcribed regions, suggesting a mechanism in which removal of NQO-DNA adducts by nucleotide excision repair exposes short single-stranded DNA tracts to enzymatic deamination. These highly enriched pairs of APOBEC signature mutations are termed didyma (Greek for twins) and are mechanistically distinct from other types of clustered mutation (omikli and kataegis). Computational analyses of lung and head & neck tumor genomes show that both APOBEC mutagenesis and didyma are elevated in cancers from smokers compared to non-smokers. APOBEC signature mutations and didyma are also elevated in normal lung tissues in smokers prior to cancer initiation. Collectively, these results indicate that DNA adducting mutagens in tobacco smoke can amplify DNA damage and mutagenesis by endogenous APOBEC enzymes and, more broadly, suggest that mutational mechanisms can interact synergistically in both cancer initiation and promotion.

Sherlock-Genome: an R Shiny application for genomic analysis and visualization

MOTIVATION

Next-generation sequencing technologies, such as whole genome sequencing (WGS), have become prominent in cancer genomics. However, managing, visualizing, and integratively analyzing WGS results across various bioinformatic pipelines remains challenging, particularly for non-bioinformaticians, hindering the usability of WGS data for biological discovery.

RESULTS

We developed Sherlock-Genome, an R Shiny app for data harmonization, visualization, and integrative analysis of WGS-based cancer genomics studies. Following FAIR principles, Sherlock-Genome provides a platform and guidelines for managing and sharing finalized sample-level WGS analysis results, enabling users to upload results, inspect analyses locally, and perform integrative analyses. It includes modules for major cancer genomic analyses, allowing interactive data visualizations and integrative analyses with other data types. Sherlock-Genome supports both local and cloud deployment, facilitating the sharing of results for related publications. This tool has the potential to be widely adopted in cancer genomics, significantly enhancing the accessibility and usability of sample-level WGS analysis results for comprehensive biological discovery and research advancements.

AVAILABILITY AND IMPLEMENTATION

The source code and installation instructions for Sherlock-Genome can be accessed via Github https://github.com/xtmgah/Sherlock-Genome . Documentation and data requirements for user project data can also be found on the same GitHub page.

Genomic and transcriptomic landscape of human gastrointestinal stromal tumors

Gastrointestinal stromal tumor (GISTs) are clinically heterogenous exhibiting varying degrees of disease aggressiveness in individual patients. We comprehensively describe the genomic and transcriptomic landscape of a cohort of 117 GISTs including 31 low-risk, 18 intermediate-risk, 29 high-risk, 34 metastatic and 5 neoadjuvant GISTs from 105 patients. GISTs have notably low tumor mutation burden but widespread copy number variations. Aggressive GISTs harbor remarkably more genomic aberrations than low-/intermediate-risk GISTs. Complex genomic alterations, chromothripsis and kataegis, occur selectively in aggressive GISTs. Despite the paucity of mutations, recurrent inactivating YLPM1 mutations are identified (10.3%, 7 of 68 patients), enriched in high-risk/metastatic GIST and functional study further demonstrates YLPM1 inactivation promotes GIST proliferation, growth and oxidative phosphorylation. Spatially and temporally separated GISTs from individual patients demonstrate complex tumor heterogeneity in metastatic GISTs. Finally, four prominent subtypes are proposed with different genomic features, expression profiles, immune characteristics, clinical characteristics and subtype-specific treatment strategies. This large-scale analysis depicts the landscape and provides further insights into GIST pathogenesis and precise treatment.

The interplay of mutagenesis and ecDNA shapes urothelial cancer evolution

Advanced urothelial cancer is a frequently lethal disease characterized by marked genetic heterogeneity1. In this study, we investigated the evolution of genomic signatures caused by endogenous and external mutagenic processes and their interplay with complex structural variants (SVs). We superimposed mutational signatures and phylogenetic analyses of matched serial tumours from patients with urothelial cancer to define the evolutionary dynamics of these processes. We show that APOBEC3-induced mutations are clonal and early, whereas chemotherapy induces mutational bursts of hundreds of late subclonal mutations. Using a genome graph computational tool2, we observed frequent high copy-number circular amplicons characteristic of extrachromosomal DNA (ecDNA)-forming SVs. We characterized the distinct temporal patterns of APOBEC3-induced and chemotherapy-induced mutations within ecDNA-forming SVs, gaining new insights into the timing of these mutagenic processes relative to ecDNA biogenesis. We discovered that most CCND1 amplifications in urothelial cancer arise within circular ecDNA-forming SVs. ecDNA-forming SVs persisted and increased in complexity, incorporating additional DNA segments and contributing to the evolution of treatment resistance. Oxford Nanopore Technologies long-read whole-genome sequencing followed by de novo assembly mapped out CCND1 ecDNA structure. Experimental modelling of CCND1 ecDNA confirmed its role as a driver of treatment resistance. Our findings define fundamental mechanisms that drive urothelial cancer evolution and have important therapeutic implications.

Mapping extrachromosomal DNA amplifications during cancer progression

To understand the role of extrachromosomal DNA (ecDNA) amplifications in cancer progression, we detected and classified focal amplifications in 8,060 newly diagnosed primary cancers, untreated metastases and heavily pretreated tumors. The ecDNAs were detected at significantly higher frequency in untreated metastatic and pretreated tumors compared to newly diagnosed cancers. Tumors from chemotherapy-pretreated patients showed significantly higher ecDNA frequency compared to untreated cancers. In particular, tubulin inhibition associated with ecDNA increases, suggesting a role for ecDNA in treatment response. In longitudinally matched tumor samples, ecDNAs were more likely to be retained compared to chromosomal amplifications. EcDNAs shared between time points, and ecDNAs in advanced cancers were more likely to harbor localized hypermutation events compared to private ecDNAs and ecDNAs in newly diagnosed tumors. Relatively high variant allele fractions of ecDNA localized hypermutations implicated early ecDNA mutagenesis. Our findings nominate ecDNAs to provide tumors with competitive advantages during cancer progression and metastasis.

Widespread mutagenesis and chromosomal instability shape somatic genomes in systemic sclerosis

Systemic sclerosis is a connective tissue disorder characterized by excessive fibrosis that primarily affects women, and can present as a multisystem pathology. Roughly 4-22% of patients with systemic sclerosis develop cancer, which drastically worsens prognosis. However, the mechanisms underlying systemic sclerosis initiation, propagation, and cancer development are poorly understood. We hypothesize that the inflammation and immune response associated with systemic sclerosis can trigger DNA damage, leading to elevated somatic mutagenesis, a hallmark of pre-cancerous tissues. To test our hypothesis, we culture clonal lineages of fibroblasts from the lung tissues of controls and systemic sclerosis patients and compare their mutation burdens and spectra. We find an overall increase in all major mutation types in systemic sclerosis samples compared to control lung samples, from small-scale events such as single base substitutions and insertions/deletions, to chromosome-level changes, including copy-number changes and structural variants. In the genomes of patients with systemic sclerosis, we find evidence of somatic hypermutation or kategis (typically only seen in cancer genomes), we identify mutation signatures closely resembling the error-prone translesion polymerase Polη activity, and observe an activation-induced deaminase-like mutation signature, which overlaps with genomic regions displaying kataegis.

Exploring the Role of Clustered Mutations in Carcinogenesis and Their Potential Clinical Implications in Cancer

Abnormal cell proliferation and growth leading to cancer primarily result from cumulative genome mutations. Single gene mutations alone do not fully explain cancer onset and progression; instead, clustered mutations-simultaneous occurrences of multiple mutations-are considered to be pivotal in cancer development and advancement. These mutations can affect different genes and pathways, resulting in cells undergoing malignant transformation with multiple functional abnormalities. Clustered mutations influence cancer growth rates, metastatic potential, and drug treatment sensitivity. This summary highlights the various types and characteristics of clustered mutations to understand their associations with carcinogenesis and discusses their potential clinical significance in cancer. As a unique mutation type, clustered mutations may involve genomic instability, DNA repair mechanism defects, and environmental exposures, potentially correlating with responsiveness to immunotherapy. Understanding the characteristics and underlying processes of clustered mutations enhances our comprehension of carcinogenesis and cancer progression, providing new diagnostic and therapeutic approaches for cancer.

Mutational impact of APOBEC3A and APOBEC3B in a human cell line and comparisons to breast cancer

A prominent source of mutation in cancer is single-stranded DNA cytosine deamination by cellular APOBEC3 enzymes, which results in signature C-to-T and C-to-G mutations in TCA and TCT motifs. Although multiple enzymes have been implicated, reports conflict and it is unclear which protein(s) are responsible. Here we report the development of a selectable system to quantify genome mutation and demonstrate its utility by comparing the mutagenic activities of three leading candidates-APOBEC3A, APOBEC3B, and APOBEC3H. The human cell line, HAP1, is engineered to express the thymidine kinase (TK) gene of HSV-1, which confers sensitivity to ganciclovir. Expression of APOBEC3A and APOBEC3B, but not catalytic mutant controls or APOBEC3H, triggers increased frequencies of TK mutation and similar TC-biased cytosine mutation profiles in the selectable TK reporter gene. Whole genome sequences from independent clones enabled an analysis of thousands of single base substitution mutations and extraction of local sequence preferences with APOBEC3A preferring YTCW motifs 70% of the time and APOBEC3B 50% of the time (Y = C/T; W = A/T). Signature comparisons with breast tumor whole genome sequences indicate that most malignancies manifest intermediate percentages of APOBEC3 signature mutations in YTCW motifs, mostly between 50 and 70%, suggesting that both enzymes contribute in a combinatorial manner to the overall mutation landscape. Although the vast majority of APOBEC3A- and APOBEC3B-induced single base substitution mutations occur outside of predicted chromosomal DNA hairpin structures, whole genome sequence analyses and supporting biochemical studies also indicate that both enzymes are capable of deaminating the single-stranded loop regions of DNA hairpins at elevated rates. These studies combine to help resolve a long-standing etiologic debate on the source of APOBEC3 signature mutations in cancer and indicate that future diagnostic and therapeutic efforts should focus on both APOBEC3A and APOBEC3B.

APOBEC3B regulates R-loops and promotes transcription-associated mutagenesis in cancer

The single-stranded DNA cytosine-to-uracil deaminase APOBEC3B is an antiviral protein implicated in cancer. However, its substrates in cells are not fully delineated. Here APOBEC3B proteomics reveal interactions with a surprising number of R-loop factors. Biochemical experiments show APOBEC3B binding to R-loops in cells and in vitro. Genetic experiments demonstrate R-loop increases in cells lacking APOBEC3B and decreases in cells overexpressing APOBEC3B. Genome-wide analyses show major changes in the overall landscape of physiological and stimulus-induced R-loops with thousands of differentially altered regions, as well as binding of APOBEC3B to many of these sites. APOBEC3 mutagenesis impacts genes overexpressed in tumors and splice factor mutant tumors preferentially, and APOBEC3-attributed kataegis are enriched in RTCW motifs consistent with APOBEC3B deamination. Taken together with the fact that APOBEC3B binds single-stranded DNA and RNA and preferentially deaminates DNA, these results support a mechanism in which APOBEC3B regulates R-loops and contributes to R-loop mutagenesis in cancer.

Topography of mutational signatures in human cancer

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

Assigning mutational signatures to individual samples and individual somatic mutations with SigProfilerAssignment

Analysis of mutational signatures is a powerful approach for understanding the mutagenic processes that have shaped the evolution of a cancer genome. Here we present SigProfilerAssignment, a desktop and an online computational framework for assigning all types of mutational signatures to individual samples. SigProfilerAssignment is the first tool that allows both analysis of copy-number signatures and probabilistic assignment of signatures to individual somatic mutations. As its computational engine, the tool uses a custom implementation of the forward stagewise algorithm for sparse regression and nonnegative least squares for numerical optimization. Analysis of 2,700 synthetic cancer genomes with and without noise demonstrates that SigProfilerAssignment outperforms four commonly used approaches for assigning mutational signatures. SigProfilerAssignment is freely available at https://github.com/AlexandrovLab/SigProfilerAssignment with a web implementation at https://cancer.sanger.ac.uk/signatures/assignment/.

Oncogenic structural aberration landscape in gastric cancer genomes

Structural variants (SVs) are responsible for driver events in gastric cancer (GC); however, their patterns and processes remain poorly understood. Here, we examine 170 GC whole genomes to unravel the oncogenic structural aberration landscape in GC genomes and identify six rearrangement signatures (RSs). Non-random combinations of RSs elucidate distinctive GC subtypes comprising one or a few dominant RS that are associated with specific driver events (BRCA1/2 defects, mismatch repair deficiency, and TP53 mutation) and epidemiological backgrounds. Twenty-seven SV hotspots are identified as GC driver candidates. SV hotspots frequently constitute complexly clustered SVs involved in driver gene amplification, such as ERBB2, CCNE1, and FGFR2. Further deconstruction of the locally clustered SVs uncovers amplicon-generating profiles characterized by super-large SVs and intensive segmental amplifications, contributing to the extensive amplification of GC oncogenes. Comprehensive analyses using adjusted SV allele frequencies indicate the significant involvement of extra-chromosomal DNA in processes linked to specific RSs.

Distinguishing Curable from Progressive Dementias for Defining Cancer Care Options

The likelihood of a diagnosis of dementia increases with a person's age, as is also the case for many cancers, including melanoma and multiple myeloma, where the median age of diagnosis is above 60 years. However, patients diagnosed with dementia are less likely to be offered invasive curative therapies for cancer. Together with analysis of diet and medication history, advanced imaging methods and genetic profiling can now indicate more about syndromes causing the neurological symptoms. Cachexia, malnutrition, dehydration, alcohol consumption, and even loneliness can all accentuate or cause the "3Ds" of dementia, delirium and depression. Many common drugs, especially in the context of polypharmacy, can cause cognitive difficulties resembling neurodegenerative disease. These syndromes may be reversed by diet, social and caregiver changes, and stopping potentially inappropriate medications (PIMs). More insidious are immune reactions to many different autoantigens, some of which are related to cancers and tumors. These can induce movement and cognitive difficulties that mimic Alzheimer's and Parkinson's diseases and other ataxias associated with aging. Paraneoplastic neurological syndromes may be reversed by directed immunotherapies if detected in their early stages but are best treated by removal of the causative tumor. A full genetic workup should be done for all individuals as soon as possible after diagnosis, to guide less invasive treatments suitable for frail individuals. While surgical interventions may be contraindicated, genetic profile guided immunotherapies, oral treatments, and radiation may be equally curative in a significant number of cancers.

DNA damage and somatic mutations in mammalian cells after irradiation with a nail polish dryer

Ultraviolet A light is commonly emitted by UV-nail polish dryers with recent reports suggesting that long-term use may increase the risk for developing skin cancer. However, no experimental evaluation has been conducted to reveal the effect of radiation emitted by UV-nail polish dryers on mammalian cells. Here, we show that irradiation by a UV-nail polish dryer causes high levels of reactive oxygen species, consistent with 8-oxo-7,8-dihydroguanine damage and mitochondrial dysfunction. Analysis of somatic mutations reveals a dose-dependent increase of C:G>A:T substitutions in irradiated samples with mutagenic patterns similar to mutational signatures previously attributed to reactive oxygen species. In summary, this study demonstrates that radiation emitted by UV-nail polish dryers can both damage DNA and permanently engrave mutations on the genomes of primary mouse embryonic fibroblasts, human foreskin fibroblasts, and human epidermal keratinocytes.

Katdetectr: an R/bioconductor package utilizing unsupervised changepoint analysis for robust kataegis detection

BACKGROUND

Kataegis refers to the occurrence of regional genomic hypermutation in cancer and is a phenomenon that has been observed in a wide range of malignancies. A kataegis locus constitutes a genomic region with a high mutation rate (i.e., a higher frequency of closely interspersed somatic variants than the overall mutational background). It has been shown that kataegis is of biological significance and possibly clinically relevant. Therefore, an accurate and robust workflow for kataegis detection is paramount.

FINDINGS

Here we present Katdetectr, an open-source R/Bioconductor-based package for the robust yet flexible and fast detection of kataegis loci in genomic data. In addition, Katdetectr houses functionalities to characterize and visualize kataegis and provides results in a standardized format useful for subsequent analysis. In brief, Katdetectr imports industry-standard formats (MAF, VCF, and VRanges), determines the intermutation distance of the genomic variants, and performs unsupervised changepoint analysis utilizing the Pruned Exact Linear Time search algorithm followed by kataegis calling according to user-defined parameters.We used synthetic data and an a priori labeled pan-cancer dataset of whole-genome sequenced malignancies for the performance evaluation of Katdetectr and 5 publicly available kataegis detection packages. Our performance evaluation shows that Katdetectr is robust regarding tumor mutational burden and shows the fastest mean computation time. Additionally, Katdetectr reveals the highest accuracy (0.99, 0.99) and normalized Matthews correlation coefficient (0.98, 0.92) of all evaluated tools for both datasets.

CONCLUSIONS

Katdetectr is a robust workflow for the detection, characterization, and visualization of kataegis and is available on Bioconductor: https://doi.org/doi:10.18129/B9.bioc.katdetectr.

Whole-exome sequencing identified mutational profiles of urothelial carcinoma post kidney transplantation

Kidney transplantation is a lifesaving option for patients with end-stage kidney disease. In Taiwan, urothelial carcinoma (UC) is the most common de novo cancer after kidney transplantation (KT). UC has a greater degree of molecular heterogeneity than do other solid tumors. Few studies have explored genomic alterations in UC after KT. We performed whole-exome sequencing to compare the genetic alterations in UC developed after kidney transplantation (UCKT) and in UC in patients on hemodialysis (UCHD). After mapping and variant calling, 18,733 and 11,093 variants were identified in patients with UCKT and UCHD, respectively. We excluded known single-nucleotide polymorphisms (SNPs) and retained genes that were annotated in the Catalogue of Somatic Mutations in Cancer (COSMIC), in the Integrative Onco Genomic cancer mutations browser (IntOGen), and in the Cancer Genome Atlas (TCGA) database of genes associated with bladder cancer. A total of 14 UCKT-specific genes with SNPs identified in more than two patients were included in further analyses. The single-base substitution (SBS) profile and signatures showed a relative high T > A pattern compared to COMSIC UC mutations. Ingenuity pathway analysis was used to explore the connections among these genes. GNAQ, IKZF1, and NTRK3 were identified as potentially involved in the signaling network of UCKT. The genetic analysis of posttransplant malignancies may elucidate a fundamental aspect of the molecular pathogenesis of UCKT.