The cancer cell map initiative: defining the hallmark networks of cancer

Revisiting the role of mesenchymal stromal cells in cancer initiation, metastasis and immunosuppression

Immune checkpoint blockade (ICB) necessitates a thorough understanding of intricate cellular interactions within the tumor microenvironment (TME). Mesenchymal stromal cells (MSCs) play a pivotal role in cancer generation, progression, and immunosuppressive tumor microenvironment. Within the TME, MSCs encompass both resident and circulating counterparts that dynamically communicate and actively participate in TME immunosurveillance and response to ICB. This review aims to reevaluate various facets of MSCs, including their potential self-transformation to function as cancer-initiating cells and contributions to the creation of a conducive environment for tumor proliferation and metastasis. Additionally, we explore the immune regulatory functions of tumor-associated MSCs (TA-MSCs) and MSC-derived extracellular vesicles (MSC-EVs) with analysis of potential connections between circulating and tissue-resident MSCs. A comprehensive understanding of the dynamics of MSC-immune cell communication and the heterogeneous cargo of tumor-educated versus naïve MSCs may unveil a new MSC-mediated immunosuppressive pathway that can be targeted to enhance cancer control by ICB.

Predicting anti-cancer drug combination responses with a temporal cell state network model

Cancer chemotherapy combines multiple drugs, but predicting the effects of drug combinations on cancer cell proliferation remains challenging, even for simple in vitro systems. We hypothesized that by combining knowledge of single drug dose responses and cell state transition network dynamics, we could predict how a population of cancer cells will respond to drug combinations. We tested this hypothesis here using three targeted inhibitors of different cell cycle states in two different cell lines in vitro. We formulated a Markov model to capture temporal cell state transitions between different cell cycle phases, with single drug data constraining how drug doses affect transition rates. This model was able to predict the landscape of all three different pairwise drug combinations across all dose ranges for both cell lines with no additional data. While further application to different cell lines, more drugs, additional cell state networks, and more complex co-culture or in vivo systems remain, this work demonstrates how currently available or attainable information could be sufficient for prediction of drug combination response for single cell lines in vitro.

Community Detection in Protein-Protein Interaction Networks and Applications

The ability to identify and characterize not only the protein-protein interactions but also their internal modular organization through network analysis is fundamental for understanding the mechanisms of biological processes at the molecular level. Indeed, the detection of the network communities can enhance our understanding of the molecular basis of disease pathology, and promote drug discovery and disease treatment in personalized medicine. This work gives an overview of recent computational methods for the detection of protein complexes and functional modules in protein-protein interaction networks, also providing a focus on some of its applications. We propose a systematic reformulation of frequently adopted taxonomies for these methods, also proposing new categories to keep up with the most recent research. We review the literature of the last five years (2017-2021) and provide links to existing data and software resources. Finally, we survey recent works exploiting module identification and analysis, in the context of a variety of disease processes for biomarker identification and therapeutic target detection. Our review provides the interested reader with an up-to-date and self-contained view of the existing research, with links to state-of-the-art literature and resources, as well as hints on open issues and future research directions in complex detection and its applications.

Discovery of novel IDH1-R132C inhibitors through structure-based virtual screening

Isocitrate dehydrogenase (IDH) belongs to a family of enzymes involved in glycometabolism. It is found in many living organisms and is one of the most mutated metabolic enzymes. In the current study, we identified novel IDH1-R132C inhibitors using docking-based virtual screening and cellular inhibition assays. A total of 100 molecules with high docking scores were obtained from docking-based virtual screening. The cellular inhibition assay demonstrated five compounds at a concentration of 10 μM could inhibit cancer cells harboring the IDH1-R132C mutation proliferation by > 50%. The compound (T001-0657) showed the most potent effect against cancer cells harboring the IDH1-R132C mutation with a half-maximal inhibitory concentration (IC₅₀) value of 1.311 μM. It also showed a cytotoxic effect against cancer cells with wild-type IDH1 and normal cells with IC₅₀ values of 49.041 μM and >50 μM, respectively. Molecular dynamics simulations were performed to investigate the stability of the kinase structure binding of allosteric inhibitor compound A and the identified compound T001-0657 binds to IDH1-R132C. Root-mean-square deviation, root-mean-square fluctuation, and binding free energy calculations showed that both compounds bind tightly to IDH1-R132C. In conclusion, the compound identified in this study had high selectivity for cancer cells harboring IDH1-R132C mutation and could be considered a promising hit compound for further development of IDH1-R132C inhibitors.

Combinatorial Inactivation of Tumor Suppressors Efficiently Initiates Lung Adenocarcinoma with Therapeutic Vulnerabilities

Lung cancer is the leading cause of cancer death worldwide, with lung adenocarcinoma being the most common subtype. Many oncogenes and tumor suppressor genes are altered in this cancer type, and the discovery of oncogene mutations has led to the development of targeted therapies that have improved clinical outcomes. However, a large fraction of lung adenocarcinomas lacks mutations in known oncogenes, and the genesis and treatment of these oncogene-negative tumors remain enigmatic. Here, we perform iterative in vivo functional screens using quantitative autochthonous mouse model systems to uncover the genetic and biochemical changes that enable efficient lung tumor initiation in the absence of oncogene alterations. Generation of hundreds of diverse combinations of tumor suppressor alterations demonstrates that inactivation of suppressors of the RAS and PI3K pathways drives the development of oncogene-negative lung adenocarcinoma. Human genomic data and histology identified RAS/MAPK and PI3K pathway activation as a common feature of an event in oncogene-negative human lung adenocarcinomas. These Onc-negativeRAS/PI3K tumors and related cell lines are vulnerable to pharmacologic inhibition of these signaling axes. These results transform our understanding of this prevalent yet understudied subtype of lung adenocarcinoma.

SIGNIFICANCE

To address the large fraction of lung adenocarcinomas lacking mutations in proto-oncogenes for which targeted therapies are unavailable, this work uncovers driver pathways of oncogene-negative lung adenocarcinomas and demonstrates their therapeutic vulnerabilities.

Artificial intelligence, molecular subtyping, biomarkers, and precision oncology

A targeted cancer therapy is only useful if there is a way to accurately identify the tumors that are susceptible to that therapy. Thus rapid expansion in the number of available targeted cancer treatments has been accompanied by a robust effort to subdivide the traditional histological and anatomical tumor classifications into molecularly defined subtypes. This review highlights the history of the paired evolution of targeted therapies and biomarkers, reviews currently used methods for subtype identification, and discusses challenges to the implementation of precision oncology as well as possible solutions.

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

[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.

A protein interaction landscape of breast cancer

[Figure: see text].

Predicting anticancer hyperfoods with graph convolutional networks

Background

Recent efforts in the field of nutritional science have allowed the discovery of disease-beating molecules within foods based on the commonality of bioactive food molecules to FDA-approved drugs. The pioneering work in this field used an unsupervised network propagation algorithm to learn the systemic-wide effect on the human interactome of 1962 FDA-approved drugs and a supervised algorithm to predict anticancer therapeutics using the learned representations. Then, a set of bioactive molecules within foods was fed into the model, which predicted molecules with cancer-beating potential.The employed methodology consisted of disjoint unsupervised feature generation and classification tasks, which can result in sub-optimal learned drug representations with respect to the classification task. Additionally, due to the disjoint nature of the tasks, the employed approach proved cumbersome to optimize, requiring testing of thousands of hyperparameter combinations and significant computational resources.To overcome the technical limitations highlighted above, we represent each drug as a graph (human interactome) with its targets as binary node features on the graph and formulate the problem as a graph classification task. To solve this task, inspired by the success of graph neural networks in graph classification problems, we use an end-to-end graph neural network model operating directly on the graphs, which learns drug representations to optimize model performance in the prediction of anticancer therapeutics.

Results

The proposed model outperforms the baseline approach in the anticancer therapeutic prediction task, achieving an F1 score of 67.99%±2.52% and an AUPR of 73.91%±3.49%. It is also shown that the model is able to capture knowledge of biological pathways to predict anticancer molecules based on the molecules’ effects on cancer-related pathways.

Conclusions

We introduce an end-to-end graph convolutional model to predict cancer-beating molecules within food. The introduced model outperforms the existing baseline approach, and shows interpretability, paving the way to the future of a personalized nutritional science approach allowing the development of nutrition strategies for cancer prevention and/or therapeutics.

Supplementary Information

The online version contains supplementary material available at (10.1186/s40246-021-00333-4).

Network-guided search for genetic heterogeneity between gene pairs

MOTIVATION

Correlating genetic loci with a disease phenotype is a common approach to improve our understanding of the genetics underlying complex diseases. Standard analyses mostly ignore two aspects, namely genetic heterogeneity and interactions between loci. Genetic heterogeneity, the phenomenon that genetic variants at different loci lead to the same phenotype, promises to increase statistical power by aggregating low-signal variants. Incorporating interactions between loci results in a computational and statistical bottleneck due to the vast amount of candidate interactions.

RESULTS

We propose a novel method SiNIMin that addresses these two aspects by finding pairs of interacting genes that are, upon combination, associated with a phenotype of interest under a model of genetic heterogeneity. We guide the interaction search using biological prior knowledge in the form of protein-protein interaction networks. Our method controls type I error and outperforms state-of-the-art methods with respect to statistical power. Additionally, we find novel associations for multiple Arabidopsis thaliana phenotypes, and, with an adapted variant of SiNIMin, for a study of rare variants in migraine patients.

AVAILABILITY AND IMPLEMENTATION

Code available at https://github.com/BorgwardtLab/SiNIMin.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Integrating 3D structural information into systems biology

Systems biology is a data-heavy field that focuses on systems-wide depictions of biological phenomena necessarily sacrificing a detailed characterization of individual components. As an example, genome-wide protein interaction networks are widely used in systems biology and continuously extended and refined as new sources of evidence become available. Despite the vast amount of information about individual protein structures and protein complexes that has accumulated in the past 50 years in the Protein Data Bank, the data, computational tools, and language of structural biology are not an integral part of systems biology. However, increasing effort has been devoted to this integration, and the related literature is reviewed here. Relationships between proteins that are detected via structural similarity offer a rich source of information not available from sequence similarity, and homology modeling can be used to leverage Protein Data Bank structures to produce 3D models for a significant fraction of many proteomes. A number of structure-informed genomic and cross-species (i.e., virus–host) interactomes will be described, and the unique information they provide will be illustrated with a number of examples. Tissue- and tumor-specific interactomes have also been developed through computational strategies that exploit patient information and through genetic interactions available from increasingly sensitive screens. Strategies to integrate structural information with these alternate data sources will be described. Finally, efforts to link protein structure space with chemical compound space offer novel sources of information in drug design, off-target identification, and the identification of targets for compounds found to be effective in phenotypic screens.

Oncoprotein-specific molecular interaction maps (SigMaps) for cancer network analyses

Tumor-specific elucidation of physical and functional oncoprotein interactions could improve tumorigenic mechanism characterization and therapeutic response prediction. Current interaction models and pathways, however, lack context specificity and are not oncoprotein specific. We introduce SigMaps as context-specific networks, comprising modulators, effectors and cognate binding-partners of a specific oncoprotein. SigMaps are reconstructed de novo by integrating diverse evidence sources-including protein structure, gene expression and mutational profiles-via the OncoSig machine learning framework. We first generated a KRAS-specific SigMap for lung adenocarcinoma, which recapitulated published KRAS biology, identified novel synthetic lethal proteins that were experimentally validated in three-dimensional spheroid models and established uncharacterized crosstalk with RAB/RHO. To show that OncoSig is generalizable, we first inferred SigMaps for the ten most mutated human oncoproteins and then for the full repertoire of 715 proteins in the COSMIC Cancer Gene Census. Taken together, these SigMaps show that the cell's regulatory and signaling architecture is highly tissue specific.

Mass spectrometry-based protein-protein interaction networks for the study of human diseases

A better understanding of the molecular mechanisms underlying disease is key for expediting the development of novel therapeutic interventions. Disease mechanisms are often mediated by interactions between proteins. Insights into the physical rewiring of protein-protein interactions in response to mutations, pathological conditions, or pathogen infection can advance our understanding of disease etiology, progression, and pathogenesis and can lead to the identification of potential druggable targets. Advances in quantitative mass spectrometry (MS)-based approaches have allowed unbiased mapping of these disease-mediated changes in protein-protein interactions on a global scale. Here, we review MS techniques that have been instrumental for the identification of protein-protein interactions at a system-level, and we discuss the challenges associated with these methodologies as well as novel MS advancements that aim to address these challenges. An overview of examples from diverse disease contexts illustrates the potential of MS-based protein-protein interaction mapping approaches for revealing disease mechanisms, pinpointing new therapeutic targets, and eventually moving toward personalized applications.

Occurrence of Inborn Errors of Metabolism in Newborns, Diagnosis and Prophylaxis

Inborn errors of metabolism (IEM) are a heterogeneous group of rare genetic disorders that are generally transmitted as autosomal or X-linked recessive disorders. These defects arise due to mutations associated with specific gene(s), especially the ones associated with key metabolic enzymes. These enzymes or their product(s) are involved in various metabolic pathways, leading to the accumulation of intermediary metabolite(s), reflecting their toxic effects upon mutations. The diagnosis of these metabolic disorders is based on the biochemical analysis of the clinical manifestations produced and their molecular mechanism. Therefore, it is imperative to devise diagnostic tests with high sensitivity and specificity for early detection of IEM. Recent advances in biochemical and polymerase chain reaction-based genetic analysis along with pedigree and prenatal diagnosis can be life-saving in nature. The latest development in exome sequencing for rapid diagnosis and enzyme replacement therapy would facilitate the successful treatment of these metabolic disorders in the future. However, the longterm clinical implications of these genetic manipulations is still a matter of debate among intellectuals and requires further research.

Metabolic Reprogramming of Fibroblasts as Therapeutic Target in Rheumatoid Arthritis and Cancer: Deciphering Key Mechanisms Using Computational Systems Biology Approaches

Fibroblasts, the most abundant cells in the connective tissue, are key modulators of the extracellular matrix (ECM) composition. These spindle-shaped cells are capable of synthesizing various extracellular matrix proteins and collagen. They also provide the structural framework (stroma) for tissues and play a pivotal role in the wound healing process. While they are maintainers of the ECM turnover and regulate several physiological processes, they can also undergo transformations responding to certain stimuli and display aggressive phenotypes that contribute to disease pathophysiology. In this review, we focus on the metabolic pathways of glucose and highlight metabolic reprogramming as a critical event that contributes to the transition of fibroblasts from quiescent to activated and aggressive cells. We also cover the emerging evidence that allows us to draw parallels between fibroblasts in autoimmune disorders and more specifically in rheumatoid arthritis and cancer. We link the metabolic changes of fibroblasts to the toxic environment created by the disease condition and discuss how targeting of metabolic reprogramming could be employed in the treatment of such diseases. Lastly, we discuss Systems Biology approaches, and more specifically, computational modeling, as a means to elucidate pathogenetic mechanisms and accelerate the identification of novel therapeutic targets.

Multi-resolution visualization and analysis of biomolecular networks through hierarchical community detection and web-based graphical tools

The visual exploration and analysis of biomolecular networks is of paramount importance for identifying hidden and complex interaction patterns among proteins. Although many tools have been proposed for this task, they are mainly focused on the query and visualization of a single protein with its neighborhood. The global exploration of the entire network and the interpretation of its underlying structure still remains difficult, mainly due to the excessively large size of the biomolecular networks. In this paper we propose a novel multi-resolution representation and exploration approach that exploits hierarchical community detection algorithms for the identification of communities occurring in biomolecular networks. The proposed graphical rendering combines two types of nodes (protein and communities) and three types of edges (protein-protein, community-community, protein-community), and displays communities at different resolutions, allowing the user to interactively zoom in and out from different levels of the hierarchy. Links among communities are shown in terms of relationships and functional correlations among the biomolecules they contain. This form of navigation can be also combined by the user with a vertex centric visualization for identifying the communities holding a target biomolecule. Since communities gather limited-size groups of correlated proteins, the visualization and exploration of complex and large networks becomes feasible on off-the-shelf computer machines. The proposed graphical exploration strategies have been implemented and integrated in UNIPred-Web, a web application that we recently introduced for combining the UNIPred algorithm, able to address both integration and protein function prediction in an imbalance-aware fashion, with an easy to use vertex-centric exploration of the integrated network. The tool has been deeply amended from different standpoints, including the prediction core algorithm. Several tests on networks of different size and connectivity have been conducted to show off the vast potential of our methodology; moreover, enrichment analyses have been performed to assess the biological meaningfulness of detected communities. Finally, a CoV-human network has been embedded in the system, and a corresponding case study presented, including the visualization and the prediction of human host proteins that potentially interact with SARS-CoV2 proteins.

Interweaving Tumor Heterogeneity into the Cancer Epigenetic/Metabolic Axis

Significance: The epigenomic/metabolic landscape in cancer has been studied extensively in the past decade and forms the basis of various drug targets. Yet, cancer treatment remains a challenge, with clinical trials exhibiting limited efficacy and high relapse rates. Patients respond differently to therapy, which is fundamentally attributed to tumor heterogeneity, both across and within tumors. This review focuses on the interactions between the heterogeneous tumor microenvironment (TME) and the epigenomic/metabolic axis in cancer, as well as the emerging technologies under development to aid heterogeneity studies. Recent Advances: Interlinks between epigenetics and metabolism in cancer have been reported. Emerging studies have unveiled interactions between the TME and cancer cells that play a critical role in regulating epigenetics and reprogramming cancer metabolism, suggesting a three-way cross talk. Critical Issues: This cross talk accentuates the multiplex nature of cancer, and the importance of considering tumor heterogeneity in various epigenomic/metabolic cancer studies. Future Directions: With the advancement in single-cell profiling, it may be possible to identify cancer subclones and their unique vulnerabilities to develop a multimodal therapy. Drugs targeting the TME are currently being studied, and a better understanding of the TME in regulating cancer epigenetics and metabolism may hold the key to identifying novel therapeutic targets.

Argonaut: A Web Platform for Collaborative Multi-omic Data Visualization and Exploration

Researchers now generate large multi-omic datasets using increasingly mature mass spectrometry techniques at an astounding pace, facing new challenges of "Big Data" dissemination, visualization, and exploration. Conveniently, web-based data portals accommodate the complexity of multi-omic experiments and the many experts involved. However, developing these tailored companion resources requires programming expertise and knowledge of web server architecture-a substantial burden for most. Here, we describe Argonaut, a simple, code-free, and user-friendly platform for creating customizable, interactive data-hosting websites. Argonaut carries out real-time statistical analyses of the data, which it organizes into easily sharable projects. Collaborating researchers worldwide can explore the results, visualized through popular plots, and modify them to streamline data interpretation. Increasing the pace and ease of access to multi-omic data, Argonaut aims to propel discovery of new biological insights. We showcase the capabilities of this tool using a published multi-omics dataset on the large mitochondrial protease deletion collection.

Crippling life support for SARS-CoV-2 and other viruses through synthetic lethality

With the rapid global spread of SARS-CoV-2, we have become acutely aware of the inadequacies of our ability to respond to viral epidemics. Although disrupting the viral life cycle is critical for limiting viral spread and disease, it has proven challenging to develop targeted and selective therapeutics. Synthetic lethality offers a promising but largely unexploited strategy against infectious viral disease; as viruses infect cells, they abnormally alter the cell state, unwittingly exposing new vulnerabilities in the infected cell. Therefore, we propose that effective therapies can be developed to selectively target the virally reconfigured host cell networks that accompany altered cellular states to cripple the host cell that has been converted into a virus factory, thus disrupting the viral life cycle.

Precision medicine - networks to the rescue

Genetic variants are often not predictive of the phenotypic outcome. Individuals carrying the same pathogenic variant, associated with Mendelian or complex disease, can manifest to different extents, from severe-to-mild to no disease. Improving the accuracy of predicted clinical manifestations of genetic variants has emerged as one of the biggest challenges in precision medicine, which can only be addressed by understanding the mechanisms underlying genotype-phenotype relationships. Efforts to understand the molecular basis of these relationships have identified complex systems of interacting biomolecules that underlie cellular function. Here, we review recent advances in how modeling cellular systems as networks of interacting proteins has fueled identification of disease-associated processes, delineation of underlying molecular mechanisms, and prediction of the pathogenicity of variants. This review is intended to be inspiring for clinicians, geneticists, and network biologists alike who aim to jointly advance our understanding of human disease and accelerate progress toward precision medicine.

HaloTag-based conjugation of proteins to barcoding-oligonucleotides

Highly sensitive protein quantification enables the detection of a small number of protein molecules that serve as markers/triggers for various biological phenomena, such as cancer. Here, we describe the development of a highly sensitive protein quantification system called HaloTag protein barcoding. The method involves covalent linking of a target protein to a unique molecule counting oligonucleotide at a 1:1 conjugation ratio based on an azido-cycloalkyne click reaction. The sensitivity of the HaloTag-based barcoding was remarkably higher than that of a conventional luciferase assay. The HaloTag system was successfully validated by analyzing a set of protein-protein interactions, with the identification rate of 44% protein interactions between positive reference pairs reported in the literature. Desmoglein 3, the target antigen of pemphigus vulgaris, an IgG-mediated autoimmune blistering disease, was used in a HaloTag protein barcode assay to detect the anti-DSG3 antibody. The dynamic range of the assay was over 104-times wider than that of a conventional enzyme-linked immunosorbent assay (ELISA). The technology was used to detect anti-DSG3 antibody in patient samples with much higher sensitivity compared to conventional ELISA. Our detection system, with its superior sensitivity, enables earlier detection of diseases possibly allowing the initiation of care/treatment at an early disease stage.

Genetics' Piece of the PI: Inferring the Origin of Complex Traits and Diseases from Proteome-Wide Protein-Protein Interaction Dynamics

How do common and rare genetic polymorphisms contribute to quantitative traits or disease risk and progression? Multiple human traits have been extensively characterized at the genomic level, revealing their complex genetic architecture. However, it is difficult to resolve the mechanisms by which specific variants contribute to a phenotype. Recently, analyses of variant effects on molecular traits have uncovered intermediate mechanisms that link sequence variation to phenotypic changes. Yet, these methods only capture a fraction of genetic contributions to phenotype. Here, in reviewing the field, it is proposed that complex traits can be understood by characterizing the dynamics of biochemical networks within living cells, and that the effects of genetic variation can be captured on these networks by using protein-protein interaction (PPI) methodologies. This synergy between PPI methodologies and the genetics of complex traits opens new avenues to investigate the molecular etiology of human diseases and to facilitate their prevention or treatment.

Differential network analysis depicts regulatory mechanisms for hepatocellular carcinoma from diverse backgrounds

Aim: To elucidate the integrative combinational gene regulatory network landscape of hepatocellular carcinoma (HCC) molecular carcinogenesis from diverse background. Materials & methods: Modified gene regulatory network analysis was used to prioritize differentially regulated genes and links. Integrative comparisons using bioinformatics methods were applied to identify potential critical molecules and pathways in HCC with different backgrounds. Results: E2F1 with its surrounding regulatory links were identified to play different key roles in the HCC risk factor dysregulation mechanisms. Hsa-mir-19a was identified as showed different effects in the three HCC differential regulation networks, and showed vital regulatory role in HBV-related HCC. Conclusion: We describe in detail the regulatory networks involved in HCC with different backgrounds. E2F1 may serve as a universal target for HCC treatment.

Mapping genetic interactions in cancer: a road to rational combination therapies

The discovery of synthetic lethal interactions between poly (ADP-ribose) polymerase (PARP) inhibitors and BRCA genes, which are involved in homologous recombination, led to the approval of PARP inhibition as a monotherapy for patients with BRCA1/2-mutated breast or ovarian cancer. Studies following the initial observation of synthetic lethality demonstrated that the reach of PARP inhibitors is well beyond just BRCA1/2 mutants. Insights into the mechanisms of action of anticancer drugs are fundamental for the development of targeted monotherapies or rational combination treatments that will synergize to promote cancer cell death and overcome mechanisms of resistance. The development of targeted therapeutic agents is premised on mapping the physical and functional dependencies of mutated genes in cancer. An important part of this effort is the systematic screening of genetic interactions in a variety of cancer types. Until recently, genetic-interaction screens have relied either on the pairwise perturbations of two genes or on the perturbation of genes of interest combined with inhibition by commonly used anticancer drugs. Here, we summarize recent advances in mapping genetic interactions using targeted, genome-wide, and high-throughput genetic screens, and we discuss the therapeutic insights obtained through such screens. We further focus on factors that should be considered in order to develop a robust analysis pipeline. Finally, we discuss the integration of functional interaction data with orthogonal methods and suggest that such approaches will increase the reach of genetic-interaction screens for the development of rational combination therapies.

Identification of hyper-rewired genomic stress non-oncogene addiction genes across 15 cancer types

Non-oncogene addiction (NOA) genes are essential for supporting the stress-burdened phenotype of tumours and thus vital for their survival. Although NOA genes are acknowledged to be potential drug targets, there has been no large-scale attempt to identify and characterise them as a group across cancer types. Here we provide the first method for the identification of conditional NOA genes and their rewired neighbours using a systems approach. Using copy number data and expression profiles from The Cancer Genome Atlas (TCGA) we performed comparative analyses between high and low genomic stress tumours for 15 cancer types. We identified 101 condition-specific differential coexpression modules, mapped to a high-confidence human interactome, comprising 133 candidate NOA rewiring hub genes. We observe that most modules lose coexpression in the high-stress state and that activated stress modules and hubs take part in homoeostasis maintenance processes such as chromosome segregation, oxireductase activity, mitotic checkpoint (PLK1 signalling), DNA replication initiation and synaptic signalling. We furthermore show that candidate NOA rewiring hubs are unique for each cancer type, but that their respective rewired neighbour genes largely are shared across cancer types.

Mapping the protein-protein and genetic interactions of cancer to guide precision medicine

Massive efforts to sequence cancer genomes have compiled an impressive catalogue of cancer mutations, revealing the recurrent exploitation of a handful of 'hallmark cancer pathways'. However, unraveling how sets of mutated proteins in these and other pathways hijack pro-proliferative signaling networks and dictate therapeutic responsiveness remains challenging. Here, we show that cancer driver protein-protein interactions are enriched for additional cancer drivers, highlighting the power of physical interaction maps to explain known, as well as uncover new, disease-promoting pathway interrelationships. We hypothesize that by systematically mapping the protein-protein and genetic interactions in cancer-thereby creating Cancer Cell Maps-we will create resources against which to contextualize a patient's mutations into perturbed pathways/complexes and thereby specify a matching targeted therapeutic cocktail.

Mapping phospho-catalytic dependencies of therapy-resistant tumours reveals actionable vulnerabilities

Phosphorylation networks intimately regulate mechanisms of response to therapies. Mapping the phospho-catalytic profile of kinases in cells or tissues remains a challenge. Here, we introduce a practical high-throughput system to measure the enzymatic activity of kinases using biological peptide targets as phospho-sensors to reveal kinase dependencies in tumour biopsies and cell lines. A 228-peptide screen was developed to detect the activity of >60 kinases, including ABLs, AKTs, CDKs and MAPKs. Focusing on BRAF^V600E tumours, we found mechanisms of intrinsic resistance to BRAF^V600E-targeted therapy in colorectal cancer, including targetable parallel activation of PDPK1 and PRKCA. Furthermore, mapping the phospho-catalytic signatures of melanoma specimens identifies RPS6KB1 and PIM1 as emerging druggable vulnerabilities predictive of poor outcome in BRAF^V600E patients. The results show that therapeutic resistance can be caused by the concerted upregulation of interdependent pathways. Our kinase activity-mapping system is a versatile strategy that innovates the exploration of actionable kinases for precision medicine.

Integrating clinical, molecular, proteomic and histopathological data within the tissue context: tissunomics

Malignant tumours show a marked degree of morphological, molecular and proteomic heterogeneity. This variability is closely related to microenvironmental factors and the location of the tumour. The activation of genetic alterations is very tissue‐dependent and only few tumours have distinct genetic alterations. Importantly, the activation state of proteins and signaling factors is heterogeneous in the primary tumour and in metastases and recurrences. The molecular diagnosis based only on genetic alterations can lead to treatments with unpredictable responses, depending on the tumour location, such as the tumour response in melanomas versus colon carcinomas with BRAF mutations. Therefore, we understand that the correct evaluation of tumours requires a system that integrates both morphological, molecular and protein information in a clinical and pathological context, where intratumoral heterogeneity can be assessed. Thus, we propose the term ‘tissunomics’, where the diagnosis will be contextualised in each tumour based on the complementation of the pathological, molecular, protein expression, environmental cells and clinical data.

Identifying Driver Interfaces Enriched for Somatic Missense Mutations in Tumors

Human cancers often harbor large numbers of somatic mutations. However, only a small proportion of these mutations are expected to contribute to tumor growth and progression. Therefore, determining causal driver mutations and the genes they target is becoming an important challenge in cancer genomics. Here we describe an approach for mapping somatic mutations onto 3D structures of human proteins in complex to identify "driver interfaces." Our strategy relies on identifying protein-interaction interfaces that are unexpectedly biased toward nonsynonymous mutations, which suggests that these interfaces are subject to positive selection during tumorigenesis, implicating the interacting proteins as candidate drivers.

Integrating molecular networks with genetic variant interpretation for precision medicine

More reliable and cheaper sequencing technologies have revealed the vast mutational landscapes characteristic of many phenotypes. The analysis of such genetic variants has led to successful identification of altered proteins underlying many Mendelian disorders. Nevertheless the simple one‐variant one‐phenotype model valid for many monogenic diseases does not capture the complexity of polygenic traits and disorders. Although experimental and computational approaches have improved detection of functionally deleterious variants and important interactions between gene products, the development of comprehensive models relating genotype and phenotypes remains a challenge in the field of genomic medicine. In this context, a new view of the pathologic state as significant perturbation of the network of interactions between biomolecules is crucial for the identification of biochemical pathways associated with complex phenotypes. Seminal studies in systems biology combined the analysis of genetic variation with protein–protein interaction networks to demonstrate that even as biological systems evolve to be robust to genetic variation, their topologies create disease vulnerabilities. More recent analyses model the impact of genetic variants as changes to the “wiring” of the interactome to better capture heterogeneity in genotype–phenotype relationships. These studies lay the foundation for using networks to predict variant effects at scale using machine‐learning or algorithmic approaches. A wealth of databases and resources for the annotation of genotype–phenotype relationships have been developed to support developments in this area. This overview describes how study of the molecular interactome has generated insights linking the organization of biological systems to disease mechanism, and how this information can enable precision medicine.

This article is categorized under:

Translational, Genomic, and Systems Medicine > Translational Medicine

Biological Mechanisms > Cell Signaling

Models of Systems Properties and Processes > Mechanistic Models

Analytical and Computational Methods > Computational Methods

Signalling maps in cancer research: construction and data analysis

Generation and usage of high-quality molecular signalling network maps can be augmented by standardizing notations, establishing curation workflows and application of computational biology methods to exploit the knowledge contained in the maps. In this manuscript, we summarize the major aims and challenges of assembling information in the form of comprehensive maps of molecular interactions. Mainly, we share our experience gained while creating the Atlas of Cancer Signalling Network. In the step-by-step procedure, we describe the map construction process and suggest solutions for map complexity management by introducing a hierarchical modular map structure. In addition, we describe the NaviCell platform, a computational technology using Google Maps API to explore comprehensive molecular maps similar to geographical maps and explain the advantages of semantic zooming principles for map navigation. We also provide the outline to prepare signalling network maps for navigation using the NaviCell platform. Finally, several examples of cancer high-throughput data analysis and visualization in the context of comprehensive signalling maps are presented.

Mutation pattern analysis reveals polygenic mini-drivers associated with relapse after surgery in lung adenocarcinoma

The genomic lesions found in malignant tumours exhibit a striking degree of heterogeneity. Many tumours lack a known driver mutation, and their genetic basis is unclear. By mapping the somatic mutations identified in primary lung adenocarcinomas onto an independent coexpression network derived from normal tissue, we identify a critical gene network enriched for metastasis-associated genes. While individual genes within this module were rarely mutated, a significant accumulation of mutations within this geneset was predictive of relapse in lung cancer patients that have undergone surgery. Since it is the density of mutations within this module that is informative, rather than the status of any individual gene, these data are in keeping with a 'mini-driver' model of tumorigenesis in which multiple mutations, each with a weak effect, combine to form a polygenic driver with sufficient power to significantly alter cell behaviour and ultimately patient outcome. These polygenic mini-drivers therefore provide a means by which heterogeneous mutation patterns can generate the consistent hallmark changes in phenotype observed across tumours.

MicroRNA-1225-5p behaves as a tumor suppressor in human glioblastoma via targeting of IRS1

Background

MicroRNAs (miRNAs) play an important role in cancer initiation, progression, and metastasis by directly regulating their target genes.

Materials and methods

In this study, we observed that the miR-1225-5p expression level in glioblastoma tissues was significantly lower as compared with that in normal brain tissues, and its low expression was significantly associated with histopathological grade and poor patient prognosis.

Results

Through establishing a miR-1225-5p overexpression glioblastoma cell line, we found that ectopic overexpression of miR-1225-5p inhibited the proliferation, migration, and invasion of glioblastoma cells in vitro. Moreover, the growth of a glioblastoma xenograft tumor was attenuated by overexpression of miR-1225-5p. Further integrative studies suggested that the insulin receptor substrate 1 (IRS1) was a direct functional target of miR-1225-5p in glioblastoma, and the mRNA and protein levels of IRS1 in six human glioblastoma cell lines (A172, SW1783, U87, LN-229, SW1088, and T98G) were significantly higher as compared with normal human astrocytes.

Conclusion

These results suggest that miR-1225-5p may be a novel candidate for glioblastoma therapy.

Synthetic Lethal Networks for Precision Oncology: Promises and Pitfalls

Synthetic lethal interactions, in which the simultaneous loss of function of two genes produces a lethal phenotype, are being explored as a means to therapeutically exploit cancer-specific vulnerabilities and expand the scope of precision oncology. Currently, three Food and Drug Administration-approved drugs work by targeting the synthetic lethal interaction between BRCA1/2 and PARP. This review examines additional efforts to discover networks of synthetic lethal interactions and discusses both challenges and opportunities regarding the translation of new synthetic lethal interactions into the clinic.

Primary brain tumours in adults

Primary CNS tumours refer to a heterogeneous group of tumours arising from cells within the CNS, and can be benign or malignant. Malignant primary brain tumours remain among the most difficult cancers to treat, with a 5 year overall survival no greater than 35%. The most common malignant primary brain tumours in adults are gliomas. Recent advances in molecular biology have improved understanding of glioma pathogenesis, and several clinically significant genetic alterations have been described. A number of these (IDH, 1p/19q codeletion, H3 Lys27Met, and RELA-fusion) are now combined with histology in the revised 2016 WHO classification of CNS tumours. It is likely that understanding such molecular alterations will contribute to the diagnosis, grading, and treatment of brain tumours. This progress in genomics, along with significant advances in cancer and CNS immunology, has defined a new era in neuro-oncology and holds promise for diagntic and therapeutic improvement. The challenge at present is to translate these advances into effective treatments. Current efforts are focused on developing molecular targeted therapies, immunotherapies, gene therapies, and novel drug-delivery technologies. Results with single-agent therapies have been disappointing so far, and combination therapies seem to be required to achieve a broad and durable antitumour response. Biomarker-targeted clinical trials could improve efficiencies of therapeutic development.

The Psychiatric Cell Map Initiative: A Convergent Systems Biological Approach to Illuminating Key Molecular Pathways in Neuropsychiatric Disorders

Although gene discovery in neuropsychiatric disorders, including autism spectrum disorder, intellectual disability, epilepsy, schizophrenia, and Tourette disorder, has accelerated, resulting in a large number of molecular clues, it has proven difficult to generate specific hypotheses without the corresponding datasets at the protein complex and functional pathway level. Here, we describe one path forward-an initiative aimed at mapping the physical and genetic interaction networks of these conditions and then using these maps to connect the genomic data to neurobiology and, ultimately, the clinic. These efforts will include a team of geneticists, structural biologists, neurobiologists, systems biologists, and clinicians, leveraging a wide array of experimental approaches and creating a collaborative infrastructure necessary for long-term investigation. This initiative will ultimately intersect with parallel studies that focus on other diseases, as there is a significant overlap with genes implicated in cancer, infectious disease, and congenital heart defects.

Cutaneous Papillomaviruses and Non-melanoma Skin Cancer: Causal Agents or Innocent Bystanders?

There is still controversy in the scientific field about whether certain types of cutaneous human papillomaviruses (HPVs) are causally involved in the development of non-melanoma skin cancer (NMSC). Deciphering the etiological role of cutaneous HPVs requires - besides tissue culture systems - appropriate preclinical models to match the obtained results with clinical data from affected patients. Clear scientific evidence about the etiology and underlying mechanisms involved in NMSC development is fundamental to provide reasonable arguments for public health institutions to classify at least certain cutaneous HPVs as group 1 carcinogens. This in turn would have implications on fundraising institutions and health care decision makers to force - similarly as for anogenital cancer - the implementation of a broad vaccination program against "high-risk" cutaneous HPVs to prevent NMSC as the most frequent cancer worldwide. Precise knowledge of the multi-step progression from normal cells to cancer is a prerequisite to understand the functional and clinical impact of cofactors that affect the individual outcome and the personalized treatment of a disease. This overview summarizes not only recent arguments that favor the acceptance of a viral etiology in NMSC development but also reflects aspects of causality in medicine, the use of empirically meaningful model systems and strategies for prevention.

The extracellular matrix as a key regulator of intracellular signalling networks

The extracellular matrix (ECM) is a salient feature of all solid tissues within the body. This complex, acellular entity is composed of hundreds of individual molecules whose assembly, architecture and biomechanical properties are critical to controlling the behaviour and phenotype of the different cell types residing within tissues. Cells are the basic unit of life and the core building block of tissues and organs. At their simplest, they follow a set of rules, governed by their genetic code and effected through the complex protein signalling networks that these genes encode. These signalling networks assimilate and process the information received by the cell to control cellular decisions that govern cell fate. The ECM is the biggest provider of external stimuli to cells and as such is responsible for influencing intracellular signalling dynamics. In this review, we discuss the inclusion of ECM as a central regulatory signalling sub-network in computational models of cellular decision making, with a focus on its role in diseases such as cancer. LINKED ARTICLES: This article is part of a themed section on Translating the Matrix. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v176.1/issuetoc.

Analysis of variability in high throughput screening data: applications to melanoma cell lines and drug responses

High-throughput screening (HTS) strategies and protocols have undergone significant development in the last decade. It is now possible to screen hundreds of thousands of compounds, each exploring multiple biological phenotypes and parameters, against various cell lines or model systems in a single setting. However, given the vast amount of data such studies generate, the fact that they use multiple reagents, and are often technician-intensive, questions have been raised about the variability, reliability and reproducibility of HTS results. Assessments of the impact of the multiple factors in HTS studies could arguably lead to more compelling insights into the robustness of the results of a particular screen, as well as the overall quality of the study. We leveraged classical, yet highly flexible, analysis of variance (ANOVA)-based linear models to explore how different factors contribute to the variation observed in a screening study of four different melanoma cell lines and 120 drugs over nine dosages studied in two independent academic laboratories. We find that factors such as plate effects, appropriate dosing ranges, and to a lesser extent, the laboratory performing the screen, are significant predictors of variation in drug responses across the cell lines. Further, we show that when sources of variation are quantified and controlled for, they contextualize claims of inconsistencies and reveal the overall quality of the HTS studies performed at each participating laboratory. In the context of the broader screening study, we show that our analysis can also elucidate the robust effects of drugs, even those within specific cell lines.

Opportunities and Challenges in Building a Spatiotemporal Multi-scale Model of the Human Pancreatic β Cell

The construction of a predictive model of an entire eukaryotic cell that describes its dynamic structure from atomic to cellular scales is a grand challenge at the intersection of biology, chemistry, physics, and computer science. Having such a model will open new dimensions in biological research and accelerate healthcare advancements. Developing the necessary experimental and modeling methods presents abundant opportunities for a community effort to realize this goal. Here, we present a vision for creation of a spatiotemporal multi-scale model of the pancreatic β-cell, a relevant target for understanding and modulating the pathogenesis of diabetes.

Looking beyond drivers and passengers in cancer genome sequencing data

Cancer arises as a result of acquired changes in the DNA sequence of the genome of somatic cells. A subset of the genetic changes, dubbed driver mutations, propels tumor growth, and the remaining changes are passengers, apparently inconsequential for neoplastic transformation. Massive genome sequencing of thousands of tumors from all major cancer types has enabled cataloging of the so-called driver and passenger mutations, and facilitated molecular classification of cancer, guiding precision medicine approach for the patients. Nonetheless, innovative analyses of cancer genomics data has led to novel, sometimes serendipitous findings that have aided to our understanding of other aspects of the biology of the disease and opened up new frontiers. For instance, emerging findings show that mutational patterns in cancer genomes can help detect signatures of known and novel DNA damage and repair processes, provide a likely chronological account of genomic changes in cancer genomes, and allow revisiting the models of cancer evolution. These findings have stimulated original approaches to identify disease etiology, stratify patients, target the disease, and monitor patient responses, complementing driver-mutation centric approaches. In this review, we discuss these emerging approaches and unexpected breakthroughs, and their implications for basic cancer research and clinical practices.

Commonly Occurring Cell Subsets in High-Grade Serous Ovarian Tumors Identified by Single-Cell Mass Cytometry

We have performed an in-depth single-cell phenotypic characterization of high-grade serous ovarian cancer (HGSOC) by multiparametric mass cytometry (CyTOF). Using a CyTOF antibody panel to interrogate features of HGSOC biology, combined with unsupervised computational analysis, we identified noteworthy cell types co-occurring across the tumors. In addition to a dominant cell subset, each tumor harbored rarer cell phenotypes. One such group co-expressed E-cadherin and vimentin (EV), suggesting their potential role in epithelial mesenchymal transition, which was substantiated by pairwise correlation analyses. Furthermore, tumors from patients with poorer outcome had an increased frequency of another rare cell type that co-expressed vimentin, HE4, and cMyc. These poorer-outcome tumors also populated more cell phenotypes, as quantified by Simpson's diversity index. Thus, despite the recognized genomic complexity of the disease, the specific cell phenotypes uncovered here offer a focus for therapeutic intervention and disease monitoring.

Protein Interactomics by Two-Hybrid Methods

Comprehensive identification of direct, physical interactions between biological macromolecules, such as protein-protein, protein-DNA, and protein-RNA interactions, is critical for our understanding of the function of gene products as well as the global organization and interworkings of various molecular machines within the cell. The accurate and comprehensive detection of direct interactions, however, remains a huge challenge due to the inherent structural complexity arising from various post-transcriptional and translational modifications coupled with huge heterogeneity in concentration, affinity, and subcellular location differences existing for any interacting molecules. This has created a need for developing multiple orthogonal and complementary assays for detecting various types of biological interactions. In this introduction, we discuss the methods developed for measuring different types of molecular interactions with an emphasis on direct protein-protein interactions, critical issues for generating high-quality interactome datasets, and the insights into biological networks and human diseases that current interaction mapping efforts provide. Further, we will discuss what future might lie ahead for the continued evolution of two-hybrid methods and the role of interactomics for expanding the advancement of biomedical science.

Computational Systems Biology Approach for the Study of Rheumatoid Arthritis: From a Molecular Map to a Dynamical Model

In this work we present a systematic effort to summarize current biological pathway knowledge concerning Rheumatoid Arthritis (RA). We are constructing a detailed molecular map based on exhaustive literature scanning, strict curation criteria, re-evaluation of previously published attempts and most importantly experts' advice. The RA map will be web-published in the coming months in the form of an interactive map, using the MINERVA platform, allowing for easy access, navigation and search of all molecular pathways implicated in RA, serving thus, as an on line knowledgebase for the disease. Moreover the map could be used as a template for Omics data visualization offering a first insight about the pathways affected in different experimental datasets. The second goal of the project is a dynamical study focused on synovial fibroblasts' behavior under different initial conditions specific to RA, as recent studies have shown that synovial fibroblasts play a crucial role in driving the persistent, destructive characteristics of the disease. Leaning on the RA knowledgebase and using the web platform Cell Collective, we are currently building a Boolean large scale dynamical model for the study of RA fibroblasts' activation.

NDEx 2.0: A Clearinghouse for Research on Cancer Pathways

We present NDEx 2.0, the latest release of the Network Data Exchange (NDEx) online data commons (www.ndexbio.org) and the ways in which it can be used to (i) improve the quality and abundance of biological networks relevant to the cancer research community; (ii) provide a medium for collaboration involving networks; and (iii) facilitate the review and dissemination of networks. We describe innovations addressing the challenges of an online data commons: scalability, data integration, data standardization, control of content and format by authors, and decentralized mechanisms for review. The practical use of NDEx is presented in the context of a novel strategy to foster network-oriented communities of interest in cancer research by adapting methods from academic publishing and social media. Cancer Res; 77(21); e58-61. ©2017 AACR.

Incorporating networks in a probabilistic graphical model to find drivers for complex human diseases

Discovering genetic mechanisms driving complex diseases is a hard problem. Existing methods often lack power to identify the set of responsible genes. Protein-protein interaction networks have been shown to boost power when detecting gene-disease associations. We introduce a Bayesian framework, Conflux, to find disease associated genes from exome sequencing data using networks as a prior. There are two main advantages to using networks within a probabilistic graphical model. First, networks are noisy and incomplete, a substantial impediment to gene discovery. Incorporating networks into the structure of a probabilistic models for gene inference has less impact on the solution than relying on the noisy network structure directly. Second, using a Bayesian framework we can keep track of the uncertainty of each gene being associated with the phenotype rather than returning a fixed list of genes. We first show that using networks clearly improves gene detection compared to individual gene testing. We then show consistently improved performance of Conflux compared to the state-of-the-art diffusion network-based method Hotnet2 and a variety of other network and variant aggregation methods, using randomly generated and literature-reported gene sets. We test Hotnet2 and Conflux on several network configurations to reveal biases and patterns of false positives and false negatives in each case. Our experiments show that our novel Bayesian framework Conflux incorporates many of the advantages of the current state-of-the-art methods, while offering more flexibility and improved power in many gene-disease association scenarios.

Networks and pathway-based methods are commonly used to improve the power of gene detection in associations with complex human diseases. Network diffusion approaches have shown their effectiveness and superior performance in cancer studies. Still, there are many problems such as noise and missingness with currently available human networks that bias the results of gene detection. We propose a novel graphical model-based method Conflux that overcomes several of the pitfalls of the existing state-of-the-art approaches while building on their successes. Conflux integrates genotype data with networks directly, using diffusion-like methods, but only as part of a structure in a probabilistic model to reduce the negative effect of the noise in the networks. This Bayesian framework allows Conflux to keep track of the uncertainty in the gene list that is being associated with the disease and consequently rank the genes with respect to our confidence in the association. It also allows for the discovery of gene sets that are not fully supported by the network if they have enough support in the data. These improvements result in a flexible approach that improves the power in many gene-disease association scenarios while reducing the number of false positives reported.