Uncovering novel mutational signatures by de novo extraction with SigProfilerExtractor

A review on trends in development and translation of omics signatures in cancer

The field of cancer genomics and transcriptomics has evolved from targeted profiling to swift sequencing of individual tumor genome and transcriptome. The steady growth in genome, epigenome, and transcriptome datasets on a genome-wide scale has significantly increased our capability in capturing signatures that represent both the intrinsic and extrinsic biological features of tumors. These biological differences can help in precise molecular subtyping of cancer, predicting tumor progression, metastatic potential, and resistance to therapeutic agents. In this review, we summarized the current development of genomic, methylomic, transcriptomic, proteomic and metabolic signatures in the field of cancer research and highlighted their potentials in clinical applications to improve diagnosis, prognosis, and treatment decision in cancer patients.

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.

Frequency and spectrum of mutations in human sperm measured using duplex sequencing correlate with trio-based de novo mutation analyses

De novo mutations (DNMs) are drivers of genetic disorders. However, the study of DNMs is hampered by technological limitations preventing accurate quantification of ultra-rare mutations. Duplex Sequencing (DS) theoretically has < 1 error/billion base-pairs (bp). To determine the DS utility to quantify and characterize DNMs, we analyzed DNA from blood and spermatozoa from six healthy, 18-year-old Swedish men using the TwinStrand DS mutagenesis panel (48 kb spanning 20 genic and intergenic loci). The mean single nucleotide variant mutation frequency (MF) was 1.2 × 10- 7 per bp in blood and 2.5 × 10- 8 per bp in sperm, with the most common base substitution being C > T. Blood MF and substitution spectrum were similar to those reported in blood cells with an orthogonal method. The sperm MF was in the same order of magnitude and had a strikingly similar spectrum to DNMs from publicly available whole genome sequencing data from human pedigrees (1.2 × 10- 8 per bp). DS revealed much larger numbers of insertions and deletions in sperm over blood, driven by an abundance of putative extra-chromosomal circular DNAs. The study indicates the strong potential of DS to characterize human DNMs to inform factors that contribute to disease susceptibility and heritable genetic risks.

A maternal germline mutator phenotype in a family affected by heritable colorectal cancer

Variation in DNA repair genes can increase cancer risk by elevating the rate of oncogenic mutation. Defects in one such gene, MUTYH, are known to elevate the incidence of colorectal cancer in a recessive Mendelian manner. Recent evidence has also linked MUTYH to a mutator phenotype affecting normal somatic cells as well as the female germline. Here, we use whole genome sequencing to measure germline de novo mutation rates in a large extended family containing both mothers and fathers who are affected by pathogenic MUTYH variation. By developing novel methodology that uses siblings as "surrogate parents" to identify de novo mutations, we were able to include mutation data from several children whose parents were unavailable for sequencing. In the children of mothers affected by the pathogenic MUTYH genotype p.Y179C/V234M, we identify an elevation of the C>A mutation rate that is weaker than mutator effects previously reported to be caused by other pathogenic MUTYH genotypes, suggesting that mutation rates in normal tissues may be useful for classifying cancer-associated variation along a continuum of severity. Surprisingly, we detect no significant elevation of the C>A mutation rate in children born to a father with the same MUTYH genotype, and we similarly find that the mutator effect of the mouse homolog Mutyh appears to be localized to embryonic development, not the spermatocytes. Our results suggest that maternal MUTYH variants can cause germline mutations by attenuating the repair of oxidative DNA damage in the early embryo.

Aristolochic acid-related renal cell carcinoma exhibits a distinct tumor-immune microenvironment favoring response to immune checkpoint blockade

Immune checkpoint blockade (ICB) is currently the standard of care for metastatic renal cell carcinoma (RCC), but treatment responses remain unpredictable. Aristolochic acid (AA), a prevalent supplement additive in Taiwan, has been associated with RCC and induces signature mutations, although its effect on the tumor-immune microenvironment (TIME) is unclear. We aimed to investigate the immune profile of AA-positive RCCs and explore its potential role as a susceptible candidate for ICB. Tissue samples from 22 patients with clear cell RCC (ccRCC) were collected for whole-exome sequencing to determine the genetic features and AA mutational signature (the discovery cohort). The corresponding RNA was sent for NanoString PanCancer IO 360 gene expression analysis to explore the immunological features. The formalin-fixed, parafilm-embedded slides of ccRCCs were sent for multiplex immunohistochemistry/immunofluorescence stain using Vectra system to evaluate the TIME. Tissues from two patients with metastatic RCC demonstrating complete response to ICB were sent for studies to validate the findings (the index patients). The results showed that AA mutational signatures with high tumor mutational burden (TMB) were present in 31.81% of the tumors in the discovery cohort. Three distinct clusters were observed through NanoString analysis. Clusters 1 and 3 were composed mainly of AA-positive RCCs. Cluster 3 RCCs exhibited higher tumor inflammation signature scores and higher immune cell type scores. Vectra analysis revealed a higher percentage of CD15+ and BATF3+ cells in cluster 1, whereas the percentage of CD8+ cells was potentially higher in cluster 3. Strong AA mutational signatures were found in the tumors of two index patients, and both were grouped to cluster 3. In conclusion, AA may induce higher TMB and alter the immune microenvironment in RCCs, which makes the tumors more susceptible to ICB. © 2024 The Author(s). The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Impaired neutrophil-mediated cell death drives Ewing's Sarcoma in the background of Down syndrome

Introduction

Ewing Sarcoma (EWS) has been reported in seven children with Down syndrome (DS). To date, a detailed assessment of this solid tumour in DS patients is yet to be made.

Methods

Here, we characterise a chemo-resistant mediastinal EWS in a 2-year-old DS child, the youngest ever reported case, by exploiting sequencing approaches.

Results

The tumour showed a neuroectodermal development driven by the EWSR1-FLI1 fusion. The inherited myeloperoxidase deficiency of the patient caused failure of neutrophil-mediated cell death and promoted genomic instability.

Discussion

In this context, the tumour underwent genome-wide near haploidisation resulting in a massive overexpression of pro-inflammatory cytokines. Recruitment of defective neutrophils fostered rapid evolution of this EWS.

Ancestry and somatic profile predict acral melanoma origin and prognosis

Acral melanoma, which is not ultraviolet (UV)-associated, is the most common type of melanoma in several low- and middle-income countries including Mexico. Latin American samples are significantly underrepresented in global cancer genomics studies, which directly affects patients in these regions as it is known that cancer risk and incidence may be influenced by ancestry and environmental exposures. To address this, here we characterise the genome and transcriptome of 128 acral melanoma tumours from 96 Mexican patients, a population notable because of its genetic admixture. Compared with other studies of melanoma, we found fewer frequent mutations in classical driver genes such as BRAF , NRAS or NF1 . While most patients had predominantly Amerindian genetic ancestry, those with higher European ancestry had increased frequency of BRAF mutations and a lower number of structural variants. These BRAF -mutated tumours have a transcriptional profile similar to cutaneous non-volar melanocytes, suggesting that acral melanomas in these patients may arise from a distinct cell of origin compared to other tumours arising in these locations. KIT mutations were found in a subset of these tumours, and transcriptional profiling defined three expression clusters; these characteristics were associated with overall survival. We highlight novel low-frequency drivers, such as SPHKAP , which correlate with a distinct genomic profile and clinical characteristics. Our study enhances knowledge of this understudied disease and underscores the importance of including samples from diverse ancestries in cancer genomics studies.

Molecular profiling of BRAF-V600E-mutant metastatic colorectal cancer in the phase 3 BEACON CRC trial

The BEACON CRC study demonstrated that encorafenib (Enco)+cetuximab (Cetux)±binimetinib (Bini) significantly improved overall survival (OS) versus Cetux + chemotherapy in previously treated patients with BRAF-V600E-mutant mCRC, providing the basis for the approval of the Enco+Cetux regimen in the United States and the European Union. A greater understanding of biomarkers predictive of response to Enco+Cetux±Bini treatment is of clinical relevance. In this prespecified, exploratory biomarker analysis of the BEACON CRC study, we characterize genomic and transcriptomic correlates of clinical outcomes and acquired resistance mechanisms through integrated clinical and molecular analysis, including whole-exome and -transcriptome tissue sequencing and circulating tumor DNA genomic profiling. Tumors with higher immune signatures showed a trend towards increased OS benefit with Enco+Bini+Cetux. RAS, MAP2K1 and MET alterations were most commonly acquired with Enco+Cetux±Bini, and more frequent in patients with a high baseline cell-cycle gene signature; baseline TP53 mutation was associated with acquired MET amplification. Acquired mutations were subclonal and polyclonal, with evidence of increased tumor mutation rate with Enco+Cetux±Bini and mutational signatures (SBS17a/b). These findings support treatment with Enco+Cetux±Bini for patients with BRAF-V600E-mutant mCRC and provide insights into the biology of response and resistance to MAPK-pathway-targeted therapy. ClinicalTrials.gov registration: NCT02928224.

Unraveling the key role of chromatin structure in cancer development through epigenetic landscape characterization of oral cancer

Epigenetic alterations, such as those in chromatin structure and DNA methylation, have been extensively studied in a number of tumor types. But oral cancer, particularly oral adenocarcinoma, has received far less attention. Here, we combined laser-capture microdissection and muti-omics mini-bulk sequencing to systematically characterize the epigenetic landscape of oral cancer, including chromatin architecture, DNA methylation, H3K27me3 modification, and gene expression. In carcinogenesis, tumor cells exhibit reorganized chromatin spatial structures, including compromised compartment structures and altered gene-gene interaction networks. Notably, some structural alterations are observed in phenotypically non-malignant paracancerous but not in normal cells. We developed transformer models to identify the cancer propensity of individual genome loci, thereby determining the carcinogenic status of each sample. Insights into cancer epigenetic landscapes provide evidence that chromatin reorganization is an important hallmark of oral cancer progression, which is also linked with genomic alterations and DNA methylation reprogramming. In particular, regions of frequent copy number alternations in cancer cells are associated with strong spatial insulation in both cancer and normal samples. Aberrant methylation reprogramming in oral squamous cell carcinomas is closely related to chromatin structure and H3K27me3 signals, which are further influenced by intrinsic sequence properties. Our findings indicate that structural changes are both significant and conserved in two distinct types of oral cancer, closely linked to transcriptomic alterations and cancer development. Notably, the structural changes remain markedly evident in oral adenocarcinoma despite the considerably lower incidence of genomic copy number alterations and lesser extent of methylation alterations compared to squamous cell carcinoma. We expect that the comprehensive analysis of epigenetic reprogramming of different types and subtypes of primary oral tumors can provide additional guidance to the design of novel detection and therapy for oral cancer.

RCC1 depletion drives protein transport defects and rupture in micronuclei

Micronuclei (MN) are a commonly used marker of chromosome instability that form when missegregated chromatin recruits its own nuclear envelope (NE) after mitosis. MN frequently rupture, which results in genome instability, upregulation of metastatic genes, and increased immune signaling. MN rupture is linked to NE defects, but the cause of these defects is poorly understood. Previous work from our lab found that chromosome identity correlates with rupture timing for small MN, i.e. MN containing a short chromosome, with more euchromatic chromosomes forming more stable MN with fewer nuclear lamina gaps. Here we demonstrate that histone methylation promotes rupture and nuclear lamina defects in small MN. This correlates with increased MN size, and we go on to find that all MN have a constitutive nuclear export defect that drives MN growth and nuclear lamina gap expansion, making the MN susceptible to rupture. We demonstrate that these export defects arise from decreased RCC1 levels in MN and that additional loss of RCC1 caused by low histone methylation in small euchromatic MN results in additional import defects that suppress nuclear lamina gaps and MN rupture. Through analysis of mutational signatures associated with early and late rupturing chromosomes in the Pan-Cancer Analysis of Whole Genomes (PCAWG) dataset, we identify an enrichment of APOBEC and DNA polymerase E hypermutation signatures in chromothripsis events on early and mid rupturing chromosomes, respectively, suggesting that MN rupture timing could determine the landscape of structural variation in chromothripsis. Our study defines a new model of MN rupture where increased MN growth, caused by defects in protein export, drives gaps in nuclear lamina organization that make the MN susceptible to membrane rupture with long-lasting effects on genome architecture.

The genomic landscape of 2,023 colorectal cancers

Colorectal carcinoma (CRC) is a common cause of mortality1, but a comprehensive description of its genomic landscape is lacking2-9. Here we perform whole-genome sequencing of 2,023 CRC samples from participants in the UK 100,000 Genomes Project, thereby providing a highly detailed somatic mutational landscape of this cancer. Integrated analyses identify more than 250 putative CRC driver genes, many not previously implicated in CRC or other cancers, including several recurrent changes outside the coding genome. We extend the molecular pathways involved in CRC development, define four new common subgroups of microsatellite-stable CRC based on genomic features and show that these groups have independent prognostic associations. We also characterize several rare molecular CRC subgroups, some with potential clinical relevance, including cancers with both microsatellite and chromosomal instability. We demonstrate a spectrum of mutational profiles across the colorectum, which reflect aetiological differences. These include the role of Escherichia colipks+ colibactin in rectal cancers10 and the importance of the SBS93 signature11-13, which suggests that diet or smoking is a risk factor. Immune-escape driver mutations14 are near-ubiquitous in hypermutant tumours and occur in about half of microsatellite-stable CRCs, often in the form of HLA copy number changes. Many driver mutations are actionable, including those associated with rare subgroups (for example, BRCA1 and IDH1), highlighting the role of whole-genome sequencing in optimizing patient care.

The Proteogenomics of Prostate Cancer Radioresistance

Prostate cancer is frequently treated with radiotherapy. Unfortunately, aggressive radioresistant relapses can arise, and the molecular underpinnings of radioresistance are unknown. Modern clinical radiotherapy is evolving to deliver higher doses of radiation in fewer fractions (hypofractionation). We therefore analyzed genomic, transcriptomic, and proteomic data to characterize prostate cancer radioresistance in cells treated with both conventionally fractionated and hypofractionated radiotherapy. Independent of fractionation schedule, resistance to radiotherapy involved massive genomic instability and abrogation of DNA mismatch repair. Specific prostate cancer driver genes were modulated at the RNA and protein levels, with distinct protein subcellular responses to radiotherapy. Conventional fractionation led to a far more aggressive biomolecular response than hypofractionation. Testing preclinical candidates identified in cell lines, we revealed POLQ (DNA Polymerase Theta) as a radiosensitizer. POLQ-modulated radioresistance in model systems and was predictive of it in large patient cohorts. The molecular response to radiation is highly multimodal and sheds light on prostate cancer lethality.

SIGNIFICANCE

Radiation is standard of care in prostate cancer. Yet, we have little understanding of its failure. We demonstrate a new paradigm that radioresistance is fractionation specific and identified POLQ as a radioresistance modulator.

Preclinical evaluation of avutometinib and defactinib in high-grade endometrioid endometrial cancer

BACKGROUND

High-grade endometrial cancers (EAC) are aggressive tumors with a high risk of progression after treatment. As EAC may harbor mutations in the RAS/MAPK pathways, we evaluated the preclinical in vitro and in vivo efficacy of avutometinib, a RAF/MEK clamp, in combination with the focal adhesion kinase (FAK) inhibitors defactinib or VS-4718, against multiple primary EAC cell lines and xenografts.

METHODS

Whole-exome sequencing (WES) was used to evaluate the genetic landscape of five primary EAC cell lines. The in vitro activity of avutometinib and defactinib as single agents and in combination was evaluated using cell viability, cell cycle, and cytotoxicity assays. Mechanistic studies were performed using Western blot assays while in vivo experiments were completed in UTE10 engrafted mice treated with either vehicle, avutometinib, VS-4718, or their combination through oral gavage.

RESULTS

WES results demonstrated multiple EAC cell lines to harbor genetic derangements in the RAS/MAPK pathway including KRAS/PTEN/PIK3CA/BRAF/ARID1A, potentially sensitizing to FAK and RAF/MEK inhibition. Five out of five of the EAC cell lines demonstrated in vitro sensitivity to FAK and/or RAF/MEK inhibition. By Western blot assays, exposure of EAC cell lines to defactinib, avutometinib, and their combination demonstrated decreased phosphorylated FAK (p-FAK) as well as decreased p-MEK and p-ERK. In vivo the combination of avutometinib/VS-4718 demonstrated superior tumor growth inhibition compared to single-agent treatment and controls starting at Day 9 (p < 0.02 and p < 0.04) in UTE10 xenografts.

CONCLUSIONS

Avutometinib, defactinib, and to a larger extent their combinations, demonstrated promising in vitro and in vivo activity against EAC cell lines and xenografts. These preclinical data support the potential clinical evaluation of this combination in high-grade EAC patients.

Prognostic genome and transcriptome signatures in colorectal cancers

Colorectal cancer is caused by a sequence of somatic genomic alterations affecting driver genes in core cancer pathways1. Here, to understand the functional and prognostic impact of cancer-causing somatic mutations, we analysed the whole genomes and transcriptomes of 1,063 primary colorectal cancers in a population-based cohort with long-term follow-up. From the 96 mutated driver genes, 9 were not previously implicated in colorectal cancer and 24 had not been linked to any cancer. Two distinct patterns of pathway co-mutations were observed, timing analyses identified nine early and three late driver gene mutations, and several signatures of colorectal-cancer-specific mutational processes were identified. Mutations in WNT, EGFR and TGFβ pathway genes, the mitochondrial CYB gene and 3 regulatory elements along with 21 copy-number variations and the COSMIC SBS44 signature correlated with survival. Gene expression classification yielded five prognostic subtypes with distinct molecular features, in part explained by underlying genomic alterations. Microsatellite-instable tumours divided into two classes with different levels of hypoxia and infiltration of immune and stromal cells. To our knowledge, this study constitutes the largest integrated genome and transcriptome analysis of colorectal cancer, and interlinks mutations, gene expression and patient outcomes. The identification of prognostic mutations and expression subtypes can guide future efforts to individualize colorectal cancer therapy.

Comprehensive molecular characterization of collecting duct carcinoma for therapeutic vulnerability

Collecting duct carcinoma (CDC) is an aggressive rare subtype of kidney cancer with unmet clinical needs. Little is known about its underlying molecular alterations and etiology, primarily due to its rarity, and lack of preclinical models. This study aims to comprehensively characterize molecular alterations in CDC and identify its therapeutic vulnerabilities. Through whole-exome and transcriptome sequencing, we identified KRAS hotspot mutations (G12A/D/V) in 3/13 (23%) of the patients, in addition to known TP53, NF2 mutations. 3/13 (23%) patients carried a mutational signature (SBS22) caused by aristolochic acid (AA) exposures, known to be more prevalent in Asia, highlighting a geologically specific disease etiology. We further discovered that cell cycle-related pathways were the most predominantly dysregulated pathways. Our drug screening with our newly established CDC preclinical models identified a CDK9 inhibitor LDC000067 that specifically inhibited CDC tumor growth and prolonged survival. Our study not only improved our understanding of oncogenic molecular alterations of Asian CDC, but also identified cell-cycle machinery as a therapeutic vulnerability, laying the foundation for clinical trials to treat patients with such aggressive cancer.

Somatic mutations in tumor-infiltrating lymphocytes impact on antitumor immunity

Immune checkpoint inhibitors (ICIs) exert clinical efficacy against various types of cancers by reinvigorating exhausted CD8+ T cells that can expand and directly attack cancer cells (cancer-specific T cells) among tumor-infiltrating lymphocytes (TILs). Although some reports have identified somatic mutations in TILs, their effect on antitumor immunity remains unclear. In this study, we successfully established 18 cancer-specific T cell clones, which have an exhaustion phenotype, from the TILs of four patients with melanoma. We conducted whole-genome sequencing for these T cell clones and identified various somatic mutations in them with high clonality. Among the somatic mutations, an SH2D2A loss-of-function frameshift mutation and TNFAIP3 deletion could activate T cell effector functions in vitro. Furthermore, we generated CD8+ T cell-specific Tnfaip3 knockout mice and showed that Tnfaip3 function loss in CD8+ T cell increased antitumor immunity, leading to remarkable response to PD-1 blockade in vivo. In addition, we analyzed bulk CD3+ T cells from TILs in additional 12 patients and identified an SH2D2A mutation in one patient through amplicon sequencing. These findings suggest that somatic mutations in TILs can affect antitumor immunity and suggest unique biomarkers and therapeutic targets.

Lynch syndrome-associated and sporadic microsatellite unstable colorectal cancers: different patterns of clonal evolution yield highly similar tumours

Microsatellite unstable colorectal cancer (MSI-CRC) can arise through germline mutations in mismatch repair (MMR) genes in individuals with Lynch syndrome (LS), or sporadically through promoter methylation of the MMR gene MLH1. Despite the different origins of hereditary and sporadic MSI tumours, their genomic features have not been extensively compared. A prominent feature of MMR-deficient genomes is the occurrence of many indels in short repeat sequences, an understudied mutation type due to the technical challenges of variant calling in these regions. In this study, we performed whole genome sequencing and RNA-sequencing on 29 sporadic and 14 hereditary MSI-CRCs. We compared the tumour groups by analysing genome-wide mutation densities, microsatellite repeat indels, recurrent protein-coding variants, signatures of single base, doublet base, and indel mutations, and changes in gene expression. We show that the mutational landscapes of hereditary and sporadic MSI-CRCs, including mutational signatures and mutation densities genome-wide and in microsatellites, are highly similar. Only a low number of differentially expressed genes were found, enriched to interferon-γ regulated immune response pathways. Analysis of the variance in allelic fractions of somatic variants in each tumour group revealed higher clonal heterogeneity in sporadic MSI-CRCs. Our results suggest that the differing molecular origins of MMR deficiency in hereditary and sporadic MSI-CRCs do not result in substantial differences in the mutational landscapes of these tumours. The divergent patterns of clonal evolution between the tumour groups may have clinical implications, as high clonal heterogeneity has been associated with decreased tumour immunosurveillance and reduced responsiveness to immunotherapy.

Divergent Evolution in Bilateral Prostate Cancer: a Case Study

Multifocal prostate cancer is a prevalent phenomenon, with most cases remaining uncharacterized from a genomic perspective. A patient presented with bilateral prostate cancer. On systematic biopsy, two indistinguishable clinicopathologic lesions were detected. Whole-genome sequencing displayed somatically unrelated tumours with distinct driver CNA regions, suggesting independent origins of the two tumors. We demonstrated that similar clinicopathologic multifocal tumours, which might be interpreted as clonal disease, can in fact represent independent cancers. Genetic prognostics can prevent mischaracterization of multifocal disease to enable optimal patient management.

Motive and opportunity: MYC rearrangements in high-grade B-cell lymphoma with MYC and BCL2 rearrangements (an LLMPP study)

ABSTRACT

Rearrangements that place the oncogenes MYC, BCL2, or BCL6 adjacent to superenhancers are common in mature B-cell lymphomas. Lymphomas with diffuse large B-cell lymphoma (DLBCL) or high-grade morphology with both MYC and BCL2 rearrangements are classified as high-grade B-cell lymphoma with MYC and BCL2 rearrangements ("double hit"; HGBCL-DH-BCL2) and are associated with aggressive disease and poor outcomes. Although it is established that MYC rearrangements involving immunoglobulin (IG) loci are associated with inferior outcomes relative to those involving other non-IG superenhancers, the frequency of and mechanisms driving IG vs non-IG MYC rearrangements have not been elucidated. Here, we used custom targeted capture and/or whole-genome sequencing to characterize oncogene rearrangements across 883 mature B-cell lymphomas including Burkitt lymphoma, follicular lymphoma, DLBCL, and HGBCL-DH-BCL2 tumors. We demonstrate that, although BCL2 rearrangement topology is consistent across entities, HGBCL-DH-BCL2 have distinct MYC rearrangement architecture relative to tumors with single MYC rearrangements or with both MYC and BCL6 rearrangements (HGBCL-DH-BCL6), including both a higher frequency of non-IG rearrangements and different architecture of MYC::IGH rearrangements. The distinct MYC rearrangement patterns in HGBCL-DH-BCL2 occur on the background of high levels of somatic hypermutation across MYC partner loci in HGBCL-DH-BCL2, creating more opportunity to form these rearrangements. Furthermore, because 1 IGH allele is already disrupted by the existing BCL2 rearrangement, the MYC rearrangement architecture in HGBCL-DH-BCL2 likely reflects selective pressure to preserve both BCL2 and B-cell receptor expression. These data provide new mechanistic explanations for the distinct patterns of MYC rearrangements observed across different lymphoma entities.

The DNA Methyltransferase Inhibitor 5-Aza-4'-thio-2'-Deoxycytidine Induces C>G Transversions and Acute Lymphoid Leukemia Development

DNA methyltransferase inhibitors (DNMTi), most commonly cytidine analogs, are compounds that decrease 5'-cytosine methylation. DNMTi are used clinically based on the hypothesis that cytosine demethylation will lead to re-expression of tumor suppressor genes. 5-Aza-4'-thio-2'-deoxycytidine (Aza-TdCyd or ATC) is a recently described thiol-substituted DNMTi that has been shown to have anti-tumor activity in solid tumor models. In this study, we investigated the therapeutic potential of ATC in a murine transplantation model of myelodysplastic syndrome. ATC treatment led to the transformation of transplanted wild-type bone marrow nucleated cells into lymphoid leukemia, and healthy mice treated with ATC also developed lymphoid leukemia. Whole-exome sequencing revealed 1,000 acquired mutations, almost all of which were C>G transversions in a specific 5'-NCG-3' context. These mutations involved dozens of genes involved in human lymphoid leukemia, such as Notch1, Pten, Pax5, Trp53, and Nf1. Human cells treated in vitro with ATC showed 1,000 acquired C>G transversions in a similar context. Deletion of Dck, the rate-limiting enzyme for the cytidine salvage pathway, eliminated C>G transversions. Taken together, these findings demonstrate a highly penetrant mutagenic and leukemogenic phenotype associated with ATC. Significance: Treatment with a DNA methyltransferase inhibitor generates a distinct mutation signature and triggers leukemic transformation, which has important implications for the research and clinical applications of these inhibitors.

Pan-cancer mutational signature analysis of 111,711 targeted sequenced tumors using SATS

Tumor mutational signatures are informative for cancer diagnosis and treatment. However, targeted sequencing, commonly used in clinical settings, lacks specialized analytical tools and a dedicated catalogue of mutational signatures. Here, we introduce SATS, a scalable mutational signature analyzer for targeted sequencing data. SATS leverages tumor mutational burdens to identify and quantify signatures in individual tumors, overcoming the challenges of sparse mutations and variable gene panels. Validations across simulated data, pseudo-targeted sequencing data, and matched whole-genome and targeted sequencing samples show that SATS can accurately detect common mutational signatures and estimate their burdens. Applying SATS to 111,711 tumors from the AACR Project GENIE, we created a pan-cancer mutational signature catalogue specific to targeted sequencing. We further validated signatures in lung, breast and colorectal cancers using an additional 16,774 independent samples. This signature catalogue is a valuable resource for estimating signature burdens in individual targeted sequenced tumors, facilitating the integration of mutational signatures with clinical data.

Dose-related Mutagenic and Clastogenic Effects of Benzo[b]fluoranthene in Mouse Somatic Tissues Detected by Duplex Sequencing and the Micronucleus Assay

Polycyclic aromatic hydrocarbons (PAHs) are common environmental pollutants that originate from the incomplete combustion of organic materials. We investigated the clastogenicity and mutagenicity of benzo[ b ]fluoranthene (BbF), one of 16 priority PAHs, in MutaMouse males after a 28-day oral exposure. BbF causes robust dose-dependent increases in micronucleus frequency in peripheral blood, indicative of chromosome damage. Duplex Sequencing (DS), an error-corrected sequencing technology, reveals that BbF induces dose-dependent increases in mutation frequencies in bone marrow (BM) and liver. Mutagenicity is increased in intergenic relative to genic regions, suggesting a role for transcription-coupled repair of BbF-induced DNA damage. At higher doses, the maximum mutagenic response to BbF is higher in liver, which has a lower mitotic index but higher metabolic capacity than BM; however, mutagenic potency is comparable between the two tissues. BbF induces primarily C:G>A:T mutations, followed by C:G>T:A and C:G>G:C, indicating that BbF metabolites mainly target guanines and cytosines. The mutation spectrum of BbF correlates with cancer mutational signatures associated with tobacco exposure, supporting its contribution to the carcinogenicity of combustion-derived PAHs in humans. Overall, BbF's mutagenic effects are similar to benzo[ a ]pyrene, a well-studied mutagenic PAH. Our work showcases the utility of DS for effective mutagenicity assessment of environmental pollutants.

Synopsis

We used Duplex Sequencing to study the mutagenicity of benzo[ b ]fluoranthene across the mouse genome. Dose-dependent changes in mutation frequency and spectrum quantify its role in PAH-induced carcinogenicity.

Genetic evolution of keratinocytes to cutaneous squamous cell carcinoma

We performed multi-omic profiling of epidermal keratinocytes, precancerous actinic keratoses, and squamous cell carcinomas to understand the molecular transitions during skin carcinogenesis. Single-cell mutational analyses showed that most keratinocytes in normal skin had lower mutation burdens than melanocytes and fibroblasts, however keratinocytes with TP53 or NOTCH1 mutations had substantially higher mutation burdens, suggesting that these mutations prime keratinocytes for transformation by increasing their mutation rate. Mutational profiling and spatial transcriptomics on squamous cell carcinomas adjacent to actinic keratoses revealed TERT promoter and CDKN2A mutations emerging in actinic keratoses, whereas additional mutations inactivating ARID2 and activating the MAPK-pathway delineated the transition to squamous cell carcinomas. Spatial variation in gene expression patterns was common in both tumor and immune cells, with high expression of checkpoint molecules at the invasive front of tumors. In conclusion, this study documents key events during the evolution of cutaneous squamous cell carcinoma.

Proteogenomic insights into the biology and treatment of pan-melanoma

Melanoma is one of the most prevalent skin cancers, with high metastatic rates and poor prognosis. Understanding its molecular pathogenesis is crucial for improving its diagnosis and treatment. Integrated analysis of multi-omics data from 207 treatment-naïve melanomas (primary-cutaneous-melanomas (CM, n = 28), primary-acral-melanomas (AM, n = 81), primary-mucosal-melanomas (MM, n = 28), metastatic-melanomas (n = 27), and nevi (n = 43)) provides insights into melanoma biology. Multivariate analysis reveals that PRKDC amplification is a prognostic molecule for melanomas. Further proteogenomic analysis combined with functional experiments reveals that the cis-effect of PRKDC amplification may lead to tumor proliferation through the activation of DNA repair and folate metabolism pathways. Proteome-based stratification of primary melanomas defines three prognosis-related subtypes, namely, the ECM subtype, angiogenesis subtype (with a high metastasis rate), and cell proliferation subtype, which provides an essential framework for the utilization of specific targeted therapies for particular melanoma subtypes. The immune classification identifies three immune subtypes. Further analysis combined with an independent anti-PD-1 treatment cohort reveals that upregulation of the MAPK7-NFKB signaling pathway may facilitate T-cell recruitment and increase the sensitivity of patients to immunotherapy. In contrast, PRKDC may reduce the sensitivity of melanoma patients to immunotherapy by promoting DNA repair in melanoma cells. These results emphasize the clinical value of multi-omics data and have the potential to improve the understanding of melanoma treatment.

Whole genome sequencing refines stratification and therapy of patients with clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is the most common form of kidney cancer, but a comprehensive description of its genomic landscape is lacking. We report the whole genome sequencing of 778 ccRCC patients enrolled in the 100,000 Genomes Project, providing for a detailed description of the somatic mutational landscape of ccRCC. We identify candidate driver genes, which as well as emphasising the major role of epigenetic regulation in ccRCC highlight additional biological pathways extending opportunities for therapeutic interventions. Genomic characterisation identified patients with divergent clinical outcome; higher number of structural copy number alterations associated with poorer prognosis, whereas VHL mutations were independently associated with a better prognosis. The observations that higher T-cell infiltration is associated with better overall survival and that genetically predicted immune evasion is not common supports the rationale for immunotherapy. These findings should inform personalised surveillance and treatment strategies for ccRCC patients.

The genomic landscape of breast and non-breast cancers from individuals with germline CHEK2 deficiency

CHEK2 is considered to be involved in homologous recombination repair (HRR). Individuals who have germline pathogenic variants (gPVs) in CHEK2 are at increased risk to develop breast cancer and likely other primary cancers. PARP inhibitors (PARPi) have been shown to be effective in the treatment of cancers that present with HRR deficiency-for example, caused by inactivation of BRCA1/2. However, clinical trials have shown little to no efficacy of PARPi in patients with CHEK2 gPVs. Here, we show that both breast and non-breast cancers from individuals who have biallelic gPVs in CHEK2 (germline CHEK2 deficiency) do not present with molecular profiles that fit with HRR deficiency. This finding provides a likely explanation why PARPi therapy is not successful in the treatment of CHEK2-deficient cancers.

Benefits for children with suspected cancer from routine whole-genome sequencing

Clinical whole-genome sequencing (WGS) has been shown to deliver potential benefits to children with cancer and to alter treatment in high-risk patient groups. It remains unknown whether offering WGS to every child with suspected cancer can change patient management. We collected WGS variant calls and clinical and diagnostic information from 281 children (282 tumors) across two English units (n = 152 from a hematology center, n = 130 from a solid tumor center) where WGS had become a routine test. Our key finding was that variants uniquely attributable to WGS changed the management in ~7% (20 out of 282) of cases while providing additional disease-relevant findings, beyond standard-of-care molecular tests, in 108 instances for 83 (29%) cases. Furthermore, WGS faithfully reproduced every standard-of-care molecular test (n = 738) and revealed several previously unknown genomic features of childhood tumors. We show that WGS can be delivered as part of routine clinical care to children with suspected cancer and can change clinical management by delivering unexpected genomic insights. Our experience portrays WGS as a clinically impactful assay for routine practice, providing opportunities for assay consolidation and for delivery of molecularly informed patient care.

Kataegis associated mutational processes linked to adverse prostate cancer presentation in African men

Kataegis, the focal hypermutation of single base substitutions (SBS) in tumour genomes, has received little attention with respect to prostate cancer (PCa) associated molecular and clinical features. Most notably, data is lacking with regards to this tumour evolutionary phenomenon and PCa racial disparities, with African men disproportionately impacted. Here through comparison between African (n = 109) and non-African (n = 79) whole genome sequenced treatment naïve primary tumours, using a single analytical workflow we assessed for shared and unique features of kataegis. Linking kataegis to aggressive presentation, structural variant burden and copy number loss, we attributed APOBEC3 activity through higher rates of SBS2 to high-risk African tumours. While kataegis positive African patients presented with elevated prostate specific antigen levels, their tumours showed evolutionary unique trajectories marked by increased subclonal and structural variant-independent kataegis. The potential to exacerbate tumour heterogeneity emphases the significance of continued exploration of biological behaviours and environmental exposures for African patients.

Human cancer genomes harbor the mutational signature of tobacco-specific nitrosamines NNN and NNK

Tobacco usage is linked to multiple cancer types and accounts for a quarter of all cancer-related deaths. Tobacco smoke contains various carcinogenic compounds, including polycyclic aromatic hydrocarbons (PAH), though the mutagenic potential of many tobacco-related chemicals remains largely unexplored. In particular, the highly carcinogenic tobacco-specific nitrosamines NNN and NNK form pre-mutagenic pyridyloxobutyl (POB) DNA adducts. In the study presented here, we identified genome-scale POB-induced mutational signatures in cell lines and rat tumors, while also investigating their role in human cancer. These signatures are characterized by T>N and C>T mutations forming from specific POB adducts damaging dT and dC residues. Analysis of 2,780 cancer genomes uncovered POB signatures in ∼180 tumors; from cancer types distinct from the ones linked to smoking-related signatures SBS4 and SBS92. This suggests that, unlike PAH compounds, the POB pathway may contribute uniquely to the mutational landscapes of certain hematological malignancies and cancers of the kidney, breast, prostate and pancreas.

Molecular effects of indoor tanning

Background

Tanning bed users have a significantly increased risk of melanoma, but it remains unclear how indoor tanning drives melanomagenesis. Tanning bed radiation is often thought of as a substitute for natural UV radiation despite differences in the maximum doses, UV content, body sites exposed, and patterns of melanoma that arise.

Methods

To better understand the epidemiologic trends and etiology of melanoma associated with tanning bed use, we described the patterns of melanoma in patients with quantifiable tanning bed usage and performed exome sequencing of 182 melanocytes from normal skin of a subset of these patients.

Results

Tanning bed users were more likely than non-users to have melanoma on body sites with low cumulative levels of sun damage and were more likely to have multiple melanomas. The melanocytes in normal appearing skin from tanning bed users had higher mutation burdens, a higher proportion of melanocytes with pathogenic mutations, and distinct mutational signatures. These differences were most prominent over body sites that experience comparatively less exposure to natural sunlight.

Conclusions

We conclude that tanning bed radiation induces melanoma by increasing the mutation burden of melanocytes and by mutagenizing a broader field of melanocytes than are typically exposed to natural sunlight. The unique signatures of mutations in skin cells of tanning users may be attributable to the distinct spectra of radiation emitted from solariums.

Molecular landscape and multi-omic measurements of heterogeneity in fetal adenocarcinoma of the lung

Fetal adenocarcinoma of the lung (FLAC) is a rare form of lung adenocarcinoma and was divided into high-grade (H-FLAC) and low-grade (L-FLAC) subtypes. Despite the existence of some small case series studies, a comprehensive multi-omics study of FLAC has yet to be undertaken. In this study, we depicted the multi-omics landscapes of this rare lung cancer type by performing multi-regional sampling on 20 FLAC cases. A comparison of multi-omics profiles revealed significant differences between H-FLAC and L-FLAC in a multi-omic landscape. Two subtypes also showed distinct relationships between multi-layer intratumor heterogeneity (ITH). We discovered that a lower genetic ITH was significantly associated with worse recurrence-free survival and overall survival in FLAC patients, whereas higher methylation ITH in H-FLAC patients suggested a short survival. Our findings highlight the complex interplay between genetic and transcriptional heterogeneity in FLAC and suggest that different types of ITH may have distinct implications for patient prognosis.

Genomic characterization and risk stratification of esophageal squamous dysplasia

Objectives

The majority of esophageal squamous dysplasia (ESD) patients progress slowly, while a subset of patients can undergo recurrence rapidly or progress to invasive cancer even after proper treatment. However, the molecular mechanisms underlying these clinical observations are still largely unknown.

Methods

By sequencing the genomic data of 160 clinical samples from 49 tumor-free ESD patients and 88 esophageal squamous cell carcinoma (ESCC) patients, we demonstrated lower somatic mutation and copy number alteration (CNA) burden in ESD compared with ESCC.

Results

Cross-species screening and functional assays identified ACSM5 as a novel driver gene for ESD progression. Furthermore, we revealed that miR-4292 promoted ESD progression and could serve as a non-invasive diagnostic marker for ESD.

Conclusions

These findings largely expanded our understanding of ESD genetics and tumorigenesis, which possessed promising significance for improving early diagnosis, reducing overtreatment, and identifying high-risk ESD patients.

Genetic landscape of interval and screen detected breast cancer

Interval breast cancers (IBCs) are cancers diagnosed between screening episodes. Understanding the biological differences between IBCs and screen-detected breast-cancers (SDBCs) has the potential to improve mammographic screening and patient management. We analysed and compared the genomic landscape of 288 IBCs and 473 SDBCs by whole genome sequencing of paired tumour-normal patient samples collected as part of the UK 100,000 Genomes Project. Compared to SDBCs, IBCs were more likely to be lobular, higher grade, and triple negative. A more aggressive clinical phenotype was reflected in IBCs displaying features of genomic instability including a higher mutation rate and number of chromosomal structural abnormalities, defective homologous recombination and TP53 mutations. We did not however, find evidence to indicate that IBCs are associated with a significantly different immune response. While IBCs do not represent a unique molecular class of invasive breast cancer they exhibit a more aggressive phenotype, which is likely to be a consequence of the timing of tumour initiation. This information is relevant both with respect to treatment as well as informing the screening interval for mammography.

Constitutional mismatch repair deficiency mimicking Lynch syndrome is associated with hypomorphic mismatch repair gene variants

Lynch syndrome (LS) and constitutional mismatch repair deficiency (CMMRD) are distinct cancer syndromes caused, respectively, by mono- and bi-allelic germline mismatch repair (MMR) variants. LS predisposes to mainly gastrointestinal and genitourinary cancers in adulthood. CMMRD predisposes to brain, haematological, and LS-spectrum cancers from childhood. Two suspected LS patients with first cancer diagnosis aged 27 or 38 years were found to be homozygous for an MMR (likely) pathogenic variant, MSH6 c.3226C>T (p.(Arg1076Cys)), or variant of uncertain significance (VUS), MLH1 c.306G>A (p.(Glu102=)). MLH1 c.306G>A was shown to cause leaky exon 3 skipping. The apparent genotype-phenotype conflict was resolved by detection of constitutional microsatellite instability in both patients, a hallmark feature of CMMRD. A hypomorphic effect of these and other variants found in additional late onset CMMRD cases, identified by literature review, likely explains a LS-like phenotype. CMMRD testing in carriers of compound heterozygous or homozygous MMR VUS may find similar cases and novel hypomorphic variants. Individualised management of mono- and bi-allelic carriers of hypomorphic MMR variants is needed until we better characterise the associated phenotypes.

Clonal evolution of hematopoietic stem cells after cancer chemotherapy

Normal hematopoietic stem and progenitor cells (HSPCs) inherently accumulate somatic mutations and lose clonal diversity with age, processes implicated in the development of myeloid malignancies 1 . The impact of exogenous stressors, such as cancer chemotherapies, on the genomic integrity and clonal dynamics of normal HSPCs is not well defined. We conducted whole-genome sequencing on 1,032 single-cell-derived HSPC colonies from 10 patients with multiple myeloma (MM), who had undergone various chemotherapy regimens. Our findings reveal that melphalan treatment distinctly increases mutational burden with a unique mutation signature, whereas other MM chemotherapies do not significantly affect the normal mutation rate of HSPCs. Among these therapy-induced mutations were several oncogenic drivers such as TET2 and PPM1D . Phylogenetic analysis showed a clonal architecture in post-treatment HSPCs characterized by extensive convergent evolution of mutations in genes such as TP53 and PPM1D . Consequently, the clonal diversity and structure of post-treatment HSPCs mirror those observed in normal elderly individuals, suggesting an accelerated clonal aging due to chemotherapy. Furthermore, analysis of matched therapy-related myeloid neoplasm (t-MN) samples, which occurred 1-8 years later, enabled us to trace the clonal origin of t-MNs to a single HSPC clone among a group of clones with competing malignant potential, indicating the critical role of secondary mutations in dictating clonal dominance and malignant transformation. Our findings suggest that cancer chemotherapy promotes an oligoclonal architecture with multiple HSPC clones possessing competing leukemic potentials, setting the stage for the selective emergence of a singular clone that evolves into t-MNs after acquiring secondary mutations. These results underscore the importance of further systematic research to elucidate the long-term hematological consequences of cancer chemotherapy.

Mutational signatures of colorectal cancers according to distinct computational workflows

Tumor mutational signatures have gained prominence in cancer research, yet the lack of standardized methods hinders reproducibility and robustness. Leveraging colorectal cancer (CRC) as a model, we explored the influence of computational parameters on mutational signature analyses across 230 CRC cell lines and 152 CRC patients. Results were validated in three independent datasets: 483 endometrial cancer patients stratified by mismatch repair (MMR) status, 35 lung cancer patients by smoking status and 12 patient-derived organoids (PDOs) annotated for colibactin exposure. Assessing various bioinformatic tools, reference datasets and input data sizes including whole genome sequencing, whole exome sequencing and a pan-cancer gene panel, we demonstrated significant variability in the results. We report that the use of distinct algorithms and references led to statistically different results, highlighting how arbitrary choices may induce variability in the mutational signature contributions. Furthermore, we found a differential contribution of mutational signatures between coding and intergenic regions and defined the minimum number of somatic variants required for reliable mutational signature assignment. To facilitate the identification of the most suitable workflows, we developed Comparative Mutational Signature analysis on Coding and Extragenic Regions (CoMSCER), a bioinformatic tool which allows researchers to easily perform comparative mutational signature analysis by coupling the results from several tools and public reference datasets and to assess mutational signature contributions in coding and non-coding genomic regions. In conclusion, our study provides a comparative framework to elucidate the impact of distinct computational workflows on mutational signatures.

The mutagenic forces shaping the genomic landscape of lung cancer in never smokers

Lung cancer in never smokers (LCINS) accounts for up to 25% of all lung cancers and has been associated with exposure to secondhand tobacco smoke and air pollution in observational studies. Here, we evaluate the mutagenic exposures in LCINS by examining deep whole-genome sequencing data from a large international cohort of 871 treatment-naïve LCINS recruited from 28 geographical locations within the Sherlock-Lung study. KRAS mutations were 3.8-fold more common in adenocarcinomas of never smokers from North America and Europe, while a 1.6-fold higher prevalence of EGFR and TP53 mutations was observed in adenocarcinomas from East Asia. Signature SBS40a, with unknown cause, was found in most samples and accounted for the largest proportion of single base substitutions in adenocarcinomas, being enriched in EGFR-mutated cases. Conversely, the aristolochic acid signature SBS22a was almost exclusively observed in patients from Taipei. Even though LCINS exposed to secondhand smoke had an 8.3% higher mutational burden and 5.4% shorter telomeres, passive smoking was not associated with driver mutations in cancer driver genes or the activities of individual mutational signatures. In contrast, patients from regions with high levels of air pollution were more likely to have TP53 mutations while exhibiting shorter telomeres and an increase in most types of somatic mutations, including a 3.9-fold elevation of signature SBS4 (q-value=3.1 × 10-5), previously linked mainly to tobacco smoking, and a 76% increase of clock-like signature SBS5 (q-value=5.0 × 10-5). A positive dose-response effect was observed with air pollution levels, which correlated with both a decrease in telomere length and an elevation in somatic mutations, notably attributed to signatures SBS4 and SBS5. Our results elucidate the diversity of mutational processes shaping the genomic landscape of lung cancer in never smokers.

Genome-wide somatic mutation analysis of sinonasal adenocarcinoma with and without wood dust exposure

BACKGROUND

Sinonasal adenocarcinoma is a rare cancer, encompassing two different entities, the intestinal-type sinonasal adenocarcinoma (ITAC) and the non-intestinal-type sinonasal adenocarcinoma (non-ITAC). Occurrence of ITAC is strongly associated with exposure to hardwood dusts. In countries with predominant exposure to softwood dust the occurrence of sinonasal adenocarcinomas is lower and the relative amount of non-ITACs to ITACs is higher. The molecular mechanisms behind the tumorigenic effects of wood dust remain largely unknown.

METHODS

We carried out whole-genome sequencing of formalin-fixed paraffin-embedded (FFPE) samples of sinonasal adenocarcinomas from ten wood dust-exposed and six non-exposed individuals, with partial tobacco exposure data. Sequences were analyzed for the presence of mutational signatures matching COSMIC database signatures. Driver mutations and CN variant regions were characterized.

RESULTS

Mutation burden was higher in samples of wood dust-exposed patients (p = 0.016). Reactive oxygen species (ROS) damage-related mutational signatures were almost exclusively identified in ITAC subtype samples (p = 0.00055). Tobacco smoke mutational signatures were observed in samples of patients with tobacco exposure or missing information, but not in samples from non-exposed patients. A tetraploidy copy number (CN) signature was enriched in ITAC subtype (p = 0.042). CN variation included recurrent gains in COSMIC Cancer Gene Census genes TERT, SDHA, RAC1, ETV1, PCM1, and MYC. Pathogenic variants were observed most frequently in TP53, NF1, CHD2, BRAF, APC, and LRP1B. Driver mutations and copy number gains did not segregate by subtype.

CONCLUSIONS

Our analysis identified distinct mutational characteristics in ITAC and non-ITAC. Mutational signature analysis may eventually become useful for documentation of occupation-related cancer, while the exact mechanisms behind wood dust-driven carcinogenesis remain elusive. The presence of homologous recombination deficiency signatures implies a novel opportunity for treatment, but further studies are needed.

MUSE-XAE: MUtational Signature Extraction with eXplainable AutoEncoder enhances tumour types classification

MOTIVATION

Mutational signatures are a critical component in deciphering the genetic alterations that underlie cancer development and have become a valuable resource to understand the genomic changes during tumorigenesis. Therefore, it is essential to employ precise and accurate methods for their extraction to ensure that the underlying patterns are reliably identified and can be effectively utilized in new strategies for diagnosis, prognosis, and treatment of cancer patients.

RESULTS

We present MUSE-XAE, a novel method for mutational signature extraction from cancer genomes using an explainable autoencoder. Our approach employs a hybrid architecture consisting of a nonlinear encoder that can capture nonlinear interactions among features, and a linear decoder which ensures the interpretability of the active signatures. We evaluated and compared MUSE-XAE with other available tools on both synthetic and real cancer datasets and demonstrated that it achieves superior performance in terms of precision and sensitivity in recovering mutational signature profiles. MUSE-XAE extracts highly discriminative mutational signature profiles by enhancing the classification of primary tumour types and subtypes in real world settings. This approach could facilitate further research in this area, with neural networks playing a critical role in advancing our understanding of cancer genomics.

AVAILABILITY AND IMPLEMENTATION

MUSE-XAE software is freely available at https://github.com/compbiomed-unito/MUSE-XAE.

3D genomic mapping reveals multifocality of human pancreatic precancers

Pancreatic intraepithelial neoplasias (PanINs) are the most common precursors of pancreatic cancer, but their small size and inaccessibility in humans make them challenging to study1. Critically, the number, dimensions and connectivity of human PanINs remain largely unknown, precluding important insights into early cancer development. Here, we provide a microanatomical survey of human PanINs by analysing 46 large samples of grossly normal human pancreas with a machine-learning pipeline for quantitative 3D histological reconstruction at single-cell resolution. To elucidate genetic relationships between and within PanINs, we developed a workflow in which 3D modelling guides multi-region microdissection and targeted and whole-exome sequencing. From these samples, we calculated a mean burden of 13 PanINs per cm3 and extrapolated that the normal intact adult pancreas harbours hundreds of PanINs, almost all with oncogenic KRAS hotspot mutations. We found that most PanINs originate as independent clones with distinct somatic mutation profiles. Some spatially continuous PanINs were found to contain multiple KRAS mutations; computational and in situ analyses demonstrated that different KRAS mutations localize to distinct cell subpopulations within these neoplasms, indicating their polyclonal origins. The extensive multifocality and genetic heterogeneity of PanINs raises important questions about mechanisms that drive precancer initiation and confer differential progression risk in the human pancreas. This detailed 3D genomic mapping of molecular alterations in human PanINs provides an empirical foundation for early detection and rational interception of pancreatic cancer.

Geographic variation of mutagenic exposures in kidney cancer genomes

International differences in the incidence of many cancer types indicate the existence of carcinogen exposures that have not yet been identified by conventional epidemiology make a substantial contribution to cancer burden1. In clear cell renal cell carcinoma, obesity, hypertension and tobacco smoking are risk factors, but they do not explain the geographical variation in its incidence2. Underlying causes can be inferred by sequencing the genomes of cancers from populations with different incidence rates and detecting differences in patterns of somatic mutations. Here we sequenced 962 clear cell renal cell carcinomas from 11 countries with varying incidence. The somatic mutation profiles differed between countries. In Romania, Serbia and Thailand, mutational signatures characteristic of aristolochic acid compounds were present in most cases, but these were rare elsewhere. In Japan, a mutational signature of unknown cause was found in more than 70% of cases but in less than 2% elsewhere. A further mutational signature of unknown cause was ubiquitous but exhibited higher mutation loads in countries with higher incidence rates of kidney cancer. Known signatures of tobacco smoking correlated with tobacco consumption, but no signature was associated with obesity or hypertension, suggesting that non-mutagenic mechanisms of action underlie these risk factors. The results of this study indicate the existence of multiple, geographically variable, mutagenic exposures that potentially affect tens of millions of people and illustrate the opportunities for new insights into cancer causation through large-scale global cancer genomics.

Contrasting somatic mutation patterns in aging human neurons and oligodendrocytes

Characterizing somatic mutations in the brain is important for disentangling the complex mechanisms of aging, yet little is known about mutational patterns in different brain cell types. Here, we performed whole-genome sequencing (WGS) of 86 single oligodendrocytes, 20 mixed glia, and 56 single neurons from neurotypical individuals spanning 0.4-104 years of age and identified >92,000 somatic single-nucleotide variants (sSNVs) and small insertions/deletions (indels). Although both cell types accumulate somatic mutations linearly with age, oligodendrocytes accumulated sSNVs 81% faster than neurons and indels 28% slower than neurons. Correlation of mutations with single-nucleus RNA profiles and chromatin accessibility from the same brains revealed that oligodendrocyte mutations are enriched in inactive genomic regions and are distributed across the genome similarly to mutations in brain cancers. In contrast, neuronal mutations are enriched in open, transcriptionally active chromatin. These stark differences suggest an assortment of active mutagenic processes in oligodendrocytes and neurons.

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

APOBEC enzymes are part of the innate immunity and are responsible for restricting viruses and retroelements by deaminating cytosine residues1,2. Most solid tumors harbor different levels of somatic mutations attributed to the off-target activities of APOBEC3A (A3A) and/or APOBEC3B (A3B)3-6. However, how APOBEC3A/B enzymes shape the tumor evolution in the presence of exogenous mutagenic processes is largely unknown. Here, by combining deep whole-genome sequencing with multi-omics profiling of 309 lung cancers from smokers with detailed tobacco smoking information, we identify two subtypes defined by low (LAS) and high (HAS) APOBEC mutagenesis. LAS are enriched for A3B-like mutagenesis and KRAS mutations, whereas HAS for A3A-like mutagenesis and TP53 mutations. Unlike APOBEC3A, APOBEC3B expression is strongly associated with an upregulation of the base excision repair pathway. Hypermutation by unrepaired A3A and tobacco smoking mutagenesis combined with TP53-induced genomic instability can trigger senescence7, apoptosis8, and cell regeneration9, as indicated by high expression of pulmonary healing signaling pathway, stemness markers and distal cell-of-origin in HAS. The expected association of tobacco smoking variables (e.g., time to first cigarette) with genomic/epigenomic changes are not observed in HAS, a plausible consequence of frequent cell senescence or apoptosis. HAS have more neoantigens, slower clonal expansion, and older age at onset compared to LAS, particularly in heavy smokers, consistent with high proportions of newly generated, unmutated cells and frequent immuno-editing. These findings show how heterogeneity in mutational burden across co-occurring mutational processes and cell types contributes to tumor development, with important clinical implications.

Diagnostic yield and clinical impact of germline sequencing in children with CNS and extracranial solid tumors-a nationwide, prospective Swedish study

Background

Childhood cancer predisposition (ChiCaP) syndromes are increasingly recognized as contributing factors to childhood cancer development. Yet, due to variable availability of germline testing, many children with ChiCaP might go undetected today. We report results from the nationwide and prospective ChiCaP study that investigated diagnostic yield and clinical impact of integrating germline whole-genome sequencing (gWGS) with tumor sequencing and systematic phenotyping in children with solid tumors.

Methods

gWGS was performed in 309 children at diagnosis of CNS (n = 123, 40%) or extracranial (n = 186, 60%) solid tumors and analyzed for disease-causing variants in 189 known cancer predisposing genes. Tumor sequencing data were available for 74% (227/309) of patients. In addition, a standardized clinical assessment for underlying predisposition was performed in 95% (293/309) of patients.

Findings

The prevalence of ChiCaP diagnoses was 11% (35/309), of which 69% (24/35) were unknown at inclusion (diagnostic yield 8%, 24/298). A second-hit and/or relevant mutational signature was observed in 19/21 (90%) tumors with informative data. ChiCaP diagnoses were more prevalent among patients with retinoblastomas (50%, 6/12) and high-grade astrocytomas (37%, 6/16), and in those with non-cancer related features (23%, 20/88), and ≥2 positive ChiCaP criteria (28%, 22/79). ChiCaP diagnoses were autosomal dominant in 80% (28/35) of patients, yet confirmed de novo in 64% (18/28). The 35 ChiCaP findings resulted in tailored surveillance (86%, 30/35) and treatment recommendations (31%, 11/35).

Interpretation

Overall, our results demonstrate that systematic phenotyping, combined with genomics-based diagnostics of ChiCaP in children with solid tumors is feasible in large-scale clinical practice and critically guides personalized care in a sizable proportion of patients.

Funding

The study was supported by the Swedish Childhood Cancer Fund and the Ministry of Health and Social Affairs.

Multiomics-Based Feature Extraction and Selection for the Prediction of Lung Cancer Survival

Lung cancer is a global health challenge, hindered by delayed diagnosis and the disease's complex molecular landscape. Accurate patient survival prediction is critical, motivating the exploration of various -omics datasets using machine learning methods. Leveraging multi-omics data, this study seeks to enhance the accuracy of survival prediction by proposing new feature extraction techniques combined with unbiased feature selection. Two lung adenocarcinoma multi-omics datasets, originating from the TCGA and CPTAC-3 projects, were employed for this purpose, emphasizing gene expression, methylation, and mutations as the most relevant data sources that provide features for the survival prediction models. Additionally, gene set aggregation was shown to be the most effective feature extraction method for mutation and copy number variation data. Using the TCGA dataset, we identified 32 molecular features that allowed the construction of a 2-year survival prediction model with an AUC of 0.839. The selected features were additionally tested on an independent CPTAC-3 dataset, achieving an AUC of 0.815 in nested cross-validation, which confirmed the robustness of the identified features.

Evolutionary trajectories of small cell lung cancer under therapy

The evolutionary processes that underlie the marked sensitivity of small cell lung cancer (SCLC) to chemotherapy and rapid relapse are unknown1-3. Here we determined tumour phylogenies at diagnosis and throughout chemotherapy and immunotherapy by multiregion sequencing of 160 tumours from 65 patients. Treatment-naive SCLC exhibited clonal homogeneity at distinct tumour sites, whereas first-line platinum-based chemotherapy led to a burst in genomic intratumour heterogeneity and spatial clonal diversity. We observed branched evolution and a shift to ancestral clones underlying tumour relapse. Effective radio- or immunotherapy induced a re-expansion of founder clones with acquired genomic damage from first-line chemotherapy. Whereas TP53 and RB1 alterations were exclusively part of the common ancestor, MYC family amplifications were frequently not constituents of the founder clone. At relapse, emerging subclonal mutations affected key genes associated with SCLC biology, and tumours harbouring clonal CREBBP/EP300 alterations underwent genome duplications. Gene-damaging TP53 alterations and co-alterations of TP53 missense mutations with TP73, CREBBP/EP300 or FMN2 were significantly associated with shorter disease relapse following chemotherapy. In summary, we uncover key processes of the genomic evolution of SCLC under therapy, identify the common ancestor as the source of clonal diversity at relapse and show central genomic patterns associated with sensitivity and resistance to chemotherapy.

Accurate and sensitive mutational signature analysis with MuSiCal

Mutational signature analysis is a recent computational approach for interpreting somatic mutations in the genome. Its application to cancer data has enhanced our understanding of mutational forces driving tumorigenesis and demonstrated its potential to inform prognosis and treatment decisions. However, methodological challenges remain for discovering new signatures and assigning proper weights to existing signatures, thereby hindering broader clinical applications. Here we present Mutational Signature Calculator (MuSiCal), a rigorous analytical framework with algorithms that solve major problems in the standard workflow. Our simulation studies demonstrate that MuSiCal outperforms state-of-the-art algorithms for both signature discovery and assignment. By reanalyzing more than 2,700 cancer genomes, we provide an improved catalog of signatures and their assignments, discover nine indel signatures absent in the current catalog, resolve long-standing issues with the ambiguous 'flat' signatures and give insights into signatures with unknown etiologies. We expect MuSiCal and the improved catalog to be a step towards establishing best practices for mutational signature analysis.

Epistasis between mutator alleles contributes to germline mutation spectrum variability in laboratory mice

Maintaining germline genome integrity is essential and enormously complex. Although many proteins are involved in DNA replication, proofreading, and repair, mutator alleles have largely eluded detection in mammals. DNA replication and repair proteins often recognize sequence motifs or excise lesions at specific nucleotides. Thus, we might expect that the spectrum of de novo mutations - the frequencies of C>T, A>G, etc. - will differ between genomes that harbor either a mutator or wild-type allele. Previously, we used quantitative trait locus mapping to discover candidate mutator alleles in the DNA repair gene Mutyh that increased the C>A germline mutation rate in a family of inbred mice known as the BXDs (Sasani et al., 2022, Ashbrook et al., 2021). In this study we developed a new method to detect alleles associated with mutation spectrum variation and applied it to mutation data from the BXDs. We discovered an additional C>A mutator locus on chromosome 6 that overlaps Ogg1, a DNA glycosylase involved in the same base-excision repair network as Mutyh (David et al., 2007). Its effect depends on the presence of a mutator allele near Mutyh, and BXDs with mutator alleles at both loci have greater numbers of C>A mutations than those with mutator alleles at either locus alone. Our new methods for analyzing mutation spectra reveal evidence of epistasis between germline mutator alleles and may be applicable to mutation data from humans and other model organisms.