In-depth characterization of the cisplatin mutational signature in human cell lines and in esophageal and liver tumors

Prenatal Exposure to Chemotherapy Increases the Mutation Burden in Human Neonatal Hematopoietic Stem Cells

SIGNIFICANCE

This study demonstrates that environmental mutagenic exposure during pregnancy can increase somatic mutation accumulation in the fetus. Given that detrimental early life exposures can adversely affect health outcomes later in life, our study highlights the need for further research into the impact of environmentally induced genomic insults during the perinatal period. See related commentary by Furudate and Takahashi, p. 870.

Error-corrected flow-based sequencing at whole-genome scale and its application to circulating cell-free DNA profiling

Differentiating sequencing errors from true variants is a central genomics challenge, calling for error suppression strategies that balance costs and sensitivity. For example, circulating cell-free DNA (ccfDNA) sequencing for cancer monitoring is limited by sparsity of circulating tumor DNA, abundance of genomic material in samples and preanalytical error rates. Whole-genome sequencing (WGS) can overcome the low abundance of ccfDNA by integrating signals across the mutation landscape, but higher costs limit its wide adoption. Here, we applied deep (~120×) lower-cost WGS (Ultima Genomics) for tumor-informed circulating tumor DNA detection within the part-per-million range. We further leveraged lower-cost sequencing by developing duplex error-corrected WGS of ccfDNA, achieving 7.7 × 10-7 error rates, allowing us to assess disease burden in individuals with melanoma and urothelial cancer without matched tumor sequencing. This error-corrected WGS approach will have broad applicability across genomics, allowing for accurate calling of low-abundance variants at efficient cost and enabling deeper mapping of somatic mosaicism as an emerging central aspect of aging and disease.

Recurrent spontaneous miscarriages from sperm after ABVD chemotherapy in a patient with Hodgkin's lymphoma: sperm DNA and methylation profiling

ABSTRACT

Lymphomas represent one of the most common malignant diseases in young men and an important issue is how treatments will affect their reproductive health. It has been hypothesized that chemotherapies, similarly to environmental chemicals, may alter the spermatogenic epigenome. Here, we report the genomic and epigenomic profiling of the sperm DNA from a 31-year-old Hodgkin lymphoma patient who faced recurrent spontaneous miscarriages in his couple 11-26 months after receiving chemotherapy with adriamycin, bleomycin, vinblastine, and dacarbazine (ABVD). In order to capture the potential deleterious impact of the ABVD treatment on mutational and methylation changes, we compared sperm DNA before and 26 months after chemotherapy with whole-genome sequencing (WGS) and reduced representation bisulfite sequencing (RRBS). The WGS analysis identified 403 variants following ABVD treatment, including 28 linked to genes crucial for embryogenesis. However, none were found in coding regions, indicating no impact of chemotherapy on protein function. The RRBS analysis identified 99 high-quality differentially methylated regions (hqDMRs) for which methylation status changed upon chemotherapy. Those hqDRMs were associated with 87 differentially methylated genes, among which 14 are known to be important or expressed during embryo development. While no variants were detected in coding regions, promoter regions of several genes potentially important for embryo development contained variants or displayed an altered methylated status. These might in turn modify the corresponding gene expression and thus affect their function during key stages of embryogenesis, leading to potential developmental disorders or miscarriages.

Structural variant and nucleosome occupancy dynamics postchemotherapy in a HER2+ breast cancer organoid model

The most common chemotherapeutics induce DNA damage to eradicate cancer cells, yet defective DNA repair can propagate mutations, instigating therapy resistance and secondary malignancies. Structural variants (SVs), arising from copy-number-imbalanced and -balanced DNA rearrangements, are a major driver of tumor evolution, yet understudied posttherapy. Here, we adapted single-cell template-strand sequencing (Strand-seq) to a HER2+ breast cancer model to investigate the formation of doxorubicin-induced de novo SVs. We coupled this approach with nucleosome occupancy (NO) measurements obtained from the same single cell to enable simultaneous SV detection and cell-type classification. Using organoids from TetO-CMYC/TetO-Neu/MMTV-rtTA mice modeling HER2+ breast cancer, we generated 459 Strand-seq libraries spanning various tumorigenesis stages, identifying a 7.4-fold increase in large chromosomal alterations post-doxorubicin. Complex DNA rearrangements, deletions, and duplications were prevalent across basal, luminal progenitor (LP), and mature luminal (ML) cells, indicating uniform susceptibility of these cell types to SV formation. Doxorubicin further elevated sister chromatid exchanges (SCEs), indicative of genomic stress persisting posttreatment. Altered nucleosome occupancy levels on distinct cancer-related genes further underscore the broad genomic impact of doxorubicin. The organoid-based system for single-cell multiomics established in this study paves the way for unraveling the most important therapy-associated SV mutational signatures, enabling systematic studies of the effect of therapy on cancer evolution.

Mutational Signatures in Colorectal Cancer: Translational Insights, Clinical Applications, and Limitations

A multitude of exogenous and endogenous processes have the potential to result in DNA damage. While the repair mechanisms are typically capable of correcting this damage, errors in the repair process can result in mutations. The findings of research conducted in 2012 indicate that mutations do not occur randomly but rather follow specific patterns that can be attributed to known or inferred mutational processes. The process of mutational signature analysis allows for the inference of the predominant mutational process for a given cancer sample, with significant potential for clinical applications. A deeper comprehension of these mutational signatures in CRC could facilitate enhanced prevention strategies, facilitate the comprehension of genotoxic drug activity, predict responses to personalized treatments, and, in the future, inform the development of targeted therapies in the context of precision oncology. The efforts of numerous researchers have led to the identification of several mutational signatures, which can be categorized into different mutational signature references. In CRC, distinct mutational signatures are identified as correlating with mismatch repair deficiency, polymerase mutations, and chemotherapy treatment. In this context, a mutational signature analysis offers considerable potential for enhancing minimal residual disease (MRD) tests in stage II (high-risk) and stage III CRC post-surgery, stratifying CRC based on the impacts of genetic and epigenetic alterations for precision oncology, identifying potential therapeutic vulnerabilities, and evaluating drug efficacy and guiding therapy, as illustrated in a proof-of-concept clinical trial.

Head and neck squamous cell carcinomas of unknown primary: Can ancillary studies help identify more primary tumor sites?

A subset of head and neck squamous cell carcinomas present solely as metastatic disease in the neck and are of unknown primary origin (SCCUP). Most primary tumors will ultimately be identified, usually in the oropharynx. In a minority of cases, the primary site remains elusive. Here, we examine the role of ancillary testing, including mutational signature analysis (MSA), to help identify likely primary sites in such cases. Twenty-two cases of SCCUP in the neck, collected over a 10-year period, were classified by morphology and viral status; including human papillomavirus (HPV) testing by p16 immunohistochemistry (IHC) and RT-qPCR, as well as Epstein-Barr virus (EBV) testing by EBER-ISH. CD5 and c-KIT (CD117) IHC was done to evaluate for possible thymic origin in all virus-negative cases. Whole exome sequencing, followed by MSA, was used to identify UV signature mutations indicative of cutaneous origin. HPV was identified in 12 of 22 tumors (54.5%), favoring an oropharyngeal origin, and closely associated with nonkeratinizing tumor morphology (Fisher's exact test; p = 0.0002). One tumor with indeterminant morphology had discordant HPV and p16 status (p16+/HPV-). All tumors were EBV-negative. Diffuse expression of CD5 and c-KIT was identified in 1 of 10 virus-negative SCCUPs (10%), suggesting a possible ectopic thymic origin rather than a metastasis. A UV mutational signature, indicating cutaneous origin, was identified in 1 of 10 (10%) virus-negative SCCUPs. A cutaneous auricular primary emerged 3 months after treatment in this patient. Primary tumors became clinically apparent in 2 others (1 hypopharynx, 1 hypopharynx/larynx). Thus, after follow-up, 6 tumors remained unclassifiable as to the possible site of origin (27%). Most SCCUPs of the neck in our series were HPV-associated and thus likely of oropharyngeal origin. UV signature mutation analysis and additional IHC for CD5 and c-KIT for possible thymic origin may aid in further classifying virus-negative unknown primaries. Close clinical inspection of hypopharyngeal mucosa may also be helpful, as a subset of primary tumors later emerged at this site.

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.

Spatial intra-tumour heterogeneity and treatment-induced genomic evolution in oesophageal adenocarcinoma: implications for prognosis and therapy

BACKGROUND

Oesophageal adenocarcinoma (OAC) is a highly heterogeneous cancer with poor survival. Standard curative treatment is chemotherapy with or without radiotherapy followed by oesophagectomy. Genomic heterogeneity is a feature of OAC and has been linked to treatment resistance.

METHODS

Whole-genome sequencing data from 59 treatment-naïve and 18 post-treatment samples from 29 OAC patients was analysed. Twenty-seven of these were enrolled in the DOCTOR trial, sponsored by the Australasian Gastro-Intestinal Trials Group. Two biopsies from each treatment-naïve tumour were assessed to define 'shared' (between both samples) and 'private' (present in one sample) mutations.

RESULTS

Mutational signatures SBS2/13 (APOBEC) and SBS3 (BRCA) were almost exclusively detected in private mutation populations of treatment-naïve tumours. Patients presenting these signatures had significantly worse disease specific survival. Furthermore, mutational signatures associated with platinum-based chemotherapy treatment as well as high platinum enrichment scores were only detected in post-treatment samples. Additionally, clones with high putative neoantigen binding scores were detected in some treatment-naïve samples suggesting immunoediting of clones.

CONCLUSIONS

This study demonstrates the high intra-tumour heterogeneity in OAC, as well as indicators for treatment-induced changes during tumour evolution. Intra-tumour heterogeneity remains a problem for successful treatment strategies in OAC.

Selective pressures of platinum compounds shape the evolution of therapy-related myeloid neoplasms

Therapy-related myeloid neoplasms (t-MN) arise as a complication of chemo- and/or radiotherapy. Although t-MN can occur both in adult and childhood cancer survivors, the mechanisms driving therapy-related leukemogenesis likely vary across different ages. Chemotherapy is thought to induce driver mutations in children, whereas in adults pre-existing mutant clones are selected by the exposure. However, selective pressures induced by chemotherapy early in life are less well studied. Here, we use single-cell whole genome sequencing and phylogenetic inference to show that the founding cell of t-MN in children starts expanding after cessation of platinum exposure. In patients with Li-Fraumeni syndrome, characterized by a germline TP53 mutation, we find that the t-MN already expands during treatment, suggesting that platinum-induced growth inhibition is TP53-dependent. Our results demonstrate that germline aberrations can interact with treatment exposures in inducing t-MN, which is important for the development of more targeted, patient-specific treatment regimens and follow-up.

Pan-cancer analysis of heterogeneity of tumor mutational burden and genomic mutation under treatment pressure

BACKGROUND

High tumor mutational burden (TMB) is one of the widely researched predictive biomarkers of immune checkpoint inhibitors and has been shown to be closely related with response to immunotherapy in multiple cancer types. However, for patients who have failed conventional therapy and are about to undergo immunotherapy, there is no consensus recommendation on the timing of tumor sampling for TMB analysis, and the effects of different therapies on TMB have not been clarified. This retrospective observational study aimed to investigate the heterogeneity of TMB and genomic mutation under the treatment pressure.

PATIENTS AND METHODS

We retrospectively collected the available genomic and therapeutic information from 8051 samples across 15 tumor types (>50 samples/tumor) found in 30 published studies and investigated the distribution and heterogeneity of TMB under treatment across diverse cohorts.

RESULTS

This integrated analysis has shown anticancer treatments increased TMB. Significant effects of treatment on TMB were more frequently observed in tumor types with lower treatment-naïve TMB, including breast, prostate, and pediatric cancers. For different cancer therapies, chemotherapy was prone to be correlated with an increased TMB in most cancer types. Meanwhile, the fraction of the TMB-high category of breast, prostate, and bladder cancers and glioma increased significantly after chemotherapy. Several actionable genes including ERS1 and NF1 in breast cancer, as well as some prognostic markers including TERT in bladder cancer and IDH1 in glioma, were significantly changed in post-chemotherapy tumors compared to treatment-naïve tumors.

CONCLUSION

Our study reveals the heterogeneity of TMB under treatment across diverse cancer types and provides evidences that chemotherapy was associated with increases in TMB as well as the fraction of TMB-high category, suggesting that resampling tumor tissues for calculating post-chemotherapy TMB could be a better option for predicting the response to immunotherapy, especially for tumors with initially low TMB.

Cisplatin in the era of PARP inhibitors and immunotherapy

Platinum compounds such as cisplatin, carboplatin and oxaliplatin are widely used in chemotherapy. Cisplatin induces cytotoxic DNA damage that blocks DNA replication and gene transcription, leading to arrest of cell proliferation. Although platinum therapy alone is effective against many tumors, cancer cells can adapt to the treatment and gain resistance. The mechanisms for cisplatin resistance are complex, including low DNA damage formation, high DNA repair capacity, changes in apoptosis signaling pathways, rewired cell metabolisms, and others. Drug resistance compromises the clinical efficacy and calls for new strategies by combining cisplatin with other therapies. Exciting progress in cancer treatment, particularly development of poly (ADP-ribose) polymerase (PARP) inhibitors and immune checkpoint inhibitors, opened a new chapter to combine cisplatin with these new cancer therapies. In this Review, we discuss how platinum synergizes with PARP inhibitors and immunotherapy to bring new hope to cancer patients.

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.

Deep whole-genome analysis of 494 hepatocellular carcinomas

Over half of hepatocellular carcinoma (HCC) cases diagnosed worldwide are in China1-3. However, whole-genome analysis of hepatitis B virus (HBV)-associated HCC in Chinese individuals is limited4-8, with current analyses of HCC mainly from non-HBV-enriched populations9,10. Here we initiated the Chinese Liver Cancer Atlas (CLCA) project and performed deep whole-genome sequencing (average depth, 120×) of 494 HCC tumours. We identified 6 coding and 28 non-coding previously undescribed driver candidates. Five previously undescribed mutational signatures were found, including aristolochic-acid-associated indel and doublet base signatures, and a single-base-substitution signature that we termed SBS_H8. Pentanucleotide context analysis and experimental validation confirmed that SBS_H8 was distinct to the aristolochic-acid-associated SBS22. Notably, HBV integrations could take the form of extrachromosomal circular DNA, resulting in elevated copy numbers and gene expression. Our high-depth data also enabled us to characterize subclonal clustered alterations, including chromothripsis, chromoplexy and kataegis, suggesting that these catastrophic events could also occur in late stages of hepatocarcinogenesis. Pathway analysis of all classes of alterations further linked non-coding mutations to dysregulation of liver metabolism. Finally, we performed in vitro and in vivo assays to show that fibrinogen alpha chain (FGA), determined as both a candidate coding and non-coding driver, regulates HCC progression and metastasis. Our CLCA study depicts a detailed genomic landscape and evolutionary history of HCC in Chinese individuals, providing important clinical implications.

Mutational signature, cancer driver genes mutations and transcriptomic subgroups predict hepatoblastoma survival

BACKGROUND

Hepatoblastoma is the most frequent pediatric liver cancer. The current treatments lead to 80% of survival rate at 5 years. In this study, we evaluated the clinical relevance of molecular features to identify patients at risk of chemoresistance, relapse and death of disease.

METHODS

All the clinical data of 86 children with hepatoblastoma were retrospectively collected. Pathological slides were reviewed, tumor DNA sequencing (by whole exome, whole genome or target) and transcriptomic profiling with RNAseq or 300-genes panel were performed. Associations between the clinical, pathological, mutational and transcriptomic data were investigated.

RESULTS

High-risk patients represented 44% of our series and the median age at diagnosis was 21.9 months (range: 0-208). Alterations of the WNT/ß-catenin pathway and of the 11p15.5 imprinted locus were identified in 98% and 74% of the tumors, respectively. Other cancer driver genes mutations were only found in less than 11% of tumors. After neoadjuvant chemotherapy, disease-specific survival and poor response to neoadjuvant chemotherapy were associated with 'Liver Progenitor' (p = 0.00049, p < 0.0001) and 'Immune Cold' (p = 0.0011, p < 0.0001) transcriptomic tumor subtypes, SBS35 cisplatin mutational signature (p = 0.018, p = 0.001), mutations in rare cancer driver genes (p = 0.0039, p = 0.0017) and embryonal predominant histological type (p = 0.0013, p = 0.0077), respectively. Integration of the clinical and molecular features revealed a cluster of molecular markers associated with resistance to chemotherapy and survival, enlightening transcriptomic 'Immune Cold' and Liver Progenitor' as a predictor of survival independent of the clinical features.

CONCLUSIONS

Response to neoadjuvant chemotherapy and survival in children treated for hepatoblastoma are associated with genomic and pathological features independently of the clinical features.

Enrichment of colibactin-associated mutational signatures in unexplained colorectal polyposis patients

BACKGROUND

Colibactin, a genotoxin produced by polyketide synthase harboring (pks+) bacteria, induces double-strand breaks and chromosome aberrations. Consequently, enrichment of pks+Escherichia coli in colorectal cancer and polyposis suggests a possible carcinogenic effect in the large intestine. Additionally, specific colibactin-associated mutational signatures; SBS88 and ID18 in the Catalogue of Somatic Mutations in Cancer database, are detected in colorectal carcinomas. Previous research showed that a recurrent APC splice variant perfectly fits SBS88.

METHODS

In this study, we explore the presence of colibactin-associated signatures and fecal pks in an unexplained polyposis cohort. Somatic targeted Next-Generation Sequencing (NGS) was performed for 379 patients. Additionally, for a subset of 29 patients, metagenomics was performed on feces and mutational signature analyses using Whole-Genome Sequencing (WGS) on Formalin-Fixed Paraffin Embedded (FFPE) colorectal tissue blocks.

RESULTS

NGS showed somatic APC variants fitting SBS88 or ID18 in at least one colorectal adenoma or carcinoma in 29% of patients. Fecal metagenomic analyses revealed enriched presence of pks genes in patients with somatic variants fitting colibactin-associated signatures compared to patients without variants fitting colibactin-associated signatures. Also, mutational signature analyses showed enrichment of SBS88 and ID18 in patients with variants fitting these signatures in NGS compared to patients without.

CONCLUSIONS

These findings further support colibactins ability to mutagenize colorectal mucosa and contribute to the development of colorectal adenomas and carcinomas explaining a relevant part of patients with unexplained polyposis.

Extensive mutational ctDNA profiles reflect High-grade serous cancer tumors and reveal emerging mutations at recurrence

OBJECTIVE

Circulating tumor DNA (ctDNA) offers a minimally-invasive alternative to study genomic changes in recurrent malignancies. With a high recurrence rate, the overall survival in high-grade serous ovarian carcinoma (HGSC) has remained low. Our objectives were to determine whether ctDNA from plasma adequately represents HGSC, and to find mutational changes at relapse suggesting therapy options that could alter patient outcome.

METHODS

We collected 152 longitudinal plasma and 92 fresh tissue samples from 29 HGSC patients, sequencing and detecting mutations with a gene panel of more than 700 cancer-related genes. Tumor content was measured using TP53 VAF. We analyzed the concordance between the mutations in tissue and plasma samples and calculated correlations to patient outcomes. We also searched for novel mutations appearing at relapse.

RESULTS

The concordance rate between mutations in plasma compared to tumor tissue was 83 % at diagnosis and 90 % at relapse. CtDNA was released similarly from the tubo-ovarian tumors, intra-abdominal metastases and ascites. CtDNA release was high when CA-125 level was elevated. The TP53 VAF in ctDNA from plasma samples before the third cycle of primary chemotherapy showed a negative correlation to patient outcome. At relapse, 19 novel, pathogenic DNA mutations appeared, suggesting possible actionable alterations and biological mechanisms related to chemoresistance.

CONCLUSION

Relapse ctDNA samples reflect tissue samples well and longitudinal sampling provides a timely source for mutational profiling. The emerging genetic mutations at recurrence propose that ctDNA accurately represents the widespread disease and provides possibilities for personalized therapy options.

5-Aza-4'-thio-2'-deoxycytidine induces C&gt;G transversions in a specific trinucleotide context and leads to acute lymphoid leukemia

DNA methyltransferase inhibitors (DNMTi), most commonly cytidine analogs, are compounds that are used clinically to decrease 5'-cytosine methylation, with the aim of re-expression of tumor suppressor genes. We used a murine pre-clinical model of myelodysplastic syndrome based on transplantation of cells expressing a NUP98::HOXD13 transgene to investigate 5-Aza-4'-thio-2'-deoxycytidine (Aza TdCyd or ATC), a thiol substituted DNMTi, as a potential therapy. We found that ATC treatment led to lymphoid leukemia in wild-type recipient cells; further study revealed that healthy mice treated with ATC also developed lymphoid leukemia. Whole exome sequencing revealed thousands of acquired mutations, almost all of which were C > G transversions in a previously unrecognized, specific 5'-NCG-3' context. These mutations involved dozens of genes well-known to be involved in human lymphoid leukemia, such as Notch1, Pten, Pax5, Trp53 , and Nf1 . Treatment of human cells in vitro showed thousands of 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 that DNMTi can be potent mutagens in human and mouse cells, both in vitro and in vivo .

Cisplatin exposure alters tRNA-derived small RNAs but does not affect epimutations in C. elegans

BACKGROUND

The individual lifestyle and environment of an organism can influence its phenotype and potentially the phenotype of its offspring. The different genetic and non-genetic components of the inheritance system and their mutual interactions are key mechanisms to generate inherited phenotypic changes. Epigenetic changes can be transmitted between generations independently from changes in DNA sequence. In Caenorhabditis elegans, epigenetic differences, i.e. epimutations, mediated by small non-coding RNAs, particularly 22G-RNAs, as well as chromatin have been identified, and their average persistence is three to five generations. In addition, previous research showed that some epimutations had a longer duration and concerned genes that were enriched for multiple components of xenobiotic response pathways. These results raise the possibility that environmental stresses might change the rate at which epimutations occur, with potential significance for adaptation.

RESULTS

In this work, we explore this question by propagating C. elegans lines either in control conditions or in moderate or high doses of cisplatin, which introduces genotoxic stress by damaging DNA. Our results show that cisplatin has a limited effect on global small non-coding RNA epimutations and epimutations in gene expression levels. However, cisplatin exposure leads to increased fluctuations in the levels of small non-coding RNAs derived from tRNA cleavage. We show that changes in tRNA-derived small RNAs may be associated with gene expression changes.

CONCLUSIONS

Our work shows that epimutations are not substantially altered by cisplatin exposure but identifies transient changes in tRNA-derived small RNAs as a potential source of variation induced by genotoxic stress.

Nucleotide excision repair deficiency is a targetable therapeutic vulnerability in clear cell renal cell carcinoma

Due to a demonstrated lack of DNA repair deficiencies, clear cell renal cell carcinoma (ccRCC) has not benefitted from targeted synthetic lethality-based therapies. We investigated whether nucleotide excision repair (NER) deficiency is present in an identifiable subset of ccRCC cases that would render those tumors sensitive to therapy targeting this specific DNA repair pathway aberration. We used functional assays that detect UV-induced 6-4 pyrimidine-pyrimidone photoproducts to quantify NER deficiency in ccRCC cell lines. We also measured sensitivity to irofulven, an experimental cancer therapeutic agent that specifically targets cells with inactivated transcription-coupled nucleotide excision repair (TC-NER). In order to detect NER deficiency in clinical biopsies, we assessed whole exome sequencing data for the presence of an NER deficiency associated mutational signature previously identified in ERCC2 mutant bladder cancer. Functional assays showed NER deficiency in ccRCC cells. Some cell lines showed irofulven sensitivity at a concentration that is well tolerated by patients. Prostaglandin reductase 1 (PTGR1), which activates irofulven, was also associated with this sensitivity. Next generation sequencing data of the cell lines showed NER deficiency-associated mutational signatures. A significant subset of ccRCC patients had the same signature and high PTGR1 expression. ccRCC cell line-based analysis showed that NER deficiency is likely present in this cancer type. Approximately 10% of ccRCC patients in the TCGA cohort showed mutational signatures consistent with ERCC2 inactivation associated NER deficiency and also substantial levels of PTGR1 expression. These patients may be responsive to irofulven, a previously abandoned anticancer agent that has minimal activity in NER-proficient cells.

Subclonal Somatic Copy-Number Alterations Emerge and Dominate in Recurrent Osteosarcoma

Multiple large-scale genomic profiling efforts have been undertaken in osteosarcoma to define the genomic drivers of tumorigenesis, therapeutic response, and disease recurrence. The spatial and temporal intratumor heterogeneity could also play a role in promoting tumor growth and treatment resistance. We conducted longitudinal whole-genome sequencing of 37 tumor samples from 8 patients with relapsed or refractory osteosarcoma. Each patient had at least one sample from a primary site and a metastatic or relapse site. Subclonal copy-number alterations were identified in all patients except one. In 5 patients, subclones from the primary tumor emerged and dominated at subsequent relapses. MYC gain/amplification was enriched in the treatment-resistant clones in 6 of 7 patients with multiple clones. Amplifications in other potential driver genes, such as CCNE1, RAD21, VEGFA, and IGF1R, were also observed in the resistant copy-number clones. A chromosomal duplication timing analysis revealed that complex genomic rearrangements typically occurred prior to diagnosis, supporting a macroevolutionary model of evolution, where a large number of genomic aberrations are acquired over a short period of time followed by clonal selection, as opposed to ongoing evolution. A mutational signature analysis of recurrent tumors revealed that homologous repair deficiency (HRD)-related SBS3 increases at each time point in patients with recurrent disease, suggesting that HRD continues to be an active mutagenic process after diagnosis. Overall, by examining the clonal relationships between temporally and spatially separated samples from patients with relapsed/refractory osteosarcoma, this study sheds light on the intratumor heterogeneity and potential drivers of treatment resistance in this disease.

SIGNIFICANCE

The chemoresistant population in recurrent osteosarcoma is subclonal at diagnosis, emerges at the time of primary resection due to selective pressure from neoadjuvant chemotherapy, and is characterized by unique oncogenic amplifications.

Chemotherapy-Induced Neoantigen Nanovaccines Enhance Checkpoint Blockade Cancer Immunotherapy

Chemotherapeutics have the potential to increase the efficacy of cancer immunotherapies by stimulating the production of damage-associated molecular patterns (DAMPs) and eliciting mutations that result in the production of neoantigens, thereby increasing the immunogenicity of cancerous lesions. However, the dose-limiting toxicity and limited immunogenicity of chemotherapeutics are not sufficient to induce a robust antitumor response. We hypothesized that cancer cells in vitro treated with ultrahigh doses of various chemotherapeutics artificially increased the abundance, variety, and specificity of DAMPs and neoantigens, thereby improving chemoimmunotherapy. The in vitro chemotherapy-induced (IVCI) nanovaccines manufactured from cell lysates comprised multiple neoantigens and DAMPs, thereby exhibiting comprehensive antigenicity and adjuvanticity. Our IVCI nanovaccines exhibited enhanced immune responses in CT26 tumor-bearing mice, with a significant increase in CD4+/CD8+ T cells in tumors in combination with immune checkpoint inhibitors. The concept of IVCI nanovaccines provides an idea for manufacturing and artificial enhancement of immunogenicity vaccines to improve chemoimmunotherapy.

Precision medicine in the era of multi-omics: can the data tsunami guide rational treatment decision?

Precision medicine for cancer is rapidly moving to an approach that integrates multiple dimensions of the biology in order to model mechanisms of cancer progression in each patient. The discovery of multiple drivers per tumor challenges medical decision that faces several treatment options. Drug sensitivity depends on the actionability of the target, its clonal or subclonal origin and coexisting genomic alterations. Sequencing has revealed a large diversity of drivers emerging at treatment failure, which are potential targets for clinical trials or drug repurposing. To effectively prioritize therapies, it is essential to rank genomic alterations based on their proven actionability. Moving beyond primary drivers, the future of precision medicine necessitates acknowledging the intricate spatial and temporal heterogeneity inherent in cancer. The advent of abundant complex biological data will make artificial intelligence algorithms indispensable for thorough analysis. Here, we will discuss the advancements brought by the use of high-throughput genomics, the advantages and limitations of precision medicine studies and future perspectives in this field.

Reliability of panel-based mutational signatures for immune-checkpoint-inhibition efficacy prediction in non-small cell lung cancer

OBJECTIVES

Mutational signatures (MS) are gaining traction for deriving therapeutic insights for immune checkpoint inhibition (ICI). We asked if MS attributions from comprehensive targeted sequencing assays are reliable enough for predicting ICI efficacy in non-small cell lung cancer (NSCLC).

METHODS

Somatic mutations of m = 126 patients were assayed using panel-based sequencing of 523 cancer-related genes. In silico simulations of MS attributions for various panels were performed on a separate dataset of m = 101 whole genome sequenced patients. Non-synonymous mutations were deconvoluted using COSMIC v3.3 signatures and used to test a previously published machine learning classifier.

RESULTS

The ICI efficacy predictor performed poorly with an accuracy of 0.51-0.09+0.09, average precision of 0.52-0.11+0.11, and an area under the receiver operating characteristic curve of 0.50-0.09+0.10. Theoretical arguments, experimental data, and in silico simulations pointed to false negative rates (FNR) related to panel size. A secondary effect was observed, where deconvolution of small ensembles of point mutations lead to reconstruction errors and misattributions.

CONCLUSION

MS attributions from current targeted panel sequencing are not reliable enough to predict ICI efficacy. We suggest that, for downstream classification tasks in NSCLC, signature attributions be based on whole exome or genome sequencing instead.

Chronic treatment with ATR and CHK1 inhibitors does not substantially increase the mutational burden of human cells

DNA replication stress (RS) entails the frequent slow down and arrest of replication forks by a variety of conditions that hinder accurate and processive genome duplication. Elevated RS leads to genome instability, replication catastrophe and eventually cell death. RS is particularly prevalent in cancer cells and its exacerbation to unsustainable levels by chemotherapeutic agents remains a cornerstone of cancer treatments. The adverse consequences of RS are normally prevented by the ATR and CHK1 checkpoint kinases that stabilize stressed forks, suppress origin firing and promote cell cycle arrest when replication is perturbed. Specific inhibitors of these kinases have been developed and shown to potentiate RS and cell death in multiple in vitro cancer settings. Ongoing clinical trials are now probing their efficacy against various cancer types, either as single agents or in combination with mainstay chemotherapeutics. Despite their promise as valuable additions to the anti-cancer pharmacopoeia, we still lack a genome-wide view of the potential mutagenicity of these new drugs. To investigate this question, we performed chronic long-term treatments of TP53-depleted human cancer cells with ATR and CHK1 inhibitors (ATRi, AZD6738/ceralasertib and CHK1i, MK8776/SCH-900776). ATR or CHK1 inhibition did not significantly increase the mutational burden of cells, nor generate specific mutational signatures. Indeed, no notable changes in the numbers of base substitutions, short insertions/deletions and larger scale rearrangements were observed despite induction of replication-associated DNA breaks during treatments. Interestingly, ATR inhibition did induce a slight increase in closely-spaced mutations, a feature previously attributed to translesion synthesis DNA polymerases. The results suggest that ATRi and CHK1i do not have substantial mutagenic effects in vitro when used as standalone agents.

Evolution of cisplatin resistance through coordinated metabolic reprogramming of the cellular reductive state

BACKGROUND

Cisplatin (CDDP) is a mainstay treatment for advanced head and neck squamous cell carcinomas (HNSCC) despite a high frequency of innate and acquired resistance. We hypothesised that tumours acquire CDDP resistance through an enhanced reductive state dependent on metabolic rewiring.

METHODS

To validate this model and understand how an adaptive metabolic programme might be imprinted, we performed an integrated analysis of CDDP-resistant HNSCC clones from multiple genomic backgrounds by whole-exome sequencing, RNA-seq, mass spectrometry, steady state and flux metabolomics.

RESULTS

Inactivating KEAP1 mutations or reductions in KEAP1 RNA correlated with Nrf2 activation in CDDP-resistant cells, which functionally contributed to resistance. Proteomics identified elevation of downstream Nrf2 targets and the enrichment of enzymes involved in generation of biomass and reducing equivalents, metabolism of glucose, glutathione, NAD(P), and oxoacids. This was accompanied by biochemical and metabolic evidence of an enhanced reductive state dependent on coordinated glucose and glutamine catabolism, associated with reduced energy production and proliferation, despite normal mitochondrial structure and function.

CONCLUSIONS

Our analysis identified coordinated metabolic changes associated with CDDP resistance that may provide new therapeutic avenues through targeting of these convergent pathways.

Using genomic scars to select immunotherapy beneficiaries in advanced non-small cell lung cancer

In advanced non-small cell lung cancer (NSCLC), response to immunotherapy is difficult to predict from pre-treatment information. Given the toxicity of immunotherapy and its financial burden on the healthcare system, we set out to identify patients for whom treatment is effective. To this end, we used mutational signatures from DNA mutations in pre-treatment tissue. Single base substitutions, doublet base substitutions, indels, and copy number alteration signatures were analysed in [Formula: see text] patients (the discovery set). We found that tobacco smoking signature (SBS4) and thiopurine chemotherapy exposure-associated signature (SBS87) were linked to durable benefit. Combining both signatures in a machine learning model separated patients with a progression-free survival hazard ratio of 0.40[Formula: see text] on the cross-validated discovery set and 0.24[Formula: see text] on an independent external validation set ([Formula: see text]). This paper demonstrates that the fingerprints of mutagenesis, codified through mutational signatures, select advanced NSCLC patients who may benefit from immunotherapy, thus potentially reducing unnecessary patient burden.

The c-MYC transcription factor conduces to resistance to cisplatin by regulating MMS19 in bladder cancer cells

Chemoresistance is one of the dominant causes for tumor progression and recurrence of bladder cancer (BC). This paper investigated the effects of transcription factor c-MYC through promoting MMS19 expression on proliferation, metastasis and cisplatin (DDP) resistance in BC cells. The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database were applied to acquire the needed BC gene data. The mRNA and protein levels of c-MYC and MMS19 were verified with q-PCR or Western blot assay. MTT and Transwell assays were utilized to detect cell viability and metastasis. Chromatin Immunoprecipitation (ChIP) assay and Luciferase reporter assay were exerted to confirm the relationship between c-MYC and MMS19. TCGA and GEO BC datasets results implied MMS19 could be an independent indicator for BC patients' prognosis. MMS19 expression was dramatically augmented in BC cell lines. Overexpression of MMS19 conduced to accelerate BC cells proliferation, metastasis and increase DDP resistance. c-MYC was positively correlated with MMS19 and acted as a transcription activator for MMS19 in BC cell lines and activated MMS19 expression. Overexpression of c-MYC facilitated BC cells proliferation, metastasis and DDP resistance. In conclusions, c-MYC gene was a transcriptional regulator of MMS19. Up-regulation of c-MYC facilitated BC cells proliferation, metastasis and DDP resistance by motivating MMS19 expression. This molecular mechanism between c-MYC and MMS19 exerts a crucial mission in BC tumorigenesis and DDP resistance, and may contribute to the diagnosis and therapy of BC for the time to come.

The evolution of non-small cell lung cancer metastases in TRACERx

Metastatic disease is responsible for the majority of cancer-related deaths1. We report the longitudinal evolutionary analysis of 126 non-small cell lung cancer (NSCLC) tumours from 421 prospectively recruited patients in TRACERx who developed metastatic disease, compared with a control cohort of 144 non-metastatic tumours. In 25% of cases, metastases diverged early, before the last clonal sweep in the primary tumour, and early divergence was enriched for patients who were smokers at the time of initial diagnosis. Simulations suggested that early metastatic divergence more frequently occurred at smaller tumour diameters (less than 8 mm). Single-region primary tumour sampling resulted in 83% of late divergence cases being misclassified as early, highlighting the importance of extensive primary tumour sampling. Polyclonal dissemination, which was associated with extrathoracic disease recurrence, was found in 32% of cases. Primary lymph node disease contributed to metastatic relapse in less than 20% of cases, representing a hallmark of metastatic potential rather than a route to subsequent recurrences/disease progression. Metastasis-seeding subclones exhibited subclonal expansions within primary tumours, probably reflecting positive selection. Our findings highlight the importance of selection in metastatic clone evolution within untreated primary tumours, the distinction between monoclonal versus polyclonal seeding in dictating site of recurrence, the limitations of current radiological screening approaches for early diverging tumours and the need to develop strategies to target metastasis-seeding subclones before relapse.

Nucleotide excision repair deficiency is a targetable therapeutic vulnerability in clear cell renal cell carcinoma

Purpose

Due to a demonstrated lack of DNA repair deficiencies, clear cell renal cell carcinoma (ccRCC) has not benefitted from targeted synthetic lethality-based therapies. We investigated whether nucleotide excision repair (NER) deficiency is present in an identifiable subset of ccRCC cases that would render those tumors sensitive to therapy targeting this specific DNA repair pathway aberration.

Experimental Design

We used functional assays that detect UV-induced 6-4 pyrimidine-pyrimidone photoproducts to quantify NER deficiency in ccRCC cell lines. We also measured sensitivity to irofulven, an experimental cancer therapeutic agent that specifically targets cells with inactivated transcription-coupled nucleotide excision repair (TC-NER). In order to detect NER deficiency in clinical biopsies, we assessed whole exome sequencing data for the presence of an NER deficiency associated mutational signature previously identified in ERCC2 mutant bladder cancer.

Results

Functional assays showed NER deficiency in ccRCC cells. Irofulven sensitivity increased in some cell lines. Prostaglandin reductase 1 (PTGR1), which activates irofulven, was also associated with this sensitivity. Next generation sequencing data of the cell lines showed NER deficiency-associated mutational signatures. A significant subset of ccRCC patients had the same signature and high PTGR1 expression.

Conclusions

ccRCC cell line based analysis showed that NER deficiency is likely present in this cancer type. Approximately 10% of ccRCC patients in the TCGA cohort showed mutational signatures consistent with ERCC2 inactivation associated NER deficiency and also substantial levels of PTGR1 expression. These patients may be responsive to irofulven, a previously abandoned anticancer agent that has minimal activity in NER-proficient cells.

The mutational impact of Illudin S on human cells

Illudin S (ILS) is a fungal sesquiterpene secondary metabolite with potent genotoxic and cytotoxic properties. Early genetic studies and more recent genome-wide CRISPR screens showed that Illudin-induced lesions are preferentially repaired by transcription-coupled nucleotide excision repair (TC-NER) with some contribution from post-replication repair pathways. In line with these results, Irofulven, a semi-synthetic ILS analog was recently shown to be particularly effective on cell lines and patient-derived xenografts with impaired NER (e.g. ERCC2/3 mutations), raising hope that ILS-derived molecules may soon enter the clinic. Despite the therapeutic potential of ILS and its analogs, we still lack a global understanding of their mutagenic potential. Here, we characterize the mutational signatures associated with chronic exposure to ILS in human cells. ILS treatment rapidly stalls DNA replication and transcription, leading to the activation of the replication stress response and the accumulation of DNA damage. Novel single and double base substitution signatures as well as a characteristic indel signature indicate that ILS treatment preferentially alkylates purine residues and induces oxidative stress, confirming prior in vitro data. Many mutation contexts exhibit a strong transcriptional strand bias, highlighting the contribution of TC-NER to the repair of ILS lesions. Finally, collateral mutations are also observed in response to ILS, suggesting a contribution of translesion synthesis pathways to ILS tolerance. Accordingly, ILS treatment led to the rapid recruitment of the Y-family DNA polymerase kappa onto chromatin, supporting its preferential use for ILS lesion bypass. Altogether, our work provides the first global assessment of the genomic impact of ILS, demonstrating the contribution of multiple DNA repair pathways to ILS resistance and mutagenicity.

Mechanistic Approach for Protective Effect of ARA290, a Specific Ligand for the Erythropoietin/CD131 Heteroreceptor, against Cisplatin-Induced Nephrotoxicity, the Involvement of Apoptosis and Inflammation Pathways

ARA 290, an 11-amino acid linear nonhematopoietic peptide derived from the three-dimensional structure of helix B of the erythropoietin (EPO), interacts selectively with the innate repair receptor (IRR) that arbitrates tissue protection. The aim of this study was to investigate the protective effects of ARA290 against cisplatin-induced nephrotoxicity. For this purpose, HEK-293 and ACHN cells were treated with ARA290 (50-400 nM) and cisplatin (2.5 μM) in pretreatment condition. Then, cytotoxicity, genotoxicity, oxidative stress parameters (ROS, GPx, SOD, and MDA), and inflammatory markers (TNFα, IL6, and IL1β) were evaluated. Furthermore, apoptotic cell death was assessed via caspase-3 activity and tunnel assay. To determine the molecular mechanisms of the possible nephroprotective effects of ARA290, gene and protein expressions of TNFα, IL1β, IL6, Caspase-3, Bax, and Bcl2 were evaluated by real-time PCR and western blot assay, respectively. The findings indicated that ARA290 significantly reduced the DNA damage parameters of comet assay and the frequency of micronuclei induced by cisplatin. Besides, ARA290 improved cisplatin-induced oxidative stress by reducing MDA/ROS levels and enhancing antioxidant enzyme levels. In addition, reduced levels of pro-inflammatory cytokines indicated that cisplatin-induced renal inflammation was mitigated upon the treatment with ARA290. Besides, ARA290 ameliorates cisplatin-induced cell injury by antagonizing apoptosis. Furthermore, the molecular findings indicated that gene and protein levels of TNFα, IL1β, IL6, Caspase-3, and Bax were significantly decreased and gene and protein levels of Bcl2 significantly increased in the ARA290 plus cisplatin group compared with the cisplatin group. These findings revealed that ARA290 as a potent chemo-preventive agent exerted a protective effect on cisplatin-induced nephrotoxicity mostly through its anti-apoptotic, anti-inflammatory, and antioxidant potentials and also suggested that ARA290 might be a new therapeutic approach for patients with acute kidney injury.

Subclonal somatic copy number alterations emerge and dominate in recurrent osteosarcoma

Multiple large-scale tumor genomic profiling efforts have been undertaken in osteosarcoma, however, little is known about the spatial and temporal intratumor heterogeneity and how it may drive treatment resistance. We performed whole-genome sequencing of 37 tumor samples from eight patients with relapsed or refractory osteosarcoma. Each patient had at least one sample from a primary site and a metastatic or relapse site. We identified subclonal copy number alterations in all but one patient. We observed that in five patients, a subclonal copy number clone from the primary tumor emerged and dominated at subsequent relapses. MYC gain/amplification was enriched in the treatment-resistant clone in 6 out of 7 patients with more than one clone. Amplifications in other potential driver genes, such as CCNE1, RAD21, VEGFA, and IGF1R, were also observed in the resistant copy number clones. Our study sheds light on intratumor heterogeneity and the potential drivers of treatment resistance in osteosarcoma.

Germline TP53 mutations undergo copy number gain years prior to tumor diagnosis

Li-Fraumeni syndrome (LFS) is a hereditary cancer predisposition syndrome associated with germline TP53 pathogenic variants. Here, we perform whole-genome sequence (WGS) analysis of tumors from 22 patients with TP53 germline pathogenic variants. We observe somatic mutations affecting Wnt, PI3K/AKT signaling, epigenetic modifiers and homologous recombination genes as well as mutational signatures associated with prior chemotherapy. We identify near-ubiquitous early loss of heterozygosity of TP53, with gain of the mutant allele. This occurs earlier in these tumors compared to tumors with somatic TP53 mutations, suggesting the timing of this mark may distinguish germline from somatic TP53 mutations. Phylogenetic trees of tumor evolution, reconstructed from bulk and multi-region WGS, reveal that LFS tumors exhibit comparatively limited heterogeneity. Overall, our study delineates early copy number gains of mutant TP53 as a characteristic mutational process in LFS tumorigenesis, likely arising years prior to tumor diagnosis.

Unravelling the instability of mutational signatures extraction via archetypal analysis

The high cosine similarity between some single-base substitution mutational signatures and their characteristic flat profiles could suggest the presence of overfitting and mathematical artefacts. The newest version (v3.3) of the signature database available in the Catalogue Of Somatic Mutations In Cancer (COSMIC) provides a collection of 79 mutational signatures, which has more than doubled with respect to previous version (30 profiles available in COSMIC signatures v2), making more critical the associations between signatures and specific mutagenic processes. This study both provides a systematic assessment of the de novo extraction task through simulation scenarios based on the latest version of the COSMIC signatures and highlights, through a novel approach using archetypal analysis, which COSMIC signatures are redundant and more likely to be considered as mathematical artefacts. 29 archetypes were able to reconstruct the profile of all the COSMIC signatures with cosine similarity > 0.8. Interestingly, these archetypes tend to group similar original signatures sharing either the same aetiology or similar biological processes. We believe that these findings will be useful to encourage the development of new de novo extraction methods avoiding the redundancy of information among the signatures while preserving the biological interpretation.

Elevated Mutational Age in Blood of Children Treated for Cancer Contributes to Therapy-Related Myeloid Neoplasms

Childhood cancer survivors are confronted with various chronic health conditions like therapy-related malignancies. However, it is unclear how exposure to chemotherapy contributes to the mutation burden and clonal composition of healthy tissues early in life. Here, we studied mutation accumulation in hematopoietic stem and progenitor cells (HSPC) before and after cancer treatment of 24 children. Of these children, 19 developed therapy-related myeloid neoplasms (t-MN). Posttreatment HSPCs had an average mutation burden increase comparable to what treatment-naïve cells accumulate during 16 years of life, with excesses up to 80 years. In most children, these additional mutations were induced by clock-like processes, which are also active during healthy aging. Other patients harbored mutations that could be directly attributed to treatments like platinum-based drugs and thiopurines. Using phylogenetic inference, we demonstrate that most t-MN in children originate after the start of treatment and that leukemic clones become dominant during or directly after chemotherapy exposure.

SIGNIFICANCE

Our study shows that chemotherapy increases the mutation burden of normal blood cells in cancer survivors. Only few drugs damage the DNA directly, whereas in most patients, chemotherapy-induced mutations are caused by processes similar to those present during normal aging. This article is highlighted in the In This Issue feature, p. 1825.

The spatio-temporal evolution of multiple myeloma from baseline to relapse-refractory states

Deciphering Multiple Myeloma evolution in the whole bone marrow is key to inform curative strategies. Here, we perform spatial-longitudinal whole-exome sequencing, including 140 samples collected from 24 Multiple Myeloma patients during up to 14 years. Applying imaging-guided sampling we observe three evolutionary patterns, including relapse driven by a single-cell expansion, competing/co-existing sub-clones, and unique sub-clones at distinct locations. While we do not find the unique relapse sub-clone in the baseline focal lesion(s), we show a close phylogenetic relationship between baseline focal lesions and relapse disease, highlighting focal lesions as hotspots of tumor evolution. In patients with ≥3 focal lesions on positron-emission-tomography at diagnosis, relapse is driven by multiple distinct sub-clones, whereas in other patients, a single-cell expansion is typically seen (p < 0.01). Notably, we observe resistant sub-clones that can be hidden over years, suggesting that a prerequisite for curative therapies would be to overcome not only tumor heterogeneity but also dormancy.

The translational challenges of precision oncology

The translational challenges in the field of precision oncology are in part related to the biological complexity and diversity of this disease. Technological advances in genomics have facilitated large sequencing efforts and discoveries that have further supported this notion. In this review, we reflect on the impact of these discoveries on our understanding of several concepts: cancer initiation, cancer prevention, early detection, adjuvant therapy and minimal residual disease monitoring, cancer drug resistance, and cancer evolution in metastasis. We discuss key areas of focus for improving cancer outcomes, from biological insights to clinical application, and suggest where the development of these technologies will lead us. Finally, we discuss practical challenges to the wider adoption of molecular profiling in the clinic and the need for robust translational infrastructure.

Functionalized Lineage Tracing Can Enable the Development of Homogenization-Based Therapeutic Strategies in Cancer

The therapeutic landscape across many cancers has dramatically improved since the introduction of potent targeted agents and immunotherapy. Nonetheless, success of these approaches is too often challenged by the emergence of therapeutic resistance, fueled by intratumoral heterogeneity and the immense evolutionary capacity inherent to cancers. To date, therapeutic strategies have attempted to outpace the evolutionary tempo of cancer but frequently fail, resulting in lack of tumor response and/or relapse. This realization motivates the development of novel therapeutic approaches which constrain evolutionary capacity by reducing the degree of intratumoral heterogeneity prior to treatment. Systematic development of such approaches first requires the ability to comprehensively characterize heterogeneous populations over the course of a perturbation, such as cancer treatment. Within this context, recent advances in functionalized lineage tracing approaches now afford the opportunity to efficiently measure multimodal features of clones within a tumor at single cell resolution, enabling the linkage of these features to clonal fitness over the course of tumor progression and treatment. Collectively, these measurements provide insights into the dynamic and heterogeneous nature of tumors and can thus guide the design of homogenization strategies which aim to funnel heterogeneous cancer cells into known, targetable phenotypic states. We anticipate the development of homogenization therapeutic strategies to better allow for cancer eradication and improved clinical outcomes.

Lost by Transcription: Fork Failures, Elevated Expression, and Clinical Consequences Related to Deletions in Metastatic Colorectal Cancer

Among the structural variants observed in metastatic colorectal cancer (mCRC), deletions (DELs) show a size preference of ~10 kb-1 Mb and are often found in common fragile sites (CFSs). To gain more insight into the biology behind the occurrence of these specific DELs in mCRC, and their possible association with outcome, we here studied them in detail in metastatic lesions of 429 CRC patients using available whole-genome sequencing and corresponding RNA-seq data. Breakpoints of DELs within CFSs are significantly more often located between two consecutive replication origins compared to DELs outside CFSs. DELs are more frequently located at the midpoint of genes inside CFSs with duplications (DUPs) at the flanks of the genes. The median expression of genes inside CFSs was significantly higher than those of similarly-sized genes outside CFSs. Patients with high numbers of these specific DELs showed a shorter progression-free survival time on platinum-containing therapy. Taken together, we propose that the observed DEL/DUP patterns in expressed genes located in CFSs are consistent with a model of transcription-dependent double-fork failure, and, importantly, that the ability to overcome the resulting stalled replication forks decreases sensitivity to platinum-containing treatment, known to induce stalled replication forks as well. Therefore, we propose that our DEL score can be used as predictive biomarker for decreased sensitivity to platinum-containing treatment, which, upon validation, may augment future therapeutic choices.

Genetic and chemotherapeutic influences on germline hypermutation

Mutations in the germline generates all evolutionary genetic variation and is a cause of genetic disease. Parental age is the primary determinant of the number of new germline mutations in an individual's genome1,2. Here we analysed the genome-wide sequences of 21,879 families with rare genetic diseases and identified 12 individuals with a hypermutated genome with between two and seven times more de novo single-nucleotide variants than expected. In most families (9 out of 12), the excess mutations came from the father. Two families had genetic drivers of germline hypermutation, with fathers carrying damaging genetic variation in DNA-repair genes. For five of the families, paternal exposure to chemotherapeutic agents before conception was probably a key driver of hypermutation. Our results suggest that the germline is well protected from mutagenic effects, hypermutation is rare, the number of excess mutations is relatively modest and most individuals with a hypermutated genome will not have a genetic disease.

APOBEC Mutational Signature and Tumor Mutational Burden as Predictors of Clinical Outcomes and Treatment Response in Patients With Advanced Urothelial Cancer

Introduction

Tumor mutational burden (TMB) and APOBEC mutational signatures are potential prognostic markers in patients with advanced urothelial carcinoma (aUC). Their utility in predicting outcomes to specific therapies in aUC warrants additional study.

Methods

We retrospectively reviewed consecutive UC cases assessed with UCSF500, an institutional assay that uses hybrid capture enrichment of target DNA to interrogate 479 common cancer genes. Hypermutated tumors (HM), defined as having TMB ≥10 mutations/Mb, were also assessed for APOBEC mutational signatures, while non-HM (NHM) tumors were not assessed due to low TMB. The logrank test was used to determine if there were differences in overall survival (OS) and progression-free survival (PFS) among patient groups of interest.

Results

Among 75 aUC patients who had UCSF500 testing, 46 patients were evaluable for TMB, of which 19 patients (41%) had HM tumors and the rest had NHM tumors (27 patients). An additional 29 patients had unknown TMB status. Among 19 HM patients, all 16 patients who were evaluable for analysis had APOBEC signatures. HM patients (N=19) were compared with NHM patients (N=27) and had improved OS from diagnosis (125.3 months vs 35.7 months, p=0.06) but inferior OS for patients treated with chemotherapy (7.0 months vs 13.1 months, p=0.04). Patients with APOBEC (N=16) were compared with remaining 56 patients, comprised of 27 NHM patients and 29 patients with unknown TMB, showing APOBEC patients to have improved OS from diagnosis (125.3 months vs 44.5 months, p=0.05) but inferior OS for patients treated with chemotherapy (7.0 months vs 13.1 months, p=0.05). Neither APOBEC nor HM status were associated with response to immunotherapy.

Conclusions

In a large, single-institution aUC cohort assessed with UCSF500, an institutional NGS panel, HM tumors were common and all such tumors that were evaluated for mutational signature analysis had APOBEC signatures. APOBEC signatures and high TMB were prognostic of improved OS from diagnosis and both analyses also predicted inferior outcomes with chemotherapy treatment.

Pancreatic cancer evolution and heterogeneity: integrating omics and clinical data

Pancreatic ductal adenocarcinoma (PDAC), already among the deadliest epithelial malignancies, is rising in both incidence and contribution to overall cancer deaths. Decades of research have improved our understanding of PDAC carcinogenesis, including characterizing germline predisposition, the cell of origin, precursor lesions, the sequence of genetic alterations, including simple and structural alterations, transcriptional changes and subtypes, tumour heterogeneity, metastatic progression and the tumour microenvironment. These fundamental advances inform contemporary translational efforts in primary prevention, screening and early detection, multidisciplinary management and survivorship, as prospective clinical trials begin to adopt molecular-based selection criteria to guide targeted therapies. Genomic and transcriptomic data on PDAC were also included in the international pan-cancer analysis of approximately 2,600 cancers, a milestone in cancer research that allows further insight through comparison with other tumour types. Thus, this is an ideal time to review our current knowledge of PDAC evolution and heterogeneity, gained from the study of preclinical models and patient biospecimens, and to propose a model of PDAC evolution that takes into consideration findings from varied sources, with a particular focus on the genomics of human PDAC.

Characterisation of the spectrum and genetic dependence of collateral mutations induced by translesion DNA synthesis

Translesion DNA synthesis (TLS) is a fundamental damage bypass pathway that utilises specialised polymerases with relaxed template specificity to achieve replication through damaged DNA. Misinsertions by low fidelity TLS polymerases may introduce additional mutations on undamaged DNA near the original lesion site, which we termed collateral mutations. In this study, we used whole genome sequencing datasets of chicken DT40 and several human cell lines to obtain evidence for collateral mutagenesis in higher eukaryotes. We found that cisplatin and UVC radiation frequently induce close mutation pairs within 25 base pairs that consist of an adduct-associated primary and a downstream collateral mutation, and genetically linked their formation to TLS activity involving PCNA ubiquitylation and polymerase κ. PCNA ubiquitylation was also indispensable for close mutation pairs observed amongst spontaneously arising base substitutions in cell lines with disrupted homologous recombination. Collateral mutation pairs were also found in melanoma genomes with evidence of UV exposure. We showed that collateral mutations frequently copy the upstream base, and extracted a base substitution signature that describes collateral mutagenesis in the presented dataset regardless of the primary mutagenic process. Using this mutation signature, we showed that collateral mutagenesis creates approximately 10–20% of non-paired substitutions as well, underscoring the importance of the process.

DNA base substitutions are the most common form of genomic mutations, formed both spontaneously and in response to environmental mutagens. One of the main mechanisms of base substitution mutagenesis is translesion synthesis, a process that relies on specialised DNA polymerases to replicate damaged DNA templates. In addition to incorrect base insertions at the site of lesions in the template, translesion polymerases may also generate ‘collateral’ mutations away from the lesion due to their lower accuracy in selecting the correct incoming nucleotide. In this study, we surveyed the whole genome sequence of experimental cell clones to examine the extent and genetic dependence of collateral mutagenesis in higher eukaryotes. Looking for close mutation pairs, we found that collateral mutations frequently occur near primary lesions generated by cisplatin or ultraviolet radiation in chicken and human cells, but are restricted to a short distance of approximately 25 base pairs. By analysing their sequence context, we showed that collateral mutations can also occur near correctly bypassed primary lesions and may be responsible for a considerable proportion of all base substitution mutations.

Mechanisms of Cisplatin Resistance in HPV Negative Head and Neck Squamous Cell Carcinomas

Head and neck squamous cell carcinomas (HNSCCs) are the eighth most common cancers worldwide. While promising new therapies are emerging, cisplatin-based chemotherapy remains the gold standard for advanced HNSCCs, although most of the patients relapse due to the development of resistance. This review aims to condense the different mechanisms involved in the development of cisplatin resistance in HNSCCs and highlight future perspectives intended to overcome its related complications. Classical resistance mechanisms include drug import and export, DNA repair and oxidative stress control. Emerging research identified the prevalence of these mechanisms in populations of cancer stem cells (CSC), which are the cells mainly contributing to cisplatin resistance. The use of old and new CSC markers has enabled the identification of the characteristics within HNSCC CSCs predisposing them to treatment resistance, such as cell quiescence, increased self-renewal capacity, low reactive oxygen species levels or the acquisition of epithelial to mesenchymal transcriptional programs. In the present review, we will discuss how cell intrinsic and extrinsic cues alter the phenotype of CSCs and how they influence resistance to cisplatin treatment. In addition, we will assess how the stromal composition and the tumor microenvironment affect drug resistance and the acquisition of CSCs' characteristics through a complex interplay between extracellular matrix content as well as immune and non-immune cell characteristics. Finally, we will describe how alterations in epigenetic modifiers or other signaling pathways can alter tumor behavior and cell plasticity to induce chemotherapy resistance. The data generated in recent years open up a wide range of promising strategies to optimize cisplatin therapy, with the potential to personalize HNSCC patient treatment strategies.

Defining genomic events involved in the evolutionary trajectories of myeloma and its precursor conditions

All patients with a diagnosis of multiple myeloma (MM) have a preceding, asymptomatic expansion of clonal plasma cells, clinically recognized as monoclonal gammopathy of undetermined significance or smoldering multiple myeloma (SMM). While most patients with monoclonal gammopathy of undetermined significance have a very small rate of progression, SMM is a widely heterogeneous condition where a fraction of patients will progress to symptomatic MM rather quickly, while others will experience an indolent clinical course. The differentiation between progressive and stable precursor condition thus represents one of the most important unmet clinical needs in the MM community. The ability to identify patients at high-risk of progression before major clonal expansion and onset of end-organ damage would enable strategies for early prevention and perhaps more effective intervention. All proposed criteria to predict the progression of myeloma precursor conditions are built around indirect markers of disease burden and, therefore, are generally able to accurately identify only a small fraction of patients in whom progression to MM is already occurring. Leveraging whole genome and exome sequencing, it has been shown that patients with stable myeloma precursor conditions are characterized by either absence or lower prevalence of distinct genomic events that are detectable in progressive precursor condition years before the progression. In this review, we discuss evolving genomic concepts and tools; and their ability to differentiate myeloma precursor conditions into two distinct entities: one benign (monoclonal gammopathy of benign significance) and another malignant (asymptomatic multiple myeloma).

Therapeutic and prognostic insights from the analysis of cancer mutational signatures

The somatic mutations in each cancer genome are caused by multiple mutational processes, each of which leaves a characteristic imprint (or 'signature'), potentially caused by specific etiologies or exposures. Deconvolution of these signatures offers a glimpse into the evolutionary history of individual tumors. Recent work has shown that mutational signatures may also yield therapeutic and prognostic insights, including the identification of cell-intrinsic signatures as biomarkers of drug response and prognosis. For example, mutational signatures indicating homologous recombination deficiency are associated with poly(ADP)-ribose polymerase (PARP) inhibitor sensitivity, whereas APOBEC-associated signatures are associated with ataxia telangiectasia and Rad3-related kinase (ATR) inhibitor sensitivity. Furthermore, therapy-induced mutational signatures implicated in cancer progression have also been uncovered, including the identification of thiopurine-induced TP53 mutations in leukemia. In this review, we explore the various ways mutational signatures can reveal new therapeutic and prognostic insights, thus extending their traditional role in identifying disease etiology.

BRCA1 deficiency specific base substitution mutagenesis is dependent on translesion synthesis and regulated by 53BP1

Defects in BRCA1, BRCA2 and other genes of the homology-dependent DNA repair (HR) pathway cause an elevated rate of mutagenesis, eliciting specific mutation patterns including COSMIC signature SBS3. Using genome sequencing of knock-out cell lines we show that Y family translesion synthesis (TLS) polymerases contribute to the spontaneous generation of base substitution and short insertion/deletion mutations in BRCA1 deficient cells, and that TLS on DNA adducts is increased in BRCA1 and BRCA2 mutants. The inactivation of 53BP1 in BRCA1 mutant cells markedly reduces TLS-specific mutagenesis, and rescues the deficiency of template switch-mediated gene conversions in the immunoglobulin V locus of BRCA1 mutant chicken DT40 cells. 53BP1 also promotes TLS in human cellular extracts in vitro. Our results show that HR deficiency-specific mutagenesis is largely caused by TLS, and suggest a function for 53BP1 in regulating the choice between TLS and error-free template switching in replicative DNA damage bypass.