Identifying driver mutations in sequenced cancer genomes: computational approaches to enable precision medicine

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.

ECD-CDGI: An efficient energy-constrained diffusion model for cancer driver gene identification

The identification of cancer driver genes (CDGs) poses challenges due to the intricate interdependencies among genes and the influence of measurement errors and noise. We propose a novel energy-constrained diffusion (ECD)-based model for identifying CDGs, termed ECD-CDGI. This model is the first to design an ECD-Attention encoder by combining the ECD technique with an attention mechanism. ECD-Attention encoder excels at generating robust gene representations that reveal the complex interdependencies among genes while reducing the impact of data noise. We concatenate topological embedding extracted from gene-gene networks through graph transformers to these gene representations. We conduct extensive experiments across three testing scenarios. Extensive experiments show that the ECD-CDGI model possesses the ability to not only be proficient in identifying known CDGs but also efficiently uncover unknown potential CDGs. Furthermore, compared to the GNN-based approach, the ECD-CDGI model exhibits fewer constraints by existing gene-gene networks, thereby enhancing its capability to identify CDGs. Additionally, ECD-CDGI is open-source and freely available. We have also launched the model as a complimentary online tool specifically crafted to expedite research efforts focused on CDGs identification.

Moving toward precision medicine to predict drug sensitivity in patients with metastatic breast cancer

Tumor heterogeneity represents a major challenge in breast cancer, being associated with disease progression and treatment resistance. Precision medicine has been extensively applied to dissect tumor heterogeneity and, through a deeper molecular understanding of the disease, to personalize therapeutic strategies. In the last years, technological advances have widely improved the understanding of breast cancer biology and several trials have been developed to translate these new insights into clinical practice, with the ultimate aim of improving patients' outcomes. In the era of molecular oncology, genomics analyses and other methodologies are shaping a new treatment algorithm in breast cancer care. In this manuscript, we review the main steps of precision medicine to predict drug sensitivity in breast cancer from a translational point of view. Genomic developments and their clinical implications are discussed, along with technological advancements that could broaden precision medicine applications. Current achievements are put into perspective to provide an overview of the state-of-art of breast cancer precision oncology as well as to identify future research directions.

Variant Allele Frequency Analysis of Circulating Tumor DNA as a Promising Tool in Assessing the Effectiveness of Treatment in Non-Small Cell Lung Carcinoma Patients

Despite the different possible paths of treatment, lung cancer remains one of the leading causes of death in oncological patients. New tools guiding the therapeutic process are under scientific investigation, and one of the promising indicators of the effectiveness of therapy in patients with NSCLC is variant allele frequency (VAF) analysis. VAF is a metric characterized as the measurement of the specific variant allele proportion within a genomic locus, and it can be determined using methods based on NGS or PCR. It can be assessed using not only tissue samples but also ctDNA (circulating tumor DNA) isolated from liquid biopsy. The non-invasive characteristic of liquid biopsy enables a more frequent collection of material and increases the potential of VAF analysis in monitoring therapy. Several studies have been performed on patients with NSCLC to evaluate the possibility of VAF usage. The research carried out so far demonstrates that the evaluation of VAF dynamics may be useful in monitoring tumor progression, remission, and recurrence during or after treatment. Moreover, the use of VAF analysis appears to be beneficial in making treatment decisions. However, several issues require better understanding and standardization before VAF testing can be implemented in clinical practice. In this review, we discuss the difficulties in the application of ctDNA VAF analysis in clinical routine, discussing the diagnostic and methodological challenges in VAF measurement in liquid biopsy. We highlight the possible applications of VAF-based measurements that are under consideration in clinical trials in the monitoring of personalized treatments for patients with NSCLC.

The Extent of Edgetic Perturbations in the Human Interactome Caused by Population-Specific Mutations

Until recently, efforts in population genetics have been focused primarily on people of European ancestry. To attenuate this bias, global population studies, such as the 1000 Genomes Project, have revealed differences in genetic variation across ethnic groups. How many of these differences can be attributed to population-specific traits? To answer this question, the mutation data must be linked with functional outcomes. A new "edgotype" concept has been proposed, which emphasizes the interaction-specific, "edgetic", perturbations caused by mutations in the interacting proteins. In this work, we performed systematic in silico edgetic profiling of ~50,000 non-synonymous SNVs (nsSNVs) from the 1000 Genomes Project by leveraging our semi-supervised learning approach SNP-IN tool on a comprehensive set of over 10,000 protein interaction complexes. We interrogated the functional roles of the variants and their impact on the human interactome and compared the results with the pathogenic variants disrupting PPIs in the same interactome. Our results demonstrated that a considerable number of nsSNVs from healthy populations could rewire the interactome. We also showed that the proteins enriched with interaction-disrupting mutations were associated with diverse functions and had implications in a broad spectrum of diseases. Further analysis indicated that distinct gene edgetic profiles among major populations could shed light on the molecular mechanisms behind the population phenotypic variances. Finally, the network analysis revealed that the disease-associated modules surprisingly harbored a higher density of interaction-disrupting mutations from healthy populations. The variation in the cumulative network damage within these modules could potentially account for the observed disparities in disease susceptibility, which are distinctly specific to certain populations. Our work demonstrates the feasibility of a large-scale in silico edgetic study, and reveals insights into the orchestrated play of population-specific mutations in the human interactome.

Variant allele frequency: a decision-making tool in precision oncology?

Precision oncology requires additional predictive biomarkers for targeted therapy selection. Variant allele frequency (VAF), measuring the proportion of variant alleles within a genomic locus, provides insights into tumor clonality in somatic genomic testing, yielding a strong rationale for targeting dominant cancer cell populations. The prognostic and predictive roles of VAF have been evaluated across different studies. Yet, the absence of validated VAF thresholds and a lack of standardization between sequencing assays currently hampers its clinical utility. Therefore, analytical and clinical validation must be further examined. This Review summarizes the evidence regarding the use of VAF as a predictive biomarker and discusses challenges and opportunities for its clinical implementation as a decision-making tool for targeted therapy selection.

Comprehensive genomic profiling for oncological advancements by precision medicine

Considerable advancements in next generation sequencing (NGS) techniques have sparked the use of comprehensive genomic profiling (CGP) as a guiding tool for precision-centered oncological treatments. The past two decades have seen the completion of the human genome project, and the consequential invention of NGS. High-throughput sequencing technologies support the discovery and commonplace use of individualized cancer treatments, specifically immune-centered checkpoint inhibitor therapies, and oncogene and tumor suppressor gene targeted therapies. Nevertheless, CGP is not commonly used in all clinical settings. This review investigates the clinically relevant applications of CGP. Studies published between the years 2000-2023 have shown substantial evidence of the benefits of integrating CGP into routine care practice, while also making important comparisons to current-standard oncological treatment strategies. Findings of a comprehensive genomic profile includes predictive, prognostic, and diagnostic biomarkers, together with somatic mutation identification which can indicate the efficacy of immunotherapies and molecularly guided therapies. This review highlights the importance of CGP in identifying driver mutations in tumors that subsequently can be effectively targeted with molecular therapeutics and lead to drug discovery, allowing for increased precision in treating tumors selectively based on their specific genetic mutations, thereby improving patient outcomes.

Adding Value to Liquid Biopsy for Brain Tumors: The Role of Imaging

Clinical management in neuro-oncology has changed to an integrative approach that incorporates molecular profiles alongside histopathology and imaging findings. While the World Health Organization (WHO) guideline recommends the genotyping of informative alterations as a routine clinical practice for central nervous system (CNS) tumors, the acquisition of tumor tissue in the CNS is invasive and not always possible. Liquid biopsy is a non-invasive approach that provides the opportunity to capture the complex molecular heterogeneity of the whole tumor through the detection of circulating tumor biomarkers in body fluids, such as blood or cerebrospinal fluid (CSF). Despite all of the advantages, the low abundance of tumor-derived biomarkers, particularly in CNS tumors, as well as their short half-life has limited the application of liquid biopsy in clinical practice. Thus, it is crucial to identify the factors associated with the presence of these biomarkers and explore possible strategies that can increase the shedding of these tumoral components into biological fluids. In this review, we first describe the clinical applications of liquid biopsy in CNS tumors, including its roles in the early detection of recurrence and monitoring of treatment response. We then discuss the utilization of imaging in identifying the factors that affect the detection of circulating biomarkers as well as how image-guided interventions such as focused ultrasound can help enhance the presence of tumor biomarkers through blood-brain barrier (BBB) disruption.

Toward the functional interpretation of somatic structural variations: bulk- and single-cell approaches

Structural variants (SVs) are genomic rearrangements that can take many different forms such as copy number alterations, inversions and translocations. During cell development and aging, somatic SVs accumulate in the genome with potentially neutral, deleterious or pathological effects. Generation of somatic SVs is a key mutational process in cancer development and progression. Despite their importance, the detection of somatic SVs is challenging, making them less studied than somatic single-nucleotide variants. In this review, we summarize recent advances in whole-genome sequencing (WGS)-based approaches for detecting somatic SVs at the tissue and single-cell levels and discuss their advantages and limitations. First, we describe the state-of-the-art computational algorithms for somatic SV calling using bulk WGS data and compare the performance of somatic SV detectors in the presence or absence of a matched-normal control. We then discuss the unique features of cutting-edge single-cell-based techniques for analyzing somatic SVs. The advantages and disadvantages of bulk and single-cell approaches are highlighted, along with a discussion of their sensitivity to copy-neutral SVs, usefulness for functional inferences and experimental and computational costs. Finally, computational approaches for linking somatic SVs to their functional readouts, such as those obtained from single-cell transcriptome and epigenome analyses, are illustrated, with a discussion of the promise of these approaches in health and diseases.

Identifying Significantly Perturbed Subnetworks in Cancer Using Multiple Protein-Protein Interaction Networks

BACKGROUND

The identification of cancer driver genes and key molecular pathways has been the focus of large-scale cancer genome studies. Network-based methods detect significantly perturbed subnetworks as putative cancer pathways by incorporating genomics data with the topological information of PPI networks. However, commonly used PPI networks have distinct topological structures, making the results of the same method vary widely when applied to different networks. Furthermore, emerging context-specific PPI networks often have incomplete topological structures, which pose serious challenges for existing subnetwork detection algorithms.

METHODS

In this paper, we propose a novel method, referred to as MultiFDRnet, to address the above issues. The basic idea is to model a set of PPI networks as a multiplex network to preserve the topological structure of individual networks, while introducing dependencies among them, and, then, to detect significantly perturbed subnetworks on the modeled multiplex network using all the structural information simultaneously.

RESULTS

To illustrate the effectiveness of the proposed approach, an extensive benchmark analysis was conducted on both simulated and real cancer data. The experimental results showed that the proposed method is able to detect significantly perturbed subnetworks jointly supported by multiple PPI networks and to identify novel modular structures in context-specific PPI networks.

The Cytotoxic Properties of Extreme Fungi's Bioactive Components-An Updated Metabolic and Omics Overview

Fungi are the most diverse living organisms on planet Earth, where their ubiquitous presence in various ecosystems offers vast potential for the research and discovery of new, naturally occurring medicinal products. Concerning human health, cancer remains one of the leading causes of mortality. While extensive research is being conducted on treatments and their efficacy in various stages of cancer, finding cytotoxic drugs that target tumor cells with no/less toxicity toward normal tissue is a significant challenge. In addition, traditional cancer treatments continue to suffer from chemical resistance. Fortunately, the cytotoxic properties of several natural products derived from various microorganisms, including fungi, are now well-established. The current review aims to extract and consolidate the findings of various scientific studies that identified fungi-derived bioactive metabolites with antitumor (anticancer) properties. The antitumor secondary metabolites identified from extremophilic and extremotolerant fungi are grouped according to their biological activity and type. It became evident that the significance of these compounds, with their medicinal properties and their potential application in cancer treatment, is tremendous. Furthermore, the utilization of omics tools, analysis, and genome mining technology to identify the novel metabolites for targeted treatments is discussed. Through this review, we tried to accentuate the invaluable importance of fungi grown in extreme environments and the necessity of innovative research in discovering naturally occurring bioactive compounds for the development of novel cancer treatments.

Role of Surgical Pathologist for the Detection of Immuno-oncologic Predictive Factors in Non-small Cell Lung Cancers

Until very recently, surgery, chemotherapy, and radiation therapy have been the mainstay of treatment in non-small cell carcinomas (NSCLCs). However, recent advances in molecular immunology have unveiled some of the complexity of the mechanisms regulating cellular immune responses and led to the successful targeting of immune checkpoints in attempts to enhance antitumor T-cell responses. Immune checkpoint molecules such as cytotoxic T-lymphocyte associated protein-4, programmed cell death protein-1, and programmed death ligand (PD-L) 1 have been shown to play central roles in evading cancer immunity. Thus, these molecules have been targeted by inhibitors for the management of cancers forming the basis of immunotherapy. Advanced NSCLC has been the paradigm for the benefits of immunotherapy in any cancer. Treatment decisions are made based on the expression of PD-L1 on the tumor cells and the presence or absence of driver mutations. Patients with high PD-L1 expression (≥50%) and no driver mutations are treated with single-agent immunotherapy whereas, for all other patients with a lower level of PD-L1 expression, a combination of chemotherapy and immunotherapy is preferred. Thus, PD-L1 blockers are the only immunotherapeutic agents approved in advanced NSCLC without any oncogenic driver mutations. PD-L1 immunohistochemistry, however, may not be the best biomarker in view of its dynamic nature in time and space, and the benefits may be seen regardless of PD -L1 expression. Each immunotherapy molecule is prescribed based on the levels of PD-L1 expression as assessed by a Food and Drug Administration-approved companion diagnostic assay. Other biomarkers that have been studied include tumor mutational burden, the T-effector signature, tumor-infiltrating lymphocytes, radiomic assays, inflammation index, presence or absence of immune-related adverse events and specific driver mutations, and gut as well as local microbiome. At the current time, none of these biomarkers are routinely used in the clinical decision-making process for immunotherapy in NSCLC. However, in individual cases, they can be useful adjuncts to conventional therapy. This review describes our current understanding of the role of biomarkers as predictors of response to immune checkpoint molecules. To begin with a brief on cancer immunology in general and in NSCLC, in particular, is discussed. In the end, recent advancements in laboratory techniques for refining biomarker assays are described.

Primary refractory plasmablastic lymphoma: A precision oncology approach

Introduction

Hematologic malignancies are currently underrepresented in multidisciplinary molecular-tumor-boards (MTB). This study assesses the potential of precision-oncology in primary-refractory plasmablastic-lymphoma (prPBL), a highly lethal blood cancer.

Methods

We evaluated clinicopathological and molecular-genetic data of 14 clinically annotated prPBL-patients from initial diagnosis. For this proof-of-concept study, we employed our certified institutional MTB-pipeline (University-Cancer-Center-Schleswig-Holstein, UCCSH) to annotate a comprehensive dataset within the scope of a virtual MTB-setting, ultimately recommending molecularly stratified therapies. Evidence-levels for MTB-recommendations were defined in accordance with the NCT/DKTK and ESCAT criteria.

Results

Median age in the cohort was 76.5 years (range 56-91), 78.6% of patients were male, 50% were HIV-positive and clinical outcome was dismal. Comprehensive genomic/transcriptomic analysis revealed potential recommendations of a molecularly stratified treatment option with evidence-levels according to NCT/DKTK of at least m2B/ESCAT of at least IIIA were detected for all 14 prPBL-cases. In addition, immunohistochemical-assessment (CD19/CD30/CD38/CD79B) revealed targeted treatment-recommendations in all 14 cases. Genetic alterations were classified by treatment-baskets proposed by Horak et al. Hereby, we identified tyrosine-kinases (TK; n=4), PI3K-MTOR-AKT-pathway (PAM; n=3), cell-cycle-alterations (CC; n=2), RAF-MEK-ERK-cascade (RME; n=2), immune-evasion (IE; n=2), B-cell-targets (BCT; n=25) and others (OTH; n=4) for targeted treatment-recommendations. The minimum requirement for consideration of a drug within the scope of the study was FDA-fast-track development.

Discussion

The presented proof-of-concept study demonstrates the clinical potential of precision-oncology, even in prPBL-patients. Due to the aggressive course of the disease, there is an urgent medical-need for personalized treatment approaches, and this population should be considered for MTB inclusion at the earliest time.

Alterations of Cytoskeleton Networks in Cell Fate Determination and Cancer Development

Cytoskeleton proteins have been long recognized as structural proteins that provide the necessary mechanical architecture for cell development and tissue homeostasis. With the completion of the cancer genome project, scientists were surprised to learn that huge numbers of mutated genes are annotated as cytoskeletal or associated proteins. Although most of these mutations are considered as passenger mutations during cancer development and evolution, some genes show high mutation rates that can even determine clinical outcomes. In addition, (phospho)proteomics study confirms that many cytoskeleton-associated proteins, e.g., β-catenin, PIK3CA, and MB21D2, are important signaling mediators, further suggesting their biofunctional roles in cancer development. With emerging evidence to indicate the involvement of mechanotransduction in stemness formation and cell differentiation, mutations in these key cytoskeleton components may change the physical/mechanical properties of the cells and determine the cell fate during cancer development. In particular, tumor microenvironment remodeling triggered by such alterations has been known to play important roles in autophagy, metabolism, cancer dormancy, and immune evasion. In this review paper, we will highlight the current understanding of how aberrant cytoskeleton networks affect cancer behaviors and cellular functions through mechanotransduction.

A New View of Activating Mutations in Cancer

A vast effort has been invested in the identification of driver mutations of cancer. However, recent studies and observations call into question whether the activating mutations or the signal strength are the major determinant of tumor development. The data argue that signal strength determines cell fate, not the mutation that initiated it. In addition to activating mutations, factors that can impact signaling strength include (i) homeostatic mechanisms that can block or enhance the signal, (ii) the types and locations of additional mutations, and (iii) the expression levels of specific isoforms of genes and regulators of proteins in the pathway. Because signal levels are largely decided by chromatin structure, they vary across cell types, states, and time windows. A strong activating mutation can be restricted by low expression, whereas a weaker mutation can be strengthened by high expression. Strong signals can be associated with cell proliferation, but too strong a signal may result in oncogene-induced senescence. Beyond cancer, moderate signal strength in embryonic neural cells may be associated with neurodevelopmental disorders, and moderate signals in aging may be associated with neurodegenerative diseases, like Alzheimer's disease. The challenge for improving patient outcomes therefore lies in determining signaling thresholds and predicting signal strength.

The Cancermuts software package for the prioritization of missense cancer variants: a case study of AMBRA1 in melanoma

Cancer genomics and cancer mutation databases have made an available wealth of information about missense mutations found in cancer patient samples. Contextualizing by means of annotation and predicting the effect of amino acid change help identify which ones are more likely to have a pathogenic impact. Those can be validated by means of experimental approaches that assess the impact of protein mutations on the cellular functions or their tumorigenic potential. Here, we propose the integrative bioinformatic approach Cancermuts, implemented as a Python package. Cancermuts is able to gather known missense cancer mutations from databases such as cBioPortal and COSMIC, and annotate them with the pathogenicity score REVEL as well as information on their source. It is also able to add annotations about the protein context these mutations are found in, such as post-translational modification sites, structured/unstructured regions, presence of short linear motifs, and more. We applied Cancermuts to the intrinsically disordered protein AMBRA1, a key regulator of many cellular processes frequently deregulated in cancer. By these means, we classified mutations of AMBRA1 in melanoma, where AMBRA1 is highly mutated and displays a tumor-suppressive role. Next, based on REVEL score, position along the sequence, and their local context, we applied cellular and molecular approaches to validate the predicted pathogenicity of a subset of mutations in an in vitro melanoma model. By doing so, we have identified two AMBRA1 mutations which show enhanced tumorigenic potential and are worth further investigation, highlighting the usefulness of the tool. Cancermuts can be used on any protein targets starting from minimal information, and it is available at https://www.github.com/ELELAB/cancermuts as free software.

Allostery: Allosteric Cancer Drivers and Innovative Allosteric Drugs

Here, we discuss the principles of allosteric activating mutations, propagation downstream of the signals that they prompt, and allosteric drugs, with examples from the Ras signaling network. We focus on Abl kinase where mutations shift the landscape toward the active, imatinib binding-incompetent conformation, likely resulting in the high affinity ATP outcompeting drug binding. Recent pharmacological innovation extends to allosteric inhibitor (GNF-5)-linked PROTAC, targeting Bcr-Abl1 myristoylation site, and broadly, allosteric heterobifunctional degraders that destroy targets, rather than inhibiting them. Designed chemical linkers in bifunctional degraders can connect the allosteric ligand that binds the target protein and the E3 ubiquitin ligase warhead anchor. The physical properties and favored conformational state of the engineered linker can precisely coordinate the distance and orientation between the target and the recruited E3. Allosteric PROTACs, noncompetitive molecular glues, and bitopic ligands, with covalent links of allosteric ligands and orthosteric warheads, increase the effective local concentration of productively oriented and placed ligands. Through covalent chemical or peptide linkers, allosteric drugs can collaborate with competitive drugs, degrader anchors, or other molecules of choice, driving innovative drug discovery.

Novel Driver Strength Index highlights important cancer genes in TCGA PanCanAtlas patients

Background

Cancer driver genes are usually ranked by mutation frequency, which does not necessarily reflect their driver strength. We hypothesize that driver strength is higher for genes preferentially mutated in patients with few driver mutations overall, because these few mutations should be strong enough to initiate cancer.

Methods

We propose formulas for the Driver Strength Index (DSI) and the Normalized Driver Strength Index (NDSI), the latter independent of gene mutation frequency. We validate them using TCGA PanCanAtlas datasets, established driver prediction algorithms and custom computational pipelines integrating SNA, CNA and aneuploidy driver contributions at the patient-level resolution.

Results

DSI and especially NDSI provide substantially different gene rankings compared to the frequency approach. E.g., NDSI prioritized members of specific protein families, including G proteins GNAQ, GNA11 and GNAS, isocitrate dehydrogenases IDH1 and IDH2, and fibroblast growth factor receptors FGFR2 and FGFR3. KEGG analysis shows that top NDSI-ranked genes comprise EGFR/FGFR2/GNAQ/GNA11-NRAS/HRAS/KRAS-BRAF pathway, AKT1-MTOR pathway, and TCEB1-VHL-HIF1A pathway.

Conclusion

Our indices are able to select for driver gene attributes not selected by frequency sorting, potentially for driver strength. Genes and pathways prioritized are likely the strongest contributors to cancer initiation and progression and should become future therapeutic targets.

SYSMut: decoding the functional significance of rare somatic mutations in cancer

Current tailored-therapy efforts in cancer are largely focused on a small number of highly recurrently mutated driver genes but therapeutic targeting of these oncogenes remains challenging. However, the vast number of genes mutated infrequently across cancers has received less attention, in part, due to a lack of understanding of their biological significance. We present SYSMut, an extendable systems biology platform that can robustly infer the biologic consequences of somatic mutations by integrating routine multiomics profiles in primary tumors. We establish SYSMut's improved performance vis-à-vis state-of-the-art driver gene identification methodologies by recapitulating the functional impact of known driver genes, while additionally identifying novel functionally impactful mutated genes across 29 cancers. Subsequent application of SYSMut on low-frequency gene mutations in head and neck squamous cell (HNSC) cancers, followed by molecular and pharmacogenetic validation, revealed the lipidogenic network as a novel therapeutic vulnerability in aggressive HNSC cancers. SYSMut is thus a robust scalable framework that enables the discovery of new targetable avenues in cancer.

DiaDeL: An Accurate Deep Learning-Based Model With Mutational Signatures for Predicting Metastasis Stage and Cancer Types

Mutational signatures help identify cancer-associated genes that are being involved in tumorigenesis pathways. Hence, these pathways guide precision medicine approaches to find appropriate drugs and treatments. The pattern of mutations varies in different cancer types. Some mutations dysregulate protein function so that their accumulation is responsible for cancer development and might be associated with different cancer types. Therefore, mutations as a feature set can be used as an informative candidate to distinguish various cancer types. There are several options for demonstrating mutations. One might employ binary values to demonstrate mutation regions. Another potential method for extracting features is utilizing mutation interpreters. In this study, we investigate the trinucleotide mutational pattern of each cancer type. Moreover, we extract salient NMF-based mutational signatures across various cancer types. Then, we identify cancer-associated genes of a target cancer based on its salient signatures. We evaluate the cancer-associated genes using survival and gene expression analysis in different stages of cancer. Furthermore, we introduce DiaDeL, which is a deep learning-based binary classifier. The DiaDeL model uses mutational signatures as input features and distinct a cancer type from the others. Our proposed model outperforms six state-of-the-art methods with 0.824 and 0.88 for accuracy and AUC, respectively. The source code is available at https://github.com/sabdollahi/DiaDeL.

Inferring Potential Cancer Driving Synonymous Variants

Synonymous single nucleotide variants (sSNVs) are often considered functionally silent, but a few cases of cancer-causing sSNVs have been reported. From available databases, we collected four categories of sSNVs: germline, somatic in normal tissues, somatic in cancerous tissues, and putative cancer drivers. We found that screening sSNVs for recurrence among patients, conservation of the affected genomic position, and synVep prediction (synVep is a machine learning-based sSNV effect predictor) recovers cancer driver variants (termed proposed drivers) and previously unknown putative cancer genes. Of the 2.9 million somatic sSNVs found in the COSMIC database, we identified 2111 proposed cancer driver sSNVs. Of these, 326 sSNVs could be further tagged for possible RNA splicing effects, RNA structural changes, and affected RBP motifs. This list of proposed cancer driver sSNVs provides computational guidance in prioritizing the experimental evaluation of synonymous mutations found in cancers. Furthermore, our list of novel potential cancer genes, galvanized by synonymous mutations, may highlight yet unexplored cancer mechanisms.

The emerging roles of NGS in clinical oncology and personalized medicine

Next-generation sequencing (NGS) has been increasingly popular in genomics studies over the last decade, as new sequencing technology has been created and improved. Recently, NGS started to be used in clinical oncology to improve cancer therapy through diverse modalities ranging from finding novel and rare cancer mutations, discovering cancer mutation carriers to reaching specific therapeutic approaches known as personalized medicine (PM). PM has the potential to minimize medical expenses by shifting the current traditional medical approach of treating cancer and other diseases to an individualized preventive and predictive approach. Currently, NGS can speed up in the early diagnosis of diseases and discover pharmacogenetic markers that help in personalizing therapies. Despite the tremendous growth in our understanding of genetics, NGS holds the added advantage of providing more comprehensive picture of cancer landscape and uncovering cancer development pathways. In this review, we provided a complete overview of potential NGS applications in scientific and clinical oncology, with a particular emphasis on pharmacogenomics in the direction of precision medicine treatment options.

The Role of Dielectrophoresis for Cancer Diagnosis and Prognosis

Liquid biopsy is emerging as a potential diagnostic tool for prostate cancer (PC) prognosis and diagnosis. Unfortunately, most circulating tumor cells (CTC) technologies, such as AdnaTest or Cellsearch®, critically rely on the epithelial cell adhesion molecule (EpCAM) marker, limiting the possibility of detecting cancer stem-like cells (CSCs) and mesenchymal-like cells (EMT-CTCs) that are present during PC progression. In this context, dielectrophoresis (DEP) is an epCAM independent, label-free enrichment system that separates rare cells simply on the basis of their specific electrical properties. As compared to other technologies, DEP may represent a superior technique in terms of running costs, cell yield and specificity. However, because of its higher complexity, it still requires further technical as well as clinical development. DEP can be improved by the use of microfluid, nanostructured materials and fluoro-imaging to increase its potential applications. In the context of cancer, the usefulness of DEP lies in its capacity to detect CTCs in the bloodstream in their epithelial, mesenchymal, or epithelial-mesenchymal phenotype forms, which should be taken into account when choosing CTC enrichment and analysis methods for PC prognosis and diagnosis.

A computational approach for the discovery of significant cancer genes by weighted mutation and asymmetric spreading strength in networks

Identifying significantly mutated genes in cancer is essential for understanding the mechanisms of tumor initiation and progression. This task is a key challenge since large-scale genomic studies have reported an endless number of genes mutated at a shallow frequency. Towards uncovering infrequently mutated genes, gene interaction networks combined with mutation data have been explored. This work proposes Discovering Significant Cancer Genes (DiSCaGe), a computational method for discovering significant genes for cancer. DiSCaGe computes a mutation score for the genes based on the type of mutations they have. The influence received for their neighbors in the network is also considered and obtained through an asymmetric spreading strength applied to a consensus gene network. DiSCaGe produces a ranking of prioritized possible cancer genes. An experimental evaluation with six types of cancer revealed the potential of DiSCaGe for discovering known and possible novel significant cancer genes.

Identification of Common Driver Gene Modules and Associations between Cancers through Integrated Network Analysis

High-throughput biological data has created an unprecedented opportunity for illuminating the mechanisms of tumor emergence and evolution. An important and challenging problem in deciphering cancers is to investigate the commonalities of driver genes and pathways and the associations between cancers. Aiming at this problem, a tool ComCovEx is developed to identify common cancer driver gene modules between two cancers by searching for the candidates in local signaling networks using an exclusivity-coverage iteration strategy and outputting those with significant coverage and exclusivity for both cancers. The associations of the cancer pairs are further evaluated by Fisher's exact test. Being applied to 11 TCGA cancer datasets, ComCovEx identifies 13 significantly associated cancer pairs with plenty of biologically significant common gene modules. The novel results of cancer relationship and common gene modules reveal the relevant pathological basis of different cancer types and provide new clues to diagnosis and drug treatment in associated cancers.

The Erlang distribution approximates the age distribution of incidence of childhood and young adulthood cancers

Background

It is widely believed that cancers develop upon acquiring a particular number of (epi) mutations in driver genes, but the law governing the kinetics of this process is not known. We have previously shown that the age distribution of incidence for the 20 most prevalent cancers of old age is best approximated by the Erlang probability distribution. The Erlang distribution describes the probability of several successive random events occurring by the given time according to the Poisson process, which allows an estimate for the number of critical driver events.

Methods

Here we employ a computational grid search method to find global parameter optima for five probability distributions on the CDC WONDER dataset of the age distribution of childhood and young adulthood cancer incidence.

Results

We show that the Erlang distribution is the only classical probability distribution we found that can adequately model the age distribution of incidence for all studied childhood and young adulthood cancers, in addition to cancers of old age.

Conclusions

This suggests that the Poisson process governs driver accumulation at any age and that the Erlang distribution can be used to determine the number of driver events for any cancer type. The Poisson process implies the fundamentally random timing of driver events and their constant average rate. As waiting times for the occurrence of the required number of driver events are counted in decades, and most cells do not live this long, it suggests that driver mutations accumulate silently in the longest-living dividing cells in the body—the stem cells.

Protein domain-based approaches for the identification and prioritization of therapeutically actionable cancer variants

The tremendous number of cancer variants that can be detected by NGS analyses has required the development of computational approaches to prioritize mutations on the basis of their biological and clinical significance. Standard strategies take a gene-centric approach to the problem, allowing exclusively the identification of highly frequent variants. On the contrary, protein domain (PD)-based approaches allow to identify functionally relevant low frequency variants by searching for mutations that recur on analogous residues across homologous proteins (i.e. containing the same PD). Such approaches enable to transfer information about the effects and druggability from one known mutation to unknown ones. Here we describe how PD-based strategies work, and discuss how they could be exploited for mutation prioritization. The principle that mutations clustered on specific residues of PDs have the same functional consequences and are therapeutically actionable in a similar manner could help the choice of patient-specific targeted drugs, eventually improving the management of cancer patients.

Identifying Driver Genes for Individual Patients through Inductive Matrix Completion

MOTIVATION

The driver genes play a key role in the evolutionary process of cancer. Effectively identifying these driver genes is crucial to cancer diagnosis and treatment. However, due to the high heterogeneity of cancers, it remains challenging to identify the driver genes for individual patients. Although some computational methods have been proposed to tackle this problem, they seldom consider the fact that the genes functionally similar to the well-established driver genes may likely play similar roles in cancer process, which potentially promotes the driver gene identification. Thus, here we developed a novel approach of IMCDriver to promote the driver gene identification both for cohorts and individual patients.

RESULTS

IMCDriver first considers the well-established driver genes as prior information, and adopts the using multi-omics data (e.g., somatic mutation, gene expression and protein-protein interaction) to compute the similarity between patients/genes. Then, IMCDriver prioritizes the personalized mutated genes according to their functional similarity to the well-established driver genes via Inductive Matrix Completion. Finally, IMCDriver identifies the highly rank-ordered genes as the personalized driver genes. The results on five cancer datasets from TCGA show that our IMCDriver outperforms other existing state-of-the-art methods both in the cohort and patient-specific driver gene identification. IMCDriver also reveals some novel driver genes that potentially drive cancer development. In addition, even for the driver genes rarely mutated among a population, IMCDriver can still identify them and prioritize them with high priorities.

AVAILABILITY

Code available at https://github.com/NWPU-903PR/IMCDriver.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A Systems Approach to Brain Tumor Treatment

Brain tumors are among the most lethal tumors. Glioblastoma, the most frequent primary brain tumor in adults, has a median survival time of approximately 15 months after diagnosis or a five-year survival rate of 10%; the recurrence rate is nearly 90%. Unfortunately, this prognosis has not improved for several decades. The lack of progress in the treatment of brain tumors has been attributed to their high rate of primary therapy resistance. Challenges such as pronounced inter-patient variability, intratumoral heterogeneity, and drug delivery across the blood-brain barrier hinder progress. A comprehensive, multiscale understanding of the disease, from the molecular to the whole tumor level, is needed to address the intratumor heterogeneity resulting from the coexistence of a diversity of neoplastic and non-neoplastic cell types in the tumor tissue. By contrast, inter-patient variability must be addressed by subtyping brain tumors to stratify patients and identify the best-matched drug(s) and therapies for a particular patient or cohort of patients. Accomplishing these diverse tasks will require a new framework, one involving a systems perspective in assessing the immense complexity of brain tumors. This would in turn entail a shift in how clinical medicine interfaces with the rapidly advancing high-throughput (HTP) technologies that have enabled the omics-scale profiling of molecular features of brain tumors from the single-cell to the tissue level. However, several gaps must be closed before such a framework can fulfill the promise of precision and personalized medicine for brain tumors. Ultimately, the goal is to integrate seamlessly multiscale systems analyses of patient tumors and clinical medicine. Accomplishing this goal would facilitate the rational design of therapeutic strategies matched to the characteristics of patients and their tumors. Here, we discuss some of the technologies, methodologies, and computational tools that will facilitate the realization of this vision to practice.

Sequence Neighborhoods Enable Reliable Prediction of Pathogenic Mutations in Cancer Genomes

Simple Summary

Cancer is caused by the accumulation of somatic mutations, some of which are responsible for the disease’s progression (drivers) while others are functionally neutral (passengers). Although several methods have been developed to distinguish between the two classes of mutations, very few have concentrated on using the neighborhood nucleotide sequences as potential discrimination features. In this study, we show that driver mutations’ neighborhood is significantly different from that of passengers. We further develop a novel machine learning tool, NBDriver, which is highly efficient at identifying pathogenic variants from multiple independent test datasets. Efficient and accurate identification of novel pathogenic variants from sequenced cancer genomes would help facilitate more effective therapies tailored to patients’ mutational profiles.

Abstract

Identifying cancer-causing mutations from sequenced cancer genomes hold much promise for targeted therapy and precision medicine. “Driver” mutations are primarily responsible for cancer progression, while “passengers” are functionally neutral. Although several computational approaches have been developed for distinguishing between driver and passenger mutations, very few have concentrated on using the raw nucleotide sequences surrounding a particular mutation as potential features for building predictive models. Using experimentally validated cancer mutation data in this study, we explored various string-based feature representation techniques to incorporate information on the neighborhood bases immediately 5′ and 3′ from each mutated position. Density estimation methods showed significant distributional differences between the neighborhood bases surrounding driver and passenger mutations. Binary classification models derived using repeated cross-validation experiments provided comparable performances across all window sizes. Integrating sequence features derived from raw nucleotide sequences with other genomic, structural, and evolutionary features resulted in the development of a pan-cancer mutation effect prediction tool, NBDriver, which was highly efficient in identifying pathogenic variants from five independent validation datasets. An ensemble predictor obtained by combining the predictions from NBDriver with three other commonly used driver prediction tools (FATHMM (cancer), CONDEL, and MutationTaster) significantly outperformed existing pan-cancer models in prioritizing a literature-curated list of driver and passenger mutations. Using the list of true positive mutation predictions derived from NBDriver, we identified a list of 138 known driver genes with functional evidence from various sources. Overall, our study underscores the efficacy of using raw nucleotide sequences as features to distinguish between driver and passenger mutations from sequenced cancer genomes.

Current cancer driver variant predictors learn to recognize driver genes instead of functional variants

BACKGROUND

Identifying variants that drive tumor progression (driver variants) and distinguishing these from variants that are a byproduct of the uncontrolled cell growth in cancer (passenger variants) is a crucial step for understanding tumorigenesis and precision oncology. Various bioinformatics methods have attempted to solve this complex task.

RESULTS

In this study, we investigate the assumptions on which these methods are based, showing that the different definitions of driver and passenger variants influence the difficulty of the prediction task. More importantly, we prove that the data sets have a construction bias which prevents the machine learning (ML) methods to actually learn variant-level functional effects, despite their excellent performance. This effect results from the fact that in these data sets, the driver variants map to a few driver genes, while the passenger variants spread across thousands of genes, and thus just learning to recognize driver genes provides almost perfect predictions.

CONCLUSIONS

To mitigate this issue, we propose a novel data set that minimizes this bias by ensuring that all genes covered by the data contain both driver and passenger variants. As a result, we show that the tested predictors experience a significant drop in performance, which should not be considered as poorer modeling, but rather as correcting unwarranted optimism. Finally, we propose a weighting procedure to completely eliminate the gene effects on such predictions, thus precisely evaluating the ability of predictors to model the functional effects of single variants, and we show that indeed this task is still open.

Biased Influences of Low Tumor Purity on Mutation Detection in Cancer

The non-cancerous components in tumor tissues, e.g., infiltrating stromal cells and immune cells, dilute tumor purity and might confound genomic mutation profile analyses and the identification of pathological biomarkers. It is necessary to systematically evaluate the influence of tumor purity. Here, using public gastric cancer samples from The Cancer Genome Atlas (TCGA), we firstly showed that numbers of mutation, separately called by four algorithms, were significant positively correlated with tumor purities (all p < 0.05, Spearman rank correlation). Similar results were also observed in other nine cancers from TCGA. Notably, the result was further confirmed by six in-house samples from two gastric cancer patients and five in-house samples from two colorectal cancer patients with different tumor purities. Furthermore, the metastasis mechanism of gastric cancer may be incorrectly characterized as numbers of mutation and tumor purities of 248 lymph node metastatic (N + M0) samples were both significantly lower than those of 121 non-metastatic (N0M0) samples (p < 0.05, Wilcoxon rank-sum test). Similar phenomena were also observed that tumor purities could confound the analysis of histological subtypes of cancer and the identification of microsatellite instability status (MSI) in both gastric and colon cancer. Finally, we suggested that the higher tumor purity, such as above 70%, rather than 60%, could be better to meet the requirement of mutation calling. In conclusion, the influence of tumor purity on the genomic mutation profile and pathological analyses should be fully considered in the further study.

An efficient and effective method to identify significantly perturbed subnetworks in cancer

The identification of key functional biological networks from high-dimensional genomics data is pivotal for cancer research. Here, we introduce FDRnet, a method for the detection of molecular subnetworks in cancer, which addresses several challenges in pathway analysis. FDRnet detects key subnetworks by solving a mixed-integer linear programming problem, using a given upper bound of false discovery rate (FDR) as a budget constraint, and minimizing a conductance score to find dense subgraphs around seed genes. A large-scale benchmark study was performed on both simulation and cancer genomics data. FDRnet outperformed other methods in the ability to detect functionally homogeneous subnetworks in a scale-free biological network, to control FDRs of the genes in detected subnetworks, to improve computational efficiency and to integrate multi-omics data. By overcoming the limitations of existing approaches, FDRnet can facilitate the detection of key functional pathways in cancer and other genetic diseases.

Correlation of genomic alterations and PD-L1 expression in thymoma

Thymic epithelial tumors (TETs) include several anterior mediastinal malignant tumours: thymomas, thymic carcinomas and thymic neuroendocrine cancers. There is significant variety in the biologic features and clinical course of TETs and many attempts have been made to identify target genes for successful therapy of TETs. Next generation sequencing (NGS) represents a huge advancement in diagnostics and these new molecular technologies revealed that thymic neoplasms have the lowest tumor mutation burden among all adult malignant tumours with a different pattern of molecular aberrations in thymomas and thymic carcinomas. As for the PD-L1 expression in tumor cells in thymoma and thymic carcinoma, it varies a lot in published studies, with findings of PD-L1 expression from 23% to 92% in thymoma and 36% to 100% in thymic carcinoma. When correlated PD-L1 expression with disease stage some controversial results were obtained, with no association with tumor stage in most studies. This is, at least in part, explained by the fact that several diverse PD-L1 immunohistochemical tests were used in each trial, with four different antibodies (SP142, SP263, 22C3, and 28-8), different definition of PD-L1 positivity and cutoff values throughout the studies as well. There is a huge interest in using genomic features to produce predictive genomic-based immunotherapy biomarkers, particularly since recent data suggest that certain tumor-specific genomic alterations, either individually or in combination, appear to influence immune checkpoint activity and better responses as the outcome, so as such in some cancer types they may complement existing biomarkers to improve the selection criteria for immunotherapy.

Discovering novel driver mutations from pan-cancer analysis of mutational and gene expression profiles

As the genomic profile across cancers varies from person to person, patient prognosis and treatment may differ based on the mutational signature of each tumour. Thus, it is critical to understand genomic drivers of cancer and identify potential mutational commonalities across tumors originating at diverse anatomical sites. Large-scale cancer genomics initiatives, such as TCGA, ICGC and GENIE have enabled the analysis of thousands of tumour genomes. Our goal was to identify new cancer-causing mutations that may be common across tumour sites using mutational and gene expression profiles. Genomic and transcriptomic data from breast, ovarian, and prostate cancers were aggregated and analysed using differential gene expression methods to identify the effect of specific mutations on the expression of multiple genes. Mutated genes associated with the most differentially expressed genes were considered to be novel candidates for driver mutations, and were validated through literature mining, pathway analysis and clinical data investigation. Our driver selection method successfully identified 116 probable novel cancer-causing genes, with 4 discovered in patients having no alterations in any known driver genes: MXRA5, OBSCN, RYR1, and TG. The candidate genes previously not officially classified as cancer-causing showed enrichment in cancer pathways and in cancer diseases. They also matched expectations pertaining to properties of cancer genes, for instance, showing larger gene and protein lengths, and having mutation patterns suggesting oncogenic or tumor suppressor properties. Our approach allows for the identification of novel putative driver genes that are common across cancer sites using an unbiased approach without any a priori knowledge on pathways or gene interactions and is therefore an agnostic approach to the identification of putative common driver genes acting at multiple cancer sites.

Comprehensive assessment of TP53 loss of function using multiple combinatorial mutagenesis libraries

The diagnosis of somatic and germline TP53 mutations in human tumors or in individuals prone to various types of cancer has now reached the clinic. To increase the accuracy of the prediction of TP53 variant pathogenicity, we gathered functional data from three independent large-scale saturation mutagenesis screening studies with experimental data for more than 10,000 TP53 variants performed in different settings (yeast or mammalian) and with different readouts (transcription, growth arrest or apoptosis). Correlation analysis and multidimensional scaling showed excellent agreement between all these variables. Furthermore, we found that some missense mutations localized in TP53 exons led to impaired TP53 splicing as shown by an analysis of the TP53 expression data from the cancer genome atlas. With the increasing availability of genomic, transcriptomic and proteomic data, it is essential to employ both protein and RNA prediction to accurately define variant pathogenicity.

Untangling a complex web: Computational analyses of tumor molecular profiles to decode driver mechanisms

Genome-scale studies focusing on molecular profiling of cancers across tissue types have revealed a plethora of aberrations across the genomic, transcriptomic, and epigenomic scales. The significant molecular heterogeneity across individual tumors even within the same tissue context complicates decoding the key etiologic mechanisms of this disease. Furthermore, it is increasingly likely that biologic mechanisms underlying the pathobiology of cancer involve multiple molecular entities interacting across functional scales. This has motivated the development of computational approaches that integrate molecular measurements with prior biological knowledge in increasingly intricate ways to enable the discovery of driver genomic aberrations across cancers. Here, we review diverse methodological approaches that have powered significant advances in our understanding of the genomic underpinnings of cancer at the cohort and at the individual tumor scales. We outline the key advances and challenges in the computational discovery of cancer mechanisms while motivating the development of systems biology approaches to comprehensively decode the biologic drivers of this complex disease.

Deep Neural Network Algorithm Feedback Model with Behavioral Intelligence and Forecast Accuracy

When attempting to apply a large-scale database that holds the behavioral intelligence training data of deep neural networks, the classification accuracy of the artificial intelligence algorithm needs to reflect the behavioral characteristics of the individual. When a change in behavior is recognized, that is, a feedback model based on a data connection model is applied, an analysis of time series data is performed by extracting feature vectors and interpolating data in a deep neural network to overcome the limitations of the existing statistical analysis. Using the results of the first feedback model as inputs to the deep neural network and, furthermore, as the input values of the second feedback model, and interpolating the behavioral intelligence data, that is, context awareness and lifelog data, including physical activities, involves applying the most appropriate conditions. The results of this study show that this method effectively improves the accuracy of the artificial intelligence results. In this paper, through an experiment, after extracting the feature vector of a deep neural network and restoring the missing value, the classification accuracy was verified to improve by about 20% on average. At the same time, by adding behavioral intelligence data to the time series data, a new data connection model, the Deep Neural Network Feedback Model, was proposed, and it was verified that the classification accuracy can be improved by about 8 to 9% on average. Based on the hypothesis, the F (X′) = X model was applied to thoroughly classify the training data set and test data set to present a symmetrical balance between the data connection model and the context-aware data. In addition, behavioral activity data were extrapolated in terms of context-aware and forecasting perspectives to prove the results of the experiment.

Reconstruction of clone- and haplotype-specific cancer genome karyotypes from bulk tumor samples

Many cancer genomes are extensively rearranged with aberrant chromosomal karyotypes. Deriving these karyotypes from high-throughput DNA sequencing of bulk tumor samples is complicated because most tumors are a heterogeneous mixture of normal cells and subpopulations of cancer cells, or clones, that harbor distinct somatic mutations. We introduce a new algorithm, Reconstructing Cancer Karyotypes (RCK), to reconstruct haplotype-specific karyotypes of one or more rearranged cancer genomes from DNA sequencing data from a bulk tumor sample. RCK leverages evolutionary constraints on the somatic mutational process in cancer to reduce ambiguity in the deconvolution of admixed sequencing data into multiple haplotype-specific cancer karyotypes. RCK models mixtures containing an arbitrary number of derived genomes and allows the incorporation of information both from short-read and long-read DNA sequencing technologies. We compare RCK to existing approaches on 17 primary and metastatic prostate cancer samples. We find that RCK infers cancer karyotypes that better explain the DNA sequencing data and conform to a reasonable evolutionary model. RCK's reconstructions of clone- and haplotype-specific karyotypes will aid further studies of the role of intra-tumor heterogeneity in cancer development and response to treatment. RCK is freely available as open source software.

Approaches for the identification of driver mutations in cancer: A tutorial from a computational perspective

Cancer is a complex disease caused by the accumulation of genetic alterations during the individual's life. Such alterations are called genetic mutations and can be divided into two groups: (1) Passenger mutations, which are not responsible for cancer and (2) Driver mutations, which are significant for cancer and responsible for its initiation and progression. Cancer cells undergo a large number of mutations, of which most are passengers, and few are drivers. The identification of driver mutations is a key point and one of the biggest challenges in Cancer Genomics. Many computational methods for such a purpose have been developed in Cancer Bioinformatics. Such computational methods are complex and are usually described in a high level of abstraction. This tutorial details some classical computational methods, from a computational perspective, with the transcription in an algorithmic format towards an easy access by researchers.

Clonal evolution driven by superdriver mutations

Background

Tumors are widely recognized to progress through clonal evolution by sequentially acquiring selectively advantageous genetic alterations that significantly contribute to tumorigenesis and thus are termned drivers. Some cancer drivers, such as TP53 point mutation or EGFR copy number gain, provide exceptional fitness gains, which, in time, can be sufficient to trigger the onset of cancer with little or no contribution from additional genetic alterations. These key alterations are called superdrivers.

Results

In this study, we employ a Wright-Fisher model to study the interplay between drivers and superdrivers in tumor progression. We demonstrate that the resulting evolutionary dynamics follow global clonal expansions of superdrivers with periodic clonal expansions of drivers. We find that the waiting time to the accumulation of a set of superdrivers and drivers in the tumor cell population can be approximated by the sum of the individual waiting times.

Conclusions

Our results suggest that superdriver dynamics dominate over driver dynamics in tumorigenesis. Furthermore, our model allows studying the interplay between superdriver and driver mutations both empirically and theoretically.

Natural Selection on Exonic SNPs Shapes Allelic Expression Imbalance (AEI) Adaptability in Lung Cancer Progression

Tumors are driven by a sequence of genetic and epigenetic alterations. Previous studies have mostly focused on the roles of somatic mutations in tumorigenesis, but how germline variants act is largely unknown. In this study, we hypothesized that allelic expression imbalance (AEI) participated in the process of germline variants on tumorigenesis. We screened single-nucleotide polymorphisms (SNPs) as representative germline variants. By using 127 patients’ RNA sequencing data from paired lung cancer and adjacent normal tissues from public databases, we analyzed the effects of the functional consequence of SNPs, function and conservativeness on genes with AEI. We found that natural selection can affect AEI. Functional adaptability of genes with a high frequency of AEI and a correlation of the incidence of AEI with conservativeness were observed in both adjacent tissues and tumor tissues. Moreover, we observed a higher incidence of AEI in genes with non-synonymous SNPs than in those with synonymous SNPs. However, we also found that AEI was affected by allele expression noise, especially in tumor tissues, which led to an increased proportion of AEI, weakened the effect of natural selection and eliminated the influence of the functional consequence of SNPs on AEI. We unveiled a previously unknown adaptive regulatory mechanism in which the effect of natural selection on SNPs can be reflected in allelic expression, which provides insight into a better understanding of cancer evolution.

Genomic profiling of thymoma using a targeted high-throughput approach

Introduction

Thymomas and thymic carcinoma (TC) are the most common neoplasms localised in the thymus. These diseases are poorly understood, but progress made in next-generation sequencing (NGS) technology has provided novel data on their molecular pathology.

Material and methods

Genomic DNA was isolated from formalin-fixed paraffin-embedded tumour tissue. We investigated somatic variants in 35 thymoma patients using amplicon-based TruSeq Amplicon Cancer Panel (TSACP) that covers 48 cancer related genes. We also analysed three samples from healthy individuals by TSACP platform and 32 healthy controls using exome sequencing.

Results

The total number of detected variants was 4447, out of which 2906 were in the coding region (median per patient 83, range: 2-300) and 1541 were in the non-coding area (median per patient 44, range: 0-172). We identified four genes, APC, ATM, ERBB4, and SMAD4, having more than 100 protein-changing variants. Additionally, more than 70% of the analysed cases harboured protein-changing variants in SMAD4, APC, ATM, PTEN, KDR, and TP53. Moreover, this study revealed 168 recurrent variants, out of which 15 were shown to be pathogenic. Comparison to controls revealed that the variants we reported in this study were somatic thymoma-specific variants. Additionally, we found that the presence of variants in SMAD4 gene predicted shorter overall survival in thymoma patients.

Conclusions

The most frequently mutated genes in thymoma samples analysed in this study belong to the EGFR, ATM, and TP53 signalling pathways, regulating cell cycle check points, gene expression, and apoptosis. The results of our study complement the knowledge of thymoma molecular pathogenesis.

Prioritizing Cancer Genes Based on an Improved Random Walk Method

Identifying driver genes that contribute to cancer progression from numerous passenger genes, although a central goal, is a major challenge. The protein-protein interaction network provides convenient and reasonable assistance for driver gene discovery. Random walk-based methods have been widely used to prioritize nodes in social or biological networks. However, most studies select the next arriving node uniformly from the random walker's neighbors. Few consider transiting preference according to the degree of random walker's neighbors. In this study, based on the random walk method, we propose a novel approach named Driver_IRW (Driver genes discovery with Improved Random Walk method), to prioritize cancer genes in cancer-related network. The key idea of Driver_IRW is to assign different transition probabilities for different edges of a constructed cancer-related network in accordance with the degree of the nodes' neighbors. Furthermore, the global centrality (here is betweenness centrality) and Katz feedback centrality are incorporated into the framework to evaluate the probability to walk to the seed nodes. Experimental results on four cancer types indicate that Driver_IRW performs more efficiently than some previously published methods for uncovering known cancer-related genes. In conclusion, our method can aid in prioritizing cancer-related genes and complement traditional frequency and network-based methods.

Non-coding somatic mutations converge on the PAX8 pathway in ovarian cancer

The functional consequences of somatic non-coding mutations in ovarian cancer (OC) are unknown. To identify regulatory elements (RE) and genes perturbed by acquired non-coding variants, here we establish epigenomic and transcriptomic landscapes of primary OCs using H3K27ac ChIP-seq and RNA-seq, and then integrate these with whole genome sequencing data from 232 OCs. We identify 25 frequently mutated regulatory elements, including an enhancer at 6p22.1 which associates with differential expression of ZSCAN16 (P = 6.6 × 10-4) and ZSCAN12 (P = 0.02). CRISPR/Cas9 knockout of this enhancer induces downregulation of both genes. Globally, there is an enrichment of single nucleotide variants in active binding sites for TEAD4 (P = 6 × 10-11) and its binding partner PAX8 (P = 2×10-10), a known lineage-specific transcription factor in OC. In addition, the collection of cis REs associated with PAX8 comprise the most frequently mutated set of enhancers in OC (P = 0.003). These data indicate that non-coding somatic mutations disrupt the PAX8 transcriptional network during OC development.

Distance measures for tumor evolutionary trees

MOTIVATION

There has been recent increased interest in using algorithmic methods to infer the evolutionary tree underlying the developmental history of a tumor. Quantitative measures that compare such trees are vital to a number of different applications including benchmarking tree inference methods and evaluating common inheritance patterns across patients. However, few appropriate distance measures exist, and those that do have low resolution for differentiating trees or do not fully account for the complex relationship between tree topology and the inheritance of the mutations labeling that topology.

RESULTS

Here, we present two novel distance measures, Common Ancestor Set distance (CASet) and Distinctly Inherited Set Comparison distance (DISC), that are specifically designed to account for the subclonal mutation inheritance patterns characteristic of tumor evolutionary trees. We apply CASet and DISC to multiple simulated datasets and two breast cancer datasets and show that our distance measures allow for more nuanced and accurate delineation between tumor evolutionary trees than existing distance measures.

AVAILABILITY AND IMPLEMENTATION

Implementations of CASet and DISC are freely available at: https://bitbucket.org/oesperlab/stereodist.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A systematic review of the applications of artificial intelligence and machine learning in autoimmune diseases

Autoimmune diseases are chronic, multifactorial conditions. Through machine learning (ML), a branch of the wider field of artificial intelligence, it is possible to extract patterns within patient data, and exploit these patterns to predict patient outcomes for improved clinical management. Here, we surveyed the use of ML methods to address clinical problems in autoimmune disease. A systematic review was conducted using MEDLINE, embase and computers and applied sciences complete databases. Relevant papers included "machine learning" or "artificial intelligence" and the autoimmune diseases search term(s) in their title, abstract or key words. Exclusion criteria: studies not written in English, no real human patient data included, publication prior to 2001, studies that were not peer reviewed, non-autoimmune disease comorbidity research and review papers. 169 (of 702) studies met the criteria for inclusion. Support vector machines and random forests were the most popular ML methods used. ML models using data on multiple sclerosis, rheumatoid arthritis and inflammatory bowel disease were most common. A small proportion of studies (7.7% or 13/169) combined different data types in the modelling process. Cross-validation, combined with a separate testing set for more robust model evaluation occurred in 8.3% of papers (14/169). The field may benefit from adopting a best practice of validation, cross-validation and independent testing of ML models. Many models achieved good predictive results in simple scenarios (e.g. classification of cases and controls). Progression to more complex predictive models may be achievable in future through integration of multiple data types.

HetRCNA: A Novel Method to Identify Recurrent Copy Number Alternations from Heterogeneous Tumor Samples Based on Matrix Decomposition Framework

A common strategy to discovering cancer associated copy number aberrations (CNAs) from a cohort of cancer samples is to detect recurrent CNAs (RCNAs). Although the previous methods can successfully identify communal RCNAs shared by nearly all tumor samples, detecting subgroup-specific RCNAs and their related subgroup samples from cancer samples with heterogeneity is still invalid for these existing approaches. In this paper, we introduce a novel integrated method called HetRCNA, which can identify statistically significant subgroup-specific RCNAs and their related subgroup samples. Based on matrix decomposition framework with weight constraint, HetRCNA can successfully measure the subgroup samples by coefficients of left vectors with weight constraint and subgroup-specific RCNAs by coefficients of the right vectors and significance test. When we evaluate HetRCNA on simulated dataset, the results show that HetRCNA gives the best performances among the competing methods and is robust to the noise factors of the simulated data. When HetRCNA is applied on a real breast cancer dataset, our approach successfully identifies a bunch of RCNA regions and the result is highly correlated with the results of the other two investigated approaches. Notably, the genomic regions identified by HetRCNA harbor many breast cancer related genes reported by previous researches.

Comprehensive Detection of Single Amino Acid Variants and Evaluation of Their Deleterious Potential in a PANC-1 Cell Line

Identifying single amino acid variants (SAAVs) in cancer is critical for precision oncology. Several advanced algorithms are now available to identify SAAVs, but attempts to combine different algorithms and optimize them on large data sets to achieve a more comprehensive coverage of SAAVs have not been implemented. Herein, we report an expanded detection of SAAVs in the PANC-1 cell line using three different strategies, which results in the identification of 540 SAAVs in the mass spectrometry data. Among the set of 540 SAAVs, 79 are evaluated as deleterious SAAVs based on analysis using the novel AssVar software in which one of the driver mutations found in each protein of KRAS, TP53, and SLC37A4 is further validated using independent selected reaction monitoring (SRM) analysis. Our study represents the most comprehensive discovery of SAAVs to date and the first large-scale detection of deleterious SAAVs in the PANC-1 cell line. This work may serve as the basis for future research in pancreatic cancer and personal immunotherapy and treatment.

Multimodal genomic features predict outcome of immune checkpoint blockade in non-small-cell lung cancer

Despite progress in immunotherapy, identifying patients that respond has remained a challenge. Through analysis of whole-exome and targeted sequence data from 5,449 tumors, we found a significant correlation between tumor mutation burden (TMB) and tumor purity, suggesting that low tumor purity tumors are likely to have inaccurate TMB estimates. We developed a new method to estimate a corrected TMB (cTMB) that was adjusted for tumor purity and more accurately predicted outcome to immune checkpoint blockade (ICB). To identify improved predictive markers together with cTMB, we performed whole-exome sequencing for 104 lung tumors treated with ICB. Through comprehensive analyses of sequence and structural alterations, we discovered a significant enrichment in activating mutations in receptor tyrosine kinase (RTK) genes in nonresponding tumors in three immunotherapy treated cohorts. An integrated multivariable model incorporating cTMB, RTK mutations, smoking-related mutational signature and human leukocyte antigen status provided an improved predictor of response to immunotherapy that was independently validated.