OncodriveCLUST: exploiting the positional clustering of somatic mutations to identify cancer genes

Advances in predicting breast cancer driver mutations: Tools for precision oncology (Review)

In the modern era of medicine, prognosis and treatment, options for a number of cancer types including breast cancer have been improved by the identification of cancer‑specific biomarkers. The availability of high‑throughput sequencing and analysis platforms, the growth of publicly available cancer databases and molecular and histological profiling facilitate the development of new drugs through a precision medicine approach. However, only a fraction of patients with breast cancer with few actionable mutations typically benefit from the precision medicine approach. In the present review, the current development in breast cancer driver gene identification, actionable breast cancer mutations, as well as the available therapeutic options, challenges and applications of breast precision oncology are systematically described. Breast cancer driver mutation‑based precision oncology helps to screen key drivers involved in disease development and progression, drug sensitivity and the genes responsible for drug resistance. Advances in precision oncology will provide more targeted therapeutic options for patients with breast cancer, improving disease‑free survival and potentially leading to significant successes in breast cancer treatment in the near future. Identification of driver mutations has allowed new targeted therapeutic approaches in combination with standard chemo‑ and immunotherapies in breast cancer. Developing new driver mutation identification strategies will help to define new therapeutic targets and improve the overall and disease‑free survival of patients with breast cancer through efficient medicine.

Protein structural context of cancer mutations reveals molecular mechanisms and candidate driver genes

Advances in protein structure determination and modeling allow us to study the structural context of human genetic variants on an unprecedented scale. Here, we analyze millions of cancer-associated missense mutations based on their structural locations and predicted perturbative effects. By considering the collective properties of mutations at the level of individual proteins, we identify distinct patterns associated with tumor suppressors and oncogenes. Tumor suppressors are enriched in structurally damaging mutations, consistent with loss-of-function mechanisms, while oncogene mutations tend to be structurally mild, reflecting selection for gain-of-function driver mutations and against loss-of-function mutations. Although oncogenes are difficult to distinguish from genes with no role in cancer using only structural damage, we find that the three-dimensional clustering of mutations is highly predictive. These observations allow us to identify candidate driver genes and speculate about their molecular roles, which we expect will have general utility in the analysis of cancer sequencing data.

Proteogenomic analysis of air-pollution-associated lung cancer reveals prevention and therapeutic opportunities

Air pollution significantly impacts lung cancer progression, but there is a lack of a comprehensive molecular characterization of clinical samples associated with air pollution. Here, we performed a proteogenomic analysis of lung adenocarcinoma (LUAD) in 169 female never-smokers from the Xuanwei area (XWLC cohort), where coal smoke is the primary contributor to the high lung cancer incidence. Genomic mutation analysis revealed XWLC as a distinct subtype of LUAD separate from cases associated with smoking or endogenous factors. Mutational signature analysis suggested that Benzo[a]pyrene (BaP) is the major risk factor in XWLC. The BaP-induced mutation hotspot, EGFR-G719X, was present in 20% of XWLC which endowed XWLC with elevated MAPK pathway activations and worse outcomes compared to common EGFR mutations. Multi-omics clustering of XWLC identified four clinically relevant subtypes. These subgroups exhibited distinct features in biological processes, genetic alterations, metabolism demands, immune landscape, and radiomic features. Finally, MAD1 and TPRN were identified as novel potential therapeutic targets in XWLC. Our study provides a valuable resource for researchers and clinicians to explore prevention and treatment strategies for air-pollution-associated lung cancers.

SSCI: Self-Supervised Deep Learning Improves Network Structure for Cancer Driver Gene Identification

The pathogenesis of cancer is complex, involving abnormalities in some genes in organisms. Accurately identifying cancer genes is crucial for the early detection of cancer and personalized treatment, among other applications. Recent studies have used graph deep learning methods to identify cancer driver genes based on biological networks. However, incompleteness and the noise of the networks will weaken the performance of models. To address this, we propose a cancer driver gene identification method based on self-supervision for graph convolutional networks, which can efficiently enhance the structure of the network and further improve predictive accuracy. The reliability of SSCI is verified by the area under the receiver operating characteristic curves (AUROC), the area under the precision-recall curves (AUPRC), and the F1 score, with respective values of 0.966, 0.964, and 0.913. The results show that our method can identify cancer driver genes with strong discriminative power and biological interpretability.

Identifying cooperating cancer driver genes in individual patients through hypergraph random walk

OBJECTIVE

Identifying cancer driver genes, especially rare or patient-specific cancer driver genes, is a primary goal in cancer therapy. Although researchers have proposed some methods to tackle this problem, these methods mostly identify cancer driver genes at single gene level, overlooking the cooperative relationship among cancer driver genes. Identifying cooperating cancer driver genes in individual patients is pivotal for understanding cancer etiology and advancing the development of personalized therapies.

METHODS

Here, we propose a novel Personalized Cooperating cancer Driver Genes (PCoDG) method by using hypergraph random walk to identify the cancer driver genes that cooperatively drive individual patient cancer progression. By leveraging the powerful ability of hypergraph in representing multi-way relationships, PCoDG first employs the personalized hypergraph to depict the complex interactions among mutated genes and differentially expressed genes of an individual patient. Then, a hypergraph random walk algorithm based on hyperedge similarity is utilized to calculate the importance scores of mutated genes, integrating these scores with signaling pathway data to identify the cooperating cancer driver genes in individual patients.

RESULTS

The experimental results on three TCGA cancer datasets (i.e., BRCA, LUAD, and COADREAD) demonstrate the effectiveness of PCoDG in identifying personalized cooperating cancer driver genes. These genes identified by PCoDG not only offer valuable insights into patient stratification correlating with clinical outcomes, but also provide an useful reference resource for tailoring personalized treatments.

CONCLUSION

We propose a novel method that can effectively identify cooperating cancer driver genes for individual patients, thereby deepening our understanding of the cooperative relationship among personalized cancer driver genes and advancing the development of precision oncology.

Drivergene.net: A Cytoscape app for the identification of driver nodes of large-scale complex networks and case studies in discovery of drug target genes

There are no tools to identify driver nodes of large-scale networks in approach of competition-based controllability. This study proposed a novel method for this computation of large-scale networks. It implemented the method in a new Cytoscape plug-in app called Drivergene.net. Experiments of the software on large-scale biomolecular networks have shown outstanding speed and computing power. Interestingly, 86.67% of the top 10 driver nodes found on these networks are anticancer drug target genes that reside mostly at the innermost K-cores of the networks. Finally, compared method with those of five other researchers and confirmed that the proposed method outperforms the other methods on identification of anticancer drug target genes. Taken together, Drivergene.net is a reliable tool that efficiently detects not only drug target genes from biomolecular networks but also driver nodes of large-scale complex networks. Drivergene.net with a user manual and example datasets are available https://github.com/tinhpd/Drivergene.git.

Variant Impact Predictor database (VIPdb), version 2: trends from three decades of genetic variant impact predictors

BACKGROUND

Variant interpretation is essential for identifying patients' disease-causing genetic variants amongst the millions detected in their genomes. Hundreds of Variant Impact Predictors (VIPs), also known as Variant Effect Predictors (VEPs), have been developed for this purpose, with a variety of methodologies and goals. To facilitate the exploration of available VIP options, we have created the Variant Impact Predictor database (VIPdb).

RESULTS

The Variant Impact Predictor database (VIPdb) version 2 presents a collection of VIPs developed over the past three decades, summarizing their characteristics, ClinGen calibrated scores, CAGI assessment results, publication details, access information, and citation patterns. We previously summarized 217 VIPs and their features in VIPdb in 2019. Building upon this foundation, we identified and categorized an additional 190 VIPs, resulting in a total of 407 VIPs in VIPdb version 2. The majority of the VIPs have the capacity to predict the impacts of single nucleotide variants and nonsynonymous variants. More VIPs tailored to predict the impacts of insertions and deletions have been developed since the 2010s. In contrast, relatively few VIPs are dedicated to the prediction of splicing, structural, synonymous, and regulatory variants. The increasing rate of citations to VIPs reflects the ongoing growth in their use, and the evolving trends in citations reveal development in the field and individual methods.

CONCLUSIONS

VIPdb version 2 summarizes 407 VIPs and their features, potentially facilitating VIP exploration for various variant interpretation applications. VIPdb is available at  https://genomeinterpretation.org/vipdb.

Somatic mutational landscape across Indian breast cancer cases by whole exome sequencing

Breast cancer (BC) has emerged as the most common malignancy among females. The genomic profile of BC is diverse in nature and complex due to heterogeneity among various geographically different ethnic groups. The primary objective of this study was to carry out a comprehensive mutational analysis of Indian BC cases by performing whole exome sequencing. The cohort included patients with a median age of 48 years. TTN, TP53, MUC16, SYNE1, and OBSCN were the frequently altered genes found in our cohort. The PIK3CA and KLC3 genes are driver genes implicated in various cellular functions and cargo transportation through microtubules, respectively. Except for CCDC168 and PIK3CA, several gene pairings were found to be significantly linked with co-occurrence. Irrespective of their hormonal receptor status, RTK/RAS was observed with frequently altered signaling pathways. Further analysis of the mutational signature revealed that SBS13, SBS6, and SBS29 were mainly observed in our cohort. This study supplements the discovery of diagnostic biomarkers and provides new therapeutic options for the improved management of BC.

Multiview representation learning for identification of novel cancer genes and their causative biological mechanisms

Tumorigenesis arises from the dysfunction of cancer genes, leading to uncontrolled cell proliferation through various mechanisms. Establishing a complete cancer gene catalogue will make precision oncology possible. Although existing methods based on graph neural networks (GNN) are effective in identifying cancer genes, they fall short in effectively integrating data from multiple views and interpreting predictive outcomes. To address these shortcomings, an interpretable representation learning framework IMVRL-GCN is proposed to capture both shared and specific representations from multiview data, offering significant insights into the identification of cancer genes. Experimental results demonstrate that IMVRL-GCN outperforms state-of-the-art cancer gene identification methods and several baselines. Furthermore, IMVRL-GCN is employed to identify a total of 74 high-confidence novel cancer genes, and multiview data analysis highlights the pivotal roles of shared, mutation-specific, and structure-specific representations in discriminating distinctive cancer genes. Exploration of the mechanisms behind their discriminative capabilities suggests that shared representations are strongly associated with gene functions, while mutation-specific and structure-specific representations are linked to mutagenic propensity and functional synergy, respectively. Finally, our in-depth analyses of these candidates suggest potential insights for individualized treatments: afatinib could counteract many mutation-driven risks, and targeting interactions with cancer gene SRC is a reasonable strategy to mitigate interaction-induced risks for NR3C1, RXRA, HNF4A, and SP1.

A Novel Affordable and Reliable Framework for Accurate Detection and Comprehensive Analysis of Somatic Mutations in Cancer

Accurate detection and analysis of somatic variants in cancer involve multiple third-party tools with complex dependencies and configurations, leading to laborious, error-prone, and time-consuming data conversions. This approach lacks accuracy, reproducibility, and portability, limiting clinical application. Musta was developed to address these issues as an end-to-end pipeline for detecting, classifying, and interpreting cancer mutations. Musta is based on a Python command-line tool designed to manage tumor-normal samples for precise somatic mutation analysis. The core is a Snakemake-based workflow that covers all key cancer genomics steps, including variant calling, mutational signature deconvolution, variant annotation, driver gene detection, pathway analysis, and tumor heterogeneity estimation. Musta is easy to install on any system via Docker, with a Makefile handling installation, configuration, and execution, allowing for full or partial pipeline runs. Musta has been validated at the CRS4-NGS Core facility and tested on large datasets from The Cancer Genome Atlas and the Beijing Institute of Genomics. Musta has proven robust and flexible for somatic variant analysis in cancer. It is user-friendly, requiring no specialized programming skills, and enables data processing with a single command line. Its reproducibility ensures consistent results across users following the same protocol.

Genomic alteration discordance in the paired primary-recurrent ovarian cancers: based on the comprehensive genomic profiling (CGP) analysis

PURPOSE

Ovarian cancer (OC) is characterized by a high recurrence rate, and homologous recombination deficiency (HRD) is an important biomarker in the clinical management of OC. We investigated the differences in clinical genomic profiles between the primary and platinum-sensitive recurrent OC (PSROC), focusing on HRD status.

MATERIALS AND METHODS

A total of 40 formalin-fixed paraffin-embedded (FFPE) tissues of primary tumors and their first platinum-sensitive recurrence from 20 OC patients were collected, and comprehensive genomic profiling (CGP) analysis of FoundationOne®CDx (F1CDx) was applied to explore the genetic (dis)similarities of the primary and recurrent tumors.

RESULTS

By comparing between paired samples, we found that genomic loss of heterozygosity (gLOH) score had a high intra-patient correlation (r2 = 0.79) and that short variants (including TP53, BRCA1/2 and NOTCH1 mutations), tumor mutational burden (TMB) and microsatellite stability status remained stable. The frequency of (likely) pathological BRCA1/2 mutations was 30% (12/40) in all samples positively correlated with gLOH scores, but the proportion of gLOH-high status (score > 16%) was 50% (10/20) and 55% (11/20) in the primary and recurrent samples, respectively. An additional 20% (4/20) of patients needed attention, a quarter of which carried the pathological BRCA1 mutation but had a gLOH-low status (gLOH < 16%), and three-quarters had different gLOH status in primary-recurrent pairs. Furthermore, we observed the PSROC samples had higher gLOH scores (16.1 ± 9.24 vs. 19.4 ± 11.1, p = 0.007), more CNVs (36.1% vs. 15.1% of discordant genomic alternations), and significant enrichment of altered genes in TGF-beta signaling and Hippo signaling pathways (p < 0.05 for all) than their paired primaries. Lastly, mutational signature and oncodrive gene analyses showed that the computed mutational signature similarity in the primary and recurrent tumors were best matched the COSMI 3 signature (Aetiology of HRD) and had consistent candidate cancer driver genes of MSH2, NOTCH1 and MSH6.

CONCLUSION

The high genetic concordance of the short variants remains stable along OC recurrence. However, the results reveal significantly higher gLOH scores in the recurrent setting than in paired primaries, supporting further clinically instantaneity HRD assay strategy.

Stratification of glioma based on stemness scores in bulk and single-cell transcriptomes

BACKGROUND

Brain tumours are known to have a high mortality and morbidity rate due to their localised and frequent invasive growth. The concept that glioma resistance could originate from the dissimilarity in the vulnerability of clonogenic glial stem cells to chemotherapeutic drugs and radiation has driven the scientific community to reexamine the comprehension of glioma growth and strategies that target these cells or modify their stemness.

METHODS

Based on the enrichment scores of 12 stemness signatures, we identified glioma subtypes in both tumour bulks and single cells by clustering analysis. Furthermore, we comprehensively compared molecular and clinical features among the glioma subtypes.

RESULTS

Consistently, in seven different datasets, hierarchical clustering uncovered three subtypes of glioma, termed Stem-H, Stem-M, and Stem-L, with high, medium, and low stemness signatures, respectively. Stem-H and Stem-L exhibited the most unfavorable and favourable overall and disease-free survival, respectively. Stem-H showed the highest enrichment scores of the EMT, invasion, proliferation, differentiation, and metastasis processes signatures, while Stem-L displayed the lowest. Stem-H harboured a greater proportion of late-stage tumours compared to Stem-L. Moreover, Stem-H manifested higher tumour mutation burden, DNA damage repair and cell cycle activity, intratumour heterogeneity, and a more frequent incidence of TP53 and EGFR mutations than Stem-L. In contrast, Stem-L had higher O6-Methylguanine-DNA Methyltransferase (MGMT) methylation levels.

CONCLUSION

The classification of glioma based on stemness may offer new insights into the biology of the tumour, as well as more accurate clinical management of the disease.

CancerGATE: Prediction of cancer-driver genes using graph attention autoencoders

Discovery of the cancer type specific-driver genes is important for understanding the molecular mechanisms of each cancer type and for providing proper treatment. Recently, graph deep learning methods became widely used in finding cancer-driver genes. However, previous methods had limited performance in individual cancer types due to a small number of cancer-driver genes used in training and biases toward the cancer-driver genes used in training the models. Here, we introduce a novel pipeline, CancerGATE that predicts the cancer-driver genes using graph attention autoencoder (GATE) to learn in a self-supervised manner and can be applied to each of the cancer types. CancerGATE utilizes biological network topology and multi-omics data from 15 types of cancer of 20,079 samples from the cancer genome atlas (TCGA). Attention coefficients calculated in the model are used to prioritize cancer-driver genes by comparing coefficients of cancer and normal contexts. CancerGATE shows a higher AUPRC with a difference ranging from 1.5 % to 36.5 % compared to the previous graph deep learning models in each cancer type. We also show that CancerGATE is free from the bias toward cancer-driver genes used in training, revealing mechanisms of the cancer-driver genes in specific cancer types. Finally, we propose novel cancer-driver gene candidates that could be therapeutic targets for specific cancer types.

Somatic mutations that affect early genetic progression and immune microenvironment in gastric carcinoma

Gastric carcinoma (GC) is a high heterogeneity and malignant tumor with a poor prognosis. The current implementation of immunotherapy in GC is limited due to the insufficient exploration of immune-related mutations and speculated early mutation events. Therefore, we performed whole-exome sequencing on 40 patients with GC to explore their genetic characteristics, shedding light on the order of genetic events, somatic mutations impacting the immune microenvironment, and potential biomarkers for immunotherapy. Regarding genetic events, TP53 disruptions were identified as frequent and early events in GC progression, often occurring alongside other gene mutations. The mutations occurring in GANS, SMAD4, and POLE were early independent events. Patients harboring CSMD3, FAT4, FLG, KMT2C, LRP1B, MUC5B, MUC16, PLEC, RNF43, SYNE1, TP53, TTN, XIRP2, and ZFHX4 mutations tended to have decreased B cells, T cells, macrophage, neutrophil, and dendritic cells infiltration, except for the ARID1A gene mutations. We also found patients with microsatellite instability-high tumors had higher homologous recombination deficiency (HRD) scores. HRD showed a positive correlation with tumor mutational burden, which might serve as indirect evidence supporting the potential of HRD as a biomarker for GC. These findings highlighted GC's high heterogeneity and complexity and provided valuable insights into the somatic mutations that affect early genetic progression and immune microenvironment.

What can we learn about acid-base transporters in cancer from studying somatic mutations in their genes?

Acidosis is a chemical signature of the tumour microenvironment that challenges intracellular pH homeostasis. The orchestrated activity of acid-base transporters of the solute-linked carrier (SLC) family is critical for removing the end-products of fermentative metabolism (lactate/H+) and maintaining a favourably alkaline cytoplasm. Given the critical role of pH homeostasis in enabling cellular activities, mutations in relevant SLC genes may impact the oncogenic process, emerging as negatively or positively selected, or as driver or passenger mutations. To address this, we performed a pan-cancer analysis of The Cancer Genome Atlas simple nucleotide variation data for acid/base-transporting SLCs (ABT-SLCs). Somatic mutation patterns of monocarboxylate transporters (MCTs) were consistent with their proposed essentiality in facilitating lactate/H+ efflux. Among all cancers, tumours of uterine corpus endometrial cancer carried more ABT-SLC somatic mutations than expected from median tumour mutation burden. Among these, somatic mutations in SLC4A3 had features consistent with meaningful consequences on cellular fitness. Definitive evidence for ABT-SLCs as 'cancer essential' or 'driver genes' will have to consider microenvironmental context in genomic sequencing because bulk approaches are insensitive to pH heterogeneity within tumours. Moreover, genomic analyses must be validated with phenotypic outcomes (i.e. SLC-carried flux) to appreciate the opportunities for targeting acid-base transport in cancers.

Personalized Driver Gene Prediction Using Graph Convolutional Networks with Conditional Random Fields

Cancer is a complex and evolutionary disease mainly driven by the accumulation of genetic variations in genes. Identifying cancer driver genes is important. However, most related studies have focused on the population level. Cancer is a disease with high heterogeneity. Thus, the discovery of driver genes at the individual level is becoming more valuable but is a great challenge. Although there have been some computational methods proposed to tackle this challenge, few can cover all patient samples well, and there is still room for performance improvement. In this study, to identify individual-level driver genes more efficiently, we propose the PDGCN method. PDGCN integrates multiple types of data features, including mutation, expression, methylation, copy number data, and system-level gene features, along with network structural features extracted using Node2vec in order to construct a sample-gene interaction network. Prediction is performed using a graphical convolutional neural network model with a conditional random field layer, which is able to better combine the network structural features with biological attribute features. Experiments on the ACC (Adrenocortical Cancer) and KICH (Kidney Chromophobe) datasets from TCGA (The Cancer Genome Atlas) demonstrated that the method performs better compared to other similar methods. It can identify not only frequently mutated driver genes, but also rare candidate driver genes and novel biomarker genes. The results of the survival and enrichment analyses of these detected genes demonstrate that the method can identify important driver genes at the individual level.

An inducible Cd79b mutation confers ibrutinib sensitivity in mouse models of Myd88-driven diffuse large B-cell lymphoma

ABSTRACT

Diffuse large B-cell lymphoma (DLBCL) is the most common aggressive lymphoma and constitutes a highly heterogenous disease. Recent comprehensive genomic profiling revealed the identity of numerous molecularly defined DLBCL subtypes, including a cluster which is characterized by recurrent aberrations in MYD88, CD79B, and BCL2, as well as various lesions promoting a block in plasma cell differentiation, including PRDM1, TBL1XR1, and SPIB. Here, we generated a series of autochthonous mouse models to mimic this DLBCL cluster and specifically focused on the impact of Cd79b mutations in this setting. We show that canonical Cd79b immunoreceptor tyrosine-based activation motif (ITAM) mutations do not accelerate Myd88- and BCL2-driven lymphomagenesis. Cd79b-mutant murine DLBCL were enriched for IgM surface expression, reminiscent of their human counterparts. Moreover, Cd79b-mutant lymphomas displayed a robust formation of cytoplasmic signaling complexes involving MYD88, CD79B, MALT1, and BTK. These complexes were disrupted upon pharmacological BTK inhibition. The BTK inhibitor-mediated disruption of these signaling complexes translated into a selective ibrutinib sensitivity of lymphomas harboring combined Cd79b and Myd88 mutations. Altogether, this in-depth cross-species comparison provides a framework for the development of molecularly targeted therapeutic intervention strategies in DLBCL.

Association of tumour mutation burden with prognosis and its clinical significance in stage III gastric cancer

Introduction

To explore the correlation between the tumour mutation burden (TMB) and prognosis and its clinical significance among patients with stage III gastric cancer (GC).

Methods

Patients with stage III GC were divided into a high TMB and low TMB group in both a study cohort of 38 patients and the Cancer Genome Atlas (TCGA) cohort of 173 patients. In the study cohort, next-generation sequencing was used to detect mutated GC genes and obtain TMB data. In the TCGA cohort, gene set enrichment analysis was performed, and the relationship between TMB, prognosis and clinicopathologic factors was analysed. Western blot and quantitative real-time polymerase chain reaction were used to detect the expression levels of both proteins and genes. Cell viability was measured using methyl thiazolyl tetrazolium and transwell cell assays.

Results

Patients in the high TMB group had better overall survival (OS) rates than patients in the low TMB group for both cohorts and TMB was associated with age, mutation signature 1 and mutation signature 17. The Cox regression analysis revealed that age, not TMB, was an independent prognosis factor. Furthermore, genes with high-frequency mutations were significantly enriched in the RTK-RAS and Notch signalling pathways. The activation of these pathways was lower in the high TMB compared with the low TMB group, and the proliferation and migration abilities of GC cells showed a similar pattern in both TMB groups.

Conclusion

Patients in the high TMB group had better OS rates than patients in the low TMB group. Genes with high-frequency mutations were significantly enriched in the RTK-RAS and Notch pathways. Hence, TMB could serve as a prognosis biomarker with potential clinical significance.

Singleton mutations in large-scale cancer genome studies: uncovering the tail of cancer genome

Singleton or low-frequency driver mutations are challenging to identify. We present a domain driver mutation estimator (DOME) to identify rare candidate driver mutations. DOME analyzes positions analogous to known statistical hotspots and resistant mutations in combination with their functional and biochemical residue context as determined by protein structures and somatic mutation propensity within conserved PFAM domains, integrating the CADD scoring scheme. Benchmarked against seven other tools, DOME exhibited superior or comparable accuracy compared to all evaluated tools in the prediction of functional cancer drivers, with the exception of one tool. DOME identified a unique set of 32 917 high-confidence predicted driver mutations from the analysis of whole proteome missense variants within domain boundaries across 1331 genes, including 1192 noncancer gene census genes, emphasizing its unique place in cancer genome analysis. Additionally, analysis of 8799 TCGA (The Cancer Genome Atlas) and in-house tumor samples revealed 847 potential driver mutations, with mutations in tyrosine kinase members forming the dominant burden, underscoring its higher significance in cancer. Overall, DOME complements current approaches for identifying novel, low-frequency drivers and resistant mutations in personalized therapy.

Genomic characteristics and immune landscape of super multiple primary lung cancer

BACKGROUND

In recent years, a growing number of patients with multiple primary lung cancer (MPLC) are being diagnosed, and a subset of these patients is found to have a large number of lesions at the time of diagnosis, which are referred to as 'super MPLC'.

METHODS

Here, we perform whole exome sequencing (WES) and immunohistochemistry (IHC) analysis of PD-L1 and CD8 on 212 tumor samples from 42 patients with super MPLC.

FINDINGS

We report the genomic alteration landscape of super MPLC. EGFR, RBM10 and TP53 mutation and TERT amplification are important molecular events in the evolution of super MPLC. We propose the conception of early intrapulmonary metastasis, which exhibits different clinical features from conventional metastasis. The IHC analyses of PD-L1 and CD8 reveal a less inflamed microenvironment of super MPLC than that of traditional non-small cell lung cancer (NSCLC). We identify the potentially susceptible germline mutations for super MPLC.

INTERPRETATION

Our study depicts the genomic characteristics and immune landscape, providing insights into the pathogenesis and possible therapeutic guidance of super MPLC.

FUNDING

A full list of funding bodies that supported this study can be found in the Acknowledgements section.

Integrating Clinical Cancer and PTM Proteomics Data Identifies a Mechanism of ACK1 Kinase Activation

Beyond the most common oncogenes activated by mutation (mut-drivers), there likely exists a variety of low-frequency mut-drivers, each of which is a possible frontier for targeted therapy. To identify new and understudied mut-drivers, we developed a machine learning (ML) model that integrates curated clinical cancer data and posttranslational modification (PTM) proteomics databases. We applied the approach to 62,746 patient cancers spanning 84 cancer types and predicted 3,964 oncogenic mutations across 1,148 genes, many of which disrupt PTMs of known and unknown function. The list of putative mut-drivers includes established drivers and others with poorly understood roles in cancer. This ML model is available as a web application. As a case study, we focused the approach on nonreceptor tyrosine kinases (NRTK) and found a recurrent mutation in activated CDC42 kinase-1 (ACK1) that disrupts the Mig6 homology region (MHR) and ubiquitin-association (UBA) domains on the ACK1 C-terminus. By studying these domains in cultured cells, we found that disruption of the MHR domain helps activate the kinase while disruption of the UBA increases kinase stability by blocking its lysosomal degradation. This ACK1 mutation is analogous to lymphoma-associated mutations in its sister kinase, TNK1, which also disrupt a C-terminal inhibitory motif and UBA domain. This study establishes a mut-driver discovery tool for the research community and identifies a mechanism of ACK1 hyperactivation shared among ACK family kinases.

IMPLICATIONS

This research identifies a potentially targetable activating mutation in ACK1 and other possible oncogenic mutations, including PTM-disrupting mutations, for further study.

Adaptation of a mutual exclusivity framework to identify driver mutations within oncogenic pathways

Distinguishing genomic alterations in cancer-associated genes that have functional impact on tumor growth and disease progression from the ones that are passengers and confer no fitness advantage have important clinical implications. Evidence-based methods for nominating drivers are limited by existing knowledge on the oncogenic effects and therapeutic benefits of specific variants from clinical trials or experimental settings. As clinical sequencing becomes a mainstay of patient care, applying computational methods to mine the rapidly growing clinical genomic data holds promise in uncovering functional candidates beyond the existing knowledge base and expanding the patient population that could potentially benefit from genetically targeted therapies. We propose a statistical and computational method (MAGPIE) that builds on a likelihood approach leveraging the mutual exclusivity pattern within an oncogenic pathway for identifying probabilistically both the specific genes within a pathway and the individual mutations within such genes that are truly the drivers. Alterations in a cancer-associated gene are assumed to be a mixture of driver and passenger mutations with the passenger rates modeled in relationship to tumor mutational burden. We use simulations to study the operating characteristics of the method and assess false-positive and false-negative rates in driver nomination. When applied to a large study of primary melanomas, the method accurately identifies the known driver genes within the RTK-RAS pathway and nominates several rare variants as prime candidates for functional validation. A comprehensive evaluation of MAGPIE against existing tools has also been conducted leveraging the Cancer Genome Atlas data.

Integrated Proteogenomic Analysis Reveals Distinct Potentially Actionable Therapeutic Vulnerabilities in Triple-Negative Breast Cancer Subtypes

Triple-negative breast cancer (TNBC) is characterized by an aggressive clinical presentation and a paucity of clinically actionable genomic alterations. Here, we utilized the Cancer Genome Atlas (TCGA) to explore the proteogenomic landscape of TNBC subtypes to see whether genomic alterations can be inferred from proteomic data. We found only 4% of the protein level changes are explained by mutations, while 21% of the protein and 35% of the transcriptomics changes were determined by copy number alterations (CNAs). We found tighter coupling between proteome and genome in some genes that are predicted to be the targets of drug inhibitors, including CDKs, PI3K, tyrosine kinase (TKI), and mTOR. The validation of our proteogenomic workflow using mass spectrometry Clinical Proteomic Tumor Analysis Consortium (MS-CPTAC) data also demonstrated the highest correlation between protein-RNA-CNA. The integrated proteogenomic approach helps to prioritize potentially actionable targets and may enable the acceleration of personalized cancer treatment.

Whole-Exome Sequencing Identifies Mutation Profile and Mutation Signature-Based Clustering Associated with Prognosis in Appendiceal Pseudomyxoma Peritonei

Pseudomyxoma peritonei (PMP) is a rare malignant clinical syndrome with little known about the global mutation profile. In this study, whole-exome sequencing (WES) was performed in 49 appendiceal PMP to investigate mutation profiles and mutation signatures. A total of 4,020 somatic mutations were detected, with a median mutation number of 56 (1-402). Tumor mutation burden (TMB) was generally low (median 1.55 mutations/Mb, 0.12-11.26 mutations/Mb). Mutations were mainly enriched in the function of cancer-related axonogenesis, extracellular matrix-related processes, calcium signaling pathway, and cAMP signaling pathway. Mutations in FCGBP, RBFOX1, SPEG, RTK-RAS, PI3K-AKT, and focal adhesion pathways were associated with high-grade mucinous carcinoma peritonei. These findings revealed distinct mutation profile in appendiceal PMP. Ten mutation signatures were identified, dividing patients into mutation signature cluster (MSC) 1 (N = 28, 57.1%) and MSC 2 (N = 21, 42.9%) groups. MSC (P = 0.007) was one of the four independent factors associated with 3-year survival. TMB (P = 0.003) and microsatellite instability (P = 0.002) were independent factors associated with MSC 2 grouping. Taken together, our findings provided a broader view in the understanding of molecular pathologic mechanism in appendiceal PMP and may be critical to developing an individualized approach to appendiceal PMP treatment.

IMPLICATIONS

This work describes exhaustive mutation profile of PMP based on WES data and derives ten mutation signatures, which divides patients into two clusters and serve as an independent prognostic factor associated with 3-year survival.

Computational Approaches: A New Frontier in Cancer Research

Cancer is a broad category of disease that can start in virtually any organ or tissue of the body when aberrant cells assault surrounding organs and proliferate uncontrollably. According to the most recent statistics, cancer will be the cause of 10 million deaths worldwide in 2020, accounting for one death out of every six worldwide. The typical approach used in anti-cancer research is highly time-consuming and expensive, and the outcomes are not particularly encouraging. Computational techniques have been employed in anti-cancer research to advance our understanding. Recent years have seen a significant and exceptional impact on anticancer research due to the rapid development of computational tools for novel drug discovery, drug design, genetic studies, genome characterization, cancer imaging and detection, radiotherapy, cancer metabolomics, and novel therapeutic approaches. In this paper, we examined the various subfields of contemporary computational techniques, including molecular docking, artificial intelligence, bioinformatics, virtual screening, and QSAR, and their applications in the study of cancer.

Advances in computational methods for identifying cancer driver genes

Cancer driver genes (CDGs) are crucial in cancer prevention, diagnosis and treatment. This study employed computational methods for identifying CDGs, categorizing them into four groups. The major frameworks for each of these four categories were summarized. Additionally, we systematically gathered data from public databases and biological networks, and we elaborated on computational methods for identifying CDGs using the aforementioned databases. Further, we summarized the algorithms, mainly involving statistics and machine learning, used for identifying CDGs. Notably, the performances of nine typical identification methods for eight types of cancer were compared to analyze the applicability areas of these methods. Finally, we discussed the challenges and prospects associated with methods for identifying CDGs. The present study revealed that the network-based algorithms and machine learning-based methods demonstrated superior performance.

Multi-omics analysis of multiple myeloma patients with differential response to first-line treatment

The genome backgrounds of multiple myeloma (MM) would affect the efficacy of specific treatment. However, the mutational and transcriptional landscapes in MM patients with differential response to first-line treatment remains unclear. We collected paired whole-exome sequencing (WES) and transcriptomic data of over 200 MM cases from MMRF-COMPASS project. R package, maftools was applied to analyze the somatic mutations and mutational signatures across MM samples. Differential expressed genes (DEG) was calculated using R package, DESeq2. The feature selection of the predictive model was determined by LASSO regression. In silico analysis revealed newly discovered recurrent mutated genes such as TTN, MUC16. TP53 mutation was observed more frequent in nonCR (complete remission) group with poor prognosis. DNA repair-associated mutational signatures were enriched in CR patients. Transcriptomic profiling showed that the activity of NF-kappa B and TGF-β pathways was suppressed in CR patients. A transcriptome-based response predictive model was constructed and showed promising predictive accuracy in MM patients receiving first-line treatment. Our study delineated distinctive mutational and transcriptional landscapes in MM patients with differential response to first-line treatment. Furthermore, we constructed a 20-gene predictive model which showed promising accuracy in predicting treatment response in newly diagnosed MM patients.

Constructing an APOBEC-related gene signature with predictive value in the overall survival and therapeutic sensitivity in lung adenocarcinoma

Background

APOBEC family play an important role in cancer mutagenesis and tumor development. The role of APOBEC family in lung adenocarcinoma (LUAD) has not been studied comprehensively.

Materials and methods

The expression data of pan-cancer as well as LUAD was obtained from public databases. The expression level of APOBEC family genes was analyzed in different normal and cancer tissues. APOBEC mutagenesis enrichment score (AMES) was utilized to evaluate the APOBEC-induced mutations and the relation of APOBEC with genomic instability. Gene set enrichment analysis was used to identify differentially enriched pathways. Univariate Cox regression and Lasso regression were applied to screen key prognostic genes. The immune cell infiltration was estimated by CIBERSORT. RT-qPCR assay, CCK-8 and Transwell assay were conducted to explore gene expression and lung cancer cell invasion.

Results

Cancer tissues had significantly altered expression of APOBEC family genes and the expression patterns of APOBEC family were different in different cancer types. APOBEC3B was the most aberrantly expressed in most cancer types. In LUAD, we observed a significantly positive correlation of AMES with intratumor heterogeneity (ITH), tumor neoantigen burden (TNB), and tumor mutation burden (TMB). High AMES group had high mutation counts of DNA damage repair pathways, and high enrichment of cell cycle and DNA repair pathways. We identified four prognostic genes (LYPD3, ANLN, MUC5B, and FOSL1) based on AMES, and constructed an AMES-related gene signature. The expressions of four genes were enhanced and accelerated the invasion ability and viability of lung cancer cells. Furthermore, we found that high group increased oxidative stress level.

Conclusions

APOBEC family was associated with genomic instability, DNA damage-related pathways, and cell cycle in LUAD. The AMES-related gene signature had a great potential to indicate the prognosis and guide immunotherapy/chemotherapy for patients suffering from LUAD.

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

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

A Graph Convolution Network-Based Model for Prioritizing Personalized Cancer Driver Genes of Individual Patients

Cancer driver genes are mutated genes that play a key role in the growth of cancer cells. Accurately identifying the cancer driver genes helps us understand cancer's pathogenesis and develop effective treatment strategies. However, cancers are highly heterogeneous diseases; patients with the same cancer type may have different genomic characteristics and clinical symptoms. Hence, it is urgent to devise effective methods to identify personalized cancer driver genes of individual patients to help determine whether a patient can be treated with a certain targeted drug. This work presents a method for predicting personalized cancer Driver genes of individual patients based on Graph Convolution Networks and Neighbor Interactions called NIGCNDriver. NIGCNDriver first constructs a gene-sample association matrix using the associations between a sample and its known driver genes. Then, it employs graph convolution models on the gene-sample network to aggregate neighbor node features, and themself features, and then combines with the element-wise level interactions between neighbors to learn new feature representations for the samples and gene nodes. Finally, a linear correlation coefficient decoder is used to reconstruct the association between the sample and the mutant gene, enabling the prediction of a personalized driver gene for the individual sample. We applied the NIGCNDriver method to predict cancer driver genes for individual samples in the TCGA and cancer cell line datasets. The results show that our method outperforms the baseline methods in cancer driver gene prediction for individual samples.

Whole exome sequencing and MicroRNA profiling of lung adenocarcinoma identified risk prediction features for tumors at stage I and its substages

OBJECTIVES

About 20% of stage I lung adenocarcinoma (LUAD) patients suffer a relapse after surgical resection. While finer substages have been defined and refined in the AJCC staging system, clinical investigations on the tumor molecular landscape are lacking.

MATERIALS AND METHODS

We performed whole exome sequencing, DNA copy number and microRNA profiling on paired tumor-normal samples from a cohort of 113 treatment-naïve stage I Taiwanese LUAD patients. We searched for molecular features associated with relapse-free survival (RFS) of stage I or its substages and validated the findings with an independent Caucasian LUAD cohort.

RESULTS

We found sixteen nonsynonymous mutations harbored at EGFR, KRAS, TP53, CTNNB1 and six other genes associated with poor RFS in a dose-dependent manner via variant allele fraction (VAF). An index, maxVAF, was constructed to quantify the overall mutation load from genes other than EGFR. High maxVAF scores discriminated a small group of high-risk LUAD at stage I (median RFS: 4.5 versus 69.5 months; HR = 10.5, 95% CI = 4.22-26.12, P < 0.001). At the substage level, higher risk was found for patients with high maxVAF or high miR-31; IA (median RFS: 32.1 versus 122.8 months, P = 0.005) and IB (median RFS: 7.1 versus 26.2, P = 0.049). MicroRNAs, miR-182, miR-183 and miR-196a were found correlated with EGFR mutation and poor RFS in stage IB patients.

CONCLUSION

Distinctive features of somatic gene mutation and microRNA expression of stage I LUAD are characterized to complement the survival prognosis by substaging. The findings open up more options for precision management of stage I LUAD patients.

Insight of novel biomarkers for papillary thyroid carcinoma through multiomics

Introduction

The overdiagnosing of papillary thyroid carcinoma (PTC) in China necessitates the development of an evidence-based diagnosis and prognosis strategy in line with precision medicine. A landscape of PTC in Chinese cohorts is needed to provide comprehensiveness.

Methods

6 paired PTC samples were employed for whole-exome sequencing, RNA sequencing, and data-dependent acquisition mass spectrum analysis. Weighted gene co-expression network analysis and protein-protein interactions networks were used to screen for hub genes. Moreover, we verified the hub genes' diagnostic and prognostic potential using online databases. Logistic regression was employed to construct a diagnostic model, and we evaluated its efficacy and specificity based on TCGA-THCA and GEO datasets.

Results

The basic multiomics landscape of PTC among local patients were drawn. The similarities and differences were compared between the Chinese cohort and TCGA-THCA cohorts, including the identification of PNPLA5 as a driver gene in addition to BRAF mutation. Besides, we found 572 differentially expressed genes and 79 differentially expressed proteins. Through integrative analysis, we identified 17 hub genes for prognosis and diagnosis of PTC. Four of these genes, ABR, AHNAK2, GPX1, and TPO, were used to construct a diagnostic model with high accuracy, explicitly targeting PTC (AUC=0.969/0.959 in training/test sets).

Discussion

Multiomics analysis of the Chinese cohort demonstrated significant distinctions compared to TCGA-THCA cohorts, highlighting the unique genetic characteristics of Chinese individuals with PTC. The novel biomarkers, holding potential for diagnosis and prognosis of PTC, were identified. Furthermore, these biomarkers provide a valuable tool for precise medicine, especially for immunotherapeutic or nanomedicine based cancer therapy.

A cancer-associated fibroblast subtypes-based signature enables the evaluation of immunotherapy response and prognosis in bladder cancer

Bladder cancer (BLCA) is one of the most prevalent and heterogeneous urinary malignant tumors. Previous researches have reported a significant association between cancer-associated fibroblasts (CAFs) and poor prognosis of tumor patients. However, uncertainty surrounds the role of CAFs in the BLCA tumor microenvironment, necessitating further investigation into the CAFs-related gene signatures in BLCA. In this study, we identified three CAF subtypes in BLCA according to single-cell RNA-seq data and constructed CAFs-related risk score (CRRS) by screening 102,714 signatures. The survival analysis, ROC curves, and nomogram suggested that CRRS was a valuable predictor in 2,042 patients from 9 available public datasets and Xiangya real-world cohort. We further revealed the significant correlation between CRRS and clinicopathological characteristics, genome alterations, and epithelial-mesenchymal transition (EMT). A high CRRS indicated a non-inflamed phenotype and a lower remission rate of immunotherapy in BLCA. In conclusion, the CRRS had the potential to predict the prognosis and immunotherapy response of BLCA patients.

Genetic alterations in peritoneal metastatic tumors predicted the outcomes for hyperthermic intraperitoneal chemotherapy

Introduction

Cytoreductive surgery (CRS) and hyperthermic intraperitoneal chemotherapy (HIPEC) are considered for patients with peritoneal metastasis (PM). However, patients selection that relies on conventional prognostic factors is not yet optimal. In this study, we performed whole exome sequencing (WES) to establish tumor molecular characteristics and expect to identify prognosis profiles for PM management.

Methods

In this study, blood and tumor samples were collected from patients with PM before HIPEC. Tumor molecular signatures were determined using WES. Patient cohort was divided into responders and non-responders according to 12-month progression-free survival (PFS). Genomic characteristics between the two cohorts were compared to study potential targets.

Results

In total, 15 patients with PM were enrolled in this study. Driver genes and enriched pathways were identified from WES results. AGAP5 mutation was found in all responders. This mutation was significantly associated with better OS (p = 0.00652).

Conclusions

We identified prognostic markers that might be useful to facilitate decision-making before CRS/HIPEC.

Systematic integration of protein-affecting mutations, gene fusions, and copy number alterations into a comprehensive somatic mutational profile

Somatic mutations occur as random genetic changes in genes through protein-affecting mutations (PAMs), gene fusions, or copy number alterations (CNAs). Mutations of different types can have a similar phenotypic effect (i.e., allelic heterogeneity) and should be integrated into a unified gene mutation profile. We developed OncoMerge to fill this niche of integrating somatic mutations to capture allelic heterogeneity, assign a function to mutations, and overcome known obstacles in cancer genetics. Application of OncoMerge to TCGA Pan-Cancer Atlas increased detection of somatically mutated genes and improved the prediction of the somatic mutation role as either activating or loss of function. Using integrated somatic mutation matrices increased the power to infer gene regulatory networks and uncovered the enrichment of switch-like feedback motifs and delay-inducing feedforward loops. These studies demonstrate that OncoMerge efficiently integrates PAMs, fusions, and CNAs and strengthens downstream analyses linking somatic mutations to cancer phenotypes.

Cancer driver mutations: predictions and reality

Cancer cells accumulate many genetic alterations throughout their lifetime, but only a few of them drive cancer progression, termed driver mutations. Driver mutations may vary between cancer types and patients, can remain latent for a long time and become drivers at particular cancer stages, or may drive oncogenesis only in conjunction with other mutations. The high mutational, biochemical, and histological tumor heterogeneity makes driver mutation identification very challenging. In this review we summarize recent efforts to identify driver mutations in cancer and annotate their effects. We underline the success of computational methods to predict driver mutations in finding novel cancer biomarkers, including in circulating tumor DNA (ctDNA). We also report on the boundaries of their applicability in clinical research.

An immunogenic and oncogenic feature-based classification for chemotherapy plus PD-1 blockade in advanced esophageal squamous cell carcinoma

Although chemotherapy plus PD-1 blockade (chemo+anti-PD-1) has become the standard first-line therapy for advanced esophageal squamous cell carcinoma (ESCC), reliable biomarkers for this regimen are lacking. Here we perform whole-exome sequencing on tumor samples from 486 patients of the JUPITER-06 study and develop a copy number alteration-corrected tumor mutational burden that depicts immunogenicity more precisely and predicts chemo+anti-PD-1 efficacy. We identify several other favorable immunogenic features (e.g., HLA-I/II diversity) and risk oncogenic alterations (e.g., PIK3CA and TET2 mutation) associated with chemo+anti-PD-1 efficacy. An esophageal cancer genome-based immuno-oncology classification (EGIC) scheme incorporating these immunogenic features and oncogenic alterations is established. Chemo+anti-PD-1 achieves significant survival improvements in EGIC1 (immunogenic feature-favorable and oncogenic alteration-negative) and EGIC2 (either immunogenic feature-favorable or oncogenic alteration-negative) subgroups, but not the EGIC3 subgroup (immunogenic feature-unfavorable and oncogenic alteration-positive). Thus, EGIC may guide future individualized treatment strategies and inform mechanistic biomarker research for chemo+anti-PD-1 treatment in patients with advanced ESCC.

LRT-CLUSTER: A New Clustering Algorithm Based on Likelihood Ratio Test to Identify Driving Genes

Somatic mutations often occur at high relapse sites in protein sequences, which indicates that the location clustering of somatic missense mutations can be used to identify driving genes. However, the traditional clustering algorithm has such problems as the background signal over-fitting, the clustering algorithm is not suitable for mutation data, and the performance of identifying low-frequency mutation genes needs to be improved. In this paper, we propose a linear clustering algorithm based on likelihood ratio test knowledge to identify driver genes. In this experiment, firstly, the polynucleotide mutation rate is calculated based on the prior knowledge of likelihood ratio test. Then, the simulation data set is obtained through the background mutation rate model. Finally, the unsupervised peak clustering algorithm is used to, respectively, evaluate the somatic mutation data and the simulation data to identify the driver genes. The experimental results show that our method achieves a better balance of precision and sensitivity. It can also identify the driver genes missed by other methods, making it an effective supplement to other methods. We also discover some potential linkages between genes and between genes and mutation sites, which is of great value to target drug therapy research. Method framework: Our proposed model framework is as follows. a. Counting mutation sites and the number of mutations in tumor gene elements. b. The nucleotide context mutation frequency is counted based on the likelihood ratio test knowledge, and the background mutation rate model is obtained. c. Based on Monte Carlo simulation method, data sets with the same number of mutations as gene elements are randomly sampled to obtain simulated mutation data, and the sampling frequency of each mutation site is related to the mutation rate of polynucleotide. d. The original mutation data and the simulated mutation data after random reconstruction are clustered by peak density, respectively, and the corresponding clustering scores are obtained. e. We can obtain the clustering information statistics in each gene segment and score of each gene segment from the original single nucleotide mutation data through step d. f. According to the observed score and the simulated clustering score, the p-value of the corresponding gene fragment is calculated. g. We can obtain the clustering information statistics in each gene segment and score of each gene segment from the simulated single nucleotide mutation data through step d.

Prognosis and Personalized Treatment Prediction in Different Mutation-Signature Hepatocellular Carcinoma

Introduction

Mutation patterns have been extensively explored to decipher the etiologies of hepatocellular carcinoma (HCC). However, the study and potential clinical role of mutation patterns to stratify high-risk patients and optimize precision therapeutic strategies remain elusive in HCC.

Methods

Using exon-sequencing data in public (n=362) and in-house (n=30) cohorts, mutation signatures were extracted to decipher relationships with the etiology and prognosis in HCC. The proteomics (n=159) and cell-line transcriptome data (n=1019) were collected to screen the implication of sensitive drugs. A novel multi-step machine-learning framework was then performed to construct a classification predictor, including recognizing stable reversed gene pairs, establishing a robust prediction model, and validating the robustness of the predictor in five independent cohorts (n=900).

Results

Two heterogeneous mutation signature clusters were identified, and a high-risk prognosis cluster was recognized for further analysis. Notably, mutation signature cluster 1 (MSC1) was featured by activated anti-tumor immune and metabolism dysfunctional states, higher genomic instability (high TMB, SNV neoantigen, indel neoantigens, and total neoantigens), and a dismal prognosis. Notably, MSC performed as an independent risk factor than clinical traits (eg, stage, vascular invasion). Additionally, afatinib and canertinib were recognized which might have potential therapeutic implications in MSC1, and the targets of these drugs presented a higher expression in both gene and protein levels in HCC.

Discussion

Our studies may provide a promising platform for improving prognosis and tailoring therapy in HCC.

Identification of three metabolic subtypes in gastric cancer and the construction of a metabolic pathway-based risk model that predicts the overall survival of GC patients

Gastric cancer (GC) is highly heterogeneous and GC patients have low overall survival rates. It is also challenging to predict the prognosis of GC patients. This is partly because little is known about the prognosis-related metabolic pathways in this disease. Hence, our objective was to identify GC subtypes and genes related to prognosis, based on changes in the activity of core metabolic pathways in GC tumor samples. Differences in the activity of metabolic pathways in GC patients were analyzed using Gene Set Variation Analysis (GSVA), leading to the identification of three clinical subtypes by non-negative matrix factorization (NMF). Based on our analysis, subtype 1 showed the best prognosis while subtype 3 exhibited the worst prognosis. Interestingly, we observed marked differences in gene expression between the three subtypes, through which we identified a new evolutionary driver gene, CNBD1. Furthermore, we used 11 metabolism-associated genes identified by LASSO and random forest algorithms to construct a prognostic model and verified our results using qRT-PCR (five matched clinical tissues of GC patients). This model was found to be both effective and robust in the GSE84437 and GSE26253 cohorts, and the results from multivariate Cox regression analyses confirmed that the 11-gene signature was an independent prognostic predictor (p < 0.0001, HR = 2.8, 95% CI 2.1-3.7). The signature was found to be relevant to the infiltration of tumor-associated immune cells. In conclusion, our work identified significant GC prognosis-related metabolic pathways in different GC subtypes and provided new insights into GC-subtype prognostic assessment.

Model performance and interpretability of semi-supervised generative adversarial networks to predict oncogenic variants with unlabeled data

BACKGROUND

It remains an important challenge to predict the functional consequences or clinical impacts of genetic variants in human diseases, such as cancer. An increasing number of genetic variants in cancer have been discovered and documented in public databases such as COSMIC, but the vast majority of them have no functional or clinical annotations. Some databases, such as CiVIC are available with manual annotation of functional mutations, but the size of the database is small due to the use of human annotation. Since the unlabeled data (millions of variants) typically outnumber labeled data (thousands of variants), computational tools that take advantage of unlabeled data may improve prediction accuracy.

RESULT

To leverage unlabeled data to predict functional importance of genetic variants, we introduced a method using semi-supervised generative adversarial networks (SGAN), incorporating features from both labeled and unlabeled data. Our SGAN model incorporated features from clinical guidelines and predictive scores from other computational tools. We also performed comparative analysis to study factors that influence prediction accuracy, such as using different algorithms, types of features, and training sample size, to provide more insights into variant prioritization. We found that SGAN can achieve competitive performances with small labeled training samples by incorporating unlabeled samples, which is a unique advantage compared to traditional machine learning methods. We also found that manually curated samples can achieve a more stable predictive performance than publicly available datasets.

CONCLUSIONS

By incorporating much larger samples of unlabeled data, the SGAN method can improve the ability to detect novel oncogenic variants, compared to other machine-learning algorithms that use only labeled datasets. SGAN can be potentially used to predict the pathogenicity of more complex variants such as structural variants or non-coding variants, with the availability of more training samples and informative features.

Integrating multi-type aberrations from DNA and RNA through dynamic mapping gene space for subtype-specific breast cancer driver discovery

Driver event discovery is a crucial demand for breast cancer diagnosis and therapy. In particular, discovering subtype-specificity of drivers can prompt the personalized biomarker discovery and precision treatment of cancer patients. Still, most of the existing computational driver discovery studies mainly exploit the information from DNA aberrations and gene interactions. Notably, cancer driver events would occur due to not only DNA aberrations but also RNA alternations, but integrating multi-type aberrations from both DNA and RNA is still a challenging task for breast cancer drivers. On the one hand, the data formats of different aberration types also differ from each other, known as data format incompatibility. On the other hand, different types of aberrations demonstrate distinct patterns across samples, known as aberration type heterogeneity. To promote the integrated analysis of subtype-specific breast cancer drivers, we design a "splicing-and-fusing" framework to address the issues of data format incompatibility and aberration type heterogeneity simultaneously. To overcome the data format incompatibility, the "splicing-step" employs a knowledge graph structure to connect multi-type aberrations from the DNA and RNA data into a unified formation. To tackle the aberration type heterogeneity, the "fusing-step" adopts a dynamic mapping gene space integration approach to represent the multi-type information by vectorized profiles. The experiments also demonstrate the advantages of our approach in both the integration of multi-type aberrations from DNA and RNA and the discovery of subtype-specific breast cancer drivers. In summary, our "splicing-and-fusing" framework with knowledge graph connection and dynamic mapping gene space fusion of multi-type aberrations data from DNA and RNA can successfully discover potential breast cancer drivers with subtype-specificity indication.

Proteogenomics of diffuse gliomas reveal molecular subtypes associated with specific therapeutic targets and immune-evasion mechanisms

Diffuse gliomas are devastating brain tumors. Here, we perform a proteogenomic profiling of 213 retrospectively collected glioma tumors. Proteogenomic analysis reveals the downstream biological events leading by EGFR-, IDH1-, TP53-mutations. The comparative analysis illustrates the distinctive features of GBMs and LGGs, indicating CDK2 inhibitor might serve as a promising drug target for GBMs. Further proteogenomic integrative analysis combined with functional experiments highlight the cis-effect of EGFR alterations might lead to glioma tumor cell proliferation through ERK5 medicates nucleotide synthesis process. Proteome-based stratification of gliomas defines 3 proteomic subgroups (S-Ne, S-Pf, S-Im), which could serve as a complement to WHO subtypes, and would provide the essential framework for the utilization of specific targeted therapies for particular glioma subtypes. Immune clustering identifies three immune subtypes with distinctive immune cell types. Further analysis reveals higher EGFR alteration frequencies accounts for elevation of immune check point protein: PD-L1 and CD70 in T-cell infiltrated tumors.

From multi-omics data to the cancer druggable gene discovery: a novel machine learning-based approach

The development of targeted drugs allows precision medicine in cancer treatment and optimal targeted therapies. Accurate identification of cancer druggable genes helps strengthen the understanding of targeted cancer therapy and promotes precise cancer treatment. However, rare cancer-druggable genes have been found due to the multi-omics data's diversity and complexity. This study proposes deep forest for cancer druggable genes discovery (DF-CAGE), a novel machine learning-based method for cancer-druggable gene discovery. DF-CAGE integrated the somatic mutations, copy number variants, DNA methylation and RNA-Seq data across ˜10 000 TCGA profiles to identify the landscape of the cancer-druggable genes. We found that DF-CAGE discovers the commonalities of currently known cancer-druggable genes from the perspective of multi-omics data and achieved excellent performance on OncoKB, Target and Drugbank data sets. Among the ˜20 000 protein-coding genes, DF-CAGE pinpointed 465 potential cancer-druggable genes. We found that the candidate cancer druggable genes (CDG) are clinically meaningful and divided the CDG into known, reliable and potential gene sets. Finally, we analyzed the omics data's contribution to identifying druggable genes. We found that DF-CAGE reports druggable genes mainly based on the copy number variations (CNVs) data, the gene rearrangements and the mutation rates in the population. These findings may enlighten the future study and development of new drugs.

Gain-of-Function Variomics and Multi-omics Network Biology for Precision Medicine

Traditionally, disease causal mutations were thought to disrupt gene function. However, it becomes more clear that many deleterious mutations could exhibit a "gain-of-function" (GOF) behavior. Systematic investigation of such mutations has been lacking and largely overlooked. Advances in next-generation sequencing have identified thousands of genomic variants that perturb the normal functions of proteins, further contributing to diverse phenotypic consequences in disease. Elucidating the functional pathways rewired by GOF mutations will be crucial for prioritizing disease-causing variants and their resultant therapeutic liabilities. In distinct cell types (with varying genotypes), precise signal transduction controls cell decision, including gene regulation and phenotypic output. When signal transduction goes awry due to GOF mutations, it would give rise to various disease types. Quantitative and molecular understanding of network perturbations by GOF mutations may provide explanations for 'missing heritability" in previous genome-wide association studies. We envision that it will be instrumental to push current paradigm toward a thorough functional and quantitative modeling of all GOF mutations and their mechanistic molecular events involved in disease development and progression. Many fundamental questions pertaining to genotype-phenotype relationships remain unresolved. For example, which GOF mutations are key for gene regulation and cellular decisions? What are the GOF mechanisms at various regulation levels? How do interaction networks undergo rewiring upon GOF mutations? Is it possible to leverage GOF mutations to reprogram signal transduction in cells, aiming to cure disease? To begin to address these questions, we will cover a wide range of topics regarding GOF disease mutations and their characterization by multi-omic networks. We highlight the fundamental function of GOF mutations and discuss the potential mechanistic effects in the context of signaling networks. We also discuss advances in bioinformatic and computational resources, which will dramatically help with studies on the functional and phenotypic consequences of GOF mutations.

Prognosis and Characterization of Microenvironment in Cervical Cancer Influenced by Fatty Acid Metabolism-Related Genes

Increasing evidence suggests that diverse activation patterns of metabolic signalling pathways may lead to molecular diversity of cervical cancer (CC). But rare research focuses on the alternation of fatty acid metabolism (FAM) in CC. Therefore, we constructed and compared models based on the expression of FAM-related genes from the Cancer Genome Atlas by different machine learning algorithms. The most reliable model was built with 14 significant genes by LASSO-Cox regression, and the CC cohort was divided into low-/high-risk groups by the median of risk score. Then, a feasible nomogram was established and validated by C-index, calibration curve, net benefit, and decision curve analysis. Furthermore, the hub genes among differential expression genes were identified and the post-transcriptional and translational regulation networks were characterized. Moreover, the somatic mutation and copy number variation landscapes were depicted. Importantly, the specific mutation drivers and signatures of the FAM phenotypes were excavated. As a result, the high-risk samples were featured by activated de novo fatty acid synthesis, epithelial to mesenchymal transition, angiogenesis, and chronic inflammation response, which might be caused by mutations of oncogenic driver genes in RTK/RAS, PI3K, and NOTCH signalling pathways. Besides the hyperactivity of cytidine deaminase and deficiency of mismatch repair, the mutations of POLE might be partially responsible for the mutations in the high-risk group. Next, the antigenome including the neoantigen and cancer germline antigens was estimated. The decreasing expression of a series of cancer germline antigens was identified to be related to reduction of CD8 T cell infiltration in the high-risk group. Then, the comprehensive evaluation of connotations between the tumour microenvironment and FAM phenotypes demonstrated that the increasing risk score was related to the suppressive immune microenvironment. Finally, the prediction of therapy targets revealed that the patients with high risk might be sensitive to the RAF inhibitor AZ628. Our findings provide a novel insight for personalized treatment in CC.

Repulsion and attraction in searching: A hybrid algorithm based on gravitational kernel and vital few for cancer driver gene prediction

By taking a new perspective to combine a machine learning method with an evolutionary algorithm, a new hybrid algorithm is developed to predict cancer driver genes. Firstly, inspired by the search strategy with the capability of global search in evolutionary algorithms, a gravitational kernel is proposed to act on the full range of gene features. Constructed by fusing PPI and mutation features, the gravitational kernel is capable to produce repulsion effects. The candidate genes with greater mutation effects and PPI have higher similarity scores. According to repulsion, the similarity score of these promising genes is larger than ordinary genes, which is beneficial to search for these promising genes. Secondly, inspired by the idea of elite populations related to evolutionary algorithms, the concept of vital few is proposed. Targeted at a local scale, it acts on the candidate genes associated with vital few genes. Under attraction effect, these vital few driver genes attract those with similar mutational effects to them, which leads to greater similarity scores. Lastly, the model and parameters are optimized by using an evolutionary algorithm, so as to obtain the optimal model and parameters for cancer driver gene prediction. Herein, a comparison is performed with six other advanced methods of cancer driver gene prediction. According to the experimental results, the method proposed in this study outperforms these six state-of-the-art algorithms on the pan-oncogene dataset.

Proteogenomic insights into the biology and treatment of pancreatic ductal adenocarcinoma

BACKGROUND

Pancreatic ductal adenocarcinoma (PDAC) is a devastating disease with poor prognosis. Proteogenomic characterization and integrative proteomic analysis provide a functional context to annotate genomic abnormalities with prognostic value.

METHODS

We performed an integrated multi-omics analysis, including whole-exome sequencing, RNA-seq, proteomic, and phosphoproteomic analysis of 217 PDAC tumors with paired non-tumor adjacent tissues. In vivo functional experiments were performed to further illustrate the biological events related to PDAC tumorigenesis and progression.

RESULTS

A comprehensive proteogenomic landscape revealed that TP53 mutations upregulated the CDK4-mediated cell proliferation process and led to poor prognosis in younger patients. Integrative multi-omics analysis illustrated the proteomic and phosphoproteomic alteration led by genomic alterations such as KRAS mutations and ADAM9 amplification of PDAC tumorigenesis. Proteogenomic analysis combined with in vivo experiments revealed that the higher amplification frequency of ADAM9 (8p11.22) could drive PDAC metastasis, though downregulating adhesion junction and upregulating WNT signaling pathway. Proteome-based stratification of PDAC revealed three subtypes (S-I, S-II, and S-III) related to different clinical and molecular features. Immune clustering defined a metabolic tumor subset that harbored FH amplicons led to better prognosis. Functional experiments revealed the role of FH in altering tumor glycolysis and in impacting PDAC tumor microenvironments. Experiments utilizing both in vivo and in vitro assay proved that loss of HOGA1 promoted the tumor growth via activating LARP7-CDK1 pathway.

CONCLUSIONS

This proteogenomic dataset provided a valuable resource for researchers and clinicians seeking for better understanding and treatment of PDAC.

Somatic mutation distribution across tumour cohorts provides a signal for positive selection in cancer

Cancer gene discovery is reliant on distinguishing driver mutations from a multitude of passenger mutations in tumour genomes. While driver genes may be revealed based on excess mutation recurrence or clustering, there is a need for orthogonal principles. Here, we take advantage of the fact that non-cancer genes, containing only passenger mutations under neutral selection, exhibit a likelihood of mutagenesis in a given tumour determined by the tumour's mutational signature and burden. This relationship can be disrupted by positive selection, leading to a difference in the distribution of mutated cases across a cohort for driver and passenger genes. We apply this principle to detect cancer drivers independently of recurrence in large pan-cancer cohorts, and show that our method (SEISMIC) performs comparably to traditional approaches and can provide resistance to known confounding mutational phenomena. Being based on a different principle, the approach provides a much-needed complement to existing methods for detecting signals of selection.

Parallel functional annotation of cancer-associated missense mutations in histone methyltransferases

Using exome sequencing for biomarker discovery and precision medicine requires connecting nucleotide-level variation with functional changes in encoded proteins. However, for functionally annotating the thousands of cancer-associated missense mutations, or variants of uncertain significance (VUS), purifying variant proteins for biochemical and functional analysis is cost-prohibitive and inefficient. We describe parallel functional annotation (PFA) of large numbers of VUS using small cultures and crude extracts in 96-well plates. Using members of a histone methyltransferase family, we demonstrate high-throughput structural and functional annotation of cancer-associated mutations. By combining functional annotation of paralogs, we discovered two phylogenetic and clustering parameters that improve the accuracy of sequence-based functional predictions to over 90%. Our results demonstrate the value of PFA for defining oncogenic/tumor suppressor functions of histone methyltransferases as well as enhancing the accuracy of sequence-based algorithms in predicting the effects of cancer-associated mutations.