The repertoire of mutational signatures in human cancer

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.

Arginine deprivation enriches lung cancer proteomes with cysteine by inducing arginine-to-cysteine substitutants

Many types of human cancers suppress the expression of argininosuccinate synthase 1 (ASS1), a rate-limiting enzyme for arginine production. Although dependency on exogenous arginine can be harnessed by arginine-deprivation therapies, the impact of ASS1 suppression on the quality of the tumor proteome is unknown. We therefore interrogated proteomes of cancer patients for arginine codon reassignments (substitutants) and surprisingly identified a strong enrichment for cysteine (R>C) in lung tumors specifically. Most R>C events did not coincide with genetically encoded R>C mutations but were likely products of tRNA misalignments. The expression of R>C substitutants was highly associated with oncogenic kelch-like epichlorohydrin (ECH)-associated protein 1 (KEAP1)-pathway mutations and suppressed by intact-KEAP1 in KEAP1-mutated cancer cells. Finally, functional interrogation indicated a key role for R>C substitutants in cell survival to cisplatin, suggesting that regulatory codon reassignments endow cancer cells with more resilience to stress. Thus, we present a mechanism for enriching lung cancer proteomes with cysteines that may affect therapeutic decisions.

Multi-omic analysis identifies hypoalbuminemia as independent biomarker of poor outcome upon PD-1 blockade in metastatic melanoma

We evaluated the prognostic value of hypoalbuminemia in context of various biomarkers at baseline, including clinical, genomic, transcriptomic, and blood-based markers, in patients with metastatic melanoma treated with anti-PD-1 monotherapy or anti-PD-1/anti-CTLA-4 combination therapy (n = 178). An independent validation cohort (n = 79) was used to validate the performance of hypoalbuminemia compared to serum LDH (lactate dehydrogenase) levels. Pre-treatment hypoalbuminemia emerged as the strongest predictor of poor outcome for both OS (HR = 4.01, 95% CI 2.10-7.67, Cox P = 2.63e-05) and PFS (HR = 3.72, 95% CI 2.06-6.73, Cox P = 1.38e-05) in univariate analysis. In multivariate analysis, the association of hypoalbuminemia with PFS was independent of serum LDH, IFN-γ signature expression, TMB, age, ECOG PS, treatment line, treatment type (combination or monotherapy), brain and liver metastasis (HR = 2.76, 95% CI 1.24-6.13, Cox P = 0.0131). Our validation cohort confirmed the prognostic power of hypoalbuminemia for OS (HR = 1.98, 95% CI 1.16-3.38; Cox P = 0.0127) and was complementary to serum LDH in analyses for both OS (LDH-adjusted HR = 2.12, 95% CI 1.2-3.72, Cox P = 0.00925) and PFS (LDH-adjusted HR = 1.91, 95% CI 1.08-3.38, Cox P = 0.0261). In conclusion, pretreatment hypoalbuminemia was a powerful predictor of outcome in ICI in melanoma and showed remarkable complementarity to previously established biomarkers, including high LDH.

ANOS1 accelerates the progression of esophageal cancer identified by multi-omic approaches

To assess the role of ANOS1 in esophageal cancer (ESCA) progression, multi-omic analysis and experimental validation were employed. It was revealed that ANOS1 expression is significantly enhanced in ESCA patients and cell lines. The expression level of ANOS1 in ESCA patients can distinguish the malignancy from normal tissue with an area under curve (AUC) >0.75. Moreover, increased expression of ANOS1 is associated with advanced T stage and worse disease-free survival of ESCA patients. Therefore, a clinically applicable nomogram with ANOS1 was established with strong predictive power. Furthermore, high expression of ANOS1 in ESCA is correlated with (i) the enrichment of epithelial-mesenchymal transition by gene set enrichment analysis, (ii) the involvement in hypoxia, angiogenesis, WNT signaling pathway, and TGFβ signaling pathway by gene set variation analysis, (iii) the presence of the small insertion and deletion mutational signature ID9, associated with chromothripsis, in the single-nucleotide polymorphism analysis, (iv) the amplification of 11q13.3 in the copy number variants analysis, (v) the enrichment of cancer-associated fibroblasts and mesenchymal stromal cells in the tumor microenvironment. All the results from multi-omic analysis indicate that ANOS1 plays a pivotal role in accelerating the progression of ESCA. Results from in vivo and in vitro experiments show that the knockdown of ANOS1 hampers the proliferation of ESCA cells, further validating the oncogenic role of ANOS1 in ESCA. Additionally, potential chemotherapeutics with sensitivity were identified in the high-ANOS1 group. In conclusion, ANOS1 accelerates the progression of ESCA.

Ataxia-Telangiectasia Mutated Loss-of-Function Displays Variant and Tissue-Specific Differences across Tumor Types

PURPOSE

Mutations in the ATM gene are common in multiple cancers, but clinical studies of therapies targeting ATM-aberrant cancers have yielded mixed results. Refinement of ATM loss of function (LOF) as a predictive biomarker of response is urgently needed.

EXPERIMENTAL DESIGN

We present the first disclosure and preclinical development of a novel, selective ATR inhibitor, ART0380, and test its antitumor activity in multiple preclinical cancer models. To refine ATM LOF as a predictive biomarker, we performed a comprehensive pan-cancer analysis of ATM variants in patient tumors and then assessed the ATM variant-to-protein relationship. Finally, we assessed a novel ATM LOF biomarker approach in retrospective clinical data sets of patients treated with platinum-based chemotherapy or ATR inhibition.

RESULTS

ART0380 had potent, selective antitumor activity in a range of preclinical cancer models with differing degrees of ATM LOF. Pan-cancer analysis identified 10,609 ATM variants in 8,587 patient tumors. Cancer lineage-specific differences were seen in the prevalence of deleterious (Tier 1) versus unknown/benign (Tier 2) variants, selective pressure for loss of heterozygosity, and concordance between a deleterious variant and ATM loss of protein (LOP). A novel ATM LOF biomarker approach that accounts for variant classification, relationship to ATM LOP, and tissue-specific penetrance significantly enriched for patients who benefited from platinum-based chemotherapy or ATR inhibition.

CONCLUSIONS

These data help to better define ATM LOF across tumor types in order to optimize patient selection and improve molecularly targeted therapeutic approaches for patients with ATM LOF cancers.

DNA lesion bypass and the stochastic dynamics of transcription-coupled repair

DNA base damage is a major source of oncogenic mutations and disruption to gene expression. The stalling of RNA polymerase II (RNAP) at sites of DNA damage and the subsequent triggering of repair processes have major roles in shaping the genome-wide distribution of mutations, clearing barriers to transcription, and minimizing the production of miscoded gene products. Despite its importance for genetic integrity, key mechanistic features of this transcription-coupled repair (TCR) process are controversial or unknown. Here, we exploited a well-powered in vivo mammalian model system to explore the mechanistic properties and parameters of TCR for alkylation damage at fine spatial resolution and with discrimination of the damaged DNA strand. For rigorous interpretation, a generalizable mathematical model of DNA damage and TCR was developed. Fitting experimental data to the model and simulation revealed that RNA polymerases frequently bypass lesions without triggering repair, indicating that small alkylation adducts are unlikely to be an efficient barrier to gene expression. Following a burst of damage, the efficiency of transcription-coupled repair gradually decays through gene bodies with implications for the occurrence and accurate inference of driver mutations in cancer. The reinitation of transcription from the repair site is not a general feature of transcription-coupled repair, and the observed data is consistent with reinitiation never taking place. Collectively, these results reveal how the directional but stochastic activity of TCR shapes the distribution of mutations following DNA damage.

G-Quadruplex Forming DNA Sequence Context Is Enriched around Points of Somatic Mutations in a Subset of Multiple Myeloma Patients

Multiple myeloma (MM) is the second most common hematological malignancy, which remains incurable despite recent advances in treatment strategies. Like other forms of cancer, MM is characterized by genomic instability, caused by defects in DNA repair. Along with mutations in DNA repair genes and genotoxic drugs used to treat MM, non-canonical secondary DNA structures (four-stranded G-quadruplex structures) can affect accumulation of somatic mutations and chromosomal abnormalities in the tumor cells of MM patients. Here, we tested the hypothesis that G-quadruplex structures may influence the distribution of somatic mutations in the tumor cells of MM patients. We sequenced exomes of normal and tumor cells of 11 MM patients and analyzed the data for the presence of G4 context around points of somatic mutations. To identify molecular mechanisms that could affect mutational profile of tumors, we also analyzed mutational signatures in tumor cells as well as germline mutations for the presence of specific SNPs in DNA repair genes or in genes regulating G-quadruplex unwinding. In several patients, we found that sites of somatic mutations are frequently located in regions with G4 context. This pattern correlated with specific germline variants found in these patients. We discuss the possible implications of these variants for mutation accumulation and specificity in MM and propose that the extent of G4 context enrichment around somatic mutation sites may be a novel metric characterizing mutational processes in tumors.

NRDE2 deficiency impairs homologous recombination repair and sensitizes hepatocellular carcinoma to PARP inhibitors

To identify novel susceptibility genes for hepatocellular carcinoma (HCC), we performed a rare-variant association study in Chinese populations consisting of 2,750 cases and 4,153 controls. We identified four HCC-associated genes, including NRDE2, RANBP17, RTEL1, and STEAP3. Using NRDE2 (index rs199890497 [p.N377I], p = 1.19 × 10-9) as an exemplary candidate, we demonstrated that it promotes homologous recombination (HR) repair and suppresses HCC. Mechanistically, NRDE2 binds to the subunits of casein kinase 2 (CK2) and facilitates the assembly and activity of the CK2 holoenzyme. This NRDE2-mediated enhancement of CK2 activity increases the phosphorylation of MDC1 and then facilitates the HR repair. These functions are eliminated almost completely by the NRDE2-p.N377I variant, which sensitizes the HCC cells to poly(ADP-ribose) polymerase (PARP) inhibitors, especially when combined with chemotherapy. Collectively, our findings highlight the relevance of the rare variants to genetic susceptibility to HCC, which would be helpful for the precise treatment of this malignancy.

Mutation rate heterogeneity at the sub-gene scale due to local DNA hypomethylation

Local mutation rates in human are highly heterogeneous, with known variability at the scale of megabase-sized chromosomal domains, and, on the other extreme, at the scale of oligonucleotides. The intermediate, kilobase-scale heterogeneity in mutation risk is less well characterized. Here, by analyzing thousands of somatic genomes, we studied mutation risk gradients along gene bodies, representing a genomic scale spanning roughly 1-10 kb, hypothesizing that different mutational mechanisms are differently distributed across gene segments. The main heterogeneity concerns several kilobases at the transcription start site and further downstream into 5' ends of gene bodies; these are commonly hypomutated with several mutational signatures, most prominently the ubiquitous C > T changes at CpG dinucleotides. The width and shape of this mutational coldspot at 5' gene ends is variable across genes, and corresponds to variable interval of lowered DNA methylation depending on gene activity level and regulation. Such hypomutated loci, at 5' gene ends or elsewhere, correspond to DNA hypomethylation that can associate with various landmarks, including intragenic enhancers, Polycomb-marked regions, or chromatin loop anchor points. Tissue-specific DNA hypomethylation begets tissue-specific local hypomutation. Of note, direction of mutation risk is inverted for AID/APOBEC3 cytosine deaminase activity, whose signatures are enriched in hypomethylated regions.

The genomic landscape of Vk*MYC myeloma highlights shared pathways of transformation between mice and humans

Multiple myeloma (MM) is a heterogeneous disease characterized by frequent MYC translocations. Sporadic MYC activation in the germinal center of genetically engineered Vk*MYC mice is sufficient to induce plasma cell tumors in which a variety of secondary mutations are spontaneously acquired and selected over time. Analysis of 119 Vk*MYC myeloma reveals recurrent copy number alterations, structural variations, chromothripsis, driver mutations, apolipoprotein B mRNA-editing enzyme, catalytic polypeptide (APOBEC) mutational activity, and a progressive decrease in immunoglobulin transcription that inversely correlates with proliferation. Moreover, we identify frequent insertional mutagenesis by endogenous retro-elements as a murine specific mechanism to activate NF-kB and IL6 signaling pathways shared with human MM. Despite the increased genomic complexity associated with progression, advanced tumors remain dependent on MYC. In summary, here we credential the Vk*MYC mouse as a unique resource to explore MM genomic evolution and describe a fully annotated collection of diverse and immortalized murine MM tumors.

The Association between Mutational Signatures and Clinical Outcomes among Patients with Early-Onset Breast Cancer

Early-onset breast cancer (EoBC), defined by a diagnosis <40 years of age, is associated with poor prognosis. This study investigated the mutational landscape of non-metastatic EoBC and the prognostic relevance of mutational signatures using 100 tumour samples from Alberta, Canada. The MutationalPatterns package in R/Bioconductor was used to extract de novo single-base substitution (SBS) and insertion-deletion (indel) mutational signatures and to fit COSMIC SBS and indel signatures. We assessed associations between these signatures and clinical characteristics of disease, in addition to recurrence-free (RFS) and overall survival (OS). Five SBS and two indel signatures were extracted. The SBS13-like signature had higher relative contributions in the HER2-enriched subtype. Patients with higher than median contribution tended to have better RFS after adjustment for other prognostic factors (HR = 0.29; 95% CI: 0.08-1.06). An unsupervised clustering algorithm based on absolute contribution revealed three clusters of fitted COSMIC SBS signatures, but cluster membership was not associated with clinical variables or survival outcomes. The results of this exploratory study reveal various SBS and indel signatures may be associated with clinical features of disease and prognosis. Future studies with larger samples are required to better understand the mechanistic underpinnings of disease progression and treatment response in EoBC.

Negative Selection Allows DNA Mismatch Repair-Deficient Mouse Fibroblasts In Vitro to Tolerate High Levels of Somatic Mutations

Substantial numbers of somatic mutations have been found to accumulate with age in different human tissues. Clonal cellular amplification of some of these mutations can cause cancer and other diseases. However, it is as yet unclear if and to what extent an increased burden of random mutations can affect cellular function without clonal amplification. We tested this in cell culture, which avoids the limitation that an increased mutation burden in vivo typically leads to cancer. We performed single-cell whole-genome sequencing of primary fibroblasts from DNA mismatch repair (MMR) deficient Msh2-/- mice and littermate control animals after long-term passaging. Apart from analyzing somatic mutation burden we analyzed clonality, mutational signatures, and hotspots in the genome, characterizing the complete landscape of somatic mutagenesis in normal and MMR-deficient mouse primary fibroblasts during passaging. While growth rate of Msh2-/- fibroblasts was not significantly different from the controls, the number of de novo single-nucleotide variants (SNVs) increased linearly up until at least 30,000 SNVs per cell, with the frequency of small insertions and deletions (INDELs) plateauing in the Msh2-/- fibroblasts to about 10,000 INDELS per cell. We provide evidence for negative selection and large-scale mutation-driven population changes, including significant clonal expansion of preexisting mutations and widespread cell-strain-specific hotspots. Overall, our results provide evidence that increased somatic mutation burden drives significant cell evolutionary changes in a dynamic cell culture system without significant effects on growth. Since similar selection processes against mutations preventing organ and tissue dysfunction during aging are difficult to envision, these results suggest that increased somatic mutation burden can play a causal role in aging and diseases other than cancer.

Non-hypermutator cancers access driver mutations through reversals in germline mutational bias

Cancer is an evolutionary disease driven by mutations in asexually-reproducing somatic cells. In asexual microbes, bias reversals in the mutation spectrum can speed adaptation by increasing access to previously undersampled beneficial mutations. By analyzing tumors from 20 tissues, along with normal tissue and the germline, we demonstrate this effect in cancer. Non-hypermutated tumors reverse the germline mutation bias and have consistent spectra across tissues. These spectra changes carry the signature of hypoxia, and they facilitate positive selection in cancer genes. Hypermutated and non-hypermutated tumors thus acquire driver mutations differently: hypermutated tumors by higher mutation rates and non-hypermutated tumors by changing the mutation spectrum to reverse the germline mutation bias.

Isocitrate Dehydrogenase Mutations in Cancer: Mechanisms of Transformation and Metabolic Liability

Isocitrate dehydrogenase 1 and 2 (IDH1 and IDH2) are metabolic enzymes that interconvert isocitrate and 2-oxoglutarate (2OG). Gain-of-function mutations in IDH1 and IDH2 occur in a number of cancers, including acute myeloid leukemia, glioma, cholangiocarcinoma, and chondrosarcoma. These mutations cripple the wild-type activity of IDH and cause the enzymes to catalyze a partial reverse reaction in which 2OG is reduced but not carboxylated, resulting in production of the (R)-enantiomer of 2-hydroxyglutarate ((R)-2HG). (R)-2HG accumulation in IDH-mutant tumors results in profound dysregulation of cellular metabolism. The most well-characterized oncogenic effects of (R)-2HG involve the dysregulation of 2OG-dependent epigenetic tumor-suppressor enzymes. However, (R)-2HG has many other effects in IDH-mutant cells, some that promote transformation and others that induce metabolic dependencies. Herein, we review how cancer-associated IDH mutations impact epigenetic regulation and cellular metabolism and discuss how these effects can potentially be leveraged to therapeutically target IDH-mutant tumors.

NeoAgDT: optimization of personal neoantigen vaccine composition by digital twin simulation of a cancer cell population

MOTIVATION

Neoantigen vaccines make use of tumor-specific mutations to enable the patient's immune system to recognize and eliminate cancer. Selecting vaccine elements, however, is a complex task which needs to take into account not only the underlying antigen presentation pathway but also tumor heterogeneity.

RESULTS

Here, we present NeoAgDT, a two-step approach consisting of: (i) simulating individual cancer cells to create a digital twin of the patient's tumor cell population and (ii) optimizing the vaccine composition by integer linear programming based on this digital twin. NeoAgDT shows improved selection of experimentally validated neoantigens over ranking-based approaches in a study of seven patients.

AVAILABILITY AND IMPLEMENTATION

The NeoAgDT code is published on Github: https://github.com/nec-research/neoagdt.

Triple Primary Cancers: An Analysis of Genetic and Environmental Factors

The occurrence of multiple primary cancers (MPC) is thought to reflect increased cancer susceptibility in patients due to a combination of genetic and environmental factors. Here we conducted a retrospective review of 2,894 consecutive patients evaluated at a single institution and identified 31 (1.14%) individuals with a history of three or more primary cancers, then analyzed the genetic and environmental influences associated with their propensity for developing malignancies. We found that 35.5% of patients had a hereditary cancer syndrome (HCS), with high penetrance HCS in 72.7% of cases, suggesting that monogenic causes underly a significant proportion of triple primary cancer risk. Analysis of cancer frequencies found that the diagnosis of breast cancer was associated with a significantly lower likelihood of HCS, while the diagnosis of colorectal, prostate, and pancreas cancer was associated with a significantly higher likelihood of HCS. Comparison of HCS-positive and HCS-negative patients revealed similar demographic characteristics, mean age at first diagnosis, and family history of cancer. Moreover, no significant differences in lifestyle behaviors, occupational exposures, chronic health conditions, or treatment with chemotherapy and radiation were observed between HCS-positive and -negative groups, though outliers in tobacco smoking, as well as systemic treatment after both first and second primary cancers were observed. These findings indicate a robust contribution of HCS to cancer susceptibility among patients with triple primary cancers while environmental influences were less evident. This emphasizes the need for larger MPC cohorts incorporating additional genetic and environmental factors to more comprehensively characterize drivers of cancer risk.

PREVENTION RELEVANCE

In patients with three or more primary cancers, genetic predisposition explained a significant proportion of cases; however, treatment history, lifestyle habits, and other exposures appeared to play a less significant role. This highlights the value of early genetic screening and the need to develop more sensitive markers of cancer susceptibility. See related Spotlight, p. 193.

Direct Measurement of 8OG Syn-Anti Flips in Mutagenic 8OG·A and Long-Range Damage-Dependent Hoogsteen Breathing Dynamics Using 1H CEST NMR

Elucidating how damage impacts DNA dynamics is essential for understanding the mechanisms of damage recognition and repair. Many DNA lesions alter their propensities to form low-populated and short-lived conformational states. However, NMR methods to measure these dynamics require isotopic enrichment, which is difficult for damaged nucleotides. Here, we demonstrate the utility of the 1H chemical exchange saturation transfer (CEST) NMR experiment in measuring the dynamics of oxidatively damaged 8-oxoguanine (8OG) in the mutagenic 8OGsyn·Aanti mismatch. Using 8OG-H7 as an NMR probe of the damaged base, we directly measured 8OG syn-anti flips to form a lowly populated (pop. ∼ 5%) and short-lived (lifetime ∼50 ms) nonmutagenic 8OGanti·Aanti. These exchange parameters were in quantitative agreement with values from 13C off-resonance R1ρ and CEST on the labeled partner adenine. The Watson-Crick-like 8OGsyn·Aanti mismatch also rescued the kinetics of Hoogsteen motions at distant A-T base pairs, which the G·A mismatch had slowed down. The results lend further support for 8OGanti·Aanti as a minor conformational state of 8OG·A, reveal that 8OG damage can impact Hoogsteen dynamics at a distance, and demonstrate the utility of 1H CEST for measuring damage-dependent dynamics in unlabeled DNA.

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

MOTIVATION

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

RESULTS

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

AVAILABILITY AND IMPLEMENTATION

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

How CRISPR Is Revolutionizing the Generation of New Models for Cancer Research

Cancers arise through acquisition of mutations in genes that regulate core biological processes like cell proliferation and cell death. Decades of cancer research have led to the identification of genes and mutations causally involved in disease development and evolution, yet defining their precise function across different cancer types and how they influence therapy responses has been challenging. Mouse models have helped define the in vivo function of cancer-associated alterations, and genome-editing approaches using CRISPR have dramatically accelerated the pace at which these models are developed and studied. Here, we highlight how CRISPR technologies have impacted the development and use of mouse models for cancer research and discuss the many ways in which these rapidly evolving platforms will continue to transform our understanding of this disease.

Frequency and Nature of Genomic Alterations in ERBB2-Altered Urothelial Bladder Cancer

BACKGROUND

Human epidermal growth factor-2 (HER2) overexpression is an oncogenic driver in many solid tumors, including urothelial bladder cancer (UBC). In addition, activating mutations in the ERBB2 gene have been shown to play an oncogenic role similar to ERBB2 amplification.

OBJECTIVE

To describe and compare the frequency and nature of genomic alterations (GA) of ERBB2-altered (mutations, amplification) and ERBB2 wild-type UBC.

PATIENTS AND METHODS

Using a hybrid capture-based comprehensive profiling assay, 9518 UBC cases were grouped by ERBB2 alteration and evaluated for all classes of genomic alterations (GA), tumor mutational burden (TMB), microsatellite instability (MSI), genome-wide loss of heterozygosity (gLOH), and genomic mutational signature. PD-L1 expression was measured by immunohistochemistry (Dako 22C3). Categorical statistical comparisons were performed using Fisher's exact tests.

RESULTS

A total of 602 (6.3%) UBC cases featured ERBB2 extracellular domain short variant (SV) GA (ECDmut+), 253 (2.7%) cases featured ERBB2 kinase domain SV GA (KDmut+), 866 (9.1%) cases had ERBB2 amplification (amp+), and 7797 (81.9%) cases were ERBB2 wild-type (wt). European genetic ancestry of ECDmut+ was higher than ERBB2wt. Numerous significant associations were observed when comparing GA by group. Notably among these, CDKN2A/MTAP loss were more frequent in ERBB2wt versus ECDmut+ and amp+. ERBB3 GA were more frequent in ECDmut+ and KDmut+ than ERBB2wt. TERT GA were more frequent in ECDmut+, KDmut+, and amp+ versus ERBB2wt. TOP2A amplification was significantly more common in ECDmut+ and amp+ versus ERBB2wt, and TP53 SV GA were significantly higher in ERBB2 amp+ versus ERBB2wt. Mean TMB levels were significantly higher in ECDmut+, KDmut+, and amp+ than in ERBB2wt. Apolipoprotein B mRNA-editing enzyme, catalytic polypeptides (APOBEC) signature was more frequent in ECDmut+, KDmut+, and amp+ versus ERBB2wt. No significant differences were observed in PD-L1 status between groups, while gLOH-high status was more common in amp+ versus ERBB2wt. MSI-high status was more frequent in KDmut+ versus ERBB2wt, and in ERBB2wt than in amp+.

CONCLUSIONS

We noted important differences in co-occurring GA in ERBB2-altered (ECDmut+, KDmut+, amp+) versus ERBB2wt UBC, as well as higher mean TMB and higher APOBEC mutational signature in the ERBB2-altered groups. Our results can help refine future clinical trial designs and elucidate possible response and resistance mechanisms for ERBB2-altered UBC.

Diffsig: Associating Risk Factors with Mutational Signatures

BACKGROUND

Somatic mutational signatures elucidate molecular vulnerabilities to therapy, and therefore detecting signatures and classifying tumors with respect to signatures has clinical value. However, identifying the etiology of the mutational signatures remains a statistical challenge, with both small sample sizes and high variability in classification algorithms posing barriers. As a result, few signatures have been strongly linked to particular risk factors.

METHODS

Here, we develop a statistical model, Diffsig, for estimating the association of one or more continuous or categorical risk factors with DNA mutational signatures. Diffsig takes into account the uncertainty associated with assigning signatures to samples as well as multiple risk factors' simultaneous effect on observed DNA mutations.

RESULTS

We applied Diffsig to breast cancer data to assess relationships between five established breast-relevant mutational signatures and etiologic variables, confirming known mechanisms of cancer development. In simulation, our model was capable of accurately estimating expected associations in a variety of contexts.

CONCLUSIONS

Diffsig allows researchers to quantify and perform inference on the associations of risk factors with mutational signatures.

IMPACT

We expect Diffsig to provide more robust associations of risk factors with signatures to lead to better understanding of the tumor development process and improved models of tumorigenesis.

APOBEC3 mutagenesis drives therapy resistance in breast cancer

Acquired genetic alterations commonly drive resistance to endocrine and targeted therapies in metastatic breast cancer 1-7 , however the underlying processes engendering these diverse alterations are largely uncharacterized. To identify the mutational processes operant in breast cancer and their impact on clinical outcomes, we utilized a well-annotated cohort of 3,880 patient samples with paired tumor-normal sequencing data. The mutational signatures associated with apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like 3 (APOBEC3) enzymes were highly prevalent and enriched in post-treatment compared to treatment-naïve hormone receptor-positive (HR+) cancers. APOBEC3 mutational signatures were independently associated with shorter progression-free survival on antiestrogen plus CDK4/6 inhibitor combination therapy in patients with HR+ metastatic breast cancer. Whole genome sequencing (WGS) of breast cancer models and selected paired primary-metastatic samples demonstrated that active APOBEC3 mutagenesis promoted resistance to both endocrine and targeted therapies through characteristic alterations such as RB1 loss-of-function mutations. Evidence of APOBEC3 activity in pre-treatment samples illustrated a pervasive role for this mutational process in breast cancer evolution. The study reveals APOBEC3 mutagenesis to be a frequent mediator of therapy resistance in breast cancer and highlights its potential as a biomarker and target for overcoming resistance.

De novo detection of somatic mutations in high-throughput single-cell profiling data sets

Characterization of somatic mutations at single-cell resolution is essential to study cancer evolution, clonal mosaicism and cell plasticity. Here, we describe SComatic, an algorithm designed for the detection of somatic mutations in single-cell transcriptomic and ATAC-seq (assay for transposase-accessible chromatin sequence) data sets directly without requiring matched bulk or single-cell DNA sequencing data. SComatic distinguishes somatic mutations from polymorphisms, RNA-editing events and artefacts using filters and statistical tests parameterized on non-neoplastic samples. Using >2.6 million single cells from 688 single-cell RNA-seq (scRNA-seq) and single-cell ATAC-seq (scATAC-seq) data sets spanning cancer and non-neoplastic samples, we show that SComatic detects mutations in single cells accurately, even in differentiated cells from polyclonal tissues that are not amenable to mutation detection using existing methods. Validated against matched genome sequencing and scRNA-seq data, SComatic achieves F1 scores between 0.6 and 0.7 across diverse data sets, in comparison to 0.2-0.4 for the second-best performing method. In summary, SComatic permits de novo mutational signature analysis, and the study of clonal heterogeneity and mutational burdens at single-cell resolution.

3D genomic mapping reveals multifocality of human pancreatic precancers

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

Single-mitosis dissection of acute and chronic DNA mutagenesis and repair

How chronic mutational processes and punctuated bursts of DNA damage drive evolution of the cancer genome is poorly understood. Here, we demonstrate a strategy to disentangle and quantify distinct mechanisms underlying genome evolution in single cells, during single mitoses and at single-strand resolution. To distinguish between chronic (reactive oxygen species (ROS)) and acute (ultraviolet light (UV)) mutagenesis, we microfluidically separate pairs of sister cells from the first mitosis following burst UV damage. Strikingly, UV mutations manifest as sister-specific events, revealing mirror-image mutation phasing genome-wide. In contrast, ROS mutagenesis in transcribed regions is reduced strand agnostically. Successive rounds of genome replication over persisting UV damage drives multiallelic variation at CC dinucleotides. Finally, we show that mutation phasing can be resolved to single strands across the entire genome of liver tumors from F1 mice. This strategy can be broadly used to distinguish the contributions of overlapping cancer relevant mutational processes.

Her2 amplification, Rel-A, and Bach1 can influence APOBEC3A expression in breast cancer cells

APOBEC-induced mutations occur in 50% of sequenced human tumors, with APOBEC3A (A3A) being a major contributor to mutagenesis in breast cancer cells. The mechanisms that cause A3A activation and mutagenesis in breast cancers are still unknown. Here, we describe factors that influence basal A3A mRNA transcript levels in breast cancer cells. We found that basal A3A mRNA correlates with A3A protein levels and predicts the amount of APOBEC signature mutations in a panel of breast cancer cell lines, indicating that increased basal transcription may be one mechanism leading to breast cancer mutagenesis. We also show that alteration of ERBB2 expression can drive A3A mRNA levels, suggesting the enrichment of the APOBEC mutation signature in Her2-enriched breast cancer could in part result from elevated A3A transcription. Hierarchical clustering of transcripts in primary breast cancers determined that A3A mRNA was co-expressed with other genes functioning in viral restriction and interferon responses. However, reduction of STAT signaling via inhibitors or shRNA in breast cancer cell lines had only minor impact on A3A abundance. Analysis of single cell RNA-seq from primary tumors indicated that A3A mRNA was highest in infiltrating immune cells within the tumor, indicating that correlations of A3A with STAT signaling in primary tumors may be result from higher immune infiltrates and are not reflective of STAT signaling controlling A3A expression in breast cancer cells. Analysis of ATAC-seq data in multiple breast cancer cell lines identified two transcription factor sites in the APOBEC3A promoter region that could promote A3A transcription. We determined that Rel-A, and Bach1, which have binding sites in these peaks, elevated basal A3A expression. Our findings highlight a complex and variable set of transcriptional activators for A3A in breast cancer cells.

Extrinsic and intrinsic drivers of natural killer cell clonality

Clonal expansion of antigen-specific lymphocytes is the fundamental mechanism enabling potent adaptive immune responses and the generation of immune memory. Accompanied by pronounced epigenetic remodeling, the massive proliferation of individual cells generates a critical mass of effectors for the control of acute infections, as well as a pool of memory cells protecting against future pathogen encounters. Classically associated with the adaptive immune system, recent work has demonstrated that innate immune memory to human cytomegalovirus (CMV) infection is stably maintained as large clonal expansions of natural killer (NK) cells, raising questions on the mechanisms for clonal selection and expansion in the absence of re-arranged antigen receptors. Here, we discuss clonal NK cell memory in the context of the mechanisms underlying clonal competition of adaptive lymphocytes and propose alternative selection mechanisms that might decide on the clonal success of their innate counterparts. We propose that the integration of external cues with cell-intrinsic sources of heterogeneity, such as variegated receptor expression, transcriptional states, and somatic variants, compose a bottleneck for clonal selection, contributing to the large size of memory NK cell clones.

Clinical and biological landscape of constitutional mismatch-repair deficiency syndrome: an International Replication Repair Deficiency Consortium cohort study

BACKGROUND

Constitutional mismatch repair deficiency (CMMRD) syndrome is a rare and aggressive cancer predisposition syndrome. Because a scarcity of data on this condition contributes to management challenges and poor outcomes, we aimed to describe the clinical spectrum, cancer biology, and impact of genetics on patient survival in CMMRD.

METHODS

In this cohort study, we collected cross-sectional and longitudinal data on all patients with CMMRD, with no age limits, registered with the International Replication Repair Deficiency Consortium (IRRDC) across more than 50 countries. Clinical data were extracted from the IRRDC database, medical records, and physician-completed case record forms. The primary objective was to describe the clinical features, cancer spectrum, and biology of the condition. Secondary objectives included estimations of cancer incidence and of the impact of the specific mismatch-repair gene and genotype on cancer onset and survival, including after cancer surveillance and immunotherapy interventions.

FINDINGS

We analysed data from 201 patients (103 males, 98 females) enrolled between June 5, 2007 and Sept 9, 2022. Median age at diagnosis of CMMRD or a related cancer was 8·9 years (IQR 5·9-12·6), and median follow-up from diagnosis was 7·2 years (3·6-14·8). Endogamy among minorities and closed communities contributed to high homozygosity within countries with low consanguinity. Frequent dermatological manifestations (117 [93%] of 126 patients with complete data) led to a clinical overlap with neurofibromatosis type 1 (35 [28%] of 126). 339 cancers were reported in 194 (97%) of 201 patients. The cumulative cancer incidence by age 18 years was 90% (95% CI 80-99). Median time between cancer diagnoses for patients with more than one cancer was 1·9 years (IQR 0·8-3·9). Neoplasms developed in 15 organs and included early-onset adult cancers. CNS tumours were the most frequent (173 [51%] cancers), followed by gastrointestinal (75 [22%]), haematological (61 [18%]), and other cancer types (30 [9%]). Patients with CNS tumours had the poorest overall survival rates (39% [95% CI 30-52] at 10 years from diagnosis; log-rank p<0·0001 across four cancer types), followed by those with haematological cancers (67% [55-82]), gastrointestinal cancers (89% [81-97]), and other solid tumours (96% [88-100]). All cancers showed high mutation and microsatellite indel burdens, and pathognomonic mutational signatures. MLH1 or MSH2 variants caused earlier cancer onset than PMS2 or MSH6 variants, and inferior survival (overall survival at age 15 years 63% [95% CI 55-73] for PMS2, 49% [35-68] for MSH6, 19% [6-66] for MLH1, and 0% for MSH2; p<0·0001). Frameshift or truncating variants within the same gene caused earlier cancers and inferior outcomes compared with missense variants (p<0·0001). The greater deleterious effects of MLH1 and MSH2 variants as compared with PMS2 and MSH6 variants persisted despite overall improvements in survival after surveillance or immune checkpoint inhibitor interventions.

INTERPRETATION

The very high cancer burden and unique genomic landscape of CMMRD highlight the benefit of comprehensive assays in timely diagnosis and precision approaches toward surveillance and immunotherapy. These data will guide the clinical management of children and patients who survive into adulthood with CMMRD.

FUNDING

The Canadian Institutes for Health Research, Stand Up to Cancer, Children's Oncology Group National Cancer Institute Community Oncology Research Program, Canadian Cancer Society, Brain Canada, The V Foundation for Cancer Research, BioCanRx, Harry and Agnieszka Hall, Meagan's Walk, BRAINchild Canada, The LivWise Foundation, St Baldrick Foundation, Hold'em for Life, and Garron Family Cancer Center.

Geographic variation of mutagenic exposures in kidney cancer genomes

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

CDK4/6 Inhibition - Therapy Sequences and the Quest to Find the Best Biomarkers - an Overview of Current Programs

In recent years, new targeted therapies have been developed to treat patients with hormone receptor-positive (HR+)/human epidermal growth factor receptor 2-negative (HER2-) breast cancer. Some of these therapies have not just become the new therapy standard but also led to significantly longer overall survival rates. The cyclin-dependent kinase 4 and 6 inhibitors (CDK4/6i) have become the therapeutic standard for first-line therapy. Around 70 - 80% of patients are treated with a CDK4/6i. In recent years, a number of biomarkers associated with progression, clonal selection or evolution have been reported for CDK4/6i and their endocrine combination partners. Understanding the mechanisms behind treatment efficacy and resistance is important. A better understanding could contribute to planning the most effective therapeutic sequences and utilizing basic molecular information to overcome endocrine resistance. One study with large numbers of patients which aims to elucidate these mechanisms is the Comprehensive Analysis of sPatial, TempORal and molecular patterns of ribociclib efficacy and resistance in advanced Breast Cancer patients (CAPTOR BC) trial. This overview summarizes the latest clinical research on resistance to endocrine therapies, focusing on CDK4/6 inhibitors and discussing current study concepts.

Breaking paradigms: Long non-coding RNAs forming gene fusions with potential implications in cancer

Long non-coding RNAs (lncRNAs) are non-coding RNAs longer than 200 nucleotides with dynamic regulatory functions. They interact with a wide range of molecules such as DNA, RNA, and proteins to modulate diverse cellular functions through several mechanisms and, if deregulated, they can lead to cancer development and progression. Recently, it has been described that lncRNAs are susceptible to form gene fusions with mRNAs or other lncRNAs, breaking the paradigm of gene fusions consisting mainly of protein-coding genes. However, their biological significance in the tumor phenotype is still uncertain. Therefore, their recent identification opens a new line of research to study their biological role in tumorigenesis, and their potential as biomarkers with clinical relevance or as therapeutic targets. The present study aimed to review the lncRNA fusions identified so far and to know which of them have been associated with a potential function. We address the current challenges to deepen their study as well as the reasons why they represent a future therapeutic window in cancer.

Chromosome evolution screens recapitulate tissue-specific tumor aneuploidy patterns

Whole chromosome and arm-level copy number alterations occur at high frequencies in tumors, but their selective advantages, if any, are poorly understood. Here, utilizing unbiased whole chromosome genetic screens combined with in vitro evolution to generate arm- and subarm-level events, we iteratively selected the fittest karyotypes from aneuploidized human renal and mammary epithelial cells. Proliferation-based karyotype selection in these epithelial lines modeled tissue-specific tumor aneuploidy patterns in patient cohorts in the absence of driver mutations. Hi-C-based translocation mapping revealed that arm-level events usually emerged in multiples of two via centromeric translocations and occurred more frequently in tetraploids than diploids, contributing to the increased diversity in evolving tetraploid populations. Isogenic clonal lineages enabled elucidation of pro-tumorigenic mechanisms associated with common copy number alterations, revealing Notch signaling potentiation as a driver of 1q gain in breast cancer. We propose that intrinsic, tissue-specific proliferative effects underlie tumor copy number patterns in cancer.

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.

Associations of combined phenotypic aging and genetic risk with incident cancer: A prospective cohort study

Background

Age is the most important risk factor for cancer, but aging rates are heterogeneous across individuals. We explored a new measure of aging-Phenotypic Age (PhenoAge)-in the risk prediction of site-specific and overall cancer.

Methods

Using Cox regression models, we examined the association of Phenotypic Age Acceleration (PhenoAgeAccel) with cancer incidence by genetic risk group among 374,463 participants from the UK Biobank. We generated PhenoAge using chronological age and nine biomarkers, PhenoAgeAccel after subtracting the effect of chronological age by regression residual, and an incidence-weighted overall cancer polygenic risk score (CPRS) based on 20 cancer site-specific polygenic risk scores (PRSs).

Results

Compared with biologically younger participants, those older had a significantly higher risk of overall cancer, with hazard ratios (HRs) of 1.22 (95% confidence interval, 1.18-1.27) in men, and 1.26 (1.22-1.31) in women, respectively. A joint effect of genetic risk and PhenoAgeAccel was observed on overall cancer risk, with HRs of 2.29 (2.10-2.51) for men and 1.94 (1.78-2.11) for women with high genetic risk and older PhenoAge compared with those with low genetic risk and younger PhenoAge. PhenoAgeAccel was negatively associated with the number of healthy lifestyle factors (Beta = -1.01 in men, p<0.001; Beta = -0.98 in women, p<0.001).

Conclusions

Within and across genetic risk groups, older PhenoAge was consistently related to an increased risk of incident cancer with adjustment for chronological age and the aging process could be retarded by adherence to a healthy lifestyle.

Funding

This work was supported by the National Natural Science Foundation of China (82230110, 82125033, 82388102 to GJ; 82273714 to MZ); and the Excellent Youth Foundation of Jiangsu Province (BK20220100 to MZ).

Integrative multi-omics characterization of hepatocellular carcinoma in Hispanic patients

Background

The incidence and mortality rates of hepatocellular carcinoma (HCC) among Hispanics in the United States are much higher than those of non-Hispanic whites. We conducted comprehensive multi-omics analyses to understand molecular alterations in HCC among Hispanic patients.

Methods

Paired tumor and adjacent non-tumor samples were collected from 31 Hispanic HCC in South Texas (STX-Hispanic) for genomic, transcriptomic, proteomic, and metabolomic profiling. Additionally, serum lipids were profiled in 40 Hispanic and non-Hispanic patients with or without clinically diagnosed HCC.

Results

Exome sequencing revealed high mutation frequencies of AXIN2 and CTNNB1 in STX Hispanic HCCs, suggesting a predominant activation of the Wnt/β-catenin pathway. The TERT promoter mutation frequency was also remarkably high in the Hispanic cohort. Cell cycles and liver functions were identified as positively- and negatively-enriched, respectively, with gene set enrichment analysis. Gene sets representing specific liver metabolic pathways were associated with dysregulation of corresponding metabolites. Negative enrichment of liver adipogenesis and lipid metabolism corroborated with a significant reduction in most lipids in the serum samples of HCC patients. Two HCC subtypes from our Hispanic cohort were identified and validated with the TCGA liver cancer cohort. The subtype with better overall survival showed higher activity of immune and angiogenesis signatures, and lower activity of liver function-related gene signatures. It also had higher levels of immune checkpoint and immune exhaustion markers.

Conclusions

Our study revealed some specific molecular features of Hispanic HCC and potential biomarkers for therapeutic management of HCC and provides a unique resource for studying Hispanic HCC.

Identifying somatic fingerprints of cancers defined by germline and environmental risk factors

Numerous studies over the past generation have identified germline variants that increase specific cancer risks. Simultaneously, a revolution in sequencing technology has permitted high-throughput annotations of somatic genomes characterizing individual tumors. However, examining the relationship between germline variants and somatic alteration patterns is hugely challenged by the large numbers of variants in a typical tumor, the rarity of most individual variants, and the heterogeneity of tumor somatic fingerprints. In this article, we propose statistical methodology that frames the investigation of germline-somatic relationships in an interpretable manner. The method uses meta-features embodying biological contexts of individual somatic alterations to implicitly group rare mutations. Our team has used this technique previously through a multilevel regression model to diagnose with high accuracy tumor site of origin. Herein, we further leverage topic models from computational linguistics to achieve interpretable lower-dimensional embeddings of the meta-features. We demonstrate how the method can identify distinctive somatic profiles linked to specific germline variants or environmental risk factors. We illustrate the method using The Cancer Genome Atlas whole-exome sequencing data to characterize somatic tumor fingerprints in breast cancer patients with germline BRCA1/2 mutations and in head and neck cancer patients exposed to human papillomavirus.

A glycolytic metabolite bypasses "two-hit" tumor suppression by BRCA2

Knudson's "two-hit" paradigm posits that carcinogenesis requires inactivation of both copies of an autosomal tumor suppressor gene. Here, we report that the glycolytic metabolite methylglyoxal (MGO) transiently bypasses Knudson's paradigm by inactivating the breast cancer suppressor protein BRCA2 to elicit a cancer-associated, mutational single-base substitution (SBS) signature in nonmalignant mammary cells or patient-derived organoids. Germline monoallelic BRCA2 mutations predispose to these changes. An analogous SBS signature, again without biallelic BRCA2 inactivation, accompanies MGO accumulation and DNA damage in Kras-driven, Brca2-mutant murine pancreatic cancers and human breast cancers. MGO triggers BRCA2 proteolysis, temporarily disabling BRCA2's tumor suppressive functions in DNA repair and replication, causing functional haploinsufficiency. Intermittent MGO exposure incites episodic SBS mutations without permanent BRCA2 inactivation. Thus, a metabolic mechanism wherein MGO-induced BRCA2 haploinsufficiency transiently bypasses Knudson's two-hit requirement could link glycolysis activation by oncogenes, metabolic disorders, or dietary challenges to mutational signatures implicated in cancer evolution.

Development of Allosteric Small Molecule APOBEC3B Inhibitors from In Silico Screening

APOBEC3B cytosine deaminase contributes to the mutational burdens of tumors, resulting in tumor progression and therapy resistance. Small molecule APOBEC3B inhibitors have potential to slow or mitigate these detrimental outcomes. Through molecular dynamics (MD) simulations and computational solvent mapping analysis, we identified a novel putative allosteric pocket on the C-terminal domain of APOBEC3B (A3Bctd), and virtually screened the ChemBridge Diversity Set (N~110,000) against both the active and potential allosteric sites. Selected high-scoring compounds were subsequently purchased, characterized for purity and composition, and tested in biochemical assays, which yielded 13 hit compounds. Orthogonal NMR assays verified binding to the target protein. Initial selectivity studies suggest these compounds preferentially target A3Bctd over related deaminase APOBEC3A (A3A), and MD simulations indicate this selectivity may be due to the steric repulsion from H56 that is unique to A3A. Taken together, our studies represent the first virtual screening effort against A3Bctd that has yielded candidate inhibitors suitable for further development.

Insight into NSCLC through novel analysis of gene interactions and characteristics

Around 80 to 85% of all lung cancers are non-small cell lung cancer (NSCLC). Previous research has aimed at exploring the genetic basis of NSCLC through individual approaches, but studies have yet to investigate the results of combining them. Here we show that analyzing NSCLC genetics through three approaches simultaneously creates unique insights into our understanding of the disease. Through a combination of previous research and bioinformatics tools, we determined 35 NSCLC candidate genes. We analyzed these genes in 3 different approaches. First, we found the gene fusions between these candidate genes. Second, we found the common superfamilies between genes. Finally, we identified mutational signatures that are possibly associated with NSCLC. Each approach has its individual, unique results. Fusion relationships identify specific gene fusion targets, common superfamilies identify possible avenues to determine novel target genes, and identifying NSCLC associated mutational signatures has diagnostic and prognostic benefits. Combining the approaches, we found that gene CD74 has significant fusion relationships, but it has no association with the other two approaches, suggesting that CD74 is associated with NSCLC mainly because of its fusion relationships. Targeting the gene fusions of CD74 may be an alternative NSCLC treatment. This genetic analysis has indeed created unique insight into NSCLC genes. Both the results from each of the approaches separately and combined allow pursuit of more effective treatment strategies for this cancer. The methodology presented can also apply to other cancers, creating insights that current analytical methods could not find.

A glycolytic metabolite that drives BRCA2 haploinsufficiency

Many types of tumor cells alter metabolic pathways to meet their energy and biosynthetic demands for proliferation or stress adaptation. In this issue of Cell, Kong et al. find that the glycolytic metabolite methylglyoxal causes cancer-associated mutant single-base substitution features by inducing BRCA2 proteolysis, leading to functional haploinsufficiency of BRCA2.

Prioritization of therapeutic targets for cancers using integrative multi-omics analysis

BACKGROUND

The integration of transcriptomic, proteomic, druggable genetic and metabolomic association studies facilitated a comprehensive investigation of molecular features and shared pathways for cancers' development and progression.

METHODS

Comprehensive approaches consisting of transcriptome-wide association studies (TWAS), proteome-wide association studies (PWAS), summary-data-based Mendelian randomization (SMR) and MR were performed to identify genes significantly associated with cancers. The results identified in above analyzes were subsequently involved in phenotype scanning and enrichment analyzes to explore the possible health effects and shared pathways. Additionally, we also conducted MR analysis   to investigate metabolic pathways related to cancers.

RESULTS

Totally 24 genes (18 transcriptomic, 1 proteomic and 5 druggable genetic) showed significant associations with cancers risk. All genes identified in multiple methods were mainly enriched in nuclear factor erythroid 2-related factor 2 (NRF2) pathway. Additionally, biosynthesis of ubiquinol and urate were found to play an important role in gastrointestinal tumors.

CONCLUSIONS

A set of putatively causal genes and pathways relevant to cancers were identified in this study, shedding light on the shared biological processes for tumorigenesis and providing compelling genetic evidence to prioritize anti-cancer drugs development.

Focus on current and emerging treatment options for glioma: A comprehensive review

This comprehensive review delves into the current updates and challenges associated with the management of low-grade gliomas (LGG), the predominant primary tumors in the central nervous system. With a general incidence rate of 5.81 per 100000, gliomas pose a significant global concern, necessitating advancements in treatment techniques to reduce mortality and morbidity. This review places a particular focus on immunotherapies, discussing promising agents such as Zotiraciclib and Lerapolturev. Zotiraciclib, a CDK9 inhibitor, has demonstrated efficacy in glioblastoma treatment in preclinical and clinical studies, showing its potential as a therapeutic breakthrough. Lerapolturev, a viral immunotherapy, induces inflammation in glioblastoma and displays positive outcomes in both adult and pediatric patients. Exploration of immunotherapy extends to Pembrolizumab, Nivolumab, and Entrectinib, revealing the challenges and variabilities in patient responses. Despite promising preclinical data, the monoclonal antibody Depatuxizumab has proven ineffective in glioblastoma treatment, emphasizing the critical need to understand resistance mechanisms. The review also covers the success of radiation therapy in pediatric LGG, with evolving techniques, such as proton therapy, showing potential improvements in patient quality of life. Surgical treatment is discussed in the context of achieving a balance between preserving the patient's quality of life and attaining gross total resection, with the extent of surgical resection significantly influencing the survival outcomes. In addition to advancements in cancer vaccine development, this review highlights the evolving landscape of LGG treatment, emphasizing a shift toward personalized and targeted therapies. Ongoing research is essential for refining strategies and enhancing outcomes in the management of LGG.

Kataegis in clinical and molecular subgroups of primary breast cancer

Kataegis is a hypermutation phenomenon characterized by localized clusters of single base pair substitution (SBS) reported in multiple cancer types. Despite a high frequency in breast cancer, large-scale analyses of kataegis patterns and associations with clinicopathological and molecular variables in established breast cancer subgroups are lacking. Therefore, WGS profiled primary breast cancers (n = 791) with associated clinical and molecular data layers, like RNA-sequencing data, were analyzed for kataegis frequency, recurrence, and associations with genomic contexts and functional elements, transcriptional patterns, driver alterations, homologous recombination deficiency (HRD), and prognosis in tumor subgroups defined by ER, PR, and HER2/ERBB2 status. Kataegis frequency was highest in the HER2-positive(p) subgroups, including both ER-negative(n)/positive(p) tumors (ERnHER2p/ERpHER2p). In TNBC, kataegis was neither associated with PAM50 nor TNBC mRNA subtypes nor with distant relapse in chemotherapy-treated patients. In ERpHER2n tumors, kataegis was associated with aggressive characteristics, including PR-negativity, molecular Luminal B subtype, higher mutational burden, higher grade, and expression of proliferation-associated genes. Recurrent kataegis loci frequently targeted regions commonly amplified in ER-positive tumors, while few recurrent loci were observed in TNBC. SBSs in kataegis loci appeared enriched in regions of open chromatin. Kataegis status was not associated with HRD in any subgroup or with distinct transcriptional patterns in unsupervised or supervised analysis. In summary, kataegis is a common hypermutation phenomenon in established breast cancer subgroups, particularly in HER2p subgroups, coinciding with an aggressive tumor phenotype in ERpHER2n disease. In TNBC, the molecular implications and associations of kataegis are less clear, including its prognostic value.

Genomic landscape of diploid and aneuploid microsatellite stable early onset colorectal cancer

Although colorectal cancer (CRC) remains the second leading cause of cancer-related death in the United States, the overall incidence and mortality from the disease have declined in recent decades. In contrast, there has been a steady increase in the incidence of CRC in individuals under 50 years of age. Hereditary syndromes contribute disproportionately to early onset CRC (EOCRC). These include microsatellite instability high (MSI+) tumors arising in patients with Lynch Syndrome. However, most EOCRCs are not associated with familial syndromes or MSI+ genotypes. Comprehensive genomic profiling has provided the basis of improved more personalized treatments for older CRC patients. However, less is known about the basis of sporadic EOCRC. To define the genomic landscape of EOCRC we used DNA content flow sorting to isolate diploid and aneuploid tumor fractions from 21 non-hereditary cases. We then generated whole exome mutational profiles for each case and whole genome copy number, telomere length, and EGFR immunohistochemistry (IHC) analyses on subsets of samples. These results discriminate the molecular features of diploid and aneuploid EOCRC and provide a basis for larger population-based studies and the development of effective strategies to monitor and treat this emerging disease.

Integrated mutational landscape analysis of poorly differentiated high-grade neuroendocrine carcinoma of the uterine cervix

High-grade neuroendocrine cervical cancers (NETc) are exceedingly rare, highly aggressive tumors. We analyzed 64 NETc tumor samples by whole-exome sequencing (WES). Human papillomavirus DNA was detected in 65.6% (42/64) of the tumors. Recurrent mutations were identified in PIK3CA, KMT2D/MLL2, K-RAS, ARID1A, NOTCH2, and RPL10. The top mutated genes included RB1, ARID1A, PTEN, KMT2D/MLL2, and WDFY3, a gene not yet implicated in NETc. Somatic CNV analysis identified two copy number gains (3q27.1 and 19q13.12) and five copy number losses (1p36.21/5q31.3/6p22.2/9q21.11/11p15.5). Also, gene fusions affecting the ACLY-CRHR1 and PVT1-MYC genes were identified in one of the eight samples subjected to RNA sequencing. To resolve evolutionary history, multiregion WES in NETc admixed with adenocarcinoma cells was performed (i.e., mixed-NETc). Phylogenetic analysis of mixed-NETc demonstrated that adenocarcinoma and neuroendocrine elements derive from a common precursor with mutations typical of adenocarcinomas. Over one-third (22/64) of NETc demonstrated a mutator phenotype of C > T at CpG consistent with deficiencies in MBD4, a member of the base excision repair (BER) pathway. Mutations in the PI3K/AMPK pathways were identified in 49/64 samples. We used two patient-derived-xenografts (PDX) (i.e., NET19 and NET21) to evaluate the activity of pan-HER (afatinib), PIK3CA (copanlisib), and ATR (elimusertib) inhibitors, alone and in combination. PDXs harboring alterations in the ERBB2/PI3K/AKT/mTOR/ATR pathway were sensitive to afatinib, copanlisib, and elimusertib (P < 0.001 vs. controls). However, combinations of copanlisib/afatinib and copanlisib/elimusertib were significantly more effective in controlling NETc tumor growth. These findings define the genetic landscape of NETc and suggest that a large subset of these highly lethal malignancies might benefit from existing targeted therapies.

Evolvability of cancer-associated genes under APOBEC3A/B selection

Evolvability is an emergent hallmark of cancer that depends on intra-tumor heterogeneity and genetic variation. Mutations generated by APOBEC3 contribute to genetic variation and tumor evolvability. However, the influence of APOBEC3 on the evolvability of the genome and its differential impact on cancer genes versus non-cancer genes remains unclear. Analyzing over 40,000 human protein-coding transcripts, we identified distinct distribution patterns of APOBEC3A/B TC motifs between cancer and non-cancer genes, suggesting unique associations with cancer. Studying a bat species with numerous APOBEC3 genes, we found distinct motif patterns in orthologs of cancer genes compared to non-cancer genes, as in humans, suggesting APOBEC3 evolution to reduce impacts on the genome rather than the converse. Simulations confirmed that APOBEC3-induced heterogeneity enhances cancer evolution through bimodal patterns of mutations in certain classes of genes. Our results suggest the bimodal distribution of APOBEC-induced mutations can significantly increase cancer heterogeneity.

Somatic mutations of esophageal adenocarcinoma: a comparison between Black and White patients

Esophageal adenocarcinoma is the most common histological subtype of esophageal cancer in Western countries and shows poor prognosis with rapid growth. EAC is characterized by a strong male predominance and racial disparity. EAC is up to fivefold more common among Whites than Blacks, yet Black patients with EAC have poorer survival rates. The racial disparity remains largely unknown, and there is limited knowledge of mutations in EAC regarding racial disparities. We used whole-exome sequencing to show somatic mutation profiles derived from tumor samples from 18 EAC male patients. We identified three molecular subgroups based on the pre-defined esophageal cancer-specific mutational signatures. Group 1 is associated with age and NTHL1 deficiency-related signatures. Group 2 occurs primarily in Black patients and is associated with signatures related to DNA damage from oxidative stress and NTHL1 deficiency-related signatures. Group 3 is associated with defective homologous recombination-based DNA often caused by BRCA mutation in White patients. We observed significantly mutated race related genes (LCE2B in Black, SDR39U1 in White) were (q-value < 0.1). Our findings underscore the possibility of distinct molecular mutation patterns in EAC among different races. Further studies are needed to validate our findings, which could contribute to precision medicine in EAC.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Tumor phylogeography reveals block-shaped spatial heterogeneity and the mode of evolution in Hepatocellular Carcinoma

Solid tumors are complex ecosystems with heterogeneous 3D structures, but the spatial intra-tumor heterogeneity (sITH) at the macroscopic (i.e., whole tumor) level is under-explored. Using a phylogeographic approach, we sequence genomes and transcriptomes from 235 spatially informed sectors across 13 hepatocellular carcinomas (HCC), generating one of the largest datasets for studying sITH. We find that tumor heterogeneity in HCC segregates into spatially variegated blocks with large genotypic and phenotypic differences. By dissecting the transcriptomic heterogeneity, we discover that 30% of patients had a "spatially competing distribution" (SCD), where different spatial blocks have distinct transcriptomic subtypes co-existing within a tumor, capturing the critical transition period in disease progression. Interestingly, the tumor regions with more advanced transcriptomic subtypes (e.g., higher cell cycle) often take clonal dominance with a wider geographic range, rejecting neutral evolution for SCD patients. Extending the statistical tests for detecting natural selection to many non-SCD patients reveal varying levels of selective signal across different tumors, implying that many evolutionary forces including natural selection and geographic isolation can influence the overall pattern of sITH. Taken together, tumor phylogeography unravels a dynamic landscape of sITH, pinpointing important evolutionary and clinical consequences of spatial heterogeneity in cancer.