The Many Faces of Gene Regulation in Cancer: A Computational Oncogenomics Outlook

Gene co-expression networks reveal sex-biased differences in musculoskeletal ageing

Aging is a universal and progressive process involving the deterioration of physiological functions and the accumulation of cellular damage. Gene regulation programs influence how phenotypes respond to environmental and intrinsic changes during aging. Although several factors, including sex, are known to impact this process, the underlying mechanisms remain incompletely understood. Here, we investigate the functional organization patterns of skeletal muscle genes across different sexes and ages using gene co-expression networks (GCNs) to explore their influence on aging. We constructed GCNs for three different age groups for male and female samples, analyzed topological similarities and differences, inferred significant associated processes for each network, and constructed null models to provide statistically robust results. We found that each network is topologically and functionally distinct, with young women having the most associated processes, likely due to reproductive tasks. The functional organization and modularity of genes decline with age, starting from middle age, potentially leading to age-related deterioration. Women maintain better gene functional organization throughout life compared to men, especially in processes like macroautophagy and sarcomere organization. The study suggests that the loss of gene co-expression could be a universal aging marker. This research offers insights into how gene organization changes with age and sex, providing a complementary method to analyze aging.

Multi-OMICS approaches in cancer biology: New era in cancer therapy

Innovative multi-omics frameworks integrate diverse datasets from the same patients to enhance our understanding of the molecular and clinical aspects of cancers. Advanced omics and multi-view clustering algorithms present unprecedented opportunities for classifying cancers into subtypes, refining survival predictions and treatment outcomes, and unravelling key pathophysiological processes across various molecular layers. However, with the increasing availability of cost-effective high-throughput technologies (HTT) that generate vast amounts of data, analyzing single layers often falls short of establishing causal relations. Integrating multi-omics data spanning genomes, epigenomes, transcriptomes, proteomes, metabolomes, and microbiomes offers unique prospects to comprehend the underlying biology of complex diseases like cancer. This discussion explores algorithmic frameworks designed to uncover cancer subtypes, disease mechanisms, and methods for identifying pivotal genomic alterations. It also underscores the significance of multi-omics in tumor classifications, diagnostics, and prognostications. Despite its unparalleled advantages, the integration of multi-omics data has been slow to find its way into everyday clinics. A major hurdle is the uneven maturity of different omics approaches and the widening gap between the generation of large datasets and the capacity to process this data. Initiatives promoting the standardization of sample processing and analytical pipelines, as well as multidisciplinary training for experts in data analysis and interpretation, are crucial for translating theoretical findings into practical applications.

Oncogenes and tumor suppressor genes: functions and roles in cancers

Cancer, being the most formidable ailment, has had a profound impact on the human health. The disease is primarily associated with genetic mutations that impact oncogenes and tumor suppressor genes (TSGs). Recently, growing evidence have shown that X-linked TSGs have specific role in cancer progression and metastasis as well. Interestingly, our genome harbors around substantial portion of genes that function as tumor suppressors, and the X chromosome alone harbors a considerable number of TSGs. The scenario becomes even more compelling as X-linked TSGs are adaptive to key epigenetic processes such as X chromosome inactivation. Therefore, delineating the new paradigm related to X-linked TSGs, for instance, their crosstalk with autosome and involvement in cancer initiation, progression, and metastasis becomes utmost importance. Considering this, herein, we present a comprehensive discussion of X-linked TSG dysregulation in various cancers as a consequence of genetic variations and epigenetic alterations. In addition, the dynamic role of X-linked TSGs in sex chromosome-autosome crosstalk in cancer genome remodeling is being explored thoroughly. Besides, the functional roles of ncRNAs, role of X-linked TSG in immunomodulation and in gender-based cancer disparities has also been highlighted. Overall, the focal idea of the present article is to recapitulate the findings on X-linked TSG regulation in the cancer landscape and to redefine their role toward improving cancer treatment strategies.

Analysis of VEGF, IGF1/2 and the Long Noncoding RNA (lncRNA) H19 Expression in Polish Women with Endometriosis

The coordinated action of VEGF, IGF1/2 and H19 factors influences the development of endometriosis. The aim of this study was to analyze the expression level of these genes in patients with endometriosis. The study group consisted of 100 patients who were diagnosed with endometriosis on laparoscopic and pathological examination. The control group consisted of 100 patients who were found to be free of endometriosis during the surgical procedure and whose eutopic endometrium wasnormal on histopathological examination. These patients were operated on for uterine fibroids. Gene expression was determined by RT-PCR. The expression of the VEGF gene was significantly higher in the samples classified as clinical stage 1-2 compared to the control material (p < 0.05). There was also a statistically significant difference between the samples studied at clinical stages 1-2 and 3-4 (p < 0.01). The expression of the VEGF gene in the group classified as 1-2 was significantly higher. IGF1 gene expression was significantly lower both in the group of samples classified as clinical stages 1-2 and 3-4 compared to the control group (p < 0.05 in both cases). The expression of the H19 gene was significantly lower in the group of samples classified as clinical stage 3-4 compared to the control group (p < 0.01). The reported studies suggest significant roles of VEGF, IGF and H19 expression in the pathogenesis of endometriosis.

Onco-Breastomics: An Eco-Evo-Devo Holistic Approach

Known as a diverse collection of neoplastic diseases, breast cancer (BC) can be hyperbolically characterized as a dynamic pseudo-organ, a living organism able to build a complex, open, hierarchically organized, self-sustainable, and self-renewable tumor system, a population, a species, a local community, a biocenosis, or an evolving dynamical ecosystem (i.e., immune or metabolic ecosystem) that emphasizes both developmental continuity and spatio-temporal change. Moreover, a cancer cell community, also known as an oncobiota, has been described as non-sexually reproducing species, as well as a migratory or invasive species that expresses intelligent behavior, or an endangered or parasite species that fights to survive, to optimize its features inside the host's ecosystem, or that is able to exploit or to disrupt its host circadian cycle for improving the own proliferation and spreading. BC tumorigenesis has also been compared with the early embryo and placenta development that may suggest new strategies for research and therapy. Furthermore, BC has also been characterized as an environmental disease or as an ecological disorder. Many mechanisms of cancer progression have been explained by principles of ecology, developmental biology, and evolutionary paradigms. Many authors have discussed ecological, developmental, and evolutionary strategies for more successful anti-cancer therapies, or for understanding the ecological, developmental, and evolutionary bases of BC exploitable vulnerabilities. Herein, we used the integrated framework of three well known ecological theories: the Bronfenbrenner's theory of human development, the Vannote's River Continuum Concept (RCC), and the Ecological Evolutionary Developmental Biology (Eco-Evo-Devo) theory, to explain and understand several eco-evo-devo-based principles that govern BC progression. Multi-omics fields, taken together as onco-breastomics, offer better opportunities to integrate, analyze, and interpret large amounts of complex heterogeneous data, such as various and big-omics data obtained by multiple investigative modalities, for understanding the eco-evo-devo-based principles that drive BC progression and treatment. These integrative eco-evo-devo theories can help clinicians better diagnose and treat BC, for example, by using non-invasive biomarkers in liquid-biopsies that have emerged from integrated omics-based data that accurately reflect the biomolecular landscape of the primary tumor in order to avoid mutilating preventive surgery, like bilateral mastectomy. From the perspective of preventive, personalized, and participatory medicine, these hypotheses may help patients to think about this disease as a process governed by natural rules, to understand the possible causes of the disease, and to gain control on their own health.

Recent progress and the emerging role of lncRNAs in cancer drug resistance; focusing on signaling pathways

It is becoming more and more apparent that many of the genetic alterations associated with cancer are located in areas that do not encode proteins. lncRNAs are a class of RNAs that do not code for proteins but play a crucial role in maintaining cell function and regulating various cellular processes. By doing this, they have recently introduced what may be a brand-new and essential layer of biological control. These have more than 200 nucleotides and are linked to several diseases; as a result, they have become potential tools for therapeutic intervention. Emerging technologies suggest the presence of mutations on genomic loci that give rise to lncRNAs rather than proteins in a disease as complex as cancer. These lncRNAs play essential parts in gene regulation, which impacts several cellular homeostasis processes, including proliferation, survival, migration, and genomic stability. The leading cause of death in the world today is cancer. Delays in diagnosis and a lack of standard and efficient treatments are the leading causes of the high death rate. Clinically, surgery is frequently used successfully to remove cancers that have not spread, but it is less successful in treating metastatic cancer, which has a drastically lower chance of survival. Chemotherapeutic drugs are a typical therapy to treat the cancer that has spread to other organs. Drug resistance to chemotherapy, however, presents a significant challenge to achieving positive outcomes and is frequently the cause of treatment failure. A substantial barrier to progress in medical oncology is cancer drug resistance. Resistance can develop clinically either before or after cancer treatment. According to this study, lncRNAs influence drug resistance through several different methods. LncRNAs often impact drug resistance by controlling the expression of a few intermediary regulatory variables rather than by directly affecting drug resistance. Additionally, lncRNAs have a variety of roles in cancer medication resistance. Most lncRNAs induce drug resistance when overexpressed; however, other lncRNAs have inhibitory effects. This study provides an overview of the current understanding of lncRNAs, relevance to cancer, and potential therapeutic applications.

The network structure of hematopoietic cancers

Hematopoietic cancers (HCs) are a heterogeneous group of malignancies that affect blood, bone marrow and lymphatic system. Here, by analyzing 1960 RNA-Seq samples from three independent datasets, we explored the co-expression landscape in HCs, by inferring gene co-expression networks (GCNs) with four cancer phenotypes (B and T-cell acute leukemia -BALL, TALL-, acute myeloid leukemia -AML-, and multiple myeloma -MM-) as well as non-cancer bone marrow. We characterized their structure (topological features) and function (enrichment analyses). We found that, as in other types of cancer, the highest co-expression interactions are intra-chromosomal, which is not the case for control GCNs. We also detected a highly co-expressed group of overexpressed pseudogenes in HC networks. The four GCNs present only a small fraction of common interactions, related to canonical functions, like immune response or erythrocyte differentiation. With this approach, we were able to reveal cancer-specific features useful for detection of disease manifestations.

Molecular mechanisms of multi-omic regulation in breast cancer

Breast cancer is a complex disease that is influenced by the concurrent influence of multiple genetic and environmental factors. Recent advances in genomics and other high throughput biomolecular techniques (-omics) have provided numerous insights into the molecular mechanisms underlying breast cancer development and progression. A number of these mechanisms involve multiple layers of regulation. In this review, we summarize the current knowledge on the role of multiple omics in the regulation of breast cancer, including the effects of DNA methylation, non-coding RNA, and other epigenomic changes. We comment on how integrating such diverse mechanisms is envisioned as key to a more comprehensive understanding of breast carcinogenesis and cancer biology with relevance to prognostics, diagnostics and therapeutics. We also discuss the potential clinical implications of these findings and highlight areas for future research. Overall, our understanding of the molecular mechanisms of multi-omic regulation in breast cancer is rapidly increasing and has the potential to inform the development of novel therapeutic approaches for this disease.

Single-cell gene regulatory network prediction by explainable AI

The molecular heterogeneity of cancer cells contributes to the often partial response to targeted therapies and relapse of disease due to the escape of resistant cell populations. While single-cell sequencing has started to improve our understanding of this heterogeneity, it offers a mostly descriptive view on cellular types and states. To obtain more functional insights, we propose scGeneRAI, an explainable deep learning approach that uses layer-wise relevance propagation (LRP) to infer gene regulatory networks from static single-cell RNA sequencing data for individual cells. We benchmark our method with synthetic data and apply it to single-cell RNA sequencing data of a cohort of human lung cancers. From the predicted single-cell networks our approach reveals characteristic network patterns for tumor cells and normal epithelial cells and identifies subnetworks that are observed only in (subgroups of) tumor cells of certain patients. While current state-of-the-art methods are limited by their ability to only predict average networks for cell populations, our approach facilitates the reconstruction of networks down to the level of single cells which can be utilized to characterize the heterogeneity of gene regulation within and across tumors.

The Role of Copy Number Variants in Gene Co-Expression Patterns for Luminal B Breast Tumors

Gene co-expression networks have become a usual approach to integrate the vast amounts of information coming from gene expression studies in cancer cohorts. The reprogramming of the gene regulatory control and the molecular pathways depending on such control are central to the characterization of the disease, aiming to unveil the consequences for cancer prognosis and therapeutics. There is, however, a multitude of factors which have been associated with anomalous control of gene expression in cancer. In the particular case of co-expression patterns, we have previously documented a phenomenon of loss of long distance co-expression in several cancer types, including breast cancer. Of the many potential factors that may contribute to this phenomenology, copy number variants (CNVs) have been often discussed. However, no systematic assessment of the role that CNVs may play in shaping gene co-expression patterns in breast cancer has been performed to date. For this reason we have decided to develop such analysis. In this study, we focus on using probabilistic modeling techniques to evaluate to what extent CNVs affect the phenomenon of long/short range co-expression in Luminal B breast tumors. We analyzed the co-expression patterns in chromosome 8, since it is known to be affected by amplifications/deletions during cancer development. We found that the CNVs pattern in chromosome 8 of Luminal B network does not alter the co-expression patterns significantly, which means that the co-expression program in this cancer phenotype is not determined by CNV structure. Additionally, we found that region 8q24.3 is highly dense in interactions, as well as region p21.3. The most connected genes in this network belong to those cytobands and are associated with several manifestations of cancer in different tissues. Interestingly, among the most connected genes, we found MAF1 and POLR3D, which may constitute an axis of regulation of gene transcription, in particular for non-coding RNA species. We believe that by advancing on our knowledge of the molecular mechanisms behind gene regulation in cancer, we will be better equipped, not only to understand tumor biology, but also to broaden the scope of diagnostic, prognostic and therapeutic interventions to ultimately benefit oncologic patients.

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.

MetaCancer: A deep learning-based pan-cancer metastasis prediction model developed using multi-omics data

Predicting metastasis in the early stages means that clinicians have more time to adjust a treatment regimen to target the primary and metastasized cancer. In this regard, several computational approaches are being developed to identify metastasis early. However, most of the approaches focus on changes on one genomic level only, and they are not being developed from a pan-cancer perspective. Thus, we here present a deep learning (DL)-based model, MetaCancer, that differentiates pan-cancer metastasis status based on three heterogeneous data layers. In particular, we built the DL-based model using 400 patients' data that includes RNA sequencing (RNA-Seq), microRNA sequencing (microRNA-Seq), and DNA methylation data from The Cancer Genome Atlas (TCGA). We quantitatively assess the proposed convolutional variational autoencoder (CVAE) and alternative feature extraction methods. We further show that integrating mRNA, microRNA, and DNA methylation data as features improves our model's performance compared to when we used mRNA data only. In addition, we show that the mRNA-related features make a more significant contribution when attempting to distinguish the primary tumors from metastatic ones computationally. Lastly, we show that our DL model significantly outperformed a machine learning (ML) ensemble method based on various metrics.

Multi-Omics Model Applied to Cancer Genetics

In this review, we focus on bioinformatic oncology as an integrative discipline that incorporates knowledge from the mathematical, physical, and computational fields to further the biomedical understanding of cancer. Before providing a deeper insight into the bioinformatics approach and utilities involved in oncology, we must understand what is a system biology framework and the genetic connection, because of the high heterogenicity of the backgrounds of people approaching precision medicine. In fact, it is essential to providing general theoretical information on genomics, epigenomics, and transcriptomics to understand the phases of multi-omics approach. We consider how to create a multi-omics model. In the last section, we describe the new frontiers and future perspectives of this field.

Exploiting the DepMap cancer dependency data using the depmap R package

The `depmap` package facilitates access in the R environment to the data from the DepMap project, a multi-year collaborative effort by the Broad Institute and Wellcome Sanger Institute, mapping genetic and chemical dependencies and other molecular biological measurements of over 1700 cancer cell lines. The 'depmap' package formats this data to simply the use of popular R data analysis and visualizing tools such as 'dplyr' and 'ggplot2'. In addition, the 'depmap' package utilizes 'ExperimentHub', storing versions of the DepMap data accessible from the Cloud, which may be selectively downloaded, providing a reproducible research framework to support exploiting this data. This paper describes a workflow demonstrating how to access and visualize the DepMap data in R using this package.

Luminal A Breast Cancer Co-expression Network: Structural and Functional Alterations

Luminal A is the most common breast cancer molecular subtype in women worldwide. These tumors have characteristic yet heterogeneous alterations at the genomic and transcriptomic level. Gene co-expression networks (GCNs) have contributed to better characterize the cancerous phenotype. We have previously shown an imbalance in the proportion of intra-chromosomal (cis-) over inter-chromosomal (trans-) interactions when comparing cancer and healthy tissue GCNs. In particular, for breast cancer molecular subtypes (Luminal A included), the majority of high co-expression interactions connect gene-pairs in the same chromosome, a phenomenon that we have called loss of trans- co-expression. Despite this phenomenon has been described, the functional implication of this specific network topology has not been studied yet. To understand the biological role that communities of co-expressed genes may have, we constructed GCNs for healthy and Luminal A phenotypes. Network modules were obtained based on their connectivity patterns and they were classified according to their chromosomal homophily (proportion of cis-/trans- interactions). A functional overrepresentation analysis was performed on communities in both networks to observe the significantly enriched processes for each community. We also investigated possible mechanisms for which the loss of trans- co-expression emerges in cancer GCN. To this end we evaluated transcription factor binding sites, CTCF binding sites, differential gene expression and copy number alterations (CNAs) in the cancer GCN. We found that trans- communities in Luminal A present more significantly enriched categories than cis- ones. Processes, such as angiogenesis, cell proliferation, or cell adhesion were found in trans- modules. The differential expression analysis showed that FOXM1, CENPA, and CIITA transcription factors, exert a major regulatory role on their communities by regulating expression of their target genes in other chromosomes. Finally, identification of CNAs, displayed a high enrichment of deletion peaks in cis- communities. With this approach, we demonstrate that network topology determine, to at certain extent, the function in Luminal A breast cancer network. Furthermore, several mechanisms seem to be acting together to avoid trans- co-expression. Since this phenomenon has been observed in other cancer tissues, a remaining question is whether the loss of long distance co-expression is a novel hallmark of cancer.

Gene Expression and Co-expression Networks Are Strongly Altered Through Stages in Clear Cell Renal Carcinoma

Clear cell renal carcinoma (ccRC) is a highly heterogeneous and progressively malignant disease. Analyzing ccRC progression in terms of modifications at the molecular and genetic level may help us to develop a broader understanding of its patho-physiology and may give us a glimpse toward improved therapeutics. In this work, by using TCGA data, we studied the molecular progression of the four main ccRC stages (i, ii, iii, iv) in two different yet complementary approaches: (a) gene expression and (b) gene co-expression. For (a) we analyzed the differential gene expression between each stage and the control non-cancer group. We compared the progression molecular signature between stages, and observed those genes that change their expression patterns through progression stages. For (b) we constructed and analyzed co-expression networks for the four ccRC progression stages, as well as for the control phenotype, to observe whether and how the co-expression landscape changes with progression. We separated genomic interactions into intra-chromosome (cis-) and inter-chromosome (trans-). Finally, we intersected those networks and performed functional enrichment analysis. All calculations were made over different network sizes, from the top 100 edges to top 1,000,000. We show that differential expression is quite similar between ccRC progression stages. However, interestingly, two genes, namely SLC6A19 and PLG show a significant progressive decrease in their expression according to ccRC stage, meanwhile two other genes, SAA2-SAA4 and CXCL13 show progressive increase. Despite the high similarity between gene expression profiles, all networks are substantially different between them in terms of their topological features. Control network has a larger proportion of trans- interactions, meanwhile for any stage, the amount of cis- interactions is higher, independent of the network cut-off. The majority of interactions in any network are phenotype-specific. Only 189 interactions are shared between the five networks, and 533 edges are ccRC-specific, independent of the stage. The small resulting connected components in both cases are formed by genes with the same differential expression trend, and are associated with important biological processes, such as cell cycle or immune system, suggesting that activity of these categories follows the differential expression trend. With this approach we have shown that, even if the expression program is similar during ccRC progression, the co-expression programs strongly differ. More research is needed to understand the delicate interplay between expression and co-expression, but this is a first approach to enclose both approaches in an integrative view aimed at a deeper understanding in gene regulation in tumor evolution.

Gene Co-expression Is Distance-Dependent in Breast Cancer

Breast carcinomas are characterized by anomalous gene regulatory programs. As is well-known, gene expression programs are able to shape phenotypes. Hence, the understanding of gene co-expression may shed light on the underlying mechanisms behind the transcriptional regulatory programs affecting tumor development and evolution. For instance, in breast cancer, there is a clear loss of inter-chromosomal (trans-) co-expression, compared with healthy tissue. At the same time cis- (intra-chromosomal) interactions are favored in breast tumors. In order to have a deeper understanding of regulatory phenomena in cancer, here, we constructed Gene Co-expression Networks by using TCGA-derived RNA-seq whole-genome samples corresponding to the four breast cancer molecular subtypes, as well as healthy tissue. We quantify the cis-/trans- co-expression imbalance in all phenotypes. Additionally, we measured the association between co-expression and physical distance between genes, and characterized the ratio of intra/inter-cytoband interactions per phenotype. We confirmed loss of trans- co-expression in all molecular subtypes. We also observed that gene cis- co-expression decays abruptly with distance in all tumors in contrast with healthy tissue. We observed co-expressed gene hotspots, that tend to be connected at cytoband regions, and coincide accurately with already known copy number altered regions, such as Chr17q12, or Chr8q24.3 for all subtypes. Our methodology recovered different alterations already reported for specific breast cancer subtypes, showing how co-expression network approaches might help to capture distinct events that modify the cell regulatory program.

Strategies for Functional Interrogation of Big Cancer Data Using Drosophila Cancer Models

Rapid development of high throughput genome analysis technologies accompanied by significant reduction in costs has led to the accumulation of an incredible amount of data during the last decade. The emergence of big data has had a particularly significant impact in biomedical research by providing unprecedented, systems-level access to many disease states including cancer, and has created promising opportunities as well as new challenges. Arguably, the most significant challenge cancer research currently faces is finding effective ways to use big data to improve our understanding of molecular mechanisms underlying tumorigenesis and developing effective new therapies. Functional exploration of these datasets and testing predictions from computational approaches using experimental models to interrogate their biological relevance is a key step towards achieving this goal. Given the daunting scale and complexity of the big data available, experimental systems like Drosophila that allow large-scale functional studies and complex genetic manipulations in a rapid, cost-effective manner will be of particular importance for this purpose. Findings from these large-scale exploratory functional studies can then be used to formulate more specific hypotheses to be explored in mammalian models. Here, I will discuss several strategies for functional exploration of big cancer data using Drosophila cancer models.

On a Class of Tensor Markov Fields

Here, we introduce a class of Tensor Markov Fields intended as probabilistic graphical models from random variables spanned over multiplexed contexts. These fields are an extension of Markov Random Fields for tensor-valued random variables. By extending the results of Dobruschin, Hammersley and Clifford to such tensor valued fields, we proved that tensor Markov fields are indeed Gibbs fields, whenever strictly positive probability measures are considered. Hence, there is a direct relationship with many results from theoretical statistical mechanics. We showed how this class of Markov fields it can be built based on a statistical dependency structures inferred on information theoretical grounds over empirical data. Thus, aside from purely theoretical interest, the Tensor Markov Fields described here may be useful for mathematical modeling and data analysis due to their intrinsic simplicity and generality.

Computational Oncology in the Multi-Omics Era: State of the Art

Cancer is the quintessential complex disease. As technologies evolve faster each day, we are able to quantify the different layers of biological elements that contribute to the emergence and development of malignancies. In this multi-omics context, the use of integrative approaches is mandatory in order to gain further insights on oncological phenomena, and to move forward toward the precision medicine paradigm. In this review, we will focus on computational oncology as an integrative discipline that incorporates knowledge from the mathematical, physical, and computational fields to further the biomedical understanding of cancer. We will discuss the current roles of computation in oncology in the context of multi-omic technologies, which include: data acquisition and processing; data management in the clinical and research settings; classification, diagnosis, and prognosis; and the development of models in the research setting, including their use for therapeutic target identification. We will discuss the machine learning and network approaches as two of the most promising emerging paradigms, in computational oncology. These approaches provide a foundation on how to integrate different layers of biological description into coherent frameworks that allow advances both in the basic and clinical settings.