The repertoire of mutational signatures in human cancer

Evolution of genome and immunogenome in esophageal squamous cell carcinomas driven by neoadjuvant chemoradiotherapy

Neoadjuvant chemoradiotherapy (NCRT) followed by surgery is a standard treatment for locally advanced esophageal squamous cell carcinomas (ESCCs). However, the evolution of genome and immunogenome in ESCCs driven by NCRT remains incompletely elucidated. We performed whole-exome sequencing of 51 ESCC tumors collected before and after NCRT, 36 of which were subjected to transcriptome sequencing. Clonal analysis identified clonal extinction in 13 ESCC patients wherein all pre-NCRT clones disappeared after NCRT, and clonal persistence in 9 patients wherein clones endured following NCRT. The clone-persistent patients showed higher pre-NCRT genomic intratumoral heterogeneity and worse prognosis than the clone-extinct ones. In contrast to the clone-extinct patients, the clone-persistent patients demonstrated a high proportion of subclonal neoantigens within pre-treatment specimens. Transcriptome analysis revealed increased immune infiltrations and up-regulated immune-related pathways after NCRT, especially in the clone-extinct patients. The number of T cell receptor-neoantigen interactions was higher in the clone-extinct patients than in the clone-persistent ones. The decrease in T cell repertoire evenness positively correlated to the decreased number of clonal neoantigens after NCRT, especially in the clone-extinct patients. In conclusion, we identified two prognosis-related clonal dynamic modes driven by NCRT in ESCCs. This study extended our knowledge of the ESCC genome and immunogenome evolutions driven by NCRT.

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.

Detection of circulating tumor DNA by tumor-informed whole-genome sequencing enables prediction of recurrence in stage III colorectal cancer patients

INTRODUCTION

Circulating tumor (ctDNA) can be used to detect residual disease after cancer treatment. Detecting low-level ctDNA is challenging, due to the limited number of recoverable ctDNA fragments at any target loci. In response, we applied tumor-informed whole-genome sequencing (WGS), leveraging thousands of mutations for ctDNA detection.

METHODS

Performance was evaluated in serial plasma samples (n = 1283) from 144 stage III colorectal cancer patients. Tumor/normal WGS was used to establish a patient-specific mutational fingerprint, which was searched for in 20x WGS plasma profiles. For reproducibility, paired aliquots of 172 plasma samples were analyzed in two independent laboratories. De novo variant calling was performed for serial plasma samples with a ctDNA level > 10 % (n = 17) to explore genomic evolution.

RESULTS

WGS-based ctDNA detection was prognostic of recurrence: post-operation (Hazard ratio [HR] 6.75, 95 %CI 3.18-14.3, p < 0.001), post-adjuvant chemotherapy (HR 28.9, 95 %CI 10.1-82.8; p < 0.001), and during surveillance (HR 22.8, 95 %CI 13.7-37.9, p < 0.0001). The 3-year cumulative incidence of ctDNA detection in recurrence patients was 95 %. ctDNA was detected a median of 8.7 months before radiological recurrence. The independently analyzed plasma aliquots showed excellent agreement (Cohens Kappa=0.9, r = 0.99). Genomic characterization of serial plasma revealed significant evolution in mutations and copy number alterations, and the timing of mutational processes, such as 5-fluorouracil-induced mutations.

CONCLUSION

Our study supports the use of WGS for sensitive ctDNA detection and demonstrates that post-treatment ctDNA detection is highly prognostic of recurrence. Furthermore, plasma WGS can identify genomic differences distinguishing the primary tumor and relapsing metastasis, and monitor treatment-induced genomic changes.

Somatic mutations in autoinflammatory and autoimmune disease

Somatic mutations (also known as acquired mutations) are emerging as common, age-related processes that occur in all cells throughout the body. Somatic mutations are canonically linked to malignant processes but over the past decade have been increasingly causally connected to benign diseases including rheumatic conditions. Here we outline the contribution of somatic mutations to complex and monogenic immunological diseases with a detailed review of unique aspects associated with such causes. Somatic mutations can cause early- or late-onset rheumatic monogenic diseases but also contribute to the pathogenesis of complex inflammatory and immune-mediated diseases, affect disease progression and define new clinical subtypes. Although even variants with a low variant allele fraction can be pathogenic, clonal dynamics could lead to changes over time in the proportion of mutant cells, with possible phenotypic consequences for the individual. Thus, somatic mutagenesis and clonal expansion have relevant implications in genetic testing and counselling. On the basis of both increased recognition of somatic diseases in clinical practice and improved technical and bioinformatic processes, we hypothesize that there will be an ever-expanding list of somatic mutations in various genes leading to inflammatory conditions, particularly in late-onset disease.

Somatic mosaicism in schizophrenia brains reveals prenatal mutational processes

Germline mutations modulate the risk of developing schizophrenia (SCZ). Much less is known about the role of mosaic somatic mutations in the context of SCZ. Deep (239×) whole-genome sequencing (WGS) of brain neurons from 61 SCZ cases and 25 controls postmortem identified mutations occurring during prenatal neurogenesis. SCZ cases showed increased somatic variants in open chromatin, with increased mosaic CpG transversions (CpG>GpG) and T>G mutations at transcription factor binding sites (TFBSs) overlapping open chromatin, a result not seen in controls. Some of these variants alter gene expression, including SCZ risk genes and genes involved in neurodevelopment. Although these mutational processes can reflect a difference in factors indirectly involved in disease, increased somatic mutations at developmental TFBSs could also potentially contribute to SCZ.

High-resolution and programmable RNA-IN and RNA-OUT genetic circuit in living mammalian cells

RNAs and their encoded proteins intricately regulate diverse cell types and states within the human body. Dysregulated RNA expressions or mutations can lead to various diseased cell states, including tumorigenesis. Detecting and manipulating these endogenous RNAs offers significant promise for restoring healthy cell states and targeting tumors both in research and clinical contexts. This study presents an RNA-IN and RNA-OUT genetic circuit capable dynamically sensing and manipulating any RNA target in a programmable manner. The RNA-IN module employes a programmable CRISPR-associated protease (CASP) complex for RNA detection, while the RNA-OUT module utilizes an engineered protease-responsive dCas9-VPR activator. Additionally, the CASP module can detect point mutations by harnessing an uncovered dual-nucleotide synergistic switching effect within the CASP complex, resulting in the amplification of point-mutation signals from initially undetectable levels (1.5-fold) to a remarkable 94-fold. We successfully showcase the circuit's ability to rewire endogenous RNA-IN signals to activate endogenous progesterone biosynthesis pathway, dynamically monitor adipogenic differentiation of mesenchymal stem cells (MSCs) and the epithelial-to-mesenchmal trans-differentiation, as well as selective killing of tumor cells. The programmable RNA-IN and RNA-OUT circuit exhibits tremendous potential for applications in gene therapy, biosensing and design of synthetic regulatory networks.

Human DNA polymerase ε is a source of C>T mutations at CpG dinucleotides

C-to-T transitions in CpG dinucleotides are the most prevalent mutations in human cancers and genetic diseases. These mutations have been attributed to deamination of 5-methylcytosine (5mC), an epigenetic modification found on CpGs. We recently linked CpG>TpG mutations to replication and hypothesized that errors introduced by polymerase ε (Pol ε) may represent an alternative source of mutations. Here we present a new method called polymerase error rate sequencing (PER-seq) to measure the error spectrum of DNA polymerases in isolation. We find that the most common human cancer-associated Pol ε mutant (P286R) produces an excess of CpG>TpG errors, phenocopying the mutation spectrum of tumors carrying this mutation and deficiencies in mismatch repair. Notably, we also discover that wild-type Pol ε has a sevenfold higher error rate when replicating 5mCpG compared to C in other contexts. Together, our results from PER-seq and human cancers demonstrate that replication errors are a major contributor to CpG>TpG mutagenesis in replicating cells, fundamentally changing our understanding of this important disease-causing mutational mechanism.

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.

High fidelity DNA ligation prevents single base insertions in the yeast genome

Finalization of eukaryotic nuclear DNA replication relies on DNA ligase 1 (LIG1) to seal DNA nicks generated during Okazaki Fragment Maturation (OFM). Using a mutational reporter in Saccharomyces cerevisiae, we previously showed that mutation of the high-fidelity magnesium binding site of LIG1Cdc9 strongly increases the rate of single-base insertions. Here we show that this rate is increased across the nuclear genome, that it is synergistically increased by concomitant loss of DNA mismatch repair (MMR), and that the additions occur in highly specific sequence contexts. These discoveries are all consistent with incorporation of an extra base into the nascent lagging DNA strand that can be corrected by MMR following mutagenic ligation by the Cdc9-EEAA variant. There is a strong preference for insertion of either dGTP or dTTP into 3-5 base pair mononucleotide sequences with stringent flanking nucleotide requirements. The results reveal unique LIG1Cdc9-dependent mutational motifs where high fidelity DNA ligation of a subset of OFs is critical for preventing mutagenesis across the genome.

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.

Topography of mutational signatures in non-small cell lung cancer: emerging concepts, clinical applications, and limitations

The genome of a cell is continuously battered by a plethora of exogenous and endogenous processes that can lead to damaged DNA. Repair mechanisms correct this damage most of the time, but failure to do so leaves mutations. Mutations do not occur in random manner, but rather typically follow a more or less specific pattern due to known or imputed mutational processes. Mutational signature analysis is the process by which the predominant mutational process can be inferred for a cancer and can be used in several contexts to study both the genesis of cancer and its response to therapy. Recent pan-cancer genomic efforts such as "The Cancer Genome Atlas" have identified numerous mutational signatures that can be categorized into single base substitutions, doublet base substitutions, or small insertions/deletions. Understanding these mutational signatures as they occur in non-small lung cancer could improve efforts at prevention, predict treatment response to personalized treatments, and guide the development of therapies targeting tumor evolution. For non-small cell lung cancer, several mutational signatures have been identified that correlate with exposures such as tobacco smoking and radon and can also reflect endogenous processes such as aging, APOBEC activity, and loss of mismatch repair. Herein, we provide an overview of the current knowledge of mutational signatures in non-small lung cancer.

TEX264 drives selective autophagy of DNA lesions to promote DNA repair and cell survival

DNA repair and autophagy are distinct biological processes vital for cell survival. Although autophagy helps maintain genome stability, there is no evidence of its direct role in the repair of DNA lesions. We discovered that lysosomes process topoisomerase 1 cleavage complexes (TOP1cc) DNA lesions in vertebrates. Selective degradation of TOP1cc by autophagy directs DNA damage repair and cell survival at clinically relevant doses of topoisomerase 1 inhibitors. TOP1cc are exported from the nucleus to lysosomes through a transient alteration of the nuclear envelope and independent of the proteasome. Mechanistically, the autophagy receptor TEX264 acts as a TOP1cc sensor at DNA replication forks, triggering TOP1cc processing by the p97 ATPase and mediating the delivery of TOP1cc to lysosomes in an MRE11-nuclease- and ATR-kinase-dependent manner. We found an evolutionarily conserved role for selective autophagy in DNA repair that enables cell survival, protects genome stability, and is clinically relevant for colorectal cancer patients.

Altered dNTP pools accelerate tumor formation in mice

Alterations in deoxyribonucleoside triphosphate (dNTP) pools have been linked to increased mutation rates and genome instability in unicellular organisms and cell cultures. However, the role of dNTP pool changes in tumor development in mammals remains unclear. In this study, we present a mouse model with a point mutation at the allosteric specificity site of ribonucleotide reductase, RRM1-Y285A. This mutation reduced ribonucleotide reductase activity, impairing the synthesis of deoxyadenosine triphosphate (dATP) and deoxyguanosine triphosphate (dGTP). Heterozygous Rrm1+/Y285A mice exhibited distinct alterations in dNTP pools across various organs, shorter lifespans and earlier tumor onset compared with wild-type controls. Mutational spectrum analysis of tumors revealed two distinct signatures, one resembling a signature extracted from a human cancer harboring a mutation of the same amino acid residue in ribonucleotide reductase, RRM1Y285C. Our findings suggest that mutations in enzymes involved in dNTP metabolism can serve as drivers of cancer development.

The evolutionary history of metastatic pancreatic neuroendocrine tumours reveals a therapy driven route to high-grade transformation

Tumour evolution with acquisition of more aggressive disease characteristics is a hallmark of disseminated cancer. Metastatic pancreatic neuroendocrine tumours (PanNETs) in particular may progress from a low/intermediate to a high-grade disease. The aim of this work was to understand the molecular mechanisms underlying metastatic progression as well as PanNET transformation from a low/intermediate to a high-grade disease. We performed multi-omics analysis (genome/exome sequencing, total RNA-sequencing and methylation array) of 32 longitudinal samples from six patients with metastatic low/intermediate grade PanNET. The clonal composition of tumour lesions and underlying phylogeny of each patient were determined with bioinformatics analyses. Findings were validated in post-alkylating chemotherapy samples from 24 patients with PanNET using targeted next generation sequencing. We validate the current PanNET evolutionary model with MEN1 inactivation that occurs very early in tumourigenesis. This was followed by pronounced genetic diversity on both spatial and temporal levels, with parallel and convergent tumour evolution involving the ATRX/DAXX and mechanistic target of the rapamycin (mTOR) pathways. Following alkylating chemotherapy treatment, some PanNETs developed mismatch repair deficiency and acquired a hypermutational phenotype. This was validated among 16 patients with PanNET who had high-grade progression after alkylating chemotherapy, of whom eight had a tumour mutational burden >50 (50%). In comparison, among the eight patients who did not show high-grade progression, 0 had a tumour mutational burden >50 (0%; odds ratio 'infinite', 95% confidence interval 1.8 to 'infinite', p = 0.02). Our findings contribute to broaden the understanding of metastatic/high-grade PanNETs and suggests that therapy driven disease evolution is an important hallmark of this disease. © 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. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

MHC Hammer reveals genetic and non-genetic HLA disruption in cancer evolution

Disruption of the class I human leukocyte antigen (HLA) molecules has important implications for immune evasion and tumor evolution. We developed major histocompatibility complex loss of heterozygosity (LOH), allele-specific mutation and measurement of expression and repression (MHC Hammer). We identified extensive variability in HLA allelic expression and pervasive HLA alternative splicing in normal lung and breast tissue. In lung TRACERx and lung and breast TCGA cohorts, 61% of lung adenocarcinoma (LUAD), 76% of lung squamous cell carcinoma (LUSC) and 35% of estrogen receptor-positive (ER+) cancers harbored class I HLA transcriptional repression, while HLA tumor-enriched alternative splicing occurred in 31%, 11% and 15% of LUAD, LUSC and ER+ cancers. Consistent with the importance of HLA dysfunction in tumor evolution, in LUADs, HLA LOH was associated with metastasis and LUAD primary tumor regions seeding a metastasis had a lower effective neoantigen burden than non-seeding regions. These data highlight the extent and importance of HLA transcriptomic disruption, including repression and alternative splicing in cancer evolution.

Establishment, characterization, and biobanking of 36 pancreatic cancer organoids: prediction of metastasis in resectable pancreatic cancer

PURPOSE

Early dissemination of primary pancreatic ductal adenocarcinoma (PDAC) is the main cause of dismal prognosis as it highly limits possible treatment options. A number of PDAC patients experience distant metastasis even after treatment due to the metastatic clones. We aimed to demonstrate the molecular architecture of borderline resectable PDAC manifests cancer dissemination of PDAC.

METHODS

Here, 36 organoids isolated from primary tumor masses of PDAC patients with diverse metastatic statues are presented. Whole-exome sequencing and RNA sequencing were performed and drug responses to clinically relevant 18 compounds were assessed.

RESULTS

Our results revealed that borderline resectable PDAC organoids exhibited distinct patterns according to their metastatic potency highlighted by multiple genetic and transcriptional factors and strong variances in drug responses.

CONCLUSIONS

These data suggest that the presence of metastatic PDAC can be identified by integrating molecular compositions and drug responses of borderline resectable PDAC.

Genomic Profiling to Contextualize the Results of Intervention for Smoldering Multiple Myeloma

PURPOSE

Early intervention for high-risk smoldering multiple myeloma (HR-SMM) achieves deep and prolonged responses. It is unclear if beneficial outcomes are due to the treatment of less complex, susceptible disease or inaccuracy in clinical definition of cases entered.

EXPERIMENTAL DESIGN

In this study, we interrogated whole-genome and whole-exome sequencing for 54 patients across two HR-SMM interventional studies (NCT01572480 and NCT02279394).

RESULTS

We reveal that the genomic landscape of treated HR-SMM is generally simple as compared with newly diagnosed multiple myeloma counterparts with less inactivation of tumor suppressor genes, RAS pathway mutations, MYC disruption, and APOBEC contribution. The absence of these events parallels that of indolent precursor conditions, possibly explaining overall excellent outcomes. However, some patients harboring genomic complexity fail to sustain response and experience resistant, progressive disease. Overall, clinical risk scores do not effectively discriminate between genomically indolent and aggressive disease.

CONCLUSIONS

Genomic profiling can contextualize the advantage of early intervention in SMM and guide personalization of therapy. See related commentary by Weinhold and Rasche, p. 4263.

Tumour mutational burden: clinical utility, challenges and emerging improvements

Tumour mutational burden (TMB), defined as the total number of somatic non-synonymous mutations present within the cancer genome, varies across and within cancer types. A first wave of retrospective and prospective research identified TMB as a predictive biomarker of response to immune-checkpoint inhibitors and culminated in the disease-agnostic approval of pembrolizumab for patients with TMB-high tumours based on data from the Keynote-158 trial. Although the applicability of outcomes from this trial to all cancer types and the optimal thresholds for TMB are yet to be ascertained, research into TMB is advancing along three principal avenues: enhancement of TMB assessments through rigorous quality control measures within the laboratory process, including the mitigation of confounding factors such as limited panel scope and low tumour purity; refinement of the traditional TMB framework through the incorporation of innovative concepts such as clonal, persistent or HLA-corrected TMB, tumour neoantigen load and mutational signatures; and integration of TMB with established and emerging biomarkers such as PD-L1 expression, microsatellite instability, immune gene expression profiles and the tumour immune contexture. Given its pivotal functions in both the pathogenesis of cancer and the ability of the immune system to recognize tumours, a profound comprehension of the foundational principles and the continued evolution of TMB are of paramount relevance for the field of oncology.

Mutations, substitutions, and selection: Linking mutagenic processes to cancer using evolutionary theory

Cancers are the product of evolutionary events, where molecular variation occurs and accumulates in tissues and tumors. Sequencing of this molecular variation informs not only which variants are driving tumorigenesis, but also the mechanisms behind what is fueling mutagenesis. Both of these details are crucial for preventing premature deaths due to cancer, whether it is by targeting the variants driving the cancer phenotype or by measures to prevent exogenous mutations from contributing to somatic evolution. Here, we review tools to determine both molecular signatures and cancer drivers, and avenues by which these metrics may be linked.

Comprehensive molecular characterization of early stage grade 3 endometrioid endometrial adenocarcinoma

OBJECTIVE

The treatment for stage IB grade 3 endometrioid endometrial adenocarcinoma is challenging with variable practice. Molecular characterization may help identify adjuvant therapy strategies beyond stage. We aimed to better understand the molecular features of these tumors by characterizing them by ProMisE classification, mutational signature, and commonly mutated genes.

METHODS

Patients with stage IB grade 3 EEC at two institutions were included. Immunohistochemistry and whole exome sequencing were performed on archival FFPE tissue sections to determine ProMisE classification. Personal Cancer Genome Reporter was used for somatic variant annotation, and mutational signatures were generated based on COSMIC single base substitution mutational signatures.

RESULTS

46 patients were included with variable adjuvant treatment. Nine patients recurred (19.6%), most with extra-abdominal disease (n = 5, or 55.6%). 10 had POLE mutations (21.7%), 18 were MMR deficient (39.1%), 6 had abnormal p53 (13.0%), and 12 were p53 wildtype (26.1%). There were no recurrences in the POLE subgroup. A dominant mutational signature was identified in 38 patients: 17 SBS5 signature (44.7%), 10 SBS15 or SBS44 signature (26.3%), 7 SBS10a or SBS10b signature (18.4%), 3 SBS14 signature (7.9%), and 1 SBS40 signature (2.6%). The six patients that recurred had a SBS5 signature. Frequently mutated genes included ARID1A (n = 30, 65%), PTEN (n = 28, 61%), MUC16 (n = 27, 59%), and PIK3CA (n = 25, 54%).

CONCLUSIONS

This comprehensive evaluation found a molecularly diverse cohort of tumors, despite the same histology, stage and grade. Mutational signature SBS5 correlated with a high risk of recurrence. Further refining of endometrial cancer classification may enable more precise patient stratification and personalized treatment approaches.

DNA mismatch repair controls the mutagenicity of Polymerase ζ-dependent translesion synthesis at methylated guanines

By replicating damaged nucleotides, error-prone DNA translesion synthesis (TLS) enables the completion of replication, albeit at the expense of fidelity. TLS of helix-distorting DNA lesions, that usually have reduced capacity of basepairing, comprises insertion opposite the lesion followed by extension, the latter in particular by polymerase ζ (Pol ζ). However, little is known about involvement of Pol ζ in TLS of non- or poorly-distorting, but miscoding, lesions such as O6-methyldeoxyguanosine (O6-medG). Using purified Pol ζ we describe that the enzyme can misincorporate thymidine opposite O6-medG and efficiently extend from terminal mismatches, suggesting its involvement in the mutagenicity of O6-medG. Surprisingly, O6-medG lesions induced by the methylating agent N-methyl-N'-nitro-N-nitrosoguanidine (MNNG) appeared more, rather than less, mutagenic in Pol ζ-deficient mouse embryonic fibroblasts (MEFs) than in wild type MEFs. This suggested that in vivo Pol ζ participates in non-mutagenic TLS of O6-medG. However, we found that the Pol ζ-dependent misinsertions at O6-medG lesions are efficiently corrected by DNA mismatch repair (MMR), which masks the error-proneness of Pol ζ. We also found that the MNNG-induced mutational signature is determined by the adduct spectrum, and modulated by MMR. The signature mimicked single base substitution signature 11 in the catalogue of somatic mutations in cancer, associated with treatment with the methylating drug temozolomide. Our results unravel the individual roles of the major contributors to methylating drug-induced mutagenesis. Moreover, these results warrant caution as to the classification of TLS as mutagenic or error-free based on in vitro data or on the analysis of mutations induced in MMR-proficient cells.

Global and local genomic features together modulate the spontaneous single nucleotide mutation rate

Spontaneous mutations are evolutionary engines as they generate variants for the evolutionary downstream processes that give rise to speciation and adaptation. Single nucleotide mutations (SNM) are the most abundant type of mutations among them. Here, we perform a meta-analysis to quantify the influence of selected global genomic parameters (genome size, genomic GC content, genomic repeat fraction, number of coding genes, gene count, and strand bias in prokaryotes) and local genomic features (local GC content, repeat content, CpG content and the number of SNM at CpG islands) on spontaneous SNM rates across the tree of life (prokaryotes, unicellular eukaryotes, multicellular eukaryotes) using wild-type sequence data in two different taxon classification systems. We find that the spontaneous SNM rates in our data are correlated with many genomic features in prokaryotes and unicellular eukaryotes irrespective of their sample sizes. On the other hand, only the number of coding genes was correlated with the spontaneous SNM rates in multicellular eukaryotes primarily contributed by vertebrates data. Considering local features, we notice that local GC content and CpG content significantly were correlated with the spontaneous SNM rates in the unicellular eukaryotes, while local repeat fraction is an important feature in prokaryotes and certain specific uni- and multi-cellular eukaryotes. Such predictive features of the spontaneous SNM rates often support non-linear models as the best fit compared to the linear model. We also observe that the strand asymmetry in prokaryotes plays an important role in determining the spontaneous SNM rates but the SNM spectrum does not.

Biphenotypic sinonasal sarcoma diagnosed by detection of PAX3-MAML3 fusion gene using integrated whole-genome and transcriptome sequencing

Biphenotypic sinonasal sarcoma (BSNS) is a double-phenotype sarcoma that shows differentiation in both the nervous and muscular systems. To date, whole-genome and transcriptome sequencing (WGTS) has not been used to analyze BSNS. We report a patient with BSNS who was diagnosed based on PAX3 rearrangement using WGTS. A 71-year-old Japanese male without remarkable symptoms showed a nasal tumor when undergoing computed tomography. Although pathological examination revealed a non-characteristic spindle cell tumor, a definitive diagnosis could not be made based on this examination. Endoscopic sinus surgery was performed for subsequent diagnosis, treatment, and WGTS. WGTS revealed a t(2; 4)(q35; q31.1) reciprocal translocation, resulting in a PAX3-MAML3 fusion gene, leading to a definitive diagnosis of BSNS. We also detected upregulation of the expression of PAX3, MAML3, and 11 known genes involved in neural and myogenic differentiation relevant to the BSNS phenotype. Hence, using WGTS in combination with conventional pathological diagnosis can contribute to a definitive diagnosis of rare cancers, including BSNS, by detecting chromosomal rearrangements or diagnostic markers.

Somatic mutations in aging and disease

Time always leaves its mark, and our genome is no exception. Mutations in the genome of somatic cells were first hypothesized to be the cause of aging in the 1950s, shortly after the molecular structure of DNA had been described. Somatic mutation theories of aging are based on the fact that mutations in DNA as the ultimate template for all cellular functions are irreversible. However, it took until the 1990s to develop the methods to test if DNA mutations accumulate with age in different organs and tissues and estimate the severity of the problem. By now, numerous studies have documented the accumulation of somatic mutations with age in normal cells and tissues of mice, humans, and other animals, showing clock-like mutational signatures that provide information on the underlying causes of the mutations. In this review, we will first briefly discuss the recent advances in next-generation sequencing that now allow quantitative analysis of somatic mutations. Second, we will provide evidence that the mutation rate differs between cell types, with a focus on differences between germline and somatic mutation rate. Third, we will discuss somatic mutational signatures as measures of aging, environmental exposure, and activities of DNA repair processes. Fourth, we will explain the concept of clonally amplified somatic mutations, with a focus on clonal hematopoiesis. Fifth, we will briefly discuss somatic mutations in the transcriptome and in our other genome, i.e., the genome of mitochondria. We will end with a brief discussion of a possible causal contribution of somatic mutations to the aging process.

Performance of somatic structural variant calling in lung cancer using Oxford Nanopore sequencing technology

BACKGROUND

Lung cancer is a heterogeneous disease and the primary cause of cancer-related mortality worldwide. Somatic mutations, including large structural variants, are important biomarkers in lung cancer for selecting targeted therapy. Genomic studies in lung cancer have been conducted using short-read sequencing. Emerging long-read sequencing technologies are a promising alternative to study somatic structural variants, however there is no current consensus on how to process data and call somatic events. In this study, we preformed whole genome sequencing of lung cancer and matched non-tumour samples using long and short read sequencing to comprehensively benchmark three sequence aligners and seven structural variant callers comprised of generic callers (SVIM, Sniffles2, DELLY in generic mode and cuteSV) and somatic callers (Severus, SAVANA, nanomonsv and DELLY in somatic modes).

RESULTS

Different combinations of aligners and variant callers influenced somatic structural variant detection. The choice of caller had a significant influence on somatic structural variant detection in terms of variant type, size, sensitivity, and accuracy. The performance of each variant caller was assessed by comparing to somatic structural variants identified by short-read sequencing. When compared to somatic structural variants detected with short-read sequencing, more events were detected with long-read sequencing. The mean recall of somatic variant events identified by long-read sequencing was higher for the somatic callers (72%) than generic callers (53%). Among the somatic callers when using the minimap2 aligner, SAVANA and Severus achieved the highest recall at 79.5% and 79.25% respectively, followed by nanomonsv with a recall of 72.5%.

CONCLUSION

Long-read sequencing can identify somatic structural variants in clincal samples. The longer reads have the potential to improve our understanding of cancer development and inform personalized cancer treatment.

Reconstructing oral cavity tumor evolution through brush biopsy

Oral potentially malignant disorders (OPMDs) with genomic alterations have a heightened risk of evolving into oral squamous cell carcinoma (OSCC). Currently, genomic data are typically obtained through invasive tissue biopsy. However, brush biopsy is a non-invasive method that has been utilized for identifying dysplastic cells in OPMD but its effectiveness in reflecting the genomic landscape of OPMDs remains uncertain. This pilot study investigates the potential of brush biopsy samples in accurately reconstructing the genomic profile and tumor evolution in a patient with both OPMD and OSCC. We analyzed single nucleotide variants (SNVs), copy number aberrations (CNAs), and subclonal architectures in paired tissue and brush biopsy samples. The results showed that brush biopsy effectively captured 90% of SNVs and had similar CNA profiles as those seen in its paired tissue biopsies in all lesions. It was specific, as normal buccal mucosa did not share these genomic alterations. Interestingly, brush biopsy revealed shared SNVs and CNAs between the distinct OPMD and OSCC lesions from the same patient, indicating a common ancestral origin. Subclonal reconstruction confirmed this shared ancestry, followed by divergent evolution of the lesions. These findings highlight the potential of brush biopsies in accurately representing the genomic profile of OPL and OSCC, proving insight into reconstructing tumor evolution.

Mutational signatures in 175 Chinese gastric cancer patients

BACKGROUND

Gastric cancer (GC), a molecularly heterogeneous disease, is the third leading cause of cancer death worldwide. The majority of GC cases worldwide occur in East Asia, predominantly China. Mutational Signature Framework offers an elegant approach to identify mutational processes present in tumors.

METHODS

To identify mutational signature patterns, we conducted whole exome sequencing (WES) analysis in Chinese patients with GC. Mutect2 and MutsigCV were used to identify significantly mutated genes in 175 Chinese GC cases using paired tumor-normal tissues. We investigated mutational signatures using Catalogue of Somatic Mutations in Cancer (COSMIC) Version 2 (V2) and Version 3 (V3).

RESULTS

We identified 104 mutated genes with P < 0.01. Seven genes (OR6B1, B2M, ELF3, RHOA, RPL22, TP53, ARIDIA) had q < 0.0001, including six previously associated with GC. Mutational signatures (COSMIC-V3) observed include 14 single base substitutions (SBS), one doublet base substitution (DBS) Signature A, and one InDel (ID2). The most frequent SBS signatures (SBS05, SBS01, SBS15, SBS20, SBS40) were also observed in 254 White GC cases from The Cancer Genome Atlas (TCGA) Project. However, SBS01 and SBS20 showed significant differences between Whites vs. All Asians (19.3% vs. 11.3% for SBS 1 (P = 0.012) and 11.4% vs. 5.9% for SBS20 (P = 0.025), respectively). Using COSMIC V2, signatures 6, 15, and 1 were the most frequent in Chinese GC cases. Further, most Chinese GC cases carried multiple signatures.

CONCLUSIONS

This effort represents the most detailed mutational signatures analysis of GC cases from China to date. Results hold promise for new insights in understanding risk and prognosis factors in GC.

A protracted war against cancer drug resistance

Currently, even the most effective anti-cancer therapies are often limited by the development of drug resistance and tumor relapse, which is a major challenge facing current cancer research. A deep understanding of the molecular and biochemical bases of drug efficacy that can help predict the clinical drug resistance, coupled with the evolution of systematic genomic and proteomic technologies, have facilitated studies identifying and elucidating the underlying mechanisms. In this review, we focus on several important issues on cancer drug resistance and provide a framework for understanding the common ways by which cancers develop resistance to therapeutic agents. With the increasing arsenal of novel anticancer agents and techniques, there are now unprecedented opportunities to understand and overcome drug resistance. The proteolysis targeting chimera (PROTAC) technology, immunotherapy, nanomedicine, and real-time monitoring of drug response all provide effective approaches for combating drug resistance. In addition to the advancement of therapeutic technologies, the revolution of treatment concept is also of great importance. We can take advantage of the interplay between drug sensitive and resistant subclones for combating cancer. However, there remains a long way to go in the protracted war against cancer drug resistance.

Morphomolecular pathology and Genomic Insights into the Cells of Origin of Cholangiocarcinoma and combined Hepatocellular-Cholangiocarcinoma

Cholangiocarcinoma is a highly heterogenous group of malignancies that, despite recent progress in the understanding of its molecular pathogenesis and clinical management, continue to pose a major challenge to public health. The traditional view posits that cholangiocarcinomas derive from the neoplastic transformation of cholangiocytes lining the biliary tree. However, increasing genetic and experimental evidence has recently pointed to a more complex - and nuanced - scenario for the potential cell of origin of cholangiocarcinomas, with hepatocytes as well as hepatic stem/progenitor cells being considered as additional potential sources, depending on microenvironmental contexts including liver injury. The hypothesis of potentially diverse cells of origin for CCA, albeit controversial, is certainly not surprising given the plasticity of the cells populating the liver as well as the existence of liver cancer subtypes with mixed histological and molecular features. This review carefully looks at the current pathological, genomic and experimental evidence supporting the existence of multiple cells of origin of liver and biliary tract cancers, with particular focus on cholangiocarcinoma and combined hepatocellular-cholangiocarcinoma.

Characterisation of APOBEC3B-Mediated RNA editing in breast cancer cells reveals regulatory roles of NEAT1 and MALAT1 lncRNAs

RNA editing is a crucial post-transcriptional process that influences gene expression and increases the diversity of the proteome as a result of amino acid substitution. Recently, the APOBEC3 family has emerged as a significant player in this mechanism, with APOBEC3A (A3A) having prominent roles in base editing during immune and stress responses. APOBEC3B (A3B), another family member, has gained attention for its potential role in generating genomic DNA mutations in breast cancer. In this study, we coupled an inducible expression cell model with a novel methodology for identifying differential variants in RNA (DVRs) to map A3B-mediated RNA editing sites in a breast cancer cell model. Our findings indicate that A3B engages in selective RNA editing including targeting NEAT1 and MALAT1 long non-coding RNAs that are often highly expressed in tumour cells. Notably, the binding of these RNAs sequesters A3B and suppresses global A3B activity against RNA and DNA. Release of A3B from NEAT1/MALAT1 resulted in increased A3B activity at the expense of A3A activity suggesting a regulatory feedback loop between the two family members. This research substantially advances our understanding of A3B's role in RNA editing, its mechanistic underpinnings, and its potential relevance in the pathogenesis of breast cancer.

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.

Clonal dynamics and somatic evolution of haematopoiesis in mouse

Haematopoietic stem cells maintain blood production throughout life. While extensively characterised using the laboratory mouse, little is known about how the population is sustained and evolves with age. We isolated stem cells and progenitors from young and old mice, identifying 221,890 somatic mutations genome-wide in 1845 single cell-derived colonies, and used phylogenetic analysis to infer the ontogeny and population dynamics of the stem cell pool. Mouse stem cells and progenitors accrue ~45 somatic mutations per year, a rate only about 2-fold greater than human progenitors despite the vastly different organismal sizes and lifespans. Phylogenetic patterns reveal that stem and multipotent progenitor cell pools are both established during embryogenesis, after which they independently self-renew in parallel over life. The stem cell pool grows steadily over the mouse lifespan to approximately 70,000 cells, self-renewing about every six weeks. Aged mice did not display the profound loss of stem cell clonal diversity characteristic of human haematopoietic ageing. However, targeted sequencing revealed small, expanded clones in the context of murine ageing, which were larger and more numerous following haematological perturbations and exhibited a selection landscape similar to humans. Our data illustrate both conserved features of population dynamics of blood and distinct patterns of age-associated somatic evolution in the short-lived mouse.

Benchmarking bulk and single-cell variant-calling approaches on Chromium scRNA-seq and scATAC-seq libraries

Single-cell sequencing methodologies such as scRNA-seq and scATAC-seq have become widespread and effective tools to interrogate tissue composition. Increasingly, variant callers are being applied to these methodologies to resolve the genetic heterogeneity of a sample, especially in the case of detecting the clonal architecture of a tumor. Typically, traditional bulk DNA variant callers are applied to the pooled reads of a single-cell library to detect candidate mutations. Recently, multiple studies have applied such callers on reads from individual cells, with some citing the ability to detect rare variants with higher sensitivity. Many studies apply these two approaches to the Chromium (10x Genomics) scRNA-seq and scATAC-seq methodologies. However, Chromium-based libraries may offer additional challenges to variant calling compared with existing single-cell methodologies, raising questions regarding the validity of variants obtained from such a workflow. To determine the merits and challenges of various variant-calling approaches on Chromium scRNA-seq and scATAC-seq libraries, we use sample libraries with matched bulk whole-genome sequencing to evaluate the performance of callers. We review caller performance, finding that bulk callers applied on pooled reads significantly outperform individual-cell approaches. We also evaluate variants unique to scRNA-seq and scATAC-seq methodologies, finding patterns of noise but also potential capture of RNA-editing events. Finally, we review the notion that variant calling at the single-cell level can detect rare somatic variants, providing empirical results that suggest resolving such variants is infeasible in single-cell Chromium libraries.

Frequent CHD1 deletions in prostate cancers of African American men is associated with rapid disease progression

We analyzed genomic data from the prostate cancer of African- and European American men to identify differences contributing to racial disparity of outcome. We also performed FISH-based studies of Chromodomain helicase DNA-binding protein 1 (CHD1) loss on prostate cancer tissue microarrays. We created CHD1-deficient prostate cancer cell lines for genomic, drug sensitivity and functional homologous recombination (HR) activity analysis. Subclonal deletion of CHD1 was nearly three times as frequent in prostate tumors of African American than in European American men and it associates with rapid disease progression. CHD1 deletion was not associated with HR deficiency associated mutational signatures or HR deficiency as detected by RAD51 foci formation. This was consistent with the moderate increase of olaparib and talazoparib sensitivity with several CHD1 deficient cell lines showing talazoparib sensitivity in the clinically relevant concentration range. CHD1 loss may contribute to worse disease outcome in African American men.

Biochemical reconstitution of heat-induced mutational processes

Non-enzymatic spontaneous deamination of 5-methylcytosine, producing thymine, is the proposed etiology of cancer mutational signature 1, which is the most predominant signature in all cancers. Here, the proposed mutational process was reconstituted using synthetic DNA and purified proteins. First, single-stranded DNA containing 5-methylcytosine at CpG context was incubated at an elevated temperature to accelerate spontaneous DNA damage. Then, the DNA was treated with uracil DNA glycosylase to remove uracil residues that were formed by deamination of cytosine. The resulting DNA was then used as a template for DNA synthesis by yeast DNA polymerase δ. The DNA products were analyzed by next-generation DNA sequencing, and mutation frequencies were quantified. The observed mutations after this process were exclusively C>T mutations at CpG context, which was very similar to signature 1. When 5-methylcytosine modification and uracil DNA glycosylase were both omitted, C>T mutations were produced on C residues in all sequence contexts, but these mutations were diminished by uracil DNA glycosylase-treatment. These results indicate that the CpG>TpG mutations were produced by the deamination of 5-methylcytosine. Additional mutations, mainly C>G, were introduced by yeast DNA polymerase ζ on the heat-damaged DNA, indicating that G residues of the templates were also damaged. However, the damage on G residues was not converted to mutations with DNA polymerase δ or ε.

Pervasive Induction of Regulatory Mutation Microclones in Sun-exposed Skin

Carcinogen-induced mutations are thought near-random, with rare cancer-driver mutations underlying clonal expansion. Using high-fidelity Duplex Sequencing to reach a mutation frequency sensitivity of 4×10 -9 per nt, we report that sun exposure creates pervasive mutations at sites with ∼100-fold UV-sensitivity in RNA-processing gene promoters - cyclobutane pyrimidine dimer (CPD) hyperhotspots - and these mutations have a mini-driver clonal expansion phenotype. Numerically, human skin harbored 10-fold more genuine mutations than previously reported, with neonatal skin containing 90,000 per cell; UV signature mutations increased 8,000-fold in sun-exposed skin, averaging 3×10 -5 per nt. Clonal expansion by neutral drift or passenger formation was nil. Tumor suppressor gene hotspots reached variant allele frequency 0.1-10% via 30-3,000 fold clonal expansion, in occasional biopsies. CPD hyperhotspots reached those frequencies in every biopsy, with modest clonal expansion. In vitro, tumor hotspot mutations arose occasionally over weeks of chronic low-dose exposure, whereas CPD hyperhotspot mutations arose in days at 1000-fold higher frequencies, growing exponentially. UV targeted mini-drivers in every skin cell.

Pan-tumor validation of a NGS fraction-based MSI analysis as a predictor of response to Pembrolizumab

Microsatellite instability high (MSI-H) and mismatch repair deficient (dMMR) tumor status have been demonstrated to predict patient response to immunotherapies. We developed and validated a next-generation sequencing (NGS)-based companion diagnostic (CDx) to detect MSI-H solid tumors via a comprehensive genomic profiling (CGP) assay, FoundationOne®CDx (F1CDx). To determine MSI status, F1CDx calculates the fraction of unstable microsatellite loci across >2000 loci using a fraction-based (FB) analysis. Across solid tumor types, F1CDx demonstrated a high analytical concordance with both PCR (n = 264) and IHC (n = 279) with an overall percent agreement (OPA) of 97.7% and 97.8%, respectively. As part of a retrospective bridging clinical study from KEYNOTE-158 Cohort K and KEYNOTE-164, patients with MSI-H tumors as determined by F1CDx demonstrated an objective response rate (ORR) of 43.0% to pembrolizumab. In real-world cancer patients from a deidentified clinicogenomic database, F1CDx was at least equivalent in assessing clinical outcome following immunotherapy compared with MMR IHC. Demonstrated analytical and clinical performance of F1CDx led to the pan-tumor FDA approval in 2022 of F1CDx to identify MSI-H solid tumor patients for treatment with pembrolizumab. F1CDx is an accurate, reliable, and FDA-approved method for the identification of MSI-H tumors for treatment with pembrolizumab.

Somatic mutations in 3929 HPV positive cervical cells associated with infection outcome and HPV type

Invasive cervical cancers (ICC), caused by HPV infections, have a heterogeneous molecular landscape. We investigate the detection, timing, and HPV type specificity of somatic mutations in 3929 HPV-positive exfoliated cervical cell samples from individuals undergoing cervical screening in the U.S. using deep targeted sequencing in ICC cases, precancers, and HPV-positive controls. We discover a subset of hotspot mutations rare in controls (2.6%) but significantly more prevalent in precancers, particularly glandular precancer lesions (10.2%), and cancers (25.7%), supporting their involvement in ICC carcinogenesis. Hotspot mutations differ by HPV type, and HPV18/45-positive ICC are more likely to have multiple hotspot mutations compared to HPV16-positive ICC. The proportion of cells containing hotspot mutations is higher (i.e., higher variant allele fraction) in ICC and mutations are detectable up to 6 years prior to cancer diagnosis. Our findings demonstrate the feasibility of using exfoliated cervical cells for detection of somatic mutations as potential diagnostic biomarkers.

Prevalence of germline variants in Brazilian pancreatic carcinoma patients

We evaluated the prevalence of pathogenic/likely pathogenic germline variants (PGV) in Brazilian pancreatic adenocarcinoma (PC) patients, that represent a multiethnic population, in a cross-sectional study. We included 192 PC patients unselected for family history of cancer. We evaluated a panel of 113 cancer genes, through genomic DNA sequencing and 46 ancestry-informative markers, through multiplex PCR. The median age was 61 years; 63.5% of the patients presented disease clinical stages III or IV; 8.3% reported personal history of cancer; 4.7% and 16.1% reported first-degree relatives with PC or breast and/or prostate cancer, respectively. Although the main ancestry was European, there was considerable genetic composition admixture. Twelve patients (6.25%) were PGV carriers in PC predisposition genes (ATM, BRCA1, BRCA2, CDKN2A, MSH2, PALB2) and another 25 (13.0%) were PGV carriers in genes with a limited association or not previously associated with PC (ACD, BLM, BRIP1, CHEK2, ERCC4, FANCA, FANCE, FANCM, GALNT12, MITF, MRE11, MUTYH, POLE, RAD51B, RAD51C, RECQL4, SDHA, TERF2IP). The most frequently affected genes were CHEK2, ATM and FANC. In tumor samples from PGV carriers in ACD, BRIP1, MRE11, POLE, SDHA, TERF2IP, which were examined through exome sequencing, the main single base substitutions (SBS) mutational signature was SBS1+5+18, probably associated with age, tobacco smoking and reactive oxygen species. SBS3 associated with homologous repair deficiency was also represented, but on a lower scale. There was no difference in the frequency of PGV carriers between: (a) patients with or without first-degree relatives with cancer; and (b) patients with admixed ancestry versus those with predominantly European ancestry. Furthermore, there was no difference in overall survival between PGV carriers and non-carriers. Therefore, genetic testing should be offered to all Brazilian pancreatic cancer patients, regardless of their ancestry. Genes with limited or previously unrecognized associations with pancreatic cancer should be further investigated to clarify their role in cancer risk.

Clinical Interest in Exome-Based Analysis of Somatic Mutational Signatures for Non-Small Cell Lung Cancer

BACKGROUND

Non-small cell lung cancer (NSCLC) remains the leading cause of cancer-related mortality. This study investigates the clinical interest of whole exome sequencing (WES) for analyzing somatic mutational signatures in patients with advanced or metastatic NSCLC treated with the current standard of care.

METHODS

Exome sequencing data and clinical characteristics from 132 patients with advanced or metastatic NSCLC were analyzed. Somatic mutational signatures including single base substitutions (SBSs), double base substitutions (DBSs), and copy number signatures were evaluated. Structural variations including tumor mutational burden (TMB), the number of neoantigens, TCR clonality, homologous recombination deficiency (HRD), copy number alterations (CNAs), and microsatellite instability (MSI) score were determined. The association between these genomic features, NSCLC subtypes, and patient outcomes (progression-free and overall survival) was evaluated.

CONCLUSIONS

Exome sequencing offers valuable insights into somatic mutational signatures in NSCLC. This study identified specific signatures associated with a poor response to immune checkpoint inhibitor (ICI) therapy and chemotherapy, potentially aiding treatment selection and identifying patients unlikely to benefit from these approaches.

MADDD-seq, a novel massively parallel sequencing tool for simultaneous detection of DNA damage and mutations

Our genome is exposed to a wide variety of DNA-damaging agents. If left unrepaired, this damage can be converted into mutations that promote carcinogenesis or the development of genetically inherited diseases. As a result, researchers and clinicians require tools that can detect DNA damage and mutations with exceptional sensitivity. In this study, we describe a massively parallel sequencing tool termed Mutation And DNA Damage Detection-seq (MADDD-seq) that is capable of detecting O6-methyl guanine lesions and mutations simultaneously, with a single assay. To illustrate the dual capabilities of MADDD-seq, we treated WT and DNA repair deficient yeast cells with the DNA-damaging agent MNNG and tracked DNA lesions and mutations over a 24-h time period. This approach allowed us to identify thousands of DNA adducts and mutations in a single sequencing run and gain deep insight into the kinetics of DNA repair and mutagenesis.

RAD52 resolves transcription-replication conflicts to mitigate R-loop induced genome instability

Collisions of the transcription and replication machineries on the same DNA strand can pose a significant threat to genomic stability. These collisions occur in part due to the formation of RNA-DNA hybrids termed R-loops, in which a newly transcribed RNA molecule hybridizes with the DNA template strand. This study investigated the role of RAD52, a known DNA repair factor, in preventing collisions by directing R-loop formation and resolution. We show that RAD52 deficiency increases R-loop accumulation, exacerbating collisions and resulting in elevated DNA damage. Furthermore, RAD52's ability to interact with the transcription machinery, coupled with its capacity to facilitate R-loop dissolution, highlights its role in preventing collisions. Lastly, we provide evidence of an increased mutational burden from double-strand breaks at conserved R-loop sites in human tumor samples, which is increased in tumors with low RAD52 expression. In summary, this study underscores the importance of RAD52 in orchestrating the balance between replication and transcription processes to prevent collisions and maintain genome stability.

Endometrial carcinomas with ambiguous histology often harbor TP53 mutations

The objective of the present study was to characterize the molecular features of endometrial carcinomas with ambiguous histology. Eighteen carcinomas that could not be conclusively typed based on morphology and immunohistochemistry underwent analysis of mismatch repair (MMR) status, microsatellite status, and whole-exome sequencing. None of the tumors had pathogenic POLE mutation. Twelve tumors (67%) were microsatellite stable, and 6 (33%) had microsatellite instability. Fourteen tumors (78%) harbored TP53 mutations, and 2 (11%) had mutations in MMR genes. Eleven carcinomas (61%) were classified as copy number high and 7 (39%) as MSI-hypermutated, the latter including 3 tumors with TP53 mutation who concomitantly had MSI or mutation in a MMR gene. Other mutations that were found in > 1 tumor affected MUC16 (7 tumors), PIK3CA (6 tumors), PPP2R1A (6 tumors), ARID1A (5 tumors), PTEN (5 tumors), FAT1 (4 tumors), FAT4 (3 tumors), BRCA2 (2 tumors), ERBB2 (2 tumors), FBXW7 (2 tumors), MET (2 tumors), MTOR (2 tumors), JAK1 (2 tumors), and CSMD3 (2 tumors). At the last follow-up (median = 68.6 months), 8 patients had no evidence of disease, 1 patient was alive with disease, 8 patients were dead of disease, and 1 patient died of other cause. In conclusion, based on this series, the molecular landscape of endometrial carcinomas with ambiguous histology is dominated by TP53 mutations and the absence of POLE mutations, with heterogeneous molecular profile with respect to other genes. A high proportion of these tumors is clinically aggressive.

Large-scale analysis of whole genome sequencing data from formalin-fixed paraffin-embedded cancer specimens demonstrates preservation of clinical utility

Whole genome sequencing (WGS) provides comprehensive, individualised cancer genomic information. However, routine tumour biopsies are formalin-fixed and paraffin-embedded (FFPE), damaging DNA, historically limiting their use in WGS. Here we analyse FFPE cancer WGS datasets from England's 100,000 Genomes Project, comparing 578 FFPE samples with 11,014 fresh frozen (FF) samples across multiple tumour types. We use an approach that characterises rather than discards artefacts. We identify three artefactual signatures, including one known (SBS57) and two previously uncharacterised (SBS FFPE, ID FFPE), and develop an "FFPEImpact" score that quantifies sample artefacts. Despite inferior sequencing quality, FFPE-derived data identifies clinically-actionable variants, mutational signatures and permits algorithmic stratification. Matched FF/FFPE validation cohorts shows good concordance while acknowledging SBS, ID and copy-number artefacts. While FF-derived WGS data remains the gold standard, FFPE-samples can be used for WGS if required, using analytical advancements developed here, potentially democratising whole cancer genomics to many.

Mutational mechanisms in multiply relapsed pediatric acute lymphoblastic leukemia

Pediatric acute lymphoblastic leukemia (ALL) is marked by low mutational load at initial diagnosis, which increases at relapse. To determine which processes are active in (relapsed) ALL and how they behave during disease progression before and after therapy, we performed whole genome sequencing on 97 tumor samples of 29 multiply relapsed ALL patients. Mutational load increased upon relapse in 28 patients and upon every subsequent relapse in 22 patients. In addition to two clock-like mutational processes, we identified UV-like damage, APOBEC activity, reactive oxygen species, thiopurine-associated damage and an unknown therapy component as drivers of mutagenesis. Mutational processes often affected patients over longer time periods, but could also occur in isolated events, suggesting the requirement of additional triggers. Thiopurine exposure was the most prominent source of new mutations in relapse, affecting over half of the studied patients in first and/or later relapse and causing potential relapse-driving mutations in multiple patients. Our data demonstrate that multiple mutational processes frequently act in parallel as prominent secondary drivers with dynamic activity during ALL development and progression.

Star wars against leukemia: attacking the clones

Leukemia, although most likely starts as a monoclonal genetic/epigenetic anomaly, is a polyclonal disease at manifestation. This polyclonal nature results from ongoing evolutionary changes in the genome/epigenome of leukemia cells to promote their survival and proliferation advantages. We discuss here how genetic and/or epigenetic aberrations alter intracellular microenvironment in individual leukemia clones and how extracellular microenvironment selects the best fitted clones. This dynamic polyclonal composition of leukemia makes designing an effective therapy a challenging task especially because individual leukemia clones often display substantial differences in response to treatment. Here, we discuss novel therapeutic approach employing single cell multiomics to identify and eradicate all individual clones in a patient.

Image-based deep learning model using DNA methylation data predicts the origin of cancer of unknown primary

Cancer of unknown primary (CUP) is a rare type of metastatic cancer in which the origin of the tumor is unknown. Since the treatment strategy for patients with metastatic tumors depends on knowing the primary site, accurate identification of the origin site is important. Here, we developed an image-based deep-learning model that utilizes a vision transformer algorithm for predicting the origin of CUP. Using DNA methylation dataset of 8,233 primary tumors from The Cancer Genome Atlas (TCGA), we categorized 29 cancer types into 18 organ classes and extracted 2,312 differentially methylated CpG sites (DMCs) from non-squamous cancer group and 420 DMCs from squamous cell cancer group. Using these DMCs, we created organ-specific DNA methylation images and used them for model training and testing. Model performance was evaluated using 394 metastatic cancer samples from TCGA (TCGA-meta) and 995 samples (693 primary and 302 metastatic cancers) obtained from 20 independent external studies. We identified that the DNA methylation image reveals a distinct pattern based on the origin of cancer. Our model achieved an overall accuracy of 96.95 % in the TCGA-meta dataset. In the external validation datasets, our classifier achieved overall accuracies of 96.39 % and 94.37 % in primary and metastatic tumors, respectively. Especially, the overall accuracies for both primary and metastatic samples of non-squamous cell cancer were exceptionally high, with 96.79 % and 96.85 %, respectively.

Ultra-sensitive molecular residual disease detection through whole genome sequencing with single-read error correction

While whole genome sequencing (WGS) of cell-free DNA (cfDNA) holds enormous promise for detection of molecular residual disease (MRD), its performance is limited by WGS error rate. Here we introduce AccuScan, an efficient cfDNA WGS technology that enables genome-wide error correction at single read-level, achieving an error rate of 4.2 × 10-7, which is about two orders of magnitude lower than a read-centric de-noising method. The application of AccuScan to MRD demonstrated analytical sensitivity down to 10-6 circulating variant allele frequency at 99% sample-level specificity. AccuScan showed 90% landmark sensitivity (within 6 weeks after surgery) and 100% specificity for predicting relapse in colorectal cancer. It also showed 67% sensitivity and 100% specificity in esophageal cancer using samples collected within one week after surgery. When AccuScan was applied to monitor immunotherapy in melanoma patients, the circulating tumor DNA (ctDNA) levels and dynamic profiles were consistent with clinical outcomes. Overall, AccuScan provides a highly accurate WGS solution for MRD detection, empowering ctDNA detection at parts per million range without requiring high sample input or personalized reagents.