MutationalPatterns: the one stop shop for the analysis of mutational processes

Heterogeneity and evolution of DNA mutation rates in microsatellite stable colorectal cancer

Historically, DNA sequence mutability has been considered relatively uniform and low in tumors with chromosomal instability (CIN), based on the assumption that high mutability would be detrimental in karyotypically aberrant contexts. Recent in silico analyses have challenged this view, suggesting some heterogeneity in mutation rates across CIN tumors; however, these predictions lack experimental validation. It also remains unclear how the intertumor variability of mutation rates compares to intratumor diversification and evolves along disease progression, whether mutation rates are functionally relevant in CIN cancers, and which mutational processes shape mutational accrual during CIN tumor onset and evolution. To address these gaps, we performed mutation accumulation experiments using clonal populations of patient-derived tumoroids from seven CIN, microsatellite-stable colorectal cancers (CRCs), and one microsatellite-unstable CRC. Each tumor exhibited a distinctive mutation rate footprint that was conserved among different clones from the same ancestor. In contrast, mutation rates diverged markedly across different tumors, with variations in magnitude within microsatellite-stable tumors as prominent as those distinguishing them from microsatellite-unstable tumors. New mutations reflected mutational processes associated with defective DNA replication and repair, which were not detected in normal tissues. Last, both mutation accumulation assays and high-depth whole-exome sequencing of subclonal variants showed higher mutation rates in metastatic lesions compared with matched primary tumors, suggesting positive selection for cells with increasing mutability during cancer dissemination. By providing an empirical assessment of mutation rates in human cancer, our data delineate heterogeneity, heritability, and progression-associated evolvability of DNA mutational instability as hallmarks of microsatellite-stable CRC.

The evolution and convergence of mutation spectra across mammals

Despite the key role genetic mutations play in shaping phenotypic differences between species, little is currently known about the evolution of germline mutation spectra across mammals. Domesticated species are likely particularly interesting case studies because of their high mutation rates and complex evolutionary histories, which can span multiple founding events and genetic bottlenecks. Here we have developed a new reusable workflow, nSPECTRa, that can undertake the key steps in characterising mutation spectra, from determining ancestral alleles to characterising multiple forms of variation. We apply nSPECTRa to seven species, including several that have undergone domestication, and highlight how nSPECTRa can provide important insights into mutation rate evolution. While mutation spectra most often show marked differences between species and even breeds, certain mutation types have risen to a high frequency in subpopulations of different species, indicative of convergent evolution in mutation rates. This includes the previously characterized TCC- > TTC change enriched among European humans, which is also enriched among East Asian cattle. We show Indicine cattle are particularly interesting examples of how different mutation spectra segregate within a population and subsequently spread across the globe. Together, this work has important implications for understanding the control and evolution of mammalian mutation rates.

Unraveling mutagenic processes influencing the tumor mutational patterns of individuals with constitutional mismatch repair deficiency

Constitutional mismatch repair deficiency (CMMRD), caused by bi-allelic germline variants in mismatch repair (MMR) genes, is associated with high cancer incidence early in life. A better understanding of mutational processes driving sequential CMMRD tumors can advance optimal treatment. Here, we describe a genomic characterization on a representative collection of CMMRD-associated tumors consisting of 41 tumors from 17 individuals. Mutational patterns in these tumors appear to be influenced by multiple factors, including the affected MMR gene and tumor type. Somatic polymerase proofreading mutations, commonly present in brain tumors, are also found in a T-cell lymphoblastic lymphoma displaying associated mutational patterns. We show prominent mutational patterns in two second primary hematological malignancies after temozolomide treatment. Furthermore, an indel signature, characterized by one-base pair cytosine insertions in cytosine homopolymers, is found in 54% of tumors. In conclusion, analysis of sequential CMMRD tumors reveals diverse mutational patterns influenced by the affected MMR gene, tumor type and treatment history.

Mercury Bioaccumulation in Female Breast Cancer Is Associated to CXCR4 Expression

The growing incidence of breast cancer over time suggests that environmental factors might contribute to the underlying causes of the disease. Mercury, a toxic metal classified as a Substance of Very High Concern, accumulates in the body through contaminated food, air, water, and soil, raising concerns about its role in tumor biology. The main aim of this study was to identify the possible associations between in situ mercury bioaccumulation and the molecular features of breast cancer. To achieve this, a total of 26 breast cancer cases were analyzed using an integrated approach that combined DNA and RNA sequencing, histological analysis, and inductively coupled plasma mass spectrometry (ICP-MS) to assess mercury bioaccumulation. Mercury was detected in 72% of the cases. A significant positive correlation was found between mercury bioaccumulation and CXCR4 expression in breast cancer tissues. Bioinformatic analysis further revealed that CXCR4 expression was significantly higher in metastatic tissues compared to primary tumors. These findings suggest that mercury accumulation may influence tumor biology through the CXCR4-CXCL12 signaling pathway, highlighting a potential mechanism by which mercury contributes to breast cancer progression.

Androgens mediate sexual dimorphism in Pilarowski-Bjornsson Syndrome

Sex-specific penetrance in autosomal dominant Mendelian conditions is largely understudied. The neurodevelopmental disorder Pilarowski-Bjornsson syndrome (PILBOS) was initially described in females. Here, we describe the clinical and genetic characteristics of the largest PILBOS cohort to date, showing that both sexes can exhibit PILBOS features, although males are overrepresented. A mouse model carrying a human-derived Chd1 missense variant (Chd1 R616Q/+) displays female-restricted phenotypes, including growth deficiency, anxiety and hypotonia. Orchiectomy unmasks a growth deficiency phenotype in male Chd1 R616Q/+ mice, while testosterone rescues the phenotype in females, implicating androgens in phenotype modulation. In the gnomAD and UK Biobank databases, rare missense variants in CHD1 are overrepresented in males, supporting a male protective effect. We identify 33 additional highly constrained autosomal genes with missense variant overrepresentation in males. Our results support androgen-regulated sexual dimorphism in PILBOS and open novel avenues to understand the mechanistic basis of sexual dimorphism in other autosomal Mendelian disorders.

Phase II Trial of Atezolizumab and Bevacizumab for Treatment of HPV-Positive Unresectable or Metastatic Squamous Cell Carcinoma of the Anal Canal

PURPOSE

Anti-PD-L1 antibodies are associated with responses in <25% of patients with metastatic human papillomavirus-associated malignancies. VEGF signaling causes immune evasion and immune suppression within the tumor. We evaluated the anti-PD-L1 antibody atezolizumab and anti-VEGF antibody bevacizumab for patients with unresectable, advanced anal cancer.

PATIENTS AND METHODS

For this phase II study, participants with previously treated, immunotherapy-naïve anal cancer received atezolizumab (1,200 mg) and bevacizumab (15 mg/kg) intravenously every 21 days. Responses were evaluated every 9 weeks (RECIST version 1.1). The primary endpoint was the best radiographic response. Median survival was estimated by Kaplan-Meier and compared for selected biomarkers (including paired pre- and on-treatment biopsies) using a log-rank test.

RESULTS

Among 20 participants, the overall response rate was 11% [95% confidence interval (CI): 1.2-32]. Median progression-free survival and overall survival were 4.1 months (95% CI, 2.6-not assessable) and 11.6 months (95% CI, 9.5-20), respectively. One grade 5 bevacizumab-related bowel perforation occurred. Analyses of 16 paired biopsies linked increases in IFN-γ (P = 0.03) and inflammatory response (P = 0.02) gene expression signatures with prolonged progression-free survival, as did increases in CD3+CD8+PD1+ (P = 0.02) cells and decreases in CD3+FoxP3+ cells (P = 0.04) from 10 paired biopsies with multiplex immunofluorescence. A subgroup of anal cancers characterized by the SBS31 "prior-platinum" signature demonstrated shorter median overall survival (HR, 6.3; 95% CI, 1.2-32; P = 0.01).

CONCLUSIONS

Atezolizumab and bevacizumab demonstrate activity similar to anti-PD-1 antibodies alone for unresectable anal cancer. Our translational data identify undescribed chromosomal and transcriptomic biomarkers associated with survival for metastatic anal cancer. These correlative findings warrant confirmation and further validation in larger, prospective immunotherapy trials for advanced anal cancer.

RNA sequence analysis of somatic mutations in aging and Parkinson's Disease

Parkinson's Disease (PD) is an age-related neurodegenerative disorder that has been associated with increased DNA damage. To test if PD is associated with increased somatic mutations, we analyzed RNA-seq data in whole blood from 5 visits of the Parkinson's Progression Markers Initiative for clonally amplified somatic variants. Comprehensive analysis of RNA-sequencing data revealed a total of 5,927 somatic variants (2.4 variants per sample on average). Mutation frequencies were significantly elevated in PD subjects as compared to age-matched controls at the time of the last visit. This was confirmed by RNA analysis of substantia nigra. By contrast, the fraction of carriers with clonal hematopoiesis, was significantly reduced in old PD patients as compared to old healthy controls. These results indicate that while the overall mutation rate is higher in PD, specific clonally amplified mutations are protective against PD, as has been found for Alzheimer's Disease.

Neuronal somatic mutations are increased in multiple sclerosis lesions

Neuroinflammation underpins neurodegeneration and clinical progression in multiple sclerosis (MS), but knowledge of processes linking these disease mechanisms remains incomplete. Here we investigated somatic single-nucleotide variants (sSNVs) in the genomes of 106 single neurons from post-mortem brain tissue of ten MS cases and 16 controls to determine whether somatic mutagenesis is involved. We observed an increase of 43.9 sSNVs per year in neurons from chronic MS lesions, a 2.5 times faster rate than in neurons from normal-appearing MS and control tissues. This difference was equivalent to 1,291 excess sSNVs in lesion neurons at 70 years of age compared to controls. We performed mutational signature analysis to investigate mechanisms underlying neuronal sSNVs and identified a signature characteristic of lesions with a strong, age-associated contribution to sSNV counts. This research suggests that neuroinflammation is mutagenic in the MS brain, potentially contributing to disease progression.

A novel mouse model of upper tract urothelial carcinoma highlights the impact of dietary intervention on gut microbiota and carcinogenesis prevention despite carcinogen exposure

Animal models of N-butyl-N-(4-hydroxy butyl) nitrosamine (BBN)-induced urothelial carcinoma (UC), particularly bladder cancer (BC), have long been established. However, the rare incidence of BBN-induced upper urinary tract UC (UTUC), which originates from the same urothelium as BC, remains elusive. The scarcity of animal models of UTUC has made it challenging to study the biology of UTUC. To address this problem, we tried to establish a novel mouse model of UTUC by treating multiple mice strains and sexes with BBN. The molecular consistency between the UTUC mouse model and human UTUC was confirmed using multi-omics analyses, including whole-exome, whole-transcriptome, and spatial transcriptome sequencing. 16S ribosomal RNA metagenome sequencing, metabolome analysis, and dietary interventions were employed to assess changes in the gut microbiome, metabolome, and carcinogenesis of UTUC. Of all treated mice, only female BALB/c mice developed UTUC over BC. Multi-omics analyses confirmed that the UTUC model reflected the molecular characteristics and heterogeneity of human UTUC with poor prognosis. Furthermore, the model exhibited increased Tnf-related inflammatory gene expression in the upper urinary tract and a low relative abundance of Parabacteroides distasonis in the gut. Dietary intervention, mainly without alanine, led to P. distasonis upregulation and successfully prevented UTUC, as well as suppressed Tnf-related inflammatory gene expression in the upper urinary tract despite the exposure to BBN. This is the first report to demonstrate a higher incidence of UTUC than BC in a non-engineered mouse model using BBN. Overall, this model could serve as a useful tool for comprehensively investigating UTUC in future studies.

Double mutant DNMT3A AML: a unique subtype experiencing increased DNA damage and poor prognosis

ABSTRACT

Mutation of DNMT3A, encoding a de novo methyltransferase essential for cytosine methylation, is a common early event in clonal hematopoiesis (CH) and adult acute myeloid leukemia (AML). Spontaneous deamination of methylated cytosines damages DNA, which is repaired by the base excision repair (BER) enzymes methyl-CpG binding domain 4 (MBD4) and thymine DNA glycosylase (TDG). Congenital MBD4 deficiency has been linked to early-onset CH and AML and is marked by exceedingly high levels of DNA damage and mutation of DNMT3A. Strikingly, wild-type (WT) DNMT3A binds TDG, thereby potentiating its repair activity. Because TDG is the only remaining BER enzyme in MBD4-deficient patients with AML capable of repairing methylation damage, we investigated whether mutant DNMT3A negatively affects the repair function of TDG. We found that, although WT DNMT3A stimulates TDG function, mutant DNMT3A impairs TDG-mediated repair of DNA damage in vitro. In light of this finding and to extrapolate our observations to the broader AML patient population, we investigate here the genetic profiles and survival outcomes of patients with AML with single mutant (SM) vs double mutant (DM) DNMT3A. Patients with DM DNMT3A AML show a characteristic driver mutation landscape and reduced overall survival compared with patients with SM DNMT3A AML. Importantly, whole-genome sequencing showed a trend for increased DNA damage in primary DM DNMT3A AML samples, especially when DNMT3A mutations are located at the DNMT3A-TDG interaction interface.

Protocol for genome-wide analysis of somatic variants at single-cell resolution using primary template-directed DNA amplification

The study of somatic mutations in single cells provides insights into aging and carcinogenesis, which is complicated by the dependency on whole-genome amplification (WGA). Here, we describe a detailed workflow starting from single-cell isolation to WGA by primary template-directed amplification (PTA), sequencing, quality control, and downstream analyses. A machine learning approach, the PTA Analysis Toolkit (PTATO), is used to filter the hundreds to thousands of artificial variants induced by WGA from true mutations at high sensitivity and accuracy. For complete details on the use and execution of this protocol, please refer to Middelkamp et al.1.

Defining the genome-wide mutagenic impact of APOBEC3 enzymes

Somatic mutations drive cancer initiation and tumor evolution. Therefore, the etiology of mutagenesis in cancer is important to preventative and treatment strategies. Somatic mutagenesis in cancer is a multifactorial process and includes both endogenous and exogenous sources of mutations. One recently recognized source of mutagenesis in cancer is the innate immune APOBEC3 family of enzymes, which catalyze cytosine deamination to restrict viral infection but can aberrantly act on the cellular genome, resulting in mutations. Single base substitution (SBS) signatures, or mutational patterns, identified in cancer genomes have demonstrated widespread mutagenesis caused by APOBEC3 enzymes throughout human tumors. To comprehensively define the consequences of APOBEC3 mutagenesis, we developed an experimental pipeline for prospective analysis of genome-wide mutations caused by APOBEC3 activity. This pipeline can be adapted to analyze additional sources of mutagenesis across a spectrum of cells.

Glioma mutational signatures associated with haloalkane exposure are enriched in firefighters

BACKGROUND

Glioma is the most common malignant primary brain tumor and is associated with significant morbidity and mortality. Modifiable risk factors remain unidentified. New advances in exposure assessment, genomic analyses, and statistical techniques permit more accurate evaluation of glioma risk associated with exogenous occupational or environmental exposures.

METHODS

By using whole-exome sequencing data from matched germline and glioma tumor samples, the authors compared tumor mutational signatures for 17 persons with glioma and a documented occupational history of firefighting with those of 18 persons with glioma without an occupational history of firefighting. All 35 individuals were participants in the University of California, San Francisco Adult Glioma Study.

RESULTS

There was a positive correlation among firefighters between the median number of sample variants attributable to single-base substitution signature 42, a single-base substitution mutational signature associated with haloalkane exposure (from the Catalogue of Somatic Mutational Signatures in Cancer) and firefighting years (p = .04; R2 = 0.29). Among nonfirefighters, the individuals with the highest number of median variants attributable to single-base substitution signature 42 also had occupations that possibly exposed them to haloalkanes, such as painting and being a mechanic.

CONCLUSIONS

In summary, the authors identified gliomas that had mutational signatures associated with haloalkane exposure that were enriched in firefighters and other occupations.

APOBEC3A drives ovarian cancer metastasis by altering epithelial-mesenchymal transition

High-grade serous ovarian cancer (HGSOC) is the most prevalent and aggressive histological subtype of ovarian cancer and often presents with metastatic disease. The drivers of metastasis in HGSOC remain enigmatic. APOBEC3A (A3A), an enzyme that generates mutations across various cancers, has been proposed as a mediator of tumor heterogeneity and disease progression. However, the role of A3A in HGSOC has not been explored. We observed an association between high levels of APOBEC3-mediated mutagenesis and poor overall survival in primary HGSOC. We experimentally addressed this correlation by modeling A3A expression in HGSOC, and this resulted in increased metastatic behavior of HGSOC cells in culture and distant metastatic spread in vivo, which was dependent on catalytic activity of A3A. A3A activity in both primary and cultured HGSOC cells yielded consistent alterations in expression of epithelial-mesenchymal transition (EMT) genes resulting in hybrid EMT and mesenchymal signatures, providing a mechanism for their increased metastatic potential. Inhibition of key EMT factors TWIST1 and IL-6 resulted in mitigation of A3A-dependent metastatic phenotypes. Our findings define the prevalence of A3A mutagenesis in HGSOC and implicate A3A as a driver of HGSOC metastasis via EMT, underscoring its clinical relevance as a potential prognostic biomarker. Our study lays the groundwork for the development of targeted therapies aimed at mitigating the deleterious effect of A3A-driven EMT in HGSOC.

UVA-light-induced mutagenesis in the exome of human nucleotide excision repair-deficient cells

Skin cancer is associated with genetic mutations caused by sunlight exposure, primarily through ultraviolet (UV) radiation that damages DNA. While UVA is less energetic, it is the predominant solar UV component reaching the Earth's surface. However, the mechanisms of UVA-induced mutagenesis and its role in skin cancer development remain poorly understood. This study employed whole exome sequencing of clones from human XP-C cells, which lack nucleotide excision repair (NER), to characterize somatic mutations induced by UVA exposure. DNA sequence analysis of UVA-irradiated XP-C cells revealed a marked increase in mutation frequency across nearly all types of base substitutions, with particular enrichment in C > T transitions within the CCN and TCN trinucleotide context-potential sites for pyrimidine dimer formation. The C > T mutation primarily occurred at the 3' base of the 5'TC dimer, and an enrichment of CC > TT tandem mutations. We also identified the SBS7b COSMIC mutational signature within irradiated cells, which has been associated with tumors in sun-exposed skin. C > A transversions, often linked to oxidized guanine, were the second most frequently induced mutation, although a specific context for this base substitution was not identified. Moreover, C > T mutations were significantly increased in unirradiated XP-C compared to NER-proficient cells, which may be caused by unrepaired spontaneous DNA damage. Thus, this study indicates that pyrimidine dimers are the primary lesions contributing to UVA-induced mutagenesis in NER-deficient human cells and demonstrates that UVA generates mutational signatures similar to those of UVB irradiation.

Tracking the Evolution of Cutaneous Melanoma by Multiparameter Flow Sorting and Genomic Profiling

Intratumoral heterogeneity and clonal evolution are pivotal in the progression and metastasis of melanoma. However, when combined with variable tumor cellularity, intratumoral heterogeneity limits the sensitivity and accuracy of uncovering a cancer's clonal evolution. In this study, we combined fluorescence-activated cell sorting (FACS) with whole-exome sequencing (WES) to investigate the clonal composition and evolutionary patterns in seven melanoma biopsies obtained from three patients, each having both primary site and metastatic samples. We employed a multiparameter ploidy sorting approach to isolate tumor populations based on DNA ploidy and melanoma biomarkers (SOX10 or S100), enabling us to investigate clonal evolution with high resolution. Our approach increased the mean tumor purity from 70% (range 19-88%) in unsorted material to 91% (range 87-96%) post-sorting. Our findings revealed significant inter- and intratumor heterogeneity, with one patient exhibiting two genomically distinct clonal tumor populations within a single primary site biopsy, each giving rise to different metastases. Our findings highlight the critical role of intratumoral heterogeneity and clonal evolution in melanoma, especially when analyzing tumor trajectories. The unique combination of multiparameter FACS and WES provides a powerful method for identifying clonal populations and reconstructing clonal evolution. This study provides valuable insights into the clonal architecture of melanoma and lays the groundwork for future research with larger patient groups.

Somatic Genomic and Transcriptomic Changes in Single Ischemic Human Heart Cardiomyocytes

Heart failure is a multifaceted syndrome contributing significantly to mortality and hospitalization rates among the global population1. One of the prevalent causes of heart failure is ischemic heart disease (IHD), often caused by a blockage in a coronary artery, ultimately leading to the loss of myocardial tissue and contractile force2. The impact of this ischemic ambiance on the cardiomyocyte genome and transcriptome has not been thoroughly studied. During normal aging, cardiomyocytes progressively accumulate somatic mutations faster than many dividing cells, suggesting that internal and external factors specific to cardiomyocytes might influence this accumulation3. In this study, we analyzed single-cell whole-genome and transcriptome data from the left ventricle of 5 individuals with IHD and 10 healthy control individuals. We found that somatic DNA alterations significantly increase in IHD cardiomyocytes, with distinct mutational patterns indicating a disrupted DNA repair system and a cytotoxic environment, potentially associated with increased inflammatory response in the myocardium and a compensatory anti-inflammatory response in IHD. An in vitro iPS-derived hypoxic cardiomyocyte mutational profile indicates similar mutational spectra. Transcriptomic analysis revealed increased expression of EGR1, FOS, and collagen genes in ischemic heart cardiomyocytes, leading to a more fibrotic heart. The aberrant accumulation of DNA alterations and changes in transcriptional patterns in the ischemic heart cardiomyocytes provide insight into the development of IHD.

Tobacco smoke carcinogens exacerbate APOBEC mutagenesis and carcinogenesis

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

The Distinct Roles of NEIL1 and XPA in Limiting Aflatoxin B1-Induced Mutagenesis in Mice

Dietary exposure to aflatoxin B1 (AFB1) is a risk factor for the development of hepatocellular carcinomas. Following metabolic activation, AFB1 reacts with guanines to form covalent DNA adducts, which induce high-frequency G > T transversions. The molecular signature associated with these mutational events aligns with the single-base substitution signature 24 (SBS24) in the Catalog of Somatic Mutations in Cancer database. Deficiencies in either base excision repair due to the absence of Nei-like DNA glycosylase 1 (NEIL1) or nucleotide excision repair due to the absence of xeroderma complementation group A protein (XPA) contribute to hepatocellular carcinomas in murine models. In the current study, ultra-low error duplex sequencing was used to characterize mutational profiles in liver DNAs of NEIL1-deficient, XPA-deficient, and DNA repair-proficient mice following neonatal injection of 1 mg/kg AFB1. Analyses of AFB1-induced mutations showed high cosine similarity to SBS24 regardless of repair proficiency status. The absence of NEIL1 resulted in an approximately 30% increase in the frequency of mutations, with the distribution suggesting preferential NEIL1-dependent repair of AFB1 lesions in open chromatin regions. A trend of increased mutagenesis was also observed in the absence of XPA. Consistent with the role of XPA in transcription-coupled repair, mutational profiles in XPA-deficient mice showed disruption of the transcriptional bias in mutations associated with SBS24. Implications: Our findings define the roles of DNA repair pathways in AFB1-induced mutagenesis and carcinogenesis in murine models, with these findings having implications in human health for those with base excision repair and nucleotide excision repair deficiencies.

Mutational signature analysis predicts bacterial hypermutation and multidrug resistance

Bacteria of clinical importance, such as Pseudomonas aeruginosa, can become hypermutators upon loss of DNA mismatch repair (MMR) and are clinically correlated with high rates of multidrug resistance (MDR). Here, we demonstrate that hypermutated MMR-deficient P. aeruginosa has a unique mutational signature and rapidly acquires MDR upon repeated exposure to first-line or last-resort antibiotics. MDR acquisition was irrespective of drug class and instead arose through common resistance mechanisms shared between the initial and secondary drugs. Rational combinations of drugs having distinct resistance mechanisms prevented MDR acquisition in hypermutated MMR-deficient P. aeruginosa. Mutational signature analysis of P. aeruginosa across different human disease contexts identified appreciable quantities of MMR-deficient clinical isolates that were already MDR or prone to future MDR acquisition. Mutational signature analysis of patient samples is a promising diagnostic tool that may predict MDR and guide precision-based medical care.

Unveiling the molecular profile of a prostate carcinoma: implications for personalized medicine

BACKGROUND

Prostate cancer is the most common diagnosed tumor and the fifth cancer related death among men in Europe. Although several genetic alterations such as ERG-TMPRSS2 fusion, MYC amplification, PTEN deletion and mutations in p53 and BRCA2 genes play a key role in the pathogenesis of prostate cancer, specific gene alteration signature that could distinguish indolent from aggressive prostate cancer or may aid in patient stratification for prognosis and/or clinical management of patients with prostate cancer is still missing. Therefore, here, by a multi-omics approach we describe a prostate cancer carrying the fusion of TMPRSS2 with ERG gene and deletion of 16q chromosome arm.

RESULTS

We have observed deletion of KDM6A gene, which may represent an additional genomic alteration to be considered for patient stratification. The cancer hallmarks gene signatures highlight intriguing molecular aspects that characterize the biology of this tumor by both a high hypoxia and immune infiltration scores. Moreover, our analysis showed a slight increase in the Tumoral Mutational Burden, as well as an over-expression of the immune checkpoints. The omics profiling integrating hypoxia, ROS and the anti-cancer immune response, optimizes therapeutic strategies and advances personalized care for prostate cancer patients.

CONCLUSION

The here data reported can lay the foundation for predicting a poor prognosis for the studied prostate cancer, as well as the possibility of targeted therapies based on the modulation of hypoxia, ROS, and the anti-cancer immune response.

Contributing factors to the oxidation-induced mutational landscape in human cells

8-oxoguanine (8-oxoG) is a common oxidative DNA lesion that causes G > T substitutions. Determinants of local and regional differences in 8-oxoG-induced mutability across genomes are currently unknown. Here, we show DNA oxidation induces G > T substitutions and insertion/deletion (INDEL) mutations in human cells and cancers. Potassium bromate (KBrO3)-induced 8-oxoGs occur with similar sequence preferences as their derived substitutions, indicating that the reactivity of specific oxidants dictates mutation sequence specificity. While 8-oxoG occurs uniformly across chromatin, 8-oxoG-induced mutations are elevated in compact genomic regions, within nucleosomes, and at inward facing guanines within strongly positioned nucleosomes. Cryo-electron microscopy structures of OGG1-nucleosome complexes indicate that these effects originate from OGG1's ability to flip outward positioned 8-oxoG lesions into the catalytic pocket while inward facing lesions are occluded by the histone octamer. Mutation spectra from human cells with DNA repair deficiencies reveals contributions of a DNA repair network limiting 8-oxoG mutagenesis, where OGG1- and MUTYH-mediated base excision repair is supplemented by the replication-associated factors Pol η and HMCES. Transcriptional asymmetry of KBrO3-induced mutations in OGG1- and Pol η-deficient cells also demonstrates transcription-coupled repair can prevent 8-oxoG-induced mutation. Thus, oxidant chemistry, chromatin structures, and DNA repair processes combine to dictate the oxidative mutational landscape in human genomes.

Integrated immunogenomic analyses of high-grade serous ovarian cancer reveal vulnerability to combination immunotherapy

Background

The efficacy of immunotherapies in high-grade serous ovarian cancer (HGSOC) is limited, but clinical trials investigating the potential of combination immunotherapy including poly-ADP-ribose polymerase inhibitors (PARPis) are ongoing. Homologous recombination repair deficiency or BRCAness and the composition of the tumor microenvironment appear to play a critical role in determining the therapeutic response.

Methods

We conducted comprehensive immunogenomic analyses of HGSOC using data from several patient cohorts. Machine learning methods were used to develop a classification model for BRCAness from gene expression data. Integrated analysis of bulk and single-cell RNA sequencing data was used to delineate the tumor immune microenvironment and was validated by immunohistochemistry. The impact of PARPi and BRCA1 mutations on the activation of immune-related pathways was studied using ovarian cancer cell lines, RNA sequencing, and immunofluorescence analysis.

Results

We identified a 24-gene signature that predicts BRCAness. Comprehensive immunogenomic analyses across patient cohorts identified samples with BRCAness and high immune infiltration. Further characterization of these samples revealed increased infiltration of immunosuppressive cells, including tumor-associated macrophages expressing TREM2, C1QA, and LILRB4, as specified by single-cell RNA sequencing data and gene expression analysis of samples from patients receiving combination therapy with PARPi and anti-PD-1. Our findings show also that genomic instability and PARPi activated the cGAS-STING signaling pathway in vitro and the downstream innate immune response in a similar manner to HGSOC patients with BRCAness status. Finally, we have developed a web application (https://ovrseq.icbi.at) and an associated R package OvRSeq, which allow for comprehensive characterization of ovarian cancer patient samples and assessment of a vulnerability score that enables stratification of patients to predict response to the combination immunotherapy.

Conclusions

Genomic instability in HGSOC affects the tumor immune environment, and TAMs play a crucial role in modulating the immune response. Based on various datasets, we have developed a diagnostic application that uses RNA sequencing data not only to comprehensively characterize HGSOC but also to predict vulnerability and response to combination immunotherapy.

Benchmarking 13 tools for mutational signature attribution, including a new and improved algorithm

Mutational signatures are characteristic patterns of mutations caused by endogenous mutational processes or by exogenous mutational exposures. There has been little benchmarking of approaches for determining which signatures are present in a sample and estimating the number of mutations due to each signature. This problem is referred to as "signature attribution." We show that there are often many combinations of signatures that can reconstruct the patterns of mutations in a sample reasonably well, even after encouraging sparse solutions. We benchmarked 13 approaches to signature attribution, including a new approach called Presence Attribute Signature Activity (PASA), on large synthetic data sets (2700 synthetic samples in total). These data sets recapitulated the single-base, insertion-deletion, and doublet-base mutational signature repertoires of nine cancer types. For single-base substitution mutations, PASA and MuSiCal outperformed other approaches on all the cancer types combined. However, the ranking of approaches varied by cancer type. For doublet-base substitutions and small insertions and deletions, while PASA outperformed the other approaches in most of the nine cancer types, the ranking of approaches again varied by cancer type. We believe that this variation reflects inherent difficulties in signature attribution. These difficulties stem from the fact that there are often many attributions that can reasonably explain the pattern of mutations in a sample and from the combinatorial search space due to the need to impose sparsity. Tables herein can provide guidance on the selection of mutational signature attribution approaches that are best suited to particular cancer types and study objectives.

Molecular and Clinicopathologic Impact of GNAS Variants Across Solid Tumors

PURPOSE

The molecular drivers underlying mucinous tumor pathogenicity are poorly understood. GNAS mutations predict metastatic burden and treatment resistance in mucinous appendiceal adenocarcinoma. We investigated the pan-cancer clinicopathologic relevance of GNAS variants.

METHODS

We assessed 58,043 patients with Memorial Sloan Kettering-Integrated Mutation Profiling of Actionable Cancer Targets (IMPACT)-sequenced solid tumors to identify oncogenic variants, including GNAS, associated with mucinous tumor phenotype. We then performed comprehensive molecular analyses to compare GNAS-mutant (mut) and wild-type tumors across cancers. Gene expression patterns associated with GNAS-mut tumors were assessed in a The Cancer Genome Atlas cohort. Associations between GNAS variant status and peritoneal metastasis, first-line systemic therapy response, progression-free survival (PFS), and overall survival (OS) were determined using a propensity-matched subcohort of patients with metastatic disease.

RESULTS

Mucinous tumors were enriched for oncogenic GNAS variants. GNAS was mutated in >1% of small bowel, cervical, colorectal, pancreatic, esophagogastric, hepatobiliary, and GI neuroendocrine cancers. Across these cancers, GNAS-mut tumors exhibited a generally conserved C-to-T mutation-high, aneuploidy-low molecular profile with co-occurring prevalent KRAS variants (65% of GNAS-mut tumors) and fewer TP53 alterations. GNAS-mut tumors exhibited recurrently comutated alternative tumor suppressors (RBM10, INPPL1) and upregulation of MAPK and cell surface modulators. GNAS-mut tumors demonstrate an increased prevalence of peritoneal metastases (odds ratio [OR], 1.7 [95% CI, 1.1 to 2.5]; P = .006), worse response to first-line systemic therapy (OR, 2.2 [95% CI, 1.3 to 3.8]; P = .003), and shorter PFS (median, 5.6 v 7.0 months; P = .047). In a multivariable analysis, GNAS mutated status was independently prognostic of worse OS (hazard ratio, 1.25 [95% CI, 1.01 to 1.56]; adjusted P = .04).

CONCLUSION

Across the assessed cancers, GNAS-mut tumors exhibit a conserved molecular and clinical phenotype defined by mucinous tumor status, increased peritoneal metastasis, poor response to first-line systemic therapy, and worse survival.

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.

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.

APOBEC3A drives metastasis of high-grade serous ovarian cancer by altering epithelial-to-mesenchymal transition

High-grade serous ovarian cancer (HGSOC) is the most prevalent and aggressive histological subtype of ovarian cancer, and often presents with metastatic disease. The drivers of metastasis in HGSOC remain enigmatic. APOBEC3A (A3A), an enzyme that generates mutations across various cancers, has been proposed as a mediator of tumor heterogeneity and disease progression. However, the role of A3A in HGSOC has not been explored. Through analysis of genome sequencing from primary HGSOC, we observed an association between high levels of APOBEC3 mutagenesis and poor overall survival. We experimentally addressed this correlation by modeling A3A activity in HGSOC cell lines and mouse models which resulted in increased metastatic behavior of HGSOC cells in culture and distant metastatic spread in vivo . A3A activity in both primary and cultured HGSOC cells yielded consistent alterations in expression of epithelial-mesenchymal-transition (EMT) genes resulting in hybrid EMT and mesenchymal signatures, and providing a mechanism for their increased metastatic potential. Our findings define the prevalence of A3A mutagenesis in HGSOC and implicate A3A as a driver of HGSOC metastasis via EMT, underscoring its clinical relevance as a potential prognostic biomarker. Our study lays the groundwork for the development of targeted therapies aimed at mitigating the deleterious impact of A3A-driven EMT in HGSOC.

Impact of cancer therapy on clonal hematopoiesis mutations and subsequent clinical outcomes

ABSTRACT

Exposure to cancer therapies is associated with an increased risk of clonal hematopoiesis (CH). The objective of our study was to investigate the genesis and evolution of CH after cancer therapy. In this prospective study, we undertook error-corrected duplex DNA sequencing in blood samples collected before and at 2 time points after chemoradiation in patients with esophageal or lung cancer recruited from 2013 to 2018. We applied a customized workflow to identify the earliest changes in CH mutation count and clone size and determine their association with clinical outcomes. Our study included 29 patients (87 samples). Their median age was 67 years, and 76% (n = 22) were male; the median follow-up period was 3.9 years. The most mutated genes were DNMT3A, TET2, TP53, and ASXL1. We observed a twofold increase in the number of mutations from before to after treatment in TP53, which differed from all other genes examined (P < .001). Among mutations detected before and after treatment, we observed an increased clone size in 38% and a decreased clone size in 5% of TP53 mutations (odds ratio, 3.7; 95% confidence interval [CI], 1.75-7.84; P < .001). Changes in mutation count and clone size were not observed in other genes. Individuals with an increase in the number of TP53 mutations after chemoradiation experienced shorter overall survival (hazard ratio, 7.07; 95% CI, 1.50-33.46; P = .014). In summary, we found an increase in the number and size of TP53 CH clones after chemoradiation that were associated with adverse clinical outcomes.

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.

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.

Comprehensive cross cancer analyses reveal mutational signature cancer specificity

Mutational signatures refer to distinct patterns of DNA mutations that occur in a specific context or under certain conditions. It is a powerful tool to describe cancer etiology. We conducted a study to show cancer heterogeneity and cancer specificity from the aspect of mutational signatures through collinearity analysis and machine learning techniques. Through thorough training and independent validation, our results show that while the majority of the mutational signatures are distinct, similarities between certain mutational signature pairs can be observed through both mutation patterns and mutational signature abundance. The observation can potentially assist to determine the etiology of yet elusive mutational signatures. Further analysis using machine learning approaches demonstrated moderate mutational signature cancer specificity. Skin cancer among all cancer types demonstrated the strongest mutational signature specificity.

A comprehensive molecular characterization of a claudin-low luminal B breast tumor

Breast cancer is the most common cause of death from cancer in women. Here, we present the case of a 43-year-old woman, who received a diagnosis of claudin-low luminal B breast cancer. The lesion revealed to be a poorly differentiated high-grade infiltrating ductal carcinoma, which was strongly estrogen receptor (ER)/progesterone receptor (PR) positive and human epidermal growth factor receptor (HER2) negative. Her tumor underwent in-depth chromosomal, mutational and gene expression analyses. We found a pathogenic protein truncating mutation in the TP53 gene, which is predicted to disrupt its transcriptional activity. The patient also harbors germline mutations in some mismatch repair (MMR) genes, and her tumor displays the presence of immune infiltrates, high tumor mutational burden (TMB) status and the apolipoprotein B mRNA editing enzyme catalytic polypeptide 3 (APOBEC3) associated signatures, which, overall, are predictive for the use of immunotherapy. Here, we propose promising prognostic indicators as well as potential therapeutic strategies based on the molecular characterization of the tumor.

Systematic review and feasibility study on pre-analytical factors and genomic analyses on archival formalin-fixed paraffin-embedded breast cancer tissue

Formalin-fixed paraffin-embedded (FFPE) tissue represents a valuable source for translational cancer research. However, the widespread application of various downstream methods remains challenging. Here, we aimed to assess the feasibility of a genomic and gene expression analysis workflow using FFPE breast cancer (BC) tissue. We conducted a systematic literature review for the assessment of concordance between FFPE and fresh-frozen matched tissue samples derived from patients with BC for DNA and RNA downstream applications. The analytical performance of three different nucleic acid extraction kits on FFPE BC clinical samples was compared. We also applied a newly developed targeted DNA Next-Generation Sequencing (NGS) 370-gene panel and the nCounter BC360® platform on simultaneously extracted DNA and RNA, respectively, using FFPE tissue from a phase II clinical trial. Of the 3701 initial search results, 40 articles were included in the systematic review. High degree of concordance was observed in various downstream application platforms. Moreover, the performance of simultaneous DNA/RNA extraction kit was demonstrated with targeted DNA NGS and gene expression profiling. Exclusion of variants below 5% variant allele frequency was essential to overcome FFPE-induced artefacts. Targeted genomic analyses were feasible in simultaneously extracted DNA/RNA from FFPE material, providing insights for their implementation in clinical trials/cohorts.

The SMC5/6 complex prevents genotoxicity upon APOBEC3A-mediated replication stress

Mutational patterns caused by APOBEC3 cytidine deaminase activity are evident throughout human cancer genomes. In particular, the APOBEC3A family member is a potent genotoxin that causes substantial DNA damage in experimental systems and human tumors. However, the mechanisms that ensure genome stability in cells with active APOBEC3A are unknown. Through an unbiased genome-wide screen, we define the Structural Maintenance of Chromosomes 5/6 (SMC5/6) complex as essential for cell viability when APOBEC3A is active. We observe an absence of APOBEC3A mutagenesis in human tumors with SMC5/6 dysfunction, consistent with synthetic lethality. Cancer cells depleted of SMC5/6 incur substantial genome damage from APOBEC3A activity during DNA replication. Further, APOBEC3A activity results in replication tract lengthening which is dependent on PrimPol, consistent with re-initiation of DNA synthesis downstream of APOBEC3A-induced lesions. Loss of SMC5/6 abrogates elongated replication tracts and increases DNA breaks upon APOBEC3A activity. Our findings indicate that replication fork lengthening reflects a DNA damage response to APOBEC3A activity that promotes genome stability in an SMC5/6-dependent manner. Therefore, SMC5/6 presents a potential therapeutic vulnerability in tumors with active APOBEC3A.

Comprehensive investigation of the mutagenic potential of six pesticides classified by IARC as probably carcinogenic to humans

Pesticides are significant environmental pollutants, and many of them possess mutagenic potential, which is closely linked to carcinogenesis. Here we tested the mutagenicity of all six pesticides classified probably carcinogenic (Group 2A) by the International Agency of Research on Cancer: 4,4'-DDT, captafol, dieldrin, diazinon, glyphosate and malathion. Whole genome sequencing of TK6 human lymphoblastoid cell clones following 30-day exposure at subtoxic concentrations revealed a clear mutagenic effect of treatment with captafol or malathion when added at 200 nM or 100 μM initial concentrations, respectively. Each pesticide induced a specific base substitution mutational signature: captafol increased C to A mutations primarily, while malathion induced mostly C to T mutations. 4,4'-DDT, dieldrin, diazinon and glyphosate were not mutagenic. Whereas captafol induced chromosomal instability, H2A.X phosphorylation and cell cycle arrest in G2/M phase, all indicating DNA damage, malathion did not induce DNA damage markers or cell cycle alterations despite its mutagenic effect. Hypersensitivity of REV1 and XPA mutant DT40 chicken cell lines suggests that captafol induces DNA adducts that are bypassed by translesion DNA synthesis and are targets for nucleotide excision repair. The experimentally identified mutational signatures of captafol and malathion could shed light on the mechanism of action of these compounds. The signatures are potentially suitable for detecting past exposure in tumour samples, but the reanalysis of large cancer genome databases did not reveal any evidence of captafol or malathion exposure.

Failure of immune checkpoint inhibitors for microsatellite instability-positive pancreatic adenocarcinoma with atypical pattern of short tandem repeat mutation

Microsatellite instability (MSI) is an important biomarker in cancer. While routine methods can detect MSI in certain tumor types, in other tumor types the results may be incorrect due to differences in the MSI loci pattern. Here, we report the case of a patient with pancreatic adenocarcinoma, with confirmed MSI by two independent next-generation sequencing tests, but not by routine methods, who had progression on pembrolizumab. Comparison of the patient's MSI loci patterns with MSI+ colorectal adenocarcinoma samples showed a lower fraction of unstable loci, low resolution of a second peak in the repeat length spectrum of unstable short tandem repeats in the patient's sample, and a lower length of indels (3.7 vs 4.5 base pairs, p < 0.01).

Proteogenomic insights into the biology and treatment of pan-melanoma

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

Research on molecular characteristics of ADME-related genes in kidney renal clear cell carcinoma

Genes involved in drug absorption, distribution, metabolism, and excretion (ADME) are named ADME genes. However, the comprehensive role of ADME genes in kidney renal clear cell carcinoma (KIRC) remains unclear. Using the clinical and gene expression data of KIRC patients downloaded from The Cancer Genome Atlas (TCGA), ArrayExpress, and the Gene Expression Omnibus (GEO) databases, we cluster patients into two patterns, and the population with a relatively poor prognosis demonstrated higher level of immunosuppressive cell infiltration and higher proportion of glycolytic subtypes. Then, 17 ADME genes combination identified through the least absolute shrinkage and selection operator algorithm (LASSO, 1000 times) was utilized to calculate the ADME score. The ADME score was found to be an independent predictor of prognosis in KIRC and to be tightly associated with the infiltration level of immune cells, metabolic properties, tumor-related signaling pathways, genetic variation, and responses to chemotherapeutics. Our work revealed the characteristics of ADME in KIRC. Assessing the ADME profiles of individual patients can deepen our comprehension of tumor microenvironment (TME) features in KIRC and can aid in developing more personalized and effective therapeutic strategies.

Non-homologous end joining shapes the genomic rearrangement landscape of chromothripsis from mitotic errors

Mitotic errors generate micronuclei entrapping mis-segregated chromosomes, which are susceptible to catastrophic fragmentation through chromothripsis. The reassembly of fragmented chromosomes by error-prone DNA double-strand break (DSB) repair generates diverse genomic rearrangements associated with human diseases. How specific repair pathways recognize and process these lesions remains poorly understood. Here we use CRISPR/Cas9 to systematically inactivate distinct DSB repair pathways and interrogate the rearrangement landscape of fragmented chromosomes. Deletion of canonical non-homologous end joining (NHEJ) components substantially reduces complex rearrangements and shifts the rearrangement landscape toward simple alterations without the characteristic patterns of chromothripsis. Following reincorporation into the nucleus, fragmented chromosomes localize within sub-nuclear micronuclei bodies (MN bodies) and undergo ligation by NHEJ within a single cell cycle. In the absence of NHEJ, chromosome fragments are rarely engaged by alternative end-joining or recombination-based mechanisms, resulting in delayed repair kinetics, persistent 53BP1-labeled MN bodies, and cell cycle arrest. Thus, we provide evidence supporting NHEJ as the exclusive DSB repair pathway generating complex rearrangements from mitotic errors.

Cancer mutational signatures identification in clinical assays using neural embedding-based representations

While mutational signatures provide a plethora of prognostic and therapeutic insights, their application in clinical-setting, targeted gene panels is extremely limited. We develop a mutational representation model (which learns and embeds specific mutation signature connections) that enables prediction of dominant signatures with only a few mutations. We predict the dominant signatures across more than 60,000 tumors with gene panels, delineating their landscape across different cancers. Dominant signature predictions in gene panels are of clinical importance. These included UV, tobacco, and apolipoprotein B mRNA editing enzyme, catalytic polypeptide (APOBEC) signatures that are associated with better survival, independently from mutational burden. Further analyses reveal gene and mutation associations with signatures, such as SBS5 with TP53 and APOBEC with FGFR3S249C. In a clinical use case, APOBEC signature is a robust and specific predictor for resistance to epidermal growth factor receptor-tyrosine kinase inhibitors (EGFR-TKIs). Our model provides an easy-to-use way to detect signatures in clinical setting assays with many possible clinical implications for an unprecedented number of cancer patients.

Rare Drivers at Low Prevalence with High Cancer Effects in T-Cell and B-Cell Pediatric Acute Lymphoblastic Leukemia

The genomic analyses of pediatric acute lymphoblastic leukemia (ALL) subtypes, particularly T-cell and B-cell lineages, have been pivotal in identifying potential therapeutic targets. Typical genomic analyses have directed attention toward the most commonly mutated genes. However, assessing the contribution of mutations to cancer phenotypes is crucial. Therefore, we estimated the cancer effects (scaled selection coefficients) for somatic substitutions in T-cell and B-cell cohorts, revealing key insights into mutation contributions. Cancer effects for well-known, frequently mutated genes like NRAS and KRAS in B-ALL were high, which underscores their importance as therapeutic targets. However, less frequently mutated genes IL7R, XBP1, and TOX also demonstrated high cancer effects, suggesting pivotal roles in the development of leukemia when present. In T-ALL, KRAS and NRAS are less frequently mutated than in B-ALL. However, their cancer effects when present are high in both subtypes. Mutations in PIK3R1 and RPL10 were not at high prevalence, yet exhibited some of the highest cancer effects in individual T-cell ALL patients. Even CDKN2A, with a low prevalence and relatively modest cancer effect, is potentially highly relevant for the epistatic effects that its mutated form exerts on other mutations. Prioritizing investigation into these moderately frequent but potentially high-impact targets not only presents novel personalized therapeutic opportunities but also enhances the understanding of disease mechanisms and advances precision therapeutics for pediatric ALL.

Mutation Rate and Effective Population Size of the Model Cooperative Bacterium Myxococcus xanthus

Intrinsic rates of genetic mutation have diverged greatly across taxa and exhibit statistical associations with several other parameters and features. These include effective population size (Ne), genome size, and gametic multicellularity, with the latter being associated with both increased mutation rates and decreased effective population sizes. However, data sufficient to test for possible relationships between microbial multicellularity and mutation rate (µ) are lacking. Here, we report estimates of two key population-genetic parameters, Ne and µ, for Myxococcus xanthus, a bacterial model organism for the study of aggregative multicellular development, predation, and social swarming. To estimate µ, we conducted an ∼400-day mutation accumulation experiment with 46 lineages subjected to regular single colony bottlenecks prior to clonal regrowth. Upon conclusion, we sequenced one clonal-isolate genome per lineage. Given collective evolution for 85,323 generations across all lines, we calculate a per base-pair mutation rate of ∼5.5 × 10-10 per site per generation, one of the highest mutation rates among free-living eubacteria. Given our estimate of µ, we derived Ne at ∼107 from neutral diversity at four-fold degenerate sites across two dozen M. xanthus natural isolates. This estimate is below average for eubacteria and strengthens an already clear negative correlation between µ and Ne in prokaryotes. The higher and lower than average mutation rate and Ne for M. xanthus, respectively, amplify the question of whether any features of its multicellular life cycle-such as group-size reduction during fruiting-body development-or its highly structured spatial distribution have significantly influenced how these parameters have evolved.

Integrated analyses of the genetic and clinicopathological features of cholangiolocarcinoma: cholangiolocarcinoma may be characterized by mismatch-repair deficiency

Cholangiolocarcinoma (CLC) is a primary liver carcinoma that resembles the canals of Hering and that has been reported to be associated with stem cell features. Due to its rarity, the nature of CLC remains unclear, and its pathological classification remains controversial. To clarify the positioning of CLC in primary liver cancers and identify characteristics that could distinguish CLC from other liver cancers, we performed integrated analyses using whole-exome sequencing (WES), immunohistochemistry, and a retrospective review of clinical information on eight CLC cases and two cases of recurrent CLC. WES demonstrated that CLC includes IDH1 and BAP1 mutations, which are characteristic of intrahepatic cholangiocarcinoma (iCCA). A mutational signature analysis showed a pattern similar to that of iCCA, which was different from that of hepatocellular carcinoma (HCC). CLC cells, including CK7, CK19, and EpCAM, were positive for cholangiocytic differentiation markers. However, the hepatocytic differentiation marker AFP and stem cell marker SALL4 were completely negative. The immunostaining patterns of CLC with CD56 and epithelial membrane antigen were similar to those of the noncancerous bile ductules. In contrast, mutational signature cluster analyses revealed that CLC formed a cluster associated with mismatch-repair deficiency (dMMR), which was separate from iCCA. Therefore, to evaluate MMR status, we performed immunostaining of four MMR proteins (PMS2, MSH6, MLH1, and MSH2) and detected dMMR in almost all CLCs. In conclusion, CLC had highly similar characteristics to iCCA but not to HCC. CLC can be categorized as a subtype of iCCA. In contrast, CLC has characteristics of dMMR tumors that are not found in iCCA, suggesting that it should be treated distinctly from iCCA. © 2024 The Pathological Society of Great Britain and Ireland.

Molecular profiling of a bladder cancer with very high tumour mutational burden

The increasing incidence of urothelial bladder cancer is a notable global concern, as evidenced by the epidemiological data in terms of frequency, distribution, as well as mortality rates. Although numerous molecular alterations have been linked to the occurrence and progression of bladder cancer, currently there is a limited knowledge on the molecular signature able of accurately predicting clinical outcomes. In this report, we present a case of a pT3b high-grade infiltrating urothelial carcinoma with areas of squamous differentiation characterized by very high tumor mutational burden (TMB), with up-regulations of immune checkpoints. The high TMB, along with elevated expressions of PD-L1, PD-L2, and PD1, underscores the rationale for developing a personalized immunotherapy focused on the use of immune-checkpoint inhibitors. Additionally, molecular analysis revealed somatic mutations in several other cancer-related genes, including TP53, TP63 and NOTCH3. Mutations of TP53 and TP63 genes provide mechanistic insights on the molecular mechanisms underlying disease development and progression. Notably, the above-mentioned mutations and the elevated hypoxia score make the targeting of p53 and/or hypoxia related pathways a plausible personalized medicine option for this bladder cancer, particularly in combination with immunotherapy. Our data suggest a requirement for molecular profiling in bladder cancer to possibly select appropriate immune-checkpoint therapy.

A pan-cancer analysis of the microbiome in metastatic cancer

Microbial communities are resident to multiple niches of the human body and are important modulators of the host immune system and responses to anticancer therapies. Recent studies have shown that complex microbial communities are present within primary tumors. To investigate the presence and relevance of the microbiome in metastases, we integrated mapping and assembly-based metagenomics, genomics, transcriptomics, and clinical data of 4,160 metastatic tumor biopsies. We identified organ-specific tropisms of microbes, enrichments of anaerobic bacteria in hypoxic tumors, associations between microbial diversity and tumor-infiltrating neutrophils, and the association of Fusobacterium with resistance to immune checkpoint blockade (ICB) in lung cancer. Furthermore, longitudinal tumor sampling revealed temporal evolution of the microbial communities and identified bacteria depleted upon ICB. Together, we generated a pan-cancer resource of the metastatic tumor microbiome that may contribute to advancing treatment strategies.

Changes in chromatin accessibility and transcriptional landscape induced by HDAC inhibitors in TP53 mutated patient-derived colon cancer organoids

Here we present the generation and characterization of patient-derived organoids (PDOs) from colorectal cancer patients. PDOs derived from two patients with TP53 mutations were tested with two different HDAC inhibitors (SAHA and NKL54). Cell death induction, transcriptome, and chromatin accessibility changes were analyzed. HDACIs promote the upregulation of low expressed genes and the downregulation of highly expressed genes. A similar differential effect is observed at the level of chromatin accessibility. Only SAHA is a potent inducer of cell death, which is characterized by the upregulation of BH3-only genes BIK and BMF. Up-regulation of BIK is associated with increased accessibility in an intronic region that has enhancer properties. SAHA, but not NKL54, also causes downregulation of BCL2L1 and decreases chromatin accessibility in three distinct regions of the BCL2L1 locus. Both inhibitors upregulate the expression of innate immunity genes and members of the MHC family. In summary, our exploratory study indicates a mechanism of action for SAHA and demonstrate the low efficacy of NKL54 as a single agent for apoptosis induction, using two PDOs. These observations need to be validated in a larger cohort of PDOs.

The SMC5/6 complex prevents genotoxicity upon APOBEC3A-mediated replication stress

Mutational patterns caused by APOBEC3 cytidine deaminase activity are evident throughout human cancer genomes. In particular, the APOBEC3A family member is a potent genotoxin that causes substantial DNA damage in experimental systems and human tumors. However, the mechanisms that ensure genome stability in cells with active APOBEC3A are unknown. Through an unbiased genome-wide screen, we define the Structural Maintenance of Chromosomes 5/6 (SMC5/6) complex as essential for cell viability when APOBEC3A is active. We observe an absence of APOBEC3A mutagenesis in human tumors with SMC5/6 dysfunction, consistent with synthetic lethality. Cancer cells depleted of SMC5/6 incur substantial genome damage from APOBEC3A activity during DNA replication. Further, APOBEC3A activity results in replication tract lengthening which is dependent on PrimPol, consistent with re-initiation of DNA synthesis downstream of APOBEC3A-induced lesions. Loss of SMC5/6 abrogates elongated replication tracts and increases DNA breaks upon APOBEC3A activity. Our findings indicate that replication fork lengthening reflects a DNA damage response to APOBEC3A activity that promotes genome stability in an SMC5/6-dependent manner. Therefore, SMC5/6 presents a potential therapeutic vulnerability in tumors with active APOBEC3A.

Genetically driven predisposition leads to an unusually genomic unstable renal cell carcinoma

Renal cell carcinoma originates from the lining of the proximal convoluted renal tubule and represents the most common type of kidney cancer. Risk factors and comorbidities might be associated to renal cell carcinoma, while a small fraction of 2-3% emerges from patients with predisposing cancer syndromes, typically associated to hereditary mutations in VHL, folliculin, fumarate hydratase or MET genes. Here, we report a case of renal cell carcinoma in patient with concurrent germline mutations in BRCA1 and RAD51 genes. This case displays an unusual high mutational burden and chromosomal aberrations compared to the typical profile of renal cell carcinoma. Mutational analysis on whole genome sequencing revealed an enrichment of the MMR2 mutational signature, which is indicative of impaired DNA repair capacity. Overall, the tumor displayed a profile of unusual high genomic instability which suggests a possible origin from germline predisposing mutations in the DNA repair genes BRCA1 and RAD51. While BRCA1 and RAD51 germline mutations are well-characterised in breast and ovarian cancer, their role in renal cell carcinoma is still largely unexplored. The genomic instability detected in this case of renal cell carcinoma, along with the presence of unusual mutations, might offer support to clinicians for the development of patient-tailored therapies.

Distinguishing preferences of human APOBEC3A and APOBEC3B for cytosines in hairpin loops, and reflection of these preferences in APOBEC-signature cancer genome mutations

The APOBEC3 enzymes convert cytosines in single-stranded DNA to uracils to protect against viruses and retrotransposons but can contribute to mutations that diversify tumors. To understand the mechanism of mutagenesis, we map the uracils resulting from expression of APOBEC3B or its catalytic carboxy-terminal domain (CTD) in Escherichia coli. Like APOBEC3A, the uracilomes of A3B and A3B-CTD show a preference to deaminate cytosines near transcription start sites and the lagging-strand replication templates and in hairpin loops. Both biochemical activities of the enzymes and genomic uracil distribution show that A3A prefers 3 nt loops the best, while A3B prefers 4 nt loops. Reanalysis of hairpin loop mutations in human tumors finds intrinsic characteristics of both the enzymes, with a much stronger contribution from A3A. We apply Hairpin Signatures 1 and 2, which define A3A and A3B preferences respectively and are orthogonal to published methods, to evaluate their contribution to human tumor mutations.