Pathway and network analysis of more than 2500 whole cancer genomes

Death-associated protein kinase 3 modulates migration and invasion of triple-negative breast cancer cells

Sixteen patient-derived xenografts (PDXs) were analyzed using a mass spectrometry (MS)-based kinase inhibitor pull-down assay (KIPA), leading to the observation that death-associated protein kinase 3 (DAPK3) is significantly and specifically overexpressed in the triple-negative breast cancer (TNBC) models. Validation studies confirmed enrichment of DAPK3 protein, in both TNBC cell lines and tumors, independent of mRNA levels. Genomic knockout of DAPK3 in TNBC cell lines inhibited in vitro migration and invasion, along with down-regulation of an epithelial-mesenchymal transition (EMT) signature, which was confirmed in vivo. The kinase and leucine-zipper domains within DAPK3 were shown by a mutational analysis to be essential for functionality. Notably, DAPK3 was found to inhibit the levels of desmoplakin (DSP), a crucial component of the desmosome complex, thereby explaining the observed migration and invasion effects. Further exploration with immunoprecipitation-mass spectrometry (IP-MS) identified that leucine-zipper protein 1 (LUZP1) is a preferential binding partner of DAPK3. LUZP1 engages in a leucine-zipper domain-mediated interaction that protects DAPK3 from proteasomal degradation. Thus, the DAPK3/LUZP1 heterodimer emerges as a newly discovered regulator of EMT/desmosome components that promote TNBC cell migration.

Mapping the Multiscale Proteomic Organization of Cellular and Disease Phenotypes

While the primary sequences of human proteins have been cataloged for over a decade, determining how these are organized into a dynamic collection of multiprotein assemblies, with structures and functions spanning biological scales, is an ongoing venture. Systematic and data-driven analyses of these higher-order structures are emerging, facilitating the discovery and understanding of cellular phenotypes. At present, knowledge of protein localization and function has been primarily derived from manual annotation and curation in resources such as the Gene Ontology, which are biased toward richly annotated genes in the literature. Here, we envision a future powered by data-driven mapping of protein assemblies. These maps can capture and decode cellular functions through the integration of protein expression, localization, and interaction data across length scales and timescales. In this review, we focus on progress toward constructing integrated cell maps that accelerate the life sciences and translational research.

ALLSTAR: Inference of ReliAble CausaL RuLes between Somatic MuTAtions and CanceR Phenotypes

MOTIVATION

Recent advances in DNA sequencing technologies have allowed the detailed characterization of genomes in large cohorts of tumors, highlighting their extreme heterogeneity, with no two tumors sharing the same complement of somatic mutations. Such heterogeneity hinders our ability to identify somatic mutations important for the disease, including mutations that determine clinically relevant phenotypes (e.g., cancer subtypes). Several tools have been developed to identify somatic mutations related to cancer phenotypes. However, such tools identify correlations between somatic mutations and cancer phenotypes, with no guarantee of highlighting causal relations.

RESULTS

We describe ALLSTAR, a novel tool to infer reliable causal relations between somatic mutations and cancer phenotypes. ALLSTAR identifies reliable causal rules highlighting combinations of somatic mutations with the highest impact in terms of average effect on the phenotype. While we prove that the underlying computational problem is NP-hard, we develop a branch-and-bound approach that employs protein-protein interaction networks and novel bounds for pruning the search space, while properly correcting for multiple hypothesis testing. Our extensive experimental evaluation on synthetic data shows that our tool is able to identify reliable causal relations in large cancer cohorts. Moreover, the reliable causal rules identified by our tool in cancer data show that our approach identifies several somatic mutations known to be relevant for cancer phenotypes as well as novel biologically meaningful relations.

AVAILABILITY AND IMPLEMENTATION

Code, data, and scripts to reproduce the experiments available at https://github.com/VandinLab/ALLSTAR.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Network analysis of driver genes in human cancers

Cancer is a heterogeneous disease that results from genetic alteration of cell cycle and proliferation controls. Identifying mutations that drive cancer, understanding cancer type specificities, and delineating how driver mutations interact with each other to establish disease is vital for identifying therapeutic vulnerabilities. Such cancer specific patterns and gene co-occurrences can be identified by studying tumor genome sequences, and networks have proven effective in uncovering relationships between sequences. We present two network-based approaches to identify driver gene patterns among tumor samples. The first approach relies on analysis using the Directed Weighted All Nearest Neighbors (DiWANN) model, which is a variant of sequence similarity network, and the second approach uses bipartite network analysis. A data reduction framework was implemented to extract the minimal relevant information for the sequence similarity network analysis, where a transformed reference sequence is generated for constructing the driver gene network. This data reduction process combined with the efficiency of the DiWANN network model, greatly lowered the computational cost (in terms of execution time and memory usage) of generating the networks enabling us to work at a much larger scale than previously possible. The DiWANN network helped us identify cancer types in which samples were more closely connected to each other suggesting they are less heterogeneous and potentially susceptible to a common drug. The bipartite network analysis provided insight into gene associations and co-occurrences. We identified genes that were broadly mutated in multiple cancer types and mutations exclusive to only a few. Additionally, weighted one-mode gene projections of the bipartite networks revealed a pattern of occurrence of driver genes in different cancers. Our study demonstrates that network-based approaches can be an effective tool in cancer genomics. The analysis identifies co-occurring and exclusive driver genes and mutations for specific cancer types, providing a better understanding of the driver genes that lead to tumor initiation and evolution.

Identifying new cancer genes based on the integration of annotated gene sets via hypergraph neural networks

MOTIVATION

Identifying cancer genes remains a significant challenge in cancer genomics research. Annotated gene sets encode functional associations among multiple genes, and cancer genes have been shown to cluster in hallmark signaling pathways and biological processes. The knowledge of annotated gene sets is critical for discovering cancer genes but remains to be fully exploited.

RESULTS

Here, we present the DIsease-Specific Hypergraph neural network (DISHyper), a hypergraph-based computational method that integrates the knowledge from multiple types of annotated gene sets to predict cancer genes. First, our benchmark results demonstrate that DISHyper outperforms the existing state-of-the-art methods and highlight the advantages of employing hypergraphs for representing annotated gene sets. Second, we validate the accuracy of DISHyper-predicted cancer genes using functional validation results and multiple independent functional genomics data. Third, our model predicts 44 novel cancer genes, and subsequent analysis shows their significant associations with multiple types of cancers. Overall, our study provides a new perspective for discovering cancer genes and reveals previously undiscovered cancer genes.

AVAILABILITY AND IMPLEMENTATION

DISHyper is freely available for download at https://github.com/genemine/DISHyper.

The intricate dance of tumor evolution: Exploring immune escape, tumor migration, drug resistance, and treatment strategies

Recent research has unveiled fascinating insights into the intricate mechanisms governing tumor evolution. These studies have illuminated how tumors adapt and proliferate by exploiting various factors, including immune evasion, resistance to therapeutic drugs, genetic mutations, and their ability to adapt to different environments. Furthermore, investigations into tumor heterogeneity and chromosomal aberrations have revealed the profound complexity that underlies the evolution of cancer. Emerging findings have also underscored the role of viral influences in the development and progression of cancer, introducing an additional layer of complexity to the field of oncology. Tumor evolution is a dynamic and complex process influenced by various factors, including immune evasion, drug resistance, tumor heterogeneity, and viral influences. Understanding these elements is indispensable for developing more effective treatments and advancing cancer therapies. A holistic approach to studying and addressing tumor evolution is crucial in the ongoing battle against cancer. The main goal of this comprehensive review is to explore the intricate relationship between tumor evolution and critical aspects of cancer biology. By delving into this complex interplay, we aim to provide a profound understanding of how tumors evolve, adapt, and respond to treatment strategies. This review underscores the pivotal importance of comprehending tumor evolution in shaping effective approaches to cancer treatment.

All paths lead to hubs in the spectroscopic networks of water isotopologues H216O and H218O

Network theory has fundamentally transformed our comprehension of complex systems, catalyzing significant advances across various domains of science and technology. In spectroscopic networks, hubs are the quantum states involved in the largest number of transitions. Here, utilizing network paths probed via precision metrology, absolute energies have been deduced, with at least 10-digit accuracy, for almost 200 hubs in the experimental spectroscopic networks of H216O and H218O. These hubs, lying on the ground vibrational states of both species and the bending fundamental of H216O, are involved in tens of thousands of observed transitions. Relying on the same hubs and other states, benchmark-quality line lists have been assembled, which supersede and improve, by three orders of magnitude, the accuracy of the massive amount of data reported in hundreds of papers dealing with Doppler-limited spectroscopy. Due to the omnipresence of water, these ultraprecise line lists could be applied to calibrate high-resolution spectra and serve ongoing and upcoming space missions.

Epigenomic Characterization of Lymphoid Neoplasms

Lymphoid neoplasms represent a heterogeneous group of disease entities and subtypes with markedly different molecular and clinical features. Beyond genetic alterations, lymphoid tumors also show widespread epigenomic changes. These severely affect the levels and distribution of DNA methylation, histone modifications, chromatin accessibility, and three-dimensional genome interactions. DNA methylation stands out as a tracer of cell identity and memory, as B cell neoplasms show epigenetic imprints of their cellular origin and proliferative history, which can be quantified by an epigenetic mitotic clock. Chromatin-associated marks are informative to uncover altered regulatory regions and transcription factor networks contributing to the development of distinct lymphoid tumors. Tumor-intrinsic epigenetic and genetic aberrations cooperate and interact with microenvironmental cells to shape the transcriptome at different phases of lymphoma evolution, and intraclonal heterogeneity can now be characterized by single-cell profiling. Finally, epigenetics offers multiple clinical applications, including powerful diagnostic and prognostic biomarkers as well as therapeutic targets.

p120 RasGAP and ZO-2 are essential for Hippo signaling and tumor-suppressor function mediated by p190A RhoGAP

ARHGAP35, which encodes p190A RhoGAP (p190A), is a major cancer gene. p190A is a tumor suppressor that activates the Hippo pathway. p190A was originally cloned via direct binding to p120 RasGAP (RasGAP). Here, we determine that interaction of p190A with the tight-junction-associated protein ZO-2 is dependent on RasGAP. We establish that both RasGAP and ZO-2 are necessary for p190A to activate large tumor-suppressor (LATS) kinases, elicit mesenchymal-to-epithelial transition, promote contact inhibition of cell proliferation, and suppress tumorigenesis. Moreover, RasGAP and ZO-2 are required for transcriptional modulation by p190A. Finally, we demonstrate that low ARHGAP35 expression is associated with shorter survival in patients with high, but not low, transcript levels of TJP2 encoding ZO-2. Hence, we define a tumor-suppressor interactome of p190A that includes ZO-2, an established constituent of the Hippo pathway, and RasGAP, which, despite strong association with Ras signaling, is essential for p190A to activate LATS kinases.

A machine learning and directed network optimization approach to uncover TP53 regulatory patterns

TP53, the Guardian of the Genome, is the most frequently mutated gene in human cancers and the functional characterization of its regulation is fundamental. To address this we employ two strategies: machine learning to predict the mutation status of TP53from transcriptomic data, and directed regulatory networks to reconstruct the effect of mutations on the transcipt levels of TP53 targets. Using data from established databases (Cancer Cell Line Encyclopedia, The Cancer Genome Atlas), machine learning could predict the mutation status, but not resolve different mutations. On the contrary, directed network optimization allowed to infer the TP53 regulatory profile across: (1) mutations, (2) irradiation in lung cancer, and (3) hypoxia in breast cancer, and we could observe differential regulatory profiles dictated by (1) mutation type, (2) deleterious consequences of the mutation, (3) known hotspots, (4) protein changes, (5) stress condition (irradiation/hypoxia). This is an important first step toward using regulatory networks for the characterization of the functional consequences of mutations, and could be extended to other perturbations, with implications for drug design and precision medicine.

Dual network analysis of transcriptome data for discovery of new therapeutic targets in non-small cell lung cancer

The drug therapy for non-small cell lung cancer (NSCLC) have always been issues of poisonous side effect, acquired drug resistance and narrow applicable population. In this study, we built a novel network analysis method (difference- correlation- enrichment- causality- node), which was based on the difference analysis, Spearman correlation network analysis, biological function analysis and Bayesian causality network analysis to discover new therapeutic target of NSCLC in the sequencing data of BEAS-2B and 7 NSCLC cell lines. Our results showed that, as a proteasome subunit coding gene in the central of cell cycle network, PSMD2 was associated with prognosis and was an independent prognostic factor for NSCLC patients. Knockout of PSMD2 inhibited the proliferation of NSCLC cells by inducing cell cycle arrest, and exhibited marked increase of cell cycle blocking protein p21, p27 and decrease of cell cycle driven protein CDK4, CDK6, CCND1 and CCNE1. IPA and molecular docking suggested bortezomib has stronger affinity to PSMD2 compared with reported targets PSMB1 and PSMB5. In vitro and In vivo experiments demonstrated the inhibitory effect of bortezomib in NSCLC with different driven mutations or with tyrosine kinase inhibitors resistance. Taken together, bortezomib could target PSMD2, PSMB1 and PSMB5 to inhibit the proteasome degradation of cell cycle check points, to block cell proliferation of NSCLC, which was potential optional drug for NSCLC patients.

Systems Biology in Cancer Diagnosis Integrating Omics Technologies and Artificial Intelligence to Support Physician Decision Making

Cancer is the second major cause of disease-related death worldwide, and its accurate early diagnosis and therapeutic intervention are fundamental for saving the patient's life. Cancer, as a complex and heterogeneous disorder, results from the disruption and alteration of a wide variety of biological entities, including genes, proteins, mRNAs, miRNAs, and metabolites, that eventually emerge as clinical symptoms. Traditionally, diagnosis is based on clinical examination, blood tests for biomarkers, the histopathology of a biopsy, and imaging (MRI, CT, PET, and US). Additionally, omics biotechnologies help to further characterize the genome, metabolome, microbiome traits of the patient that could have an impact on the prognosis and patient's response to the therapy. The integration of all these data relies on gathering of several experts and may require considerable time, and, unfortunately, it is not without the risk of error in the interpretation and therefore in the decision. Systems biology algorithms exploit Artificial Intelligence (AI) combined with omics technologies to perform a rapid and accurate analysis and integration of patient's big data, and support the physician in making diagnosis and tailoring the most appropriate therapeutic intervention. However, AI is not free from possible diagnostic and prognostic errors in the interpretation of images or biochemical-clinical data. Here, we first describe the methods used by systems biology for combining AI with omics and then discuss the potential, challenges, limitations, and critical issues in using AI in cancer research.

Extending protein interaction networks using proteoforms and small molecules

MOTIVATION

Biological network analysis for high-throughput biomedical data interpretation relies heavily on topological characteristics. Networks are commonly composed of nodes representing genes or proteins that are connected by edges when interacting. In this study, we use the rich information available in the Reactome pathway database to build biological networks accounting for small molecules and proteoforms modeled using protein isoforms and post-translational modifications to study the topological changes induced by this refinement of the network representation.

RESULTS

We find that improving the interactome modeling increases the number of nodes and interactions, but that isoform and post-translational modification annotation is still limited compared to what can be expected biologically. We also note that small molecule information can distort the topology of the network due to the high connectedness of these molecules, which does not necessarily represent the reality of biology. However, by restricting the connections of small molecules to the context of biochemical reactions, we find that these improve the overall connectedness of the network and reduce the prevalence of isolated components and nodes. Overall, changing the representation of the network alters the prevalence of articulation points and bridges globally but also within and across pathways. Hence, some molecules can gain or lose in biological importance depending on the level of detail of the representation of the biological system, which might in turn impact network-based studies of diseases or druggability.

AVAILABILITY AND IMPLEMENTATION

Networks are constructed based on data publicly available in the Reactome Pathway knowledgebase: reactome.org.

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.

p120 RasGAP and ZO-2 are essential for Hippo signaling and tumor suppressor function mediated by p190A RhoGAP

ARHGAP35 , which encodes p190A RhoGAP (p190A), is a major cancer gene. p190A is a tumor suppressor that activates the Hippo pathway. p190A was originally cloned via direct binding to p120 RasGAP (RasGAP). Here, we determine that a novel interaction of p190A with the tight junction-associated protein ZO-2 is dependent on RasGAP. We establish that both RasGAP and ZO-2 are necessary for p190A to activate LATS kinases, elicit mesenchymal-to-epithelial transition, promote contact inhibition of cell proliferation and suppress tumorigenesis. Moreover, RasGAP and ZO-2 are required for transcriptional modulation by p190A. Finally, we demonstrate that low ARHGAP35 expression is associated with shorter survival in patients with high, but not low, transcript levels of TJP2 encoding ZO-2. Hence, we define a tumor suppressor interactome of p190A that includes ZO-2, an established constituent of the Hippo pathway, and RasGAP, which despite strong association with Ras signaling, is essential for p190A to activate LATS kinases.

Molecular and phenotypic characteristics influencing the degree of cytoreduction in high-grade serous ovarian carcinomas

BACKGROUND

High-grade serous ovarian carcinoma (HGSOC) is the deadliest ovarian cancer subtype, and survival relates to initial cytoreductive surgical treatment. The existing tools for surgical outcome prediction remain inadequate for anticipating the outcomes of the complex relationship between tumour biology, clinical phenotypes, co-morbidity and surgical skills. In this genotype-phenotype association study, we combine phenotypic markers with targeted DNA sequencing to discover novel biomarkers to guide the surgical management of primary HGSOC.

METHODS

Primary tumour tissue samples (n = 97) and matched blood from a phenotypically well-characterised treatment-naïve HGSOC patient cohort were analysed by targeted massive parallel DNA sequencing (next generation sequencing [NGS]) of a panel of 360 cancer-related genes. Association analyses were performed on phenotypic traits related to complete cytoreductive surgery, while logistic regression analysis was applied for the predictive model.

RESULTS

The positive influence of complete cytoreductive surgery (R0) on overall survival was confirmed (p = 0.003). Before surgery, low volumes of ascitic fluid, lower CA125 levels, higher platelet counts and relatively lower clinical stage at diagnosis were all indicators, alone and combined, for complete cytoreduction (R0). Mutations in either the chromatin remodelling SWI_SNF (p = 0.036) pathway or the histone H3K4 methylation pathway (p = 0.034) correlated with R0. The R0 group also demonstrated higher tumour mutational burden levels (p = 0.028). A predictive model was developed by combining two phenotypes and the mutational status of five genes and one genetic pathway, enabling the prediction of surgical outcomes in 87.6% of the cases in this cohort.

CONCLUSION

Inclusion of molecular biomarkers adds value to the pre-operative stratification of HGSOC patients. A potential preoperative risk stratification model combining phenotypic traits and single-gene mutational status is suggested, but the set-up needs to be validated in larger cohorts.

Computational Methods Summarizing Mutational Patterns in Cancer: Promise and Limitations for Clinical Applications

Since the rise of next-generation sequencing technologies, the catalogue of mutations in cancer has been continuously expanding. To address the complexity of the cancer-genomic landscape and extract meaningful insights, numerous computational approaches have been developed over the last two decades. In this review, we survey the current leading computational methods to derive intricate mutational patterns in the context of clinical relevance. We begin with mutation signatures, explaining first how mutation signatures were developed and then examining the utility of studies using mutation signatures to correlate environmental effects on the cancer genome. Next, we examine current clinical research that employs mutation signatures and discuss the potential use cases and challenges of mutation signatures in clinical decision-making. We then examine computational studies developing tools to investigate complex patterns of mutations beyond the context of mutational signatures. We survey methods to identify cancer-driver genes, from single-driver studies to pathway and network analyses. In addition, we review methods inferring complex combinations of mutations for clinical tasks and using mutations integrated with multi-omics data to better predict cancer phenotypes. We examine the use of these tools for either discovery or prediction, including prediction of tumor origin, treatment outcomes, prognosis, and cancer typing. We further discuss the main limitations preventing widespread clinical integration of computational tools for the diagnosis and treatment of cancer. We end by proposing solutions to address these challenges using recent advances in machine learning.

On the gene expression landscape of cancer

Kauffman picture of normal and tumor states as attractors in an abstract state space is used in order to interpret gene expression data for 15 cancer localizations obtained from The Cancer Genome Atlas. A principal component analysis of this data unveils the following qualitative aspects about tumors: 1) The state of a tissue in gene expression space can be described by a few variables. In particular, there is a single variable describing the progression from a normal tissue to a tumor. 2) Each cancer localization is characterized by a gene expression profile, in which genes have specific weights in the definition of the cancer state. There are no less than 2500 differentially-expressed genes, which lead to power-like tails in the expression distribution functions. 3) Tumors in different localizations share hundreds or even thousands of differentially expressed genes. There are 6 genes common to the 15 studied tumor localizations. 4) The tumor region is a kind of attractor. Tumors in advanced stages converge to this region independently of patient age or genetic characteristics. 5) There is a landscape of cancer in gene expression space with an approximate border separating normal tissues from tumors.

Cancer Segmentation by Entropic Analysis of Ordered Gene Expression Profiles

The availability of massive gene expression data has been challenging in terms of how to cure, process, and extract useful information. Here, we describe the use of entropic measures as discriminating criteria in cancer using the whole data set of gene expression levels. These methods were applied in classifying samples between tumor and normal type for 13 types of tumors with a high success ratio. Using gene expression, ordered by pathways, results in complexity-entropy diagrams. The map allows the clustering of the tumor and normal types samples, with a high success rate for nine of the thirteen, studied cancer types. Further analysis using information distance also shows good discriminating behavior, but, more importantly, allows for discriminating between cancer types. Together, our results allow the classification of tissues without the need to identify relevant genes or impose a particular cancer model. The used procedure can be extended to classification problems beyond the reported results.

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.

Multi-omics Uncovering Different Faces of Clear Cell Ovarian Cancer

The diagnosis of clear cell ovarian cancer relies on expert histopathology review. Further characterization from deep genomic and transcriptomic analyses can identify different subgroups. International collaboration is required to define the clinical impact and therapy opportunities in these specific subclassifications. See related article by Bolton et al., p. 4947.

Big data in basic and translational cancer research

Historically, the primary focus of cancer research has been molecular and clinical studies of a few essential pathways and genes. Recent years have seen the rapid accumulation of large-scale cancer omics data catalysed by breakthroughs in high-throughput technologies. This fast data growth has given rise to an evolving concept of 'big data' in cancer, whose analysis demands large computational resources and can potentially bring novel insights into essential questions. Indeed, the combination of big data, bioinformatics and artificial intelligence has led to notable advances in our basic understanding of cancer biology and to translational advancements. Further advances will require a concerted effort among data scientists, clinicians, biologists and policymakers. Here, we review the current state of the art and future challenges for harnessing big data to advance cancer research and treatment.

Chromatin accessibility of primary human cancers ties regional mutational processes and signatures with tissues of origin

Somatic mutations in cancer genomes are associated with DNA replication timing (RT) and chromatin accessibility (CA), however these observations are based on normal tissues and cell lines while primary cancer epigenomes remain uncharacterised. Here we use machine learning to model megabase-scale mutation burden in 2,500 whole cancer genomes and 17 cancer types via a compendium of 900 CA and RT profiles covering primary cancers, normal tissues, and cell lines. CA profiles of primary cancers, rather than those of normal tissues, are most predictive of regional mutagenesis in most cancer types. Feature prioritisation shows that the epigenomes of matching cancer types and organ systems are often the strongest predictors of regional mutation burden, highlighting disease-specific associations of mutational processes. The genomic distributions of mutational signatures are also shaped by the epigenomes of matched cancer and tissue types, with SBS5/40, carcinogenic and unknown signatures most accurately predicted by our models. In contrast, fewer associations of RT and regional mutagenesis are found. Lastly, the models highlight genomic regions with overrepresented mutations that dramatically exceed epigenome-derived expectations and show a pan-cancer convergence to genes and pathways involved in development and oncogenesis, indicating the potential of this approach for coding and non-coding driver discovery. The association of regional mutational processes with the epigenomes of primary cancers suggests that the landscape of passenger mutations is predominantly shaped by the epigenomes of cancer cells after oncogenic transformation.

Individualized discovery of rare cancer drivers in global network context

Late advances in genome sequencing expanded the space of known cancer driver genes several-fold. However, most of this surge was based on computational analysis of somatic mutation frequencies and/or their impact on the protein function. On the contrary, experimental research necessarily accounted for functional context of mutations interacting with other genes and conferring cancer phenotypes. Eventually, just such results become ‘hard currency’ of cancer biology. The new method, NEAdriver employs knowledge accumulated thus far in the form of global interaction network and functionally annotated pathways in order to recover known and predict novel driver genes. The driver discovery was individualized by accounting for mutations’ co-occurrence in each tumour genome – as an alternative to summarizing information over the whole cancer patient cohorts. For each somatic genome change, probabilistic estimates from two lanes of network analysis were combined into joint likelihoods of being a driver. Thus, ability to detect previously unnoticed candidate driver events emerged from combining individual genomic context with network perspective. The procedure was applied to 10 largest cancer cohorts followed by evaluating error rates against previous cancer gene sets. The discovered driver combinations were shown to be informative on cancer outcome. This revealed driver genes with individually sparse mutation patterns that would not be detectable by other computational methods and related to cancer biology domains poorly covered by previous analyses. In particular, recurrent mutations of collagen, laminin, and integrin genes were observed in the adenocarcinoma and glioblastoma cancers. Considering constellation patterns of candidate drivers in individual cancer genomes opens a novel avenue for personalized cancer medicine.

Reconstruction and analysis of a large-scale binary Ras-effector signaling network

BACKGROUND

Ras is a key cellular signaling hub that controls numerous cell fates via multiple downstream effector pathways. While pathways downstream of effectors such as Raf, PI3K and RalGDS are extensively described in the literature, how other effectors signal downstream of Ras is often still enigmatic.

METHODS

A comprehensive and unbiased Ras-effector network was reconstructed downstream of 43 effector proteins (converging onto 12 effector classes) using public pathway and protein-protein interaction (PPI) databases. The output is an oriented graph of pairwise interactions defining a 3-layer signaling network downstream of Ras. The 2290 proteins comprising the network were studied for their implication in signaling crosstalk and feedbacks, their subcellular localizations, and their cellular functions.

RESULTS

The final Ras-effector network consists of 2290 proteins that are connected via 19,080 binary PPIs, increasingly distributed across the downstream layers, with 441 PPIs in layer 1, 1660 in layer 2, and 16,979 in layer 3. We identified a high level of crosstalk among proteins of the 12 effector classes. A class-specific Ras sub-network was generated in CellDesigner (.xml file) and a functional enrichment analysis thereof shows that 58% of the processes have previously been associated to a respective effector pathway, with the remaining providing insights into novel and unexplored functions of specific effector pathways.

CONCLUSIONS

Our large-scale and cell general Ras-effector network is a crucial steppingstone towards defining the network boundaries. It constitutes a 'reference interactome' and can be contextualized for specific conditions, e.g. different cell types or biopsy material obtained from cancer patients. Further, it can serve as a basis for elucidating systems properties, such as input-output relationships, crosstalk, and pathway redundancy. Video Abstract.

Extracting functional insights from loss-of-function screens using deep link prediction

We present deep link prediction (DLP), a method for the interpretation of loss-of-function screens. Our approach uses representation-based link prediction to reprioritize phenotypic readouts by integrating screening experiments with gene-gene interaction networks. We validate on 2 different loss-of-function technologies, RNAi and CRISPR, using datasets obtained from DepMap. Extensive benchmarking shows that DLP-DeepWalk outperforms other methods in recovering cell-specific dependencies, achieving an average precision well above 90% across 7 different cancer types and on both RNAi and CRISPR data. We show that the genes ranked highest by DLP-DeepWalk are appreciably more enriched in drug targets compared to the ranking based on original screening scores. Interestingly, this enrichment is more pronounced on RNAi data compared to CRISPR data, consistent with the greater inherent noise of RNAi screens. Finally, we demonstrate how DLP-DeepWalk can infer the molecular mechanism through which putative targets trigger cell line mortality.

Leveraging biochemical reactions to unravel functional impacts of cancer somatic variants affecting protein interaction interfaces

Background: Considering protein mutations in their biological context is essential for understanding their functional impact, interpretation of high-dimensional datasets and development of effective targeted therapies in personalized medicine. Methods: We combined the curated knowledge of biochemical reactions from Reactome with the analysis of interaction-mediating 3D interfaces from Mechismo. In addition, we provided a software tool for users to explore and browse the analysis results in a multi-scale perspective starting from pathways and reactions to protein-protein interactions and protein 3D structures. Results: We analyzed somatic mutations from TCGA, revealing several significantly impacted reactions and pathways in specific cancer types. We found examples of genes not yet listed as oncodrivers, whose rare mutations were predicted to affect cancer processes similarly to known oncodrivers. Some identified processes lack any known oncodrivers, which suggests potentially new cancer-related processes (e.g. complement cascade reactions). Furthermore, we found that mutations perturbing certain processes are significantly associated with distinct phenotypes (i.e. survival time) in specific cancer types (e.g. PIK3CA centered pathways in LGG and UCEC cancer types), suggesting the translational potential of our approach for patient stratification. Our analysis also uncovered several druggable processes (e.g. GPCR signalling pathways) containing enriched reactions, providing support for new off-label therapeutic options. Conclusions: In summary, we have established a multi-scale approach to study genetic variants based on protein-protein interaction 3D structures. Our approach is different from previously published studies in its focus on biochemical reactions and can be applied to other data types (e.g. post-translational modifications) collected for many types of disease.

Leveraging biochemical reactions to unravel functional impacts of cancer somatic variants affecting protein interaction interfaces

Background: Considering protein mutations in their biological context is essential for understanding their functional impact, interpretation of high-dimensional datasets and development of effective targeted therapies in personalized medicine. Methods: We combined the curated knowledge of biochemical reactions from Reactome with the analysis of interaction-mediating 3D interfaces from Mechismo. In addition, we provided a software tool for users to explore and browse the analysis results in a multi-scale perspective starting from pathways and reactions to protein-protein interactions and protein 3D structures. Results: We analyzed somatic mutations from TCGA, revealing several significantly impacted reactions and pathways in specific cancer types. We found examples of genes not yet listed as oncodrivers, whose rare mutations were predicted to affect cancer processes similarly to known oncodrivers. Some identified processes lack any known oncodrivers, which suggests potentially new cancer-related processes (e.g. complement cascade reactions). Furthermore, we found that mutations perturbing certain processes are significantly associated with distinct phenotypes (i.e. survival time) in specific cancer types (e.g. PIK3CA centered pathways in LGG and UCEC cancer types), suggesting the translational potential of our approach for patient stratification. Our analysis also uncovered several druggable processes (e.g. GPCR signalling pathways) containing enriched reactions, providing support for new off-label therapeutic options. Conclusions: In summary, we have established a multi-scale approach to study genetic variants based on protein-protein interaction 3D structures. Our approach is different from previously published studies in its focus on biochemical reactions and can be applied to other data types (e.g. post-translational modifications) collected for many types of disease.

Leveraging biochemical reactions to unravel functional impacts of cancer somatic variants affecting protein interaction interfaces

Background: Considering protein mutations in their biological context is essential for understanding their functional impact, interpretation of high-dimensional datasets and development of effective targeted therapies in personalized medicine. Methods: We combined the curated knowledge of biochemical reactions from Reactome with the analysis of interaction-mediating 3D interfaces from Mechismo. In addition, we provided a software tool for users to explore and browse the analysis results in a multi-scale perspective starting from pathways and reactions to protein-protein interactions and protein 3D structures. Results: We analyzed somatic mutations from TCGA, revealing several significantly impacted reactions and pathways in specific cancer types. We found examples of genes not yet listed as oncodrivers, whose rare mutations were predicted to affect cancer processes similarly to known oncodrivers. Some identified processes lack any known oncodrivers, which suggests potentially new cancer-related processes (e.g. complement cascade reactions). Furthermore, we found that mutations perturbing certain processes are significantly associated with distinct phenotypes (i.e. survival time) in specific cancer types (e.g. PIK3CA centered pathways in LGG and UCEC cancer types), suggesting the translational potential of our approach for patient stratification. Our analysis also uncovered several druggable processes (e.g. GPCR signalling pathways) containing enriched reactions, providing support for new off-label therapeutic options. Conclusions: In summary, we have established a multi-scale approach to study genetic variants based on protein-protein interaction 3D structures. Our approach is different from previously published studies in its focus on biochemical reactions and can be applied to other data types (e.g. post-translational modifications) collected for many types of disease.

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.

Machine learning applications for therapeutic tasks with genomics data

Thanks to the increasing availability of genomics and other biomedical data, many machine learning algorithms have been proposed for a wide range of therapeutic discovery and development tasks. In this survey, we review the literature on machine learning applications for genomics through the lens of therapeutic development. We investigate the interplay among genomics, compounds, proteins, electronic health records, cellular images, and clinical texts. We identify 22 machine learning in genomics applications that span the whole therapeutics pipeline, from discovering novel targets, personalizing medicine, developing gene-editing tools, all the way to facilitating clinical trials and post-market studies. We also pinpoint seven key challenges in this field with potentials for expansion and impact. This survey examines recent research at the intersection of machine learning, genomics, and therapeutic development.

Modularity in Biological Networks

Network modeling, from the ecological to the molecular scale has become an essential tool for studying the structure, dynamics and complex behavior of living systems. Graph representations of the relationships between biological components open up a wide variety of methods for discovering the mechanistic and functional properties of biological systems. Many biological networks are organized into a modular structure, so methods to discover such modules are essential if we are to understand the biological system as a whole. However, most of the methods used in biology to this end, have a limited applicability, as they are very specific to the system they were developed for. Conversely, from the statistical physics and network science perspective, graph modularity has been theoretically studied and several methods of a very general nature have been developed. It is our perspective that in particular for the modularity detection problem, biology and theoretical physics/network science are less connected than they should. The central goal of this review is to provide the necessary background and present the most applicable and pertinent methods for community detection in a way that motivates their further usage in biological research.

A protein network map of head and neck cancer reveals PIK3CA mutant drug sensitivity

[Figure: see text].

A protein interaction landscape of breast cancer

[Figure: see text].

Interpretation of cancer mutations using a multiscale map of protein systems

A major goal of cancer research is to understand how mutations distributed across diverse genes affect common cellular systems, including multiprotein complexes and assemblies. Two challenges—how to comprehensively map such systems and how to identify which are under mutational selection—have hindered this understanding. Accordingly, we created a comprehensive map of cancer protein systems integrating both new and published multi-omic interaction data at multiple scales of analysis. We then developed a unified statistical model that pinpoints 395 specific systems under mutational selection across 13 cancer types. This map, called NeST (Nested Systems in Tumors), incorporates canonical processes and notable discoveries, including a PIK3CA-actomyosin complex that inhibits phosphatidylinositol 3-kinase signaling and recurrent mutations in collagen complexes that promote tumor proliferation. These systems can be used as clinical biomarkers and implicate a total of 548 genes in cancer evolution and progression. This work shows how disparate tumor mutations converge on protein assemblies at different scales.

Mapping genomic and epigenomic evolution in cancer ecosystems

[Figure: see text].

Transcriptomic landscape of male and female reproductive cancers: Similar pathways and molecular signatures predicting response to endocrine therapy

Reproductive cancers in both genders represent serious health problems, whose incidence has significantly risen over the past decades. Although considerable differences among reproductive cancers exist, we aimed to identify similar signaling pathways and key molecular oncomarkers shared among six human reproductive cancers that can advance the current knowledge of cancer biology to propose new strategies for more effective therapies. Using a computational analysis approach, here we uncover aberrant miRNAs-mRNAs networks shared in six reproductive tumor types, and identify common molecular mechanisms strictly associated with cancer promotion and aggressiveness. Based on the fact that estrogenic and androgenic signaling pathways were most active in prostate and breast cancers, we further demonstrated that both androgen and estrogen deprivation therapy are capable of regulating the expression of the same key molecular sensors associated with endoplasmic reticulum dysfunction and cell cycle in these cancers. Overall, our data reveal a potential mechanistic framework of cellular processes that are shared among reproductive cancers, and particularly, highlight the importance of hormonal deprivation in breast and prostate cancers and potentially new biomarkers of response to these therapeutic approaches.

Bridging Tumorigenesis and Therapy Resistance With a Non-Darwinian and Non-Lamarckian Mechanism of Adaptive Evolution

Although neo-Darwinian (and less often Lamarckian) dynamics are regularly invoked to interpret cancer's multifarious molecular profiles, they shine little light on how tumorigenesis unfolds and often fail to fully capture the frequency and breadth of resistance mechanisms. This uncertainty frames one of the most problematic gaps between science and practice in modern times. Here, we offer a theory of adaptive cancer evolution, which builds on a molecular mechanism that lies outside neo-Darwinian and Lamarckian schemes. This mechanism coherently integrates non-genetic and genetic changes, ecological and evolutionary time scales, and shifts the spotlight away from positive selection towards purifying selection, genetic drift, and the creative-disruptive power of environmental change. The surprisingly simple use-it or lose-it rationale of the proposed theory can help predict molecular dynamics during tumorigenesis. It also provides simple rules of thumb that should help improve therapeutic approaches in cancer.

Moving pan-cancer studies from basic research toward the clinic

Although all cancers share common hallmarks, we have long realized that there is no silver-bullet treatment for the disease. Many clinical oncologists specialize in a single cancer type, based predominantly on the tissue of origin. With advances brought by genetics and cancer genomic research, we now know that cancers are profoundly different, both in origins and in genetic alterations. At the same time, commonalities such as key driver mutations, altered pathways, mutational, immune and microbial signatures and other areas (many revealed by pan-cancer studies) point to the intriguing possibility of targeting common traits across diverse cancer types with the same therapeutic strategies. Studies designed to delineate differences and similarities across cancer types are thus critical in discerning the basic dynamics of oncogenesis, as well as informing diagnoses, prognoses and therapies. We anticipate growing emphases on the development and application of therapies targeting underlying commonalities of different cancer types, while tailoring to the unique tissue environment and intrinsic molecular fingerprints of each cancer type and subtype. Here we summarize the facets of pan-cancer research and how they are pushing progress toward personalized medicine.

Decoding human cancer with whole genome sequencing: a review of PCAWG Project studies published in February 2020

Cancer is underlined by genetic changes. In an unprecedented international effort, the Pan-Cancer Analysis of Whole Genomes (PCAWG) of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA) sequenced the tumors of over two thousand five hundred patients across 38 different cancer types, as well as the corresponding healthy tissue, with the aim of identifying genome-wide mutations exclusively found in cancer and uncovering new genetic changes that drive tumor formation. What set this project apart from earlier efforts is the use of whole genome sequencing (WGS) that enabled to explore alterations beyond the coding DNA, into cancer's non-coding genome. WGS of the entire cohort allowed to tease apart driving mutations that initiate and support carcinogenesis from passenger mutations that do not play an overt role in the disease. At least one causative mutation was found in 95% of all cancers, with many tumors showing an average of 5 driver mutations. The PCAWG Project also assessed the transcriptional output altered in cancer and rebuilt the evolutionary history of each tumor showing that initial driver mutations can occur years if not decades prior to a diagnosis. Here, I provide a concise review of the Pan-Cancer Project papers published on February 2020, along with key computational tools and the digital framework generated as part of the project. This represents an historic effort by hundreds of international collaborators, which provides a comprehensive understanding of cancer genetics, with publicly available data and resources representing a treasure trove of information to advance cancer research for years to come.

The current issues and future perspective of artificial intelligence for developing new treatment strategy in non-small cell lung cancer: harmonization of molecular cancer biology and artificial intelligence

Comprehensive analysis of omics data, such as genome, transcriptome, proteome, metabolome, and interactome, is a crucial technique for elucidating the complex mechanism of cancer onset and progression. Recently, a variety of new findings have been reported based on multi-omics analysis in combination with various clinical information. However, integrated analysis of multi-omics data is extremely labor intensive, making the development of new analysis technology indispensable. Artificial intelligence (AI), which has been under development in recent years, is quickly becoming an effective approach to reduce the labor involved in analyzing large amounts of complex data and to obtain valuable information that is often overlooked in manual analysis and experiments. The use of AI, such as machine learning approaches and deep learning systems, allows for the efficient analysis of massive omics data combined with accurate clinical information and can lead to comprehensive predictive models that will be desirable for further developing individual treatment strategies of immunotherapy and molecular target therapy. Here, we aim to review the potential of AI in the integrated analysis of omics data and clinical information with a special focus on recent advances in the discovery of new biomarkers and the future direction of personalized medicine in non-small lung cancer.

Identification and Functional Characterization of Novel MYC-Regulated Long Noncoding RNAs in Group 3 Medulloblastoma

Simple Summary

Medulloblastoma is the most common malignant pediatric brain tumor, which accounts for approximately 20% of all childhood brain tumors. To date, no pharmacological approaches are decisive in the treatment of this cancer, while the secondary effects of conventional therapies as chemotherapy, radiotherapy or surgical interventions heavily affect the quality of life of patients. This requires the rapid development of alternative molecular therapies, which are the future challenge of personalized medicine. In this context, we addressed our research towards the most aggressive form of Medulloblastoma to identify novel genes responsible for its onset and/or progression. We discovered three newly implicated genes, for which we highlighted a contribution in the control of cancer cell features. Deepening into the Medulloblastoma biology, this study represents a further step forward for the development of molecular therapies in the era of precision oncology.

Abstract

The impact of protein-coding genes on cancer onset and progression is a well-established paradigm in molecular oncology. Nevertheless, unveiling the contribution of the noncoding genes—including long noncoding RNAs (lncRNAs)—to tumorigenesis represents a great challenge for personalized medicine, since they (i) constitute the majority of the human genome, (ii) are essential and flexible regulators of gene expression and (iii) present all types of genomic alterations described for protein-coding genes. LncRNAs have been increasingly associated with cancer, their highly tissue- and cancer type-specific expression making them attractive candidates as both biomarkers and therapeutic targets. Medulloblastoma is one of the most common malignant pediatric brain tumors. Group 3 is the most aggressive subgroup, showing the highest rate of metastasis at diagnosis. Transcriptomics and reverse genetics approaches were combined to identify lncRNAs implicated in Group 3 Medulloblastoma biology. Here we present the first collection of lncRNAs dependent on the activity of the MYC oncogene, the major driver gene of Group 3 Medulloblastoma. We assessed the expression profile of selected lncRNAs in Group 3 primary tumors and functionally characterized these species. Overall, our data demonstrate the direct involvement of three lncRNAs in Medulloblastoma cancer cell phenotypes.

Integrative Analysis of Next-Generation Sequencing for Next-Generation Cancer Research toward Artificial Intelligence

The rapid improvement of next-generation sequencing (NGS) technologies and their application in large-scale cohorts in cancer research led to common challenges of big data. It opened a new research area incorporating systems biology and machine learning. As large-scale NGS data accumulated, sophisticated data analysis methods became indispensable. In addition, NGS data have been integrated with systems biology to build better predictive models to determine the characteristics of tumors and tumor subtypes. Therefore, various machine learning algorithms were introduced to identify underlying biological mechanisms. In this work, we review novel technologies developed for NGS data analysis, and we describe how these computational methodologies integrate systems biology and omics data. Subsequently, we discuss how deep neural networks outperform other approaches, the potential of graph neural networks (GNN) in systems biology, and the limitations in NGS biomedical research. To reflect on the various challenges and corresponding computational solutions, we will discuss the following three topics: (i) molecular characteristics, (ii) tumor heterogeneity, and (iii) drug discovery. We conclude that machine learning and network-based approaches can add valuable insights and build highly accurate models. However, a well-informed choice of learning algorithm and biological network information is crucial for the success of each specific research question.

Diagnostic and prognostic potential of the proteomic profiling of serum-derived extracellular vesicles in prostate cancer

Extracellular vesicles (EVs) and their cargo represent an intriguing source of cancer biomarkers for developing robust and sensitive molecular tests by liquid biopsy. Prostate cancer (PCa) is still one of the most frequent and deadly tumor in men and analysis of EVs from biological fluids of PCa patients has proven the feasibility and the unprecedented potential of such an approach. Here, we exploited an antibody-based proteomic technology, i.e. the Reverse-Phase Protein microArrays (RPPA), to measure key antigens and activated signaling in EVs isolated from sera of PCa patients. Notably, we found tumor-specific protein profiles associated with clinical settings as well as candidate markers for EV-based tumor diagnosis. Among others, PD-L1, ERG, Integrin-β5, Survivin, TGF-β, phosphorylated-TSC2 as well as partners of the MAP-kinase and mTOR pathways emerged as differentially expressed endpoints in tumor-derived EVs. In addition, the retrospective analysis of EVs from a 15-year follow-up cohort generated a protein signature with prognostic significance. Our results confirm that serum-derived EV cargo may be exploited to improve the current diagnostic procedures while providing potential prognostic and predictive information. The approach proposed here has been already applied to tumor entities other than PCa, thus proving its value in translational medicine and paving the way to innovative, clinically meaningful tools.

Construction of a prognostic model with histone modification-related genes and identification of potential drugs in pancreatic cancer

Background

Pancreatic cancer (PC) is a highly fatal and aggressive disease with its incidence and mortality quite discouraging. An effective prediction model is urgently needed for the accurate assessment of patients’ prognosis to assist clinical decision-making.

Methods

Gene expression data and clinicopathological data of the samples were acquired from The Cancer Genome Atlas (TCGA), Genotype-Tissue Expression (GTEx), and Gene Expression Omnibus (GEO) databases. Differential expressed genes (DEGs) analysis, univariate Cox regression analysis, least absolute shrinkage and selection operator (LASSO) regression analysis, random forest screening and multivariate Cox regression analysis were applied to construct the risk signature. The effectiveness and independence of the model were validated by time-dependent receiver operating characteristic (ROC) curve, Kaplan–Meier (KM) survival analysis and survival point graph in training set, test set, TCGA entire set and GSE57495 set. The validity of the core gene was verified by immunohistochemistry and our own independent cohort. Meanwhile, functional enrichment analysis of DEGs between the high and low risk groups revealed the potential biological pathways. Finally, CMap database and drug sensitivity assay were utilized to identify potential small molecular drugs as the risk model-related treatments for PC patients.

Results

Four histone modification-related genes were identified to establish the risk signature, including CBX8, CENPT, DPY30 and PADI1. The predictive performance of risk signature was validated in training set, test set, TCGA entire set and GSE57495 set, with the areas under ROC curve (AUCs) for 3-year survival were 0.773, 0.729, 0.775 and 0.770 respectively. Furthermore, KM survival analysis, univariate and multivariate Cox regression analysis proved it as an independent prognostic factor. Mechanically, functional enrichment analysis showed that the poor prognosis of high-risk population was related to the metabolic disorders caused by inadequate insulin secretion, which was fueled by neuroendocrine aberration. Lastly, a cluster of small molecule drugs were identified with significant potentiality in treating PC patients.

Conclusions

Based on a histone modification-related gene signature, our model can serve as a reliable prognosis assessment tool and help to optimize the treatment for PC patients. Meanwhile, a cluster of small molecule drugs were also identified with significant potentiality in treating PC patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-01928-6.

Functional and genetic determinants of mutation rate variability in regulatory elements of cancer genomes

Background

Cancer genomes are shaped by mutational processes with complex spatial variation at multiple scales. Entire classes of regulatory elements are affected by local variations in mutation frequency. However, the underlying mechanisms with functional and genetic determinants remain poorly understood.

Results

We characterise the mutational landscape of 1.3 million gene-regulatory and chromatin architectural elements in 2419 whole cancer genomes with transcriptional and pathway activity, functional conservation and recurrent driver events. We develop RM2, a statistical model that quantifies mutational enrichment or depletion in classes of genomic elements through genetic, trinucleotide and megabase-scale effects. We report a map of localised mutational processes affecting CTCF binding sites, transcription start sites (TSS) and tissue-specific open-chromatin regions. Increased mutation frequency in TSSs associates with mRNA abundance in most cancer types, while open-chromatin regions are generally enriched in mutations. We identify ~ 10,000 CTCF binding sites with core DNA motifs and constitutive binding in 66 cell types that represent focal points of mutagenesis. We detect site-specific mutational signature enrichments, such as SBS40 in open-chromatin regions in prostate cancer and SBS17b in CTCF binding sites in gastrointestinal cancers. Candidate drivers of localised mutagenesis are also apparent: BRAF mutations associate with mutational enrichments at CTCF binding sites in melanoma, and ARID1A mutations with TSS-specific mutagenesis in pancreatic cancer.

Conclusions

Our method and catalogue of localised mutational processes provide novel perspectives to cancer genome evolution, mutagenesis, DNA repair and driver gene discovery. The functional and genetic correlates of mutational processes suggest mechanistic hypotheses for future studies.

ActiveDriverDB: Interpreting Genetic Variation in Human and Cancer Genomes Using Post-translational Modification Sites and Signaling Networks (2021 Update)

Deciphering the functional impact of genetic variation is required to understand phenotypic diversity and the molecular mechanisms of inherited disease and cancer. While millions of genetic variants are now mapped in genome sequencing projects, distinguishing functional variants remains a major challenge. Protein-coding variation can be interpreted using post-translational modification (PTM) sites that are core components of cellular signaling networks controlling molecular processes and pathways. ActiveDriverDB is an interactive proteo-genomics database that uses more than 260,000 experimentally detected PTM sites to predict the functional impact of genetic variation in disease, cancer and the human population. Using machine learning tools, we prioritize proteins and pathways with enriched PTM-specific amino acid substitutions that potentially rewire signaling networks via induced or disrupted short linear motifs of kinase binding. We then map these effects to site-specific protein interaction networks and drug targets. In the 2021 update, we increased the PTM datasets by nearly 50%, included glycosylation, sumoylation and succinylation as new types of PTMs, and updated the workflows to interpret inherited disease mutations. We added a recent phosphoproteomics dataset reflecting the cellular response to SARS-CoV-2 to predict the impact of human genetic variation on COVID-19 infection and disease course. Overall, we estimate that 16-21% of known amino acid substitutions affect PTM sites among pathogenic disease mutations, somatic mutations in cancer genomes and germline variants in the human population. These data underline the potential of interpreting genetic variation through the lens of PTMs and signaling networks. The open-source database is freely available at www.ActiveDriverDB.org.

Molecular characterization of a marine turtle tumor epizootic, profiling external, internal and postsurgical regrowth tumors

Sea turtle populations are under threat from an epizootic tumor disease (animal epidemic) known as fibropapillomatosis. Fibropapillomatosis continues to spread geographically, with prevalence of the disease also growing at many longer-affected sites globally. However, we do not yet understand the precise environmental, mutational and viral events driving fibropapillomatosis tumor formation and progression.Here we perform transcriptomic and immunohistochemical profiling of five fibropapillomatosis tumor types: external new, established and postsurgical regrowth tumors, and internal lung and kidney tumors. We reveal that internal tumors are molecularly distinct from the more common external tumors. However, they have a small number of conserved potentially therapeutically targetable molecular vulnerabilities in common, such as the MAPK, Wnt, TGFβ and TNF oncogenic signaling pathways. These conserved oncogenic drivers recapitulate remarkably well the core pan-cancer drivers responsible for human cancers. Fibropapillomatosis has been considered benign, but metastatic-related transcriptional signatures are strongly activated in kidney and established external tumors. Tumors in turtles with poor outcomes (died/euthanized) have genes associated with apoptosis and immune function suppressed, with these genes providing putative predictive biomarkers.Together, these results offer an improved understanding of fibropapillomatosis tumorigenesis and provide insights into the origins, inter-tumor relationships, and therapeutic treatment for this wildlife epizootic.

Prolactin synergizes with canonical Wnt signals to drive development of ER+ mammary tumors via activation of the Notch pathway

Prolactin (PRL) cooperates with other factors to orchestrate mammary development and lactation, and is epidemiologically linked to higher risk for breast cancer. However, how PRL collaborates with oncogenes to foster tumorigenesis and influence breast cancer phenotype is not well understood. To understand its interactions with canonical Wnt signals, which elevate mammary stem cell activity, we crossed heterozygous NRL-PRL mice with Apc^Min/+ mice and treated pubertal females with a single dose of mutagen. PRL in the context of Apc^Min/+ fueled a dramatic increase in tumor incidence in nulliparous mice, compared to Apc^Min/+ alone. Although carcinomas in both NRL-PRL/Apc^Min/+ and Apc^Min/+ females acquired a mutation in the remaining wildtype Apc allele and expressed abundant β-catenin, PRL-promoted tumors displayed higher levels of Notch-driven target genes and Notch-dependent cancer stem cell activity, compared to β-catenin-driven activity in Apc^Min/+ tumors. This PRL-induced shift to dominant Notch signals was evident in preneoplastic epithelial hyperplasias at 120 days of age. In NRL-PRL/Apc^Min/+ females, rapidly proliferating hyperplasias, characterized by β-catenin at cell junctions and high NOTCH1 expression, contrasted with slower growing lesions with nuclear β-catenin in Apc^Min/+ females. These studies demonstrate that PRL can powerfully modulate the incidence and phenotype of mammary tumors, shedding light on mechanisms whereby PRL elevates risk of breast cancer.

Context Matters: NOTCH Signatures and Pathway in Cancer Progression and Metastasis

The Notch signaling pathway is a critical player in embryogenesis but also plays various roles in tumorigenesis, with both tumor suppressor and oncogenic activities. Mutations, deletions, amplifications, or over-expression of Notch receptors, ligands, and a growing list of downstream Notch-activated genes have by now been described for most human cancer types. Yet, it often remains unclear what may be the functional impact of these changes for tumor biology, initiation, and progression, for cancer therapy, and for personalized medicine. Emerging data indicate that Notch signaling can also contribute to increased aggressive properties such as invasion, tumor heterogeneity, angiogenesis, or tumor cell dormancy within solid cancer tissues; especially in epithelial cancers, which are in the center of this review. Notch further supports the “stemness” of cancer cells and helps define the stem cell niche for their long-term survival, by integrating the interaction between cancer cells and the cells of the tumor microenvironment (TME). The complexity of Notch crosstalk with other signaling pathways and its roles in cell fate and trans-differentiation processes such as epithelial-to-mesenchymal transition (EMT) point to this pathway as a decisive player that may tip the balance between tumor suppression and promotion, differentiation and invasion. Here we not only review the literature, but also explore genomic databases with a specific focus on Notch signatures, and how they relate to different stages in tumor development. Altered Notch signaling hereby plays a key role for tumor cell survival and coping with a broad spectrum of vital issues, contributing to failed therapies, poor patient outcome, and loss of lives.