Functional impact bias reveals cancer drivers

The integrated molecular and histological analysis defines subtypes of esophageal squamous cell carcinoma

Esophageal squamous cell carcinoma (ESCC) is highly heterogeneous. Our understanding of full molecular and immune landscape of ESCC remains limited, hindering the development of personalised therapeutic strategies. To address this, we perform genomic-transcriptomic characterizations and AI-aided histopathological image analysis of 120 Chinese ESCC patients. Here we show that ESCC can be categorized into differentiated, metabolic, immunogenic and stemness subtypes based on bulk and single-cell RNA-seq, each exhibiting specific molecular and histopathological features based on an amalgamated deep-learning model. The stemness subgroup with signature genes, such as WFDC2, SFRP1, LGR6 and VWA2, has the poorest prognosis and is associated with downregulated immune activities, a high frequency of EP300 mutation/activation, functional mutation enrichment in Wnt signalling and the highest level of intratumoural heterogeneity. The immune profiling by transcriptomics and immunohistochemistry reveals ESCC cells overexpress natural killer cell markers XCL1 and CD160 as immune evasion. Strikingly, XCL1 expression also affects the sensitivity of ESCC cells to common chemotherapy drugs. This study opens avenues for ESCC treatment and provides a valuable public resource to better understand ESCC.

Bioinformatics tools and resources for cancer and application

ABSTRACT

Tumor bioinformatics plays an important role in cancer research and precision medicine. The primary focus of traditional cancer research has been molecular and clinical studies of a number of fundamental pathways and genes. In recent years, driven by breakthroughs in high-throughput technologies, large-scale cancer omics data have accumulated rapidly. How to effectively utilize and share these data is particularly important. To address this crucial task, many computational tools and databases have been developed over the past few years. To help researchers quickly learn and understand the functions of these tools, in this review, we summarize publicly available bioinformatics tools and resources for pan-cancer multi-omics analysis, regulatory analysis of tumorigenesis, tumor treatment and prognosis, immune infiltration analysis, immune repertoire analysis, cancer driver gene and driver mutation analysis, and cancer single-cell analysis, which may further help researchers find more suitable tools for their research.

CIBRA identifies genomic alterations with a system-wide impact on tumor biology

MOTIVATION

Genomic instability is a hallmark of cancer, leading to many somatic alterations. Identifying which alterations have a system-wide impact is a challenging task. Nevertheless, this is an essential first step for prioritizing potential biomarkers. We developed CIBRA (Computational Identification of Biologically Relevant Alterations), a method that determines the system-wide impact of genomic alterations on tumor biology by integrating two distinct omics data types: one indicating genomic alterations (e.g. genomics), and another defining a system-wide expression response (e.g. transcriptomics). CIBRA was evaluated with genome-wide screens in 33 cancer types using primary and metastatic cancer data from the Cancer Genome Atlas and Hartwig Medical Foundation.

RESULTS

We demonstrate the capability of CIBRA by successfully confirming the impact of point mutations in experimentally validated oncogenes and tumor suppressor genes (0.79 AUC). Surprisingly, many genes affected by structural variants were identified to have a strong system-wide impact (30.3%), suggesting that their role in cancer development has thus far been largely under-reported. Additionally, CIBRA can identify impact with only 10 cases and controls, providing a novel way to prioritize genomic alterations with a prominent role in cancer biology. Our findings demonstrate that CIBRA can identify cancer drivers by combining genomics and transcriptomics data. Moreover, our work shows an unexpected substantial system-wide impact of structural variants in cancer. Hence, CIBRA has the potential to preselect and refine current definitions of genomic alterations to derive more nuanced biomarkers for diagnostics, disease progression, and treatment response.

AVAILABILITY AND IMPLEMENTATION

The R package CIBRA is available at https://github.com/AIT4LIFE-UU/CIBRA.

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.

Somatic mutations associate with clonal expansion of CD8+ T cells

Somatic mutations in T cells can cause cancer but also have implications for immunological diseases and cell therapies. The mutation spectrum in nonmalignant T cells is unclear. Here, we examined somatic mutations in CD4+ and CD8+ T cells from 90 patients with hematological and immunological disorders and used T cell receptor (TCR) and single-cell sequencing to link mutations with T cell expansions and phenotypes. CD8+ cells had a higher mutation burden than CD4+ cells. Notably, the biggest variant allele frequency (VAF) of non-synonymous variants was higher than synonymous variants in CD8+ T cells, indicating non-random occurrence. The non-synonymous VAF in CD8+ T cells strongly correlated with the TCR frequency, but not age. We identified mutations in pathways essential for T cell function and often affected lymphoid neoplasia. Single-cell sequencing revealed cytotoxic TEMRA phenotypes of mutated T cells. Our findings suggest that somatic mutations contribute to CD8+ T cell expansions without malignant transformation.

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.

A latent variable model for evaluating mutual exclusivity and co-occurrence between driver mutations in cancer

A key challenge in cancer genomics is understanding the functional relationships and dependencies between combinations of somatic mutations that drive cancer development. Such driver mutations frequently exhibit patterns of mutual exclusivity or co-occurrence across tumors, and many methods have been developed to identify such dependency patterns from bulk DNA sequencing data of a cohort of patients. However, while mutual exclusivity and co-occurrence are described as properties of driver mutations, existing methods do not explicitly disentangle functional, driver mutations from neutral, passenger mutations. In particular, nearly all existing methods evaluate mutual exclusivity or co-occurrence at the gene level, marking a gene as mutated if any mutation - driver or passenger - is present. Since some genes have a large number of passenger mutations, existing methods either restrict their analyses to a small subset of suspected driver genes - limiting their ability to identify novel dependencies - or make spurious inferences of mutual exclusivity and co-occurrence involving genes with many passenger mutations. We introduce DIALECT, an algorithm to identify dependencies between pairs of driver mutations from somatic mutation counts. We derive a latent variable mixture model for drivers and passengers that combines existing probabilistic models of passenger mutation rates with a latent variable describing the unknown status of a mutation as a driver or passenger. We use an expectation maximization (EM) algorithm to estimate the parameters of our model, including the rates of mutually exclusivity and co-occurrence between drivers. We demonstrate that DIALECT more accurately infers mutual exclusivity and co-occurrence between driver mutations compared to existing methods on both simulated mutation data and somatic mutation data from 5 cancer types in The Cancer Genome Atlas (TCGA).

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.

Haplotype-specific assembly of shattered chromosomes in esophageal adenocarcinomas

The epigenetic landscape of cancer is regulated by many factors, but primarily it derives from the underlying genome sequence. Chromothripsis is a catastrophic localized genome shattering event that drives, and often initiates, cancer evolution. We characterized five esophageal adenocarcinoma organoids with chromothripsis using long-read sequencing and transcriptome and epigenome profiling. Complex structural variation and subclonal variants meant that haplotype-aware de novo methods were required to generate contiguous cancer genome assemblies. Chromosomes were assembled separately and scaffolded using haplotype-resolved Hi-C reads, producing accurate assemblies even with up to 900 structural rearrangements. There were widespread differences between the chromothriptic and wild-type copies of chromosomes in topologically associated domains, chromatin accessibility, histone modifications, and gene expression. Differential epigenome peaks were most enriched within 10 kb of chromothriptic structural variants. Alterations in transcriptome and higher-order chromosome organization frequently occurred near differential epigenetic marks. Overall, chromothripsis reshapes gene regulation, causing coordinated changes in epigenetic landscape, transcription, and chromosome conformation.

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.

Exploring new roles for RNA-binding proteins in epigenetic and gene regulation

A significant portion of the human proteome comprises RNA-binding proteins (RBPs) that play fundamental roles in numerous biological processes. In the last decade, there has been a staggering increase in RBP identification and classification, which has fueled interest in the evolving roles of RBPs and RBP-driven molecular mechanisms. Here, we focus on recent insights into RBP-dependent regulation of the epigenetic and transcriptional landscape. We describe advances in methodologies that define the RNA-protein interactome and machine-learning algorithms that are streamlining RBP discovery and predicting new RNA-binding regions. Finally, we present how RBP dysregulation leads to alterations in tumor-promoting gene expression and discuss the potential for targeting these RBPs for the development of new cancer therapeutics.

Exploring gene-patient association to identify personalized cancer driver genes by linear neighborhood propagation

BACKGROUND

Driver genes play a vital role in the development of cancer. Identifying driver genes is critical for diagnosing and understanding cancer. However, challenges remain in identifying personalized driver genes due to tumor heterogeneity of cancer. Although many computational methods have been developed to solve this problem, few efforts have been undertaken to explore gene-patient associations to identify personalized driver genes.

RESULTS

Here we propose a method called LPDriver to identify personalized cancer driver genes by employing linear neighborhood propagation model on individual genetic data. LPDriver builds personalized gene network based on the genetic data of individual patients, extracts the gene-patient associations from the bipartite graph of the personalized gene network and utilizes a linear neighborhood propagation model to mine gene-patient associations to detect personalized driver genes. The experimental results demonstrate that as compared to the existing methods, our method shows competitive performance and can predict cancer driver genes in a more accurate way. Furthermore, these results also show that besides revealing novel driver genes that have been reported to be related with cancer, LPDriver is also able to identify personalized cancer driver genes for individual patients by their network characteristics even if the mutation data of genes are hidden.

CONCLUSIONS

LPDriver can provide an effective approach to predict personalized cancer driver genes, which could promote the diagnosis and treatment of cancer. The source code and data are freely available at https://github.com/hyr0771/LPDriver .

Prediction of cancer driver genes and mutations: the potential of integrative computational frameworks

The vast amount of available sequencing data allows the scientific community to explore different genetic alterations that may drive cancer or favor cancer progression. Software developers have proposed a myriad of predictive tools, allowing researchers and clinicians to compare and prioritize driver genes and mutations and their relative pathogenicity. However, there is little consensus on the computational approach or a golden standard for comparison. Hence, benchmarking the different tools depends highly on the input data, indicating that overfitting is still a massive problem. One of the solutions is to limit the scope and usage of specific tools. However, such limitations force researchers to walk on a tightrope between creating and using high-quality tools for a specific purpose and describing the complex alterations driving cancer. While the knowledge of cancer development increases daily, many bioinformatic pipelines rely on single nucleotide variants or alterations in a vacuum without accounting for cellular compartments, mutational burden or disease progression. Even within bioinformatics and computational cancer biology, the research fields work in silos, risking overlooking potential synergies or breakthroughs. Here, we provide an overview of databases and datasets for building or testing predictive cancer driver tools. Furthermore, we introduce predictive tools for driver genes, driver mutations, and the impact of these based on structural analysis. Additionally, we suggest and recommend directions in the field to avoid silo-research, moving towards integrative frameworks.

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.

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.

Identification of Driver Epistatic Gene Pairs Combining Germline and Somatic Mutations in Cancer

Cancer arises from the complex interplay of various factors. Traditionally, the identification of driver genes focuses primarily on the analysis of somatic mutations. We describe a new method for the detection of driver gene pairs based on an epistasis analysis that considers both germline and somatic variations. Specifically, the identification of significantly mutated gene pairs entails the calculation of a contingency table, wherein one of the co-mutated genes can exhibit a germline variant. By adopting this approach, it is possible to select gene pairs in which the individual genes do not exhibit significant associations with cancer. Finally, a survival analysis is used to select clinically relevant gene pairs. To test the efficacy of the new algorithm, we analyzed the colon adenocarcinoma (COAD) and lung adenocarcinoma (LUAD) samples available at The Cancer Genome Atlas (TCGA). In the analysis of the COAD and LUAD samples, we identify epistatic gene pairs significantly mutated in tumor tissue with respect to normal tissue. We believe that further analysis of the gene pairs detected by our method will unveil new biological insights, enhancing a better description of the cancer mechanism.

Noninvasive genomic profiling of somatic mutations in oral cavity cancers

OBJECTIVES

Somatic mutations may predict prognosis, therapeutic response, or cancer progression. We evaluated targeted sequencing of oral rinse samples (ORS) for non-invasive mutational profiling of oral squamous cell carcinomas (OSCC).

MATERIALS AND METHODS

A custom hybrid capture panel targeting 42 frequently mutated genes in OSCC was used to identify DNA sequence variants in matched ORS and fresh-frozen tumors from 120 newly-diagnosed patients. Receiver operating characteristic (ROC) curves determined the optimal variant allele fraction (VAF) cutoff for variant discrimination in ORS. Behavioral, clinical, and analytical factors were evaluated for impacts on assay performance.

RESULTS

Half of tumors involved oral tongue (50 %), and a majority were T1-T2 tumor stage (55 %). Median depth of sequencing coverage was 260X for OSCC and 1,563X for ORS. Frequencies of single nucleotide variants (SNVs) at highly mutated genes (including TP53, FAT1, HRAS, NOTCH1, CDKN2A, CASP8, NFE2L2, and PIK3CA) in OSCC were highly correlated with TCGA data (R = 0.96, p = 2.5E-22). An ROC curve with area-under-the-curve (AUC) of 0.80 showed that, at an optimal VAF cutoff of 0.10 %, ORS provided 76 % sensitivity, 96 % specificity, but precision of only 2.6E-4. At this VAF cutoff, 206 of 270 SNVs in OSCC were detected in matched ORS. Sensitivity varied by patient, T stage and target gene. Neither downsampled ORS as matched control nor a naïve Bayesian classifier adjusting for sequencing bias appreciably improved assay performance.

CONCLUSION

Targeted sequencing of ORS provides moderate assay performance for noninvasive detection of SNVs in OSCC. Our findings strongly rationalize further clinical and laboratory optimization of this assay, including strategies to improve precision.

Oncogenic CDK13 mutations impede nuclear RNA surveillance

RNA surveillance pathways detect and degrade defective transcripts to ensure RNA fidelity. We found that disrupted nuclear RNA surveillance is oncogenic. Cyclin-dependent kinase 13 (CDK13) is mutated in melanoma, and patient-mutated CDK13 accelerates zebrafish melanoma. CDK13 mutation causes aberrant RNA stabilization. CDK13 is required for ZC3H14 phosphorylation, which is necessary and sufficient to promote nuclear RNA degradation. Mutant CDK13 fails to activate nuclear RNA surveillance, causing aberrant protein-coding transcripts to be stabilized and translated. Forced aberrant RNA expression accelerates melanoma in zebrafish. We found recurrent mutations in genes encoding nuclear RNA surveillance components in many malignancies, establishing nuclear RNA surveillance as a tumor-suppressive pathway. Activating nuclear RNA surveillance is crucial to avoid accumulation of aberrant RNAs and their ensuing consequences in development and disease.

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.

WMDS.net: a network control framework for identifying key players in transcriptome programs

MOTIVATION

Mammalian cells can be transcriptionally reprogramed to other cellular phenotypes. Controllability of such complex transitions in transcriptional networks underlying cellular phenotypes is an inherent biological characteristic. This network controllability can be interpreted by operating a few key regulators to guide the transcriptional program from one state to another. Finding the key regulators in the transcriptional program can provide key insights into the network state transition underlying cellular phenotypes.

RESULTS

To address this challenge, here, we proposed to identify the key regulators in the transcriptional co-expression network as a minimum dominating set (MDS) of driver nodes that can fully control the network state transition. Based on the theory of structural controllability, we developed a weighted MDS network model (WMDS.net) to find the driver nodes of differential gene co-expression networks. The weight of WMDS.net integrates the degree of nodes in the network and the significance of gene co-expression difference between two physiological states into the measurement of node controllability of the transcriptional network. To confirm its validity, we applied WMDS.net to the discovery of cancer driver genes in RNA-seq datasets from The Cancer Genome Atlas. WMDS.net is powerful among various cancer datasets and outperformed the other top-tier tools with a better balance between precision and recall.

AVAILABILITY AND IMPLEMENTATION

https://github.com/chaofen123/WMDS.net.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A network-based method for identifying cancer driver genes based on node control centrality

Cancer is one of the major contributors to human mortality and has a serious influence on human survival and health. In biomedical research, the identification of cancer driver genes (cancer drivers for short) is an important task; cancer drivers can promote the progression and generation of cancer. To identify cancer drivers, many methods have been developed. These computational models only identify coding cancer drivers; however, non-coding drivers likewise play significant roles in the progression of cancer. Hence, we propose a Network-based Method for identifying cancer Driver Genes based on node Control Centrality (NMDGCC), which can identify coding and non-coding cancer driver genes. The process of NMDGCC for identifying driver genes mainly includes the following two steps. In the first step, we construct a gene interaction network by using mRNAs and miRNAs expression data in the cancer state. In the second step, the control centrality of the node is used to identify cancer drivers in the constructed network. We use the breast cancer dataset from The Cancer Genome Atlas (TCGA) to verify the effectiveness of NMDGCC. Compared with the existing methods of cancer driver genes identification, NMDGCC has a better performance. NMDGCC also identifies 295 miRNAs as non-coding cancer drivers, of which 158 are related to tumorigenesis of BRCA. We also apply NMDGCC to identify driver genes related to the different breast cancer subtypes. The result shows that NMDGCC detects many cancer drivers of specific cancer subtypes.

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.

Genome interpretation using in silico predictors of variant impact

Estimating the effects of variants found in disease driver genes opens the door to personalized therapeutic opportunities. Clinical associations and laboratory experiments can only characterize a tiny fraction of all the available variants, leaving the majority as variants of unknown significance (VUS). In silico methods bridge this gap by providing instant estimates on a large scale, most often based on the numerous genetic differences between species. Despite concerns that these methods may lack reliability in individual subjects, their numerous practical applications over cohorts suggest they are already helpful and have a role to play in genome interpretation when used at the proper scale and context. In this review, we aim to gain insights into the training and validation of these variant effect predicting methods and illustrate representative types of experimental and clinical applications. Objective performance assessments using various datasets that are not yet published indicate the strengths and limitations of each method. These show that cautious use of in silico variant impact predictors is essential for addressing genome interpretation challenges.

Identification of key somatic oncogenic mutation based on a confounder-free causal inference model

Abnormal cell proliferation and epithelial-mesenchymal transition (EMT) are the essential events that induce cancer initiation and progression. A fundamental goal in cancer research is to develop an efficient method to detect mutational genes capable of driving cancer. Although several computational methods have been proposed to identify these key mutations, many of them focus on the association between genetic mutations and functional changes in relevant biological processes, but not their real causality. Causal effect inference provides a way to estimate the real induce effect of a certain mutation on vital biological processes of cancer initiation and progression, through addressing the confounder bias due to neutral mutations and unobserved latent variables. In this study, integrating genomic and transcriptomic data, we construct a novel causal inference model based on a deep variational autoencoder to identify key oncogenic somatic mutations. Applied to 10 cancer types, our method quantifies the causal effect of genetic mutations on cell proliferation and EMT by reducing both observed and unobserved confounding biases. The experimental results indicate that genes with higher mutation frequency do not necessarily mean they are more potent in inducing cancer and promoting cancer development. Moreover, our study fills a gap in the use of machine learning for causal inference to identify oncogenic mutations.

Identifying key mutations of cancers is helpful to better understand the mechanisms of cancer cell transformation and is critical for therapeutic approaches. Besides sequence and structure-based computational approaches, some functional impact-based methods which consider the association between mutation events and the activity of cancer-related biological processes have also been developed to detect key mutations. However, these methods mainly consider the correlation but ignore that the correlation is far from causality due to the existence of observed and unobserved confounding factors. We develop a confounder-free machine learning-based causal inference framework to estimate the causal effect of mutations on abnormal cell proliferation and epithelial-mesenchymal transition (EMT). It fills a gap in the use of causal mechanisms to discover potential driver mutations in cancer biological systems. Applying our method to 10 cancer types, the identified key mutations are highly consistent with public well-verified ones. Additionally, some new key mutations have also been discovered.

Uncovering Oncogenic Mechanisms of Tumor Suppressor Genes in Breast Cancer Multi-Omics Data

Tumor suppressor genes (TSGs) are essential genes in the development of cancer. While they have many roles in normal cells, mutation and dysregulation of the TSGs result in aberrant molecular processes in cancer cells. Therefore, understanding TSGs and their roles in the oncogenic process is crucial for prevention and treatment of cancer. In this research, multi-omics breast cancer data were used to identify molecular mechanisms of TSGs in breast cancer. Differentially expressed genes and differentially coexpressed genes were identified in four large-scale transcriptomics data from public repositories and multi-omics data analyses of copy number, methylation and gene expression were performed. The results of the analyses were integrated using enrichment analysis and meta-analysis of a p-value summation method. The integrative analysis revealed that TSGs have a significant relationship with genes of gene ontology terms that are related to cell cycle, genome stability, RNA processing and metastasis, indicating the regulatory mechanisms of TSGs on cancer cells. The analysis frame and research results will provide valuable information for the further identification of TSGs in different types of cancers.

Discovering the drivers of clonal hematopoiesis

Mutations in genes that confer a selective advantage to hematopoietic stem cells (HSCs) drive clonal hematopoiesis (CH). While some CH drivers have been identified, the compendium of all genes able to drive CH upon mutations in HSCs remains incomplete. Exploiting signals of positive selection in blood somatic mutations may be an effective way to identify CH driver genes, analogously to cancer. Using the tumor sample in blood/tumor pairs as reference, we identify blood somatic mutations across more than 12,000 donors from two large cancer genomics cohorts. The application of IntOGen, a driver discovery pipeline, to both cohorts, and more than 24,000 targeted sequenced samples yields a list of close to 70 genes with signals of positive selection in CH, available at http://www.intogen.org/ch . This approach recovers known CH genes, and discovers other candidates.

DriverRWH: discovering cancer driver genes by random walk on a gene mutation hypergraph

Background

Recent advances in next-generation sequencing technologies have helped investigators generate massive amounts of cancer genomic data. A critical challenge in cancer genomics is identification of a few cancer driver genes whose mutations cause tumor growth. However, the majority of existing computational approaches underuse the co-occurrence mutation information of the individuals, which are deemed to be important in tumorigenesis and tumor progression, resulting in high rate of false positive.

Results

To make full use of co-mutation information, we present a random walk algorithm referred to as DriverRWH on a weighted gene mutation hypergraph model, using somatic mutation data and molecular interaction network data to prioritize candidate driver genes. Applied to tumor samples of different cancer types from The Cancer Genome Atlas, DriverRWH shows significantly better performance than state-of-art prioritization methods in terms of the area under the curve scores and the cumulative number of known driver genes recovered in top-ranked candidate genes. Besides, DriverRWH discovers several potential drivers, which are enriched in cancer-related pathways. DriverRWH recovers approximately 50% known driver genes in the top 30 ranked candidate genes for more than half of the cancer types. In addition, DriverRWH is also highly robust to perturbations in the mutation data and gene functional network data.

Conclusion

DriverRWH is effective among various cancer types in prioritizes cancer driver genes and provides considerable improvement over other tools with a better balance of precision and sensitivity. It can be a useful tool for detecting potential driver genes and facilitate targeted cancer therapies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04788-7.

The genomic landscape of cholangiocarcinoma reveals the disruption of post-transcriptional modifiers

Molecular variation between geographical populations and subtypes indicate potential genomic heterogeneity and novel genomic features within CCA. Here, we analyze exome-sequencing data of 87 perihilar cholangiocarcinoma (pCCA) and 261 intrahepatic cholangiocarcinoma (iCCA) cases from 3 Asian centers (including 43 pCCAs and 24 iCCAs from our center). iCCA tumours demonstrate a higher tumor mutation burden and copy number alteration burden (CNAB) than pCCA tumours, and high CNAB indicates a poorer pCCA prognosis. We identify 12 significantly mutated genes and 5 focal CNA regions, and demonstrate common mutations in post-transcriptional modification-related potential driver genes METTL14 and RBM10 in pCCA tumours. Finally we demonstrate the tumour-suppressive role of METTL14, a major RNA N6-adenosine methyltransferase (m6A), and illustrate that its loss-of-function mutation R298H may act through m6A modification on potential driver gene MACF1. Our results may be valuable for better understanding of how post-transcriptional modification can affect CCA development, and highlight both similarities and differences between pCCA and iCCA.

Cholangiocarcinoma is a heterogenous group of cancers, with large genetic variation seen within subtypes. Here, the authors find 12 significantly mutated genes and 5 focal CNA regions were found in perihilar cholangiocarcinoma, and identified METTL14 to have a potential tumour suppressive role.

OMEN: Network-based Driver Gene Identification using Mutual Exclusivity

MOTIVATION

Network-based driver identification methods that can exploit mutual exclusivity typically fail to detect rare drivers because of their statistical rigor. Propagation-based methods in contrast allow recovering rare driver genes, but the interplay between network topology and high-scoring nodes often results in spurious predictions. The specificity of driver gene detection can be improved by taking into account both gene-specific and gene-set properties. Combining these requires a formalism that can adjust gene-set properties depending on the exact network context within which a gene is analyzed.

RESULTS

We developed OMEN: a logic programming framework based on random walk semantics. OMEN presents a number of novel concepts. In particular, its design is unique in that it presents an effective approach to combine both gene-specific driver properties and gene-set properties, and includes a novel method to avoid restrictive, a priori filtering of genes by exploiting the gene-set property of mutual exclusivity, expressed in terms of the functional impact scores of mutations, rather than in terms of simple binary mutation calls. Applying OMEN to a benchmark data set derived from TCGA illustrates how OMEN is able to robustly identify driver genes and modules of driver genes as proxies of driver pathways.

AVAILABILITY

The source code is freely available for download at www.github.com/DriesVanDaele/OMEN The data set is archived at https://doi.org/10.5281/zenodo.6419097 and the code at https://doi.org/10.5281/zenodo.6419764.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A general calculus of fitness landscapes finds genes under selection in cancers

Genetic variants drive the evolution of traits and diseases. We previously modeled these variants as small displacements in fitness landscapes and estimated their functional impact by differentiating the evolutionary relationship between genotype and phenotype. Conversely, here we integrate these derivatives to identify genes steering specific traits. Over cancer cohorts, integration identified 460 likely tumor-driving genes. Many have literature and experimental support but had eluded prior genomic searches for positive selection in tumors. Beyond providing cancer insights, these results introduce a general calculus of evolution to quantify the genotype-phenotype relationship and discover genes associated with complex traits and diseases.

Comprehensive profiling of 1015 patients' exomes reveals genomic-clinical associations in colorectal cancer

The genetic basis of colorectal cancer (CRC) and its clinical associations remain poorly understood due to limited samples or targeted genes in current studies. Here, we perform ultradeep whole-exome sequencing on 1015 patients with CRC as part of the ChangKang Project. We identify 46 high-confident significantly mutated genes, 8 of which mutate in 14.9% of patients: LYST, DAPK1, CR2, KIF16B, NPIPB15, SYTL2, ZNF91, and KIAA0586. With an unsupervised clustering algorithm, we propose a subtyping strategy that classisfies CRC patients into four genomic subtypes with distinct clinical characteristics, including hypermutated, chromosome instability with high risk, chromosome instability with low risk, and genome stability. Analysis of immunogenicity uncover the association of immunogenicity reduction with genomic subtypes and poor prognosis in CRC. Moreover, we find that mitochondrial DNA copy number is an independent factor for predicting the survival outcome of CRCs. Overall, our results provide CRC-related molecular features for clinical practice and a valuable resource for translational research.

Effects of Multi-Omics Characteristics on Identification of Driver Genes Using Machine Learning Algorithms

Cancer is a complex disease caused by genomic and epigenetic alterations; hence, identifying meaningful cancer drivers is an important and challenging task. Most studies have detected cancer drivers with mutated traits, while few studies consider multiple omics characteristics as important factors. In this study, we present a framework to analyze the effects of multi-omics characteristics on the identification of driver genes. We utilize four machine learning algorithms within this framework to detect cancer driver genes in pan-cancer data, including 75 characteristics among 19,636 genes. The 75 features are divided into four types and analyzed using Kullback–Leibler divergence based on CGC genes and non-CGC genes. We detect cancer driver genes in two different ways. One is to detect driver genes from a single feature type, while the other is from the top N features. The first analysis denotes that the mutational features are the best characteristics. The second analysis reveals that the top 45 features are the most effective feature combinations and superior to the mutational features. The top 45 features not only contain mutational features but also three other types of features. Therefore, our study extends the detection of cancer driver genes and provides a more comprehensive understanding of cancer mechanisms.

EPIMUTESTR: a nearest neighbor machine learning approach to predict cancer driver genes from the evolutionary action of coding variants

Discovering rare cancer driver genes is difficult because their mutational frequency is too low for statistical detection by computational methods. EPIMUTESTR is an integrative nearest-neighbor machine learning algorithm that identifies such marginal genes by modeling the fitness of their mutations with the phylogenetic Evolutionary Action (EA) score. Over cohorts of sequenced patients from The Cancer Genome Atlas representing 33 tumor types, EPIMUTESTR detected 214 previously inferred cancer driver genes and 137 new candidates never identified computationally before of which seven genes are supported in the COSMIC Cancer Gene Census. EPIMUTESTR achieved better robustness and specificity than existing methods in a number of benchmark methods and datasets.

Prognostic Pathways Guide Drug Indications in Pan-Cancers

Pathway-level analysis is a powerful approach enabling the interpretation of post-genomic data at a higher level than that of individual molecules. Molecular-targeted therapy focusing on cascade signaling pathways has become a new paradigm in anticancer therapy, instead of a single protein. However, the approaches to narrowing down the long list of biological pathways are limited. Here, we proposed a strategy for in silico Drug Prescription on biological pathways across pan-Cancers (CDP), by connecting drugs to candidate pathways. Applying on a list of 120 traditional Chinese medicines (TCM), we especially identified the “TCM–pathways–cancers” triplet and constructed it into a heterogeneous network across pan-cancers. Applying them into TCMs, the computational prescribing methods deepened the understanding of the efficacy of TCM at the molecular level. Further applying them into Western medicines, CDP could promote drug reposition avoiding time-consuming developments of new drugs.

DGRanker: Cancer Driver Gene Detection in Human Transcriptional Regulatory Network

BACKGROUND

Cancer is a group of diseases that have received much attention in biological research because of its high mortality rate and the lack of accurate identification of its root causes. In such studies, researchers usually try to identify cancer driver genes (CDGs) that start cancer in a cell. The majority of the methods that have ever been proposed for the identification of CDGs are based on gene expression data and the concept of mutation in genomic data. Recently, using networking techniques and the concept of influence maximization, some models have been proposed to identify these genes.

OBJECTIVES

We aimed to construct the cancer transcriptional regulatory network and identify cancer driver genes using a network science approach without the use of mutation and genomic data.

MATERIALS AND METHODS

In this study, we will employ the social influence network theory to identify CDGs in the human gene regulatory network (GRN) that is based on the concept of influence and power of webpages. First, we will create GRN Networks using gene expression data and Existing nodes and edges. Next, we will implement the modified algorithm on GRN networks being studied by weighting the regulatory interaction edges using the influence spread concept. Nodes with the highest ratings will be selected as the CDGs.

RESULTS

The results show our proposed method outperforms most of the other computational and network-based methods and show its superiority in identifying CDGs compared to many other methods. In addition, the proposed method can identify many CDGs that are overlooked by all previously published methods.

CONCLUSIONS

Our study demonstrated that the Google's PageRank algorithm can be utilized and modified as a network-based method for identifying cancer driver gene in transcriptional regulatory network. Furthermore, the proposed method can be considered as a complementary method to the computational-based cancer driver gene identification tools.

Machine learning methods for prediction of cancer driver genes: a survey paper

Identifying the genes and mutations that drive the emergence of tumors is a critical step to improving our understanding of cancer and identifying new directions for disease diagnosis and treatment. Despite the large volume of genomics data, the precise detection of driver mutations and their carrying genes, known as cancer driver genes, from the millions of possible somatic mutations remains a challenge. Computational methods play an increasingly important role in discovering genomic patterns associated with cancer drivers and developing predictive models to identify these elements. Machine learning (ML), including deep learning, has been the engine behind many of these efforts and provides excellent opportunities for tackling remaining gaps in the field. Thus, this survey aims to perform a comprehensive analysis of ML-based computational approaches to identify cancer driver mutations and genes, providing an integrated, panoramic view of the broad data and algorithmic landscape within this scientific problem. We discuss how the interactions among data types and ML algorithms have been explored in previous solutions and outline current analytical limitations that deserve further attention from the scientific community. We hope that by helping readers become more familiar with significant developments in the field brought by ML, we may inspire new researchers to address open problems and advance our knowledge towards cancer driver discovery.

Comprehensive evaluation of computational methods for predicting cancer driver genes

Optimal methods could effectively improve the accuracy of predicting and identifying candidate driver genes. Various computational methods based on mutational frequency, network and function approaches have been developed to identify mutation driver genes in cancer genomes. However, a comprehensive evaluation of the performance levels of network-, function- and frequency-based methods is lacking. In the present study, we assessed and compared eight performance criteria for eight network-based, one function-based and three frequency-based algorithms using eight benchmark datasets. Under different conditions, the performance of approaches varied in terms of network, measurement and sample size. The frequency-based driverMAPS and network-based HotNet2 methods showed the best overall performance. Network-based algorithms using protein-protein interaction networks outperformed the function- and the frequency-based approaches. Precision, F1 score and Matthews correlation coefficient were low for most approaches. Thus, most of these algorithms require stringent cutoffs to correctly distinguish driver and non-driver genes. We constructed a website named Cancer Driver Catalog (http://159.226.67.237/sun/cancer_driver/), wherein we integrated the gene scores predicted by the foregoing software programs. This resource provides valuable guidance for cancer researchers and clinical oncologists prioritizing cancer driver gene candidates by using an optimal tool.

Systematic illumination of druggable genes in cancer genomes

By combining 6 druggable genome resources, we identify 6,083 genes as potential druggable genes (PDGs). We characterize their expression, recurrent genomic alterations, cancer dependencies, and therapeutic potentials by integrating genome, functionome, and druggome profiles across cancers. 81.5% of PDGs are reliably expressed in major adult cancers, 46.9% show selective expression patterns, and 39.1% exhibit at least one recurrent genomic alteration. We annotate a total of 784 PDGs as dependent genes for cancer cell growth. We further quantify 16 cancer-related features and estimate a PDG cancer drug target score (PCDT score). PDGs with higher PCDT scores are significantly enriched for genes encoding kinases and histone modification enzymes. Importantly, we find that a considerable portion of high PCDT score PDGs are understudied genes, providing unexplored opportunities for drug development in oncology. By integrating the druggable genome and the cancer genome, our study thus generates a comprehensive blueprint of potential druggable genes across cancers.

Jiang et al. generate a comprehensive blueprint of potential druggable genes (PDGs) across cancers by a systematic integration of the druggable genome and the cancer genome. This resource is publicly available to the cancer research community in The Cancer Druggable Gene Atlas (TCDA) through the Functional Cancer Genome data portal.

A network-based method for detecting cancer driver gene in transcriptional regulatory networks using the structure analysis of weighted regulatory interactions

Background:

The identification of genes that instigate cell anomalies and cause cancer in humans is an important field in oncology research. Abnormalities in these genes are transferred to other genes in the cell, disrupting its normal functionality. Such genes are known as cancer driver genes (CDGs). Various methods have been proposed for predicting CDGs, most of which are based on genomic data and computational methods. Some novel bioinformatic approaches have been developed.

Objective:

In this article, we propose a network-based algorithm, SalsaDriver (Stochastic approach for link-structure analysis to driver detection), which can calculate the receiving and influencing power of each gene using the stochastic analysis of regulatory interaction structures in gene regulatory networks.

Method:

First, regulatory networks related to breast, colon, and lung cancers were constructed using gene expression data and a list of regulatory interactions, the weights of which were then calculated using biological and topological features of the network. After that, the weighted regulatory interactions were used in the structure analysis of interactions achieved using two separate Markov chains on the bipartite graph taken from the main graph of the gene network and implementing the stochastic approach for link-structure analysis. The proposed algorithm categorizes higher-ranked genes as driver genes.

Results:

The proposed algorithm was compared with 24 other computational and network tools based on the F-measure value and the number of detected CDGs. The results were validated using four valid databases. The findings of this study show that SalsaDriver outperforms other methods and can identify a significant number of driver genes not identified using other methods.

Conclusion:

The SalsaDriver network-based approach is suitable for predicting CDGs and can be used as a complementary method along with other computational tools.

Using CRISPR-Cas9 to Dissect Cancer Mutations in Cell Lines

The CRISPR-Cas9 technology has revolutionized the scope and pace of biomedical research, enabling the targeting of specific genomic sequences for a wide spectrum of applications. Here we describe assays to functionally interrogate mutations identified in cancer cells utilizing both CRISPR-Cas9 nuclease and base editors. We provide guidelines to interrogate known cancer driver mutations or functionally screen for novel vulnerability mutations with these systems in characterized human cancer cell lines. The proposed platform should be transferable to primary cancer cells, opening up a path for precision oncology on a functional level.

Network-Based Analysis to Identify Drivers of Metastatic Prostate Cancer Using GoNetic

Most known driver genes of metastatic prostate cancer are frequently mutated. To dig into the long tail of rarely mutated drivers, we performed network-based driver identification on the Hartwig Medical Foundation metastatic prostate cancer data set (HMF cohort). Hereto, we developed GoNetic, a method based on probabilistic pathfinding, to identify recurrently mutated subnetworks. In contrast to most state-of-the-art network-based methods, GoNetic can leverage sample-specific mutational information and the weights of the underlying prior network. When applied to the HMF cohort, GoNetic successfully recovered known primary and metastatic drivers of prostate cancer that are frequently mutated in the HMF cohort (TP53, RB1, and CTNNB1). In addition, the identified subnetworks contain frequently mutated genes, reflect processes related to metastatic prostate cancer, and contain rarely mutated driver candidates. To further validate these rarely mutated genes, we assessed whether the identified genes were more mutated in metastatic than in primary samples using an independent cohort. Then we evaluated their association with tumor evolution and with the lymph node status of the patients. This resulted in forwarding several novel putative driver genes for metastatic prostate cancer, some of which might be prognostic for disease evolution.

Genomic Analysis Reveals Heterogeneity Between Lesions in Synchronous Primary Right-Sided and Left-Sided Colon Cancer

Background: The synchronous primary right-sided and left-sided colon cancer (sRL-CC) is a peculiar subtype of colorectal cancer. However, the genomic landscape of sRL-CC remains elusive. Methods: Twenty-eight paired tumor samples and their corresponding normal mucosa samples from 14 patients were collected from the Second Affiliated Hospital of Harbin Medical University from 2011 to 2018. The clinical-pathological data were obtained, and whole-exome sequencing was performed based on formalin-fixed and paraffin-embedded samples of these patients, and then, comprehensive bioinformatic analyses were conducted. Results: Both the lesions of sRL-CC presented dissimilar histological grade and differentiation. Based on sequencing data, few overlapping SNV signatures, onco-driver gene mutations, and SMGs were identified. Moreover, the paired lesions harbored a different distribution of copy number variants (CNVs) and loss of heterozygosity. The clonal architecture analysis demonstrated the polyclonal origin of sRL-CC and inter-cancerous heterogeneity between two lesions. Conclusion: Our work provides evidence that lesions of sRL-CC share few overlapping mutational signatures and CNVs, and may originate from different clones.

The impact of gene-body H3K36me3 patterns on gene expression level changes in chronic myelogenous leukemia

Chronic myelogenous leukemia (CML) is a malignant clonal disease of hematopoietic stem cells. Researches have exhibited that the progression of CML is related to histone modifications. Here, we perform the systematic analyses of H3K36me3 patterns and gene expression level changes. We observe that the genes with higher gene-body H3K36me3 levels in normal cells show fewer expression changes during leukemogenesis, while the genes with lower gene-body H3K36me3 levels in normal cells yield obvious expression changes during leukemogenesis (ρ = -0.98, P = 9.30 × 10^-8). These findings are conserved in human lung/breast cancers and mouse CML, regardless of gene expression levels and gene lengths. Regulatory element analysis and Random Forest regression display that Hoxd13, Rara, Scl, Smad3, Smad4 and Tgif1 induce the up-regulation of genes with lower H3K36me3 levels (ρ = 0.97, P = 2.35 × 10^-56). Enrichment analysis shows that the differentially expressed genes with lower H3K36me3 levels are involved in leukemia-related pathways, such as leukocyte migration and regulation of leukocyte activation. Finally, six driver genes (Tp53, Wt1, Dnmt3a, Cacna1b, Phactr1 and Gbp4) with lower H3K36me3 levels are identified. Our analyses indicate that lower gene-body H3K36me3 levels may serve as a biomarker for the progression of CML.

A systematic view of computational methods for identifying driver genes based on somatic mutation data

Abnormal changes of driver genes are serious for human health and biomedical research. Identifying driver genes, exactly from enormous genes with mutations, promotes accurate diagnosis and treatment of cancer. A lot of works about uncovering driver genes have been developed over the past decades. By analyzing previous works, we find that computational methods are more efficient than traditional biological experiments when distinguishing driver genes from massive data. In this study, we summarize eight common computational algorithms only using somatic mutation data. We first group these methods into three categories according to mutation features they apply. Then, we conclude a general process of nominating candidate cancer driver genes. Finally, we evaluate three representative methods on 10 kinds of cancer derived from The Cancer Genome Atlas Program and five Chinese projects from the International Cancer Genome Consortium. In addition, we compare results of methods with various parameters. Evaluation is performed from four perspectives, including CGC, OG/TSG, Q-value and QQQuantile-Quantileplot. To sum up, we present algorithms using somatic mutation data in order to offer a systematic view of various mutation features and lay the foundation of methods based on integration of mutation information and other types of data.

Prospective Validation of the Italian Alliance Against Cancer Lung Panel in Patients With Advanced Non-Small-Cell Lung Cancer

BACKGROUND

The deeper knowledge of non-small-cell lung cancer (NSCLC) biology and the discovery of driver molecular alterations have opened the era of precision medicine in lung oncology, thus significantly revolutionizing the diagnostic and therapeutic approach to NSCLC. In Italy, however, molecular assessment remains heterogeneous across the country, and numbers of patients accessing personalized treatments remain relatively low. Nationwide programs have demonstrated that the creation of consortia represent a successful strategy to increase the number of patients with a molecular classification.

PATIENTS AND METHODS

The Alliance Against Cancer (ACC), a network of 25 Italian Research Institutes, has developed a targeted sequencing panel for the detection of genomic alterations in 182 genes in patients with a diagnosis of NSCLC (ACC lung panel). One thousand metastatic NSCLC patients will be enrolled onto a prospective trial designed to measure the sensitivity and specificity of the ACC lung panel as a tool for molecular screening compared to standard methods.

RESULTS AND CONCLUSION

The ongoing trial is part of a nationwide strategy of ACC to develop infrastructures and improve competences to make the Italian research institutes independent for genomic profiling of cancer patients.

EARN: an ensemble machine learning algorithm to predict driver genes in metastatic breast cancer

BACKGROUND

Today, there are a lot of markers on the prognosis and diagnosis of complex diseases such as primary breast cancer. However, our understanding of the drivers that influence cancer aggression is limited.

METHODS

In this work, we study somatic mutation data consists of 450 metastatic breast tumor samples from cBio Cancer Genomics Portal. We use four software tools to extract features from this data. Then, an ensemble classifier (EC) learning algorithm called EARN (Ensemble of Artificial Neural Network, Random Forest, and non-linear Support Vector Machine) is proposed to evaluate plausible driver genes for metastatic breast cancer (MBCA). The decision-making strategy for the proposed ensemble machine is based on the aggregation of the predicted scores obtained from individual learning classifiers to be prioritized homo sapiens genes annotated as protein-coding from NCBI.

RESULTS

This study is an attempt to focus on the findings in several aspects of MBCA prognosis and diagnosis. First, drivers and passengers predicted by SVM, ANN, RF, and EARN are introduced. Second, biological inferences of predictions are discussed based on gene set enrichment analysis. Third, statistical validation and comparison of all learning methods are performed by some evaluation metrics. Finally, the pathway enrichment analysis (PEA) using ReactomeFIVIz tool (FDR < 0.03) for the top 100 genes predicted by EARN leads us to propose a new gene set panel for MBCA. It includes HDAC3, ABAT, GRIN1, PLCB1, and KPNA2 as well as NCOR1, TBL1XR1, SIRT4, KRAS, CACNA1E, PRKCG, GPS2, SIN3A, ACTB, KDM6B, and PRMT1. Furthermore, we compare results for MBCA to other outputs regarding 983 primary tumor samples of breast invasive carcinoma (BRCA) obtained from the Cancer Genome Atlas (TCGA). The comparison between outputs shows that ROC-AUC reaches 99.24% using EARN for MBCA and 99.79% for BRCA. This statistical result is better than three individual classifiers in each case.

CONCLUSIONS

This research using an integrative approach assists precision oncologists to design compact targeted panels that eliminate the need for whole-genome/exome sequencing. The schematic representation of the proposed model is presented as the Graphic abstract.

Systematic Characterization of Recurrent Genomic Alterations in Cyclin-Dependent Kinases Reveals Potential Therapeutic Strategies for Cancer Treatment

Recurrent copy-number alterations, mutations, and transcript fusions of the genes encoding CDKs/cyclins are characterized in >10,000 tumors. Genomic alterations of CDKs/cyclins are dominantly driven by copy number aberrations. In contrast to cell-cycle-related CDKs/cyclins, which are globally amplified, transcriptional CDKs/cyclins recurrently lose copy numbers across cancers. Although mutations and transcript fusions are relatively rare events, CDK12 exhibits recurrent mutations in multiple cancers. Among the transcriptional CDKs, CDK7 and CDK12 show the most significant copy number loss and mutation, respectively. Their genomic alterations are correlated with increased sensitivities to DNA-damaging drugs. Inhibition of CDK7 preferentially represses the expression of genes in the DNA-damage-repair pathways and impairs the activity of homologous recombination. Low-dose CDK7 inhibitor treatment sensitizes cancer cells to PARP inhibitor-induced DNA damage and cell death. Our analysis provides genomic information for identification and prioritization of drug targets for CDKs and reveals rationales for treatment strategies.

pDriver : A novel method for unravelling personalised coding and miRNA cancer drivers

MOTIVATION

Unravelling cancer driver genes is important in cancer research. Although computational methods have been developed to identify cancer drivers, most of them detect cancer drivers at population level. However, two patients who have the same cancer type and receive the same treatment may have different outcomes because each patient has a different genome and their disease might be driven by different driver genes. Therefore new methods are being developed for discovering cancer drivers at individual level, but existing personalised methods only focus on coding drivers while microRNAs (miRNAs) have been shown to drive cancer progression as well. Thus, novel methods are required to discover both coding and miRNA cancer drivers at individual level.

RESULTS

We propose the novel method, pDriver, to discover personalised cancer drivers. pDriver includes two stages: (1) Constructing gene networks for each cancer patient and (2) Discovering cancer drivers for each patient based on the constructed gene networks. To demonstrate the effectiveness of pDriver, we have applied it to five TCGA cancer datasets and compared it with the state-of-the-art methods. The result indicates that pDriver is more effective than other methods. Furthermore, pDriver can also detect miRNA cancer drivers and most of them have been confirmed to be associated with cancer by literature. We further analyse the predicted personalised drivers for breast cancer patients and the result shows that they are significantly enriched in many GO processes and KEGG pathways involved in breast cancer.

AVAILABILITY AND IMPLEMENTATION

pDriver is available at https://github.com/pvvhoang/pDriver.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.