Characterizing genetic intra-tumor heterogeneity across 2,658 human cancer genomes.

Intra-tumor heterogeneity (ITH) is a mechanism of therapeutic resistance and therefore an important clinical challenge. However, the extent, origin, and drivers of ITH across cancer types are poorly understood. To address this, we extensively characterize ITH across whole-genome sequences of 2,658 cancer samples spanning 38 cancer types. Nearly all informative samples (95.1%) contain evidence of distinct subclonal expansions with frequent branching relationships between subclones. We observe positive selection of subclonal driver mutations across most cancer types and identify cancer type-specific subclonal patterns of driver gene mutations, fusions, structural variants, and copy number alterations as well as dynamic changes in mutational processes between subclonal expansions. Our results underline the importance of ITH and its drivers in tumor evolution and provide a pan-cancer resource of comprehensively annotated subclonal events from whole-genome sequencing data.

Autonomous rhythmic activity in glioma networks drives brain tumour growth.

Diffuse gliomas, particularly glioblastomas, are incurable brain tumours1. They are characterized by networks of interconnected brain tumour cells that communicate via Ca2+ transients2-6. However, the networks' architecture and communication strategy and how these influence tumour biology remain unknown. Here we describe how glioblastoma cell networks include a small, plastic population of highly active glioblastoma cells that display rhythmic Ca2+ oscillations and are particularly connected to others. Their autonomous periodic Ca2+ transients preceded Ca2+ transients of other network-connected cells, activating the frequency-dependent MAPK and NF-κB pathways. Mathematical network analysis revealed that glioblastoma network topology follows scale-free and small-world properties, with periodic tumour cells frequently located in network hubs. This network design enabled resistance against random damage but was vulnerable to losing its key hubs. Targeting of autonomous rhythmic activity by selective physical ablation of periodic tumour cells or by genetic or pharmacological interference with the potassium channel KCa3.1 (also known as IK1, SK4 or KCNN4) strongly compromised global network communication. This led to a marked reduction of tumour cell viability within the entire network, reduced tumour growth in mice and extended animal survival. The dependency of glioblastoma networks on periodic Ca2+ activity generates a vulnerability7 that can be exploited for the development of novel therapies, such as with KCa3.1-inhibiting drugs.

Hypermutation of the inactive X chromosome is a frequent event in cancer.

Mutation is a fundamental process in tumorigenesis. However, the degree to which the rate of somatic mutation varies across the human genome and the mechanistic basis underlying this variation remain to be fully elucidated. Here, we performed a cross-cancer comparison of 402 whole genomes comprising a diverse set of childhood and adult tumors, including both solid and hematopoietic malignancies. Surprisingly, we found that the inactive X chromosome of many female cancer genomes accumulates on average twice and up to four times as many somatic mutations per megabase, as compared to the individual autosomes. Whole-genome sequencing of clonally expanded hematopoietic stem/progenitor cells (HSPCs) from healthy individuals and a premalignant myelodysplastic syndrome (MDS) sample revealed no X chromosome hypermutation. Our data suggest that hypermutation of the inactive X chromosome is an early and frequent feature of tumorigenesis resulting from DNA replication stress in aberrantly proliferating cells.

The evolutionary history of 2,658 cancers.

Cancer develops through a process of somatic evolution1,2. Sequencing data from a single biopsy represent a snapshot of this process that can reveal the timing of specific genomic aberrations and the changing influence of mutational processes3. Here, by whole-genome sequencing analysis of 2,658 cancers as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA)4, we reconstruct the life history and evolution of mutational processes and driver mutation sequences of 38 types of cancer. Early oncogenesis is characterized by mutations in a constrained set of driver genes, and specific copy number gains, such as trisomy 7 in glioblastoma and isochromosome 17q in medulloblastoma. The mutational spectrum changes significantly throughout tumour evolution in 40% of samples. A nearly fourfold diversification of driver genes and increased genomic instability are features of later stages. Copy number alterations often occur in mitotic crises, and lead to simultaneous gains of chromosomal segments. Timing analyses suggest that driver mutations often precede diagnosis by many years, if not decades. Together, these results determine the evolutionary trajectories of cancer, and highlight opportunities for early cancer detection.

Explainable artificial intelligence for omics data: a systematic mapping study.

Researchers increasingly turn to explainable artificial intelligence (XAI) to analyze omics data and gain insights into the underlying biological processes. Yet, given the interdisciplinary nature of the field, many findings have only been shared in their respective research community. An overview of XAI for omics data is needed to highlight promising approaches and help detect common issues. Toward this end, we conducted a systematic mapping study. To identify relevant literature, we queried Scopus, PubMed, Web of Science, BioRxiv, MedRxiv and arXiv. Based on keywording, we developed a coding scheme with 10 facets regarding the studies' AI methods, explainability methods and omics data. Our mapping study resulted in 405 included papers published between 2010 and 2023. The inspected papers analyze DNA-based (mostly genomic), transcriptomic, proteomic or metabolomic data by means of neural networks, tree-based methods, statistical methods and further AI methods. The preferred post-hoc explainability methods are feature relevance (n = 166) and visual explanation (n = 52), while papers using interpretable approaches often resort to the use of transparent models (n = 83) or architecture modifications (n = 72). With many research gaps still apparent for XAI for omics data, we deduced eight research directions and discuss their potential for the field. We also provide exemplary research questions for each direction. Many problems with the adoption of XAI for omics data in clinical practice are yet to be resolved. This systematic mapping study outlines extant research on the topic and provides research directions for researchers and practitioners.

Synergy of retinoic acid and BH3 mimetics in MYC(N)-driven embryonal nervous system tumours.

BACKGROUND: Certain paediatric nervous system malignancies have dismal prognoses. Retinoic acid (RA) is used in neuroblastoma treatment, and preclinical data indicate potential benefit in selected paediatric brain tumour entities. However, limited single-agent efficacy necessitates combination treatment approaches. METHODS: We performed drug sensitivity profiling of 76 clinically relevant drugs in combination with RA in 16 models (including patient-derived tumouroids) of the most common paediatric nervous system tumours. Drug responses were assessed by viability assays, high-content imaging, and apoptosis assays and RA relevant pathways by RNAseq from treated models and patient samples obtained through the precision oncology programme INFORM (n = 2288). Immunoprecipitation detected BCL-2 family interactions, and zebrafish embryo xenografts were used for in vivo efficacy testing. RESULTS: Group 3 medulloblastoma (MBG3) and neuroblastoma models were highly sensitive to RA treatment. RA induced differentiation and regulated apoptotic genes. RNAseq analysis revealed high expression of BCL2L1 in MBG3 and BCL2 in neuroblastomas. Co-treatments with RA and BCL-2/XL inhibitor navitoclax synergistically decreased viability at clinically achievable concentrations. The combination of RA with navitoclax disrupted the binding of BIM to BCL-XL in MBG3 and to BCL-2 in neuroblastoma, inducing apoptosis in vitro and in vivo. CONCLUSIONS: RA treatment primes MBG3 and NB cells for apoptosis, triggered by navitoclax cotreatment.

The chromosomal translocation t(1;6)(p35.3;p25.2), recurrent in chronic lymphocytic leukaemia, leads to RCC1::IRF4 fusion.

The chromosomal translocation t(1;6)(p35.3;p25.2) is a rare but recurrent aberration in chronic lymphocytic leukaemia (CLL). We report molecular characterization of 10 cases and show that this translocation juxtaposes interferon regulatory factor 4 (IRF4) on 6p25 with regulator of chromosome condensation 1 (RCC1) on 1p35. The breakpoints fell within the 5' untranslated regions of both genes, resulting in RCC1::IRF4 fusion transcripts without alterations of the protein-coding sequences. Levels of expression of both RCC1 and IRF4 proteins were not obviously deregulated. The cases showed other mutations typical of CLL and we confirm previously reported skewing towards the IGHV-unmutated subtype. RCC1::IRF4 fusion characterizes a rare subset of CLL.

PerSurge (NOA-30) phase II trial of perampanel treatment around surgery in patients with progressive glioblastoma.

BACKGROUND: Glioblastoma is the most frequent and a particularly malignant primary brain tumor with no efficacy-proven standard therapy for recurrence. It has recently been discovered that excitatory synapses of the AMPA-receptor subtype form between non-malignant brain neurons and tumor cells. This neuron-tumor network connectivity contributed to glioma progression and could be efficiently targeted with the EMA/FDA approved antiepileptic AMPA receptor inhibitor perampanel in preclinical studies. The PerSurge trial was designed to test the clinical potential of perampanel to reduce tumor cell network connectivity and tumor growth with an extended window-of-opportunity concept. METHODS: PerSurge is a phase IIa clinical and translational treatment study around surgical resection of progressive or recurrent glioblastoma. In this multicenter, 2-arm parallel-group, double-blind superiority trial, patients are 1:1 randomized to either receive placebo or perampanel (n = 66 in total). It consists of a treatment and observation period of 60 days per patient, starting 30 days before a planned surgical resection, which itself is not part of the study interventions. Only patients with an expected safe waiting interval are included, and a safety MRI is performed. Tumor cell network connectivity from resected tumor tissue on single cell transcriptome level as well as AI-based assessment of tumor growth dynamics in T2/FLAIR MRI scans before resection will be analyzed as the co-primary endpoints. Secondary endpoints will include further imaging parameters such as pre- and postsurgical contrast enhanced MRI scans, postsurgical T2/FLAIR MRI scans, quality of life, cognitive testing, overall and progression-free survival as well as frequency of epileptic seizures. Further translational research will focus on additional biological aspects of neuron-tumor connectivity. DISCUSSION: This trial is set up to assess first indications of clinical efficacy and tolerability of perampanel in recurrent glioblastoma, a repurposed drug which inhibits neuron-glioma synapses and thereby glioblastoma growth in preclinical models. If perampanel proved to be successful in the clinical setting, it would provide the first evidence that interference with neuron-cancer interactions may indeed lead to a benefit for patients, which would lay the foundation for a larger confirmatory trial in the future. TRIAL REGISTRATION: EU-CT number: 2023-503938-52-00 30.11.2023.

Stroma AReactive Invasion Front Areas (SARIFA) proves prognostic relevance in gastric carcinoma and is based on a tumor-adipocyte interaction indicating an altered immune response.

BACKGROUND: Recently, we presented Stroma AReactive Invasion Front Areas (SARIFA) as a new histomorphologic negative prognostic biomarker in gastric cancer. It is defined as direct contact between tumor cells and fat cells. The aim of this study was to further elucidate the underlying genomic, transcriptional, and immunological mechanisms of the SARIFA phenomenon. METHODS: To address these questions, SARIFA was classified on H&E-stained tissue sections of three cohorts: an external cohort (n = 489, prognostic validation), the TCGA-STAD cohort (n = 194, genomic and transcriptomic analysis), and a local cohort (n = 60, digital spatial profiling (whole transcriptome) and double RNA in situ hybridization/immunostaining of cytokines). RESULTS: SARIFA status proved to be an independent negative prognostic factor for overall survival in an external cohort of gastric carcinomas. In TCGA-STAD cohort, SARIFA is not driven by distinct genomic alterations, whereas the gene expression analyses showed an upregulation of FABP4 in SARIFA-positive tumors. In addition, the transcriptional regulations of white adipocyte differentiation, triglyceride metabolism, and catabolism were upregulated in pathway analyses. In the DSP analysis of SARIFA-positive tumors, FABP4 and the transcriptional regulation of white adipocyte differentiation were upregulated in macrophages. Additionally, a significantly lower expression of the cytokines IL6 and TNFα was observed at the invasion front. CONCLUSIONS: SARIFA proves to be a strong negative prognostic biomarker in advanced gastric cancer, implicating an interaction of tumor cells with tumor-promoting adipocytes with crucial changes in tumor cell metabolism. SARIFA is not driven by tumor genetics but is very likely driven by an altered immune response as a causative mechanism.

DNMT1 Deficiency Impacts on Plasmacytoid Dendritic Cells in Homeostasis and Autoimmune Disease.

Dendritic cells (DCs) are heterogeneous immune regulators involved in autoimmune diseases. Epigenomic mechanisms orchestrating DC development and DC subset diversification remain insufficiently understood but could be important to modulate DC fate for clinical purposes. By combining whole-genome methylation assessment with the analysis of mice expressing reduced DNA methyltransferase 1 levels, we show that distinct DNA methylation levels and patterns are required for the development of plasmacytoid DC and conventional DC subsets. We provide clonal in vivo evidence for DC lineage establishment at the stem cell level, and we show that a high DNA methylation threshold level is essential for Flt3-dependent survival of DC precursors. Importantly, reducing methylation predominantly depletes plasmacytoid DC and alleviates systemic lupus erythematosus in an autoimmunity mouse model. This study shows how DNA methylation regulates the production of DC subsets and provides a potential rationale for targeting autoimmune disease using hypomethylating agents.

Focal structural variants revealed by whole genome sequencing disrupt the histone demethylase KDM4C in B-cell lymphomas.

Histone methylation-modifiers, such as EZH2 and KMT2D, are recurrently altered in B-cell lymphomas. To comprehensively describe the landscape of alterations affecting genes encoding histone methylation-modifiers in lymphomagenesis we investigated whole genome and transcriptome data of 186 mature B-cell lymphomas sequenced in the ICGC MMML-Seq project. Besides confirming common alterations of KMT2D (47% of cases), EZH2 (17%), SETD1B (5%), PRDM9 (4%), KMT2C (4%), and SETD2 (4%), also identified by prior exome or RNA-sequencing studies, we here found recurrent alterations to KDM4C in chromosome 9p24, encoding a histone demethylase. Focal structural variation was the main mechanism of KDM4C alterations, and was independent from 9p24 amplification. We also identified KDM4C alterations in lymphoma cell lines including a focal homozygous deletion in a classical Hodgkin lymphoma cell line. By integrating RNA-sequencing and genome sequencing data we predict that KDM4C structural variants result in loss-offunction. By functional reconstitution studies in cell lines, we provide evidence that KDM4C can act as a tumor suppressor. Thus, we show that identification of structural variants in whole genome sequencing data adds to the comprehensive description of the mutational landscape of lymphomas and, moreover, establish KDM4C as a putative tumor suppressive gene recurrently altered in subsets of B-cell derived lymphomas.

cola: an R/Bioconductor package for consensus partitioning through a general framework.

Classification of high-throughput genomic data is a powerful method to assign samples to subgroups with specific molecular profiles. Consensus partitioning is the most widely applied approach to reveal subgroups by summarizing a consensus classification from a list of individual classifications generated by repeatedly executing clustering on random subsets of the data. It is able to evaluate the stability of the classification. We implemented a new R/Bioconductor package, cola, that provides a general framework for consensus partitioning. With cola, various parameters and methods can be user-defined and easily integrated into different steps of an analysis, e.g., feature selection, sample classification or defining signatures. cola provides a new method named ATC (ability to correlate to other rows) to extract features and recommends spherical k-means clustering (skmeans) for subgroup classification. We show that ATC and skmeans have better performance than other commonly used methods by a comprehensive benchmark on public datasets. We also benchmark key parameters in the consensus partitioning procedure, which helps users to select optimal parameter values. Moreover, cola provides rich functionalities to apply multiple partitioning methods in parallel and directly compare their results, as well as rich visualizations. cola can automate the complete analysis and generates a comprehensive HTML report.

T-cell prolymphocytic leukemia is associated with deregulation of oncogenic microRNAs on transcriptional and epigenetic level.

Deregulation of micro(mi)-RNAs is a common mechanism in tumorigenesis. We investigated the expression of 2083 miRNAs in T-cell prolymphocytic leukemia (T-PLL). Compared to physiologic CD4+ and CD8+ T-cell subsets, 111 miRNAs were differentially expressed in T-PLL. Of these, 33 belonged to miRNA gene clusters linked to cancer. Genomic variants affecting miRNAs were infrequent with the notable exception of copy number aberrations. Remarkably, we found strong upregulation of the miR-200c/-141 cluster in T-PLL to be associated with DNA hypomethylation and active promoter marks. Our findings suggest that copy number aberrations and epigenetic changes could contribute to miRNA deregulation in T-PLL.

Whole exome sequencing identifies novel germline variants of SLC15A4 gene as potentially cancer predisposing in familial colorectal cancer.

About 15% of colorectal cancer (CRC) patients have first-degree relatives affected by the same malignancy. However, for most families the cause of familial aggregation of CRC is unknown. To identify novel high-to-moderate-penetrance germline variants underlying CRC susceptibility, we performed whole exome sequencing (WES) on four CRC cases and two unaffected members of a Polish family without any mutation in known CRC predisposition genes. After WES, we used our in-house developed Familial Cancer Variant Prioritization Pipeline and identified two novel variants in the solute carrier family 15 member 4 (SLC15A4) gene. The heterozygous missense variant, p. Y444C, was predicted to affect the phylogenetically conserved PTR2/POT domain and to have a deleterious effect on the function of the encoded peptide/histidine transporter. The other variant was located in the upstream region of the same gene (GRCh37.p13, 12_129308531_C_T; 43 bp upstream of transcription start site, ENST00000266771.5) and it was annotated to affect the promoter region of SLC15A4 as well as binding sites of 17 different transcription factors. Our findings of two distinct variants in the same gene may indicate a synergistic up-regulation of SLC15A4 as the underlying genetic cause and implicate this gene for the first time in genetic inheritance of familial CRC.

Comprehensive genomic analysis of refractory multiple myeloma reveals a complex mutational landscape associated with drug resistance and novel therapeutic vulnerabilities.

The outcomes of patients with multiple myeloma (MM) refractory to immunomodulatory agents (IMiDs) and proteasome inhibitors (PIs) remain poor. In this study, we performed whole genome and transcriptome sequencing of 39 heavily pretreated relapsed/refractory MM (RRMM) patients to identify mechanisms of resistance and potential therapeutic targets. We observed a high mutational load and indications of increased genomic instability. Recurrently mutated genes in RRMM, which had not been previously reported or only observed at a lower frequency in newly diagnosed MM, included NRAS, BRAF, TP53, SLC4A7, MLLT4, EWSR1, HCFC2, and COPS3. We found multiple genomic regions with bi-allelic events affecting tumor suppressor genes and demonstrated a significant adverse impact of bi-allelic TP53 alterations on survival. With regard to potentially resistance conferring mutations, recurrently mutated gene networks included genes with relevance for PI and IMiD activity; the latter particularly affecting members of the Cereblon and the COP9 signalosome complex. We observed a major impact of signatures associated with exposure to melphalan or impaired DNA double-strand break homologous recombination repair in RRMM. The latter coincided with mutations in genes associated with PARP inhibitor sensitivity in 49% of RRMM patients; a finding with potential therapeutic implications. In conclusion, this comprehensive genomic characterization revealed a complex mutational and structural landscape in RRMM and highlights potential implications for therapeutic strategies.

iTReX: Interactive exploration of mono- and combination therapy dose response profiling data.

High throughput screening methods, measuring the sensitivity and resistance of tumor cells to drug treatments have been rapidly evolving. Not only do these screens allow correlating response profiles to tumor genomic features for developing novel predictors of treatment response, but they can also add evidence for therapy decision making in precision oncology. Recent analysis methods developed for either assessing single agents or combination drug efficacies enable quantification of dose-response curves with restricted symmetric fit settings. Here, we introduce iTReX, a user-friendly and interactive Shiny/R application, for both the analysis of mono- and combination therapy responses. The application features an extended version of the drug sensitivity score (DSS) based on the integral of an advanced five-parameter dose-response curve model and a differential DSS for combination therapy profiling. Additionally, iTReX includes modules that visualize drug target interaction networks and support the detection of matches between top therapy hits and the sample omics features to enable the identification of druggable targets and biomarkers. iTReX enables the analysis of various quantitative drug or therapy response readouts (e.g. luminescence, fluorescence microscopy) and multiple treatment strategies (drug treatments, radiation). Using iTReX we validate a cost-effective drug combination screening approach and reveal the application's ability to identify potential sample-specific biomarkers based on drug target interaction networks. The iTReX web application is accessible at https://itrex.kitz-heidelberg.de.

Whole-Exome Sequencing Identifies a Novel Germline Variant in PTK7 Gene in Familial Colorectal Cancer.

Colorectal cancer (CRC) is the third most frequently diagnosed malignancy worldwide. Only 5% of all CRC cases are due to germline mutations in known predisposition genes, and the remaining genetic burden still has to be discovered. In this study, we performed whole-exome sequencing on six members of a Polish family diagnosed with CRC and identified a novel germline variant in the protein tyrosine kinase 7 (inactive) gene (PTK7, ENST00000230419, V354M). Targeted screening of the variant in 1705 familial CRC cases and 1674 healthy elderly individuals identified the variant in an additional familial CRC case. Introduction of this variant in HT-29 cells resulted in increased cell proliferation, migration, and invasion; it also caused down-regulation of CREB, p21 and p53 mRNA and protein levels, and increased AKT phosphorylation. These changes indicated inhibition of apoptosis pathways and activation of AKT signaling. Our study confirmed the oncogenic function of PTK7 and supported its role in genetic predisposition of familial CRC.

Analysis of mutational signatures with yet another package for signature analysis.

Different mutational processes leave characteristic patterns of somatic mutations in the genome that can be identified as mutational signatures. Determining the contributions of mutational signatures to cancer genomes allows not only to reconstruct the etiology of somatic mutations, but can also be used for improved tumor classification and support therapeutic decisions. We here present the R package yet another package for signature analysis (YAPSA) to deconvolute the contributions of mutational signatures to tumor genomes. YAPSA provides in-built collections from the COSMIC and PCAWG SNV signature sets as well as the PCAWG Indel signatures and employs signature-specific cutoffs to increase sensitivity and specificity. Furthermore, YAPSA allows to determine 95% confidence intervals for signature exposures, to perform constrained stratified signature analyses to obtain enrichment and depletion patterns of the identified signatures and, when applied to whole exome sequencing data, to correct for the triplet content of individual target capture kits. With this functionality, YAPSA has proved to be a valuable tool for analysis of mutational signatures in molecular tumor boards in a precision oncology context. YAPSA is available at R/Bioconductor (http://bioconductor.org/packages/3.12/bioc/html/YAPSA.html).

The mutational landscape of Burkitt-like lymphoma with 11q aberration is distinct from that of Burkitt lymphoma.

The new recently described provisional lymphoma category Burkitt-like lymphoma with 11q aberration comprises cases similar to Burkitt lymphoma (BL) on morphological, immunophenotypic and gene-expression levels but lacking the IG-MYC translocation. They are characterized by a peculiar imbalance pattern on chromosome 11, but the landscape of mutations is not yet described. Thus, we investigated 15 MYC-negative Burkitt-like lymphoma with 11q aberration (mnBLL,11q,) cases by copy-number analysis and whole-exome sequencing. We refined the regions of 11q imbalance and identified the INO80 complex-associated gene NFRKB as a positional candidate in 11q24.3. Next to recurrent gains in 12q13.11-q24.32 and 7q34-qter as well as losses in 13q32.3-q34, we identified 47 genes recurrently affected by protein-changing mutations (each ≥3 of 15 cases). Strikingly, we did not detect recurrent mutations in genes of the ID3-TCF3 axis or the SWI/SNF complex that are frequently altered in BL, or in genes frequently mutated in germinal center-derived B-cell lymphomas like KMT2D or CREBBP An exception is GNA13, which was mutated in 7 of 15 cases. We conclude that the genomic landscape of mnBLL,11q, differs from that of BL both at the chromosomal and mutational levels. Our findings implicate that mnBLL,11q, is a lymphoma category distinct from BL at the molecular level.