@article {403, title = {A comparison of mechanistic signaling pathway activity analysis methods.}, journal = {Brief Bioinform}, volume = {20}, year = {2019}, month = {2019 09 27}, pages = {1655-1668}, abstract = {

Understanding the aspects of cell functionality that account for disease mechanisms or drug modes of action is a main challenge for precision medicine. Classical gene-based approaches ignore the modular nature of most human traits, whereas conventional pathway enrichment approaches produce only illustrative results of limited practical utility. Recently, a family of new methods has emerged that change the focus from the whole pathways to the definition of elementary subpathways within them that have any mechanistic significance and to the study of their activities. Thus, mechanistic pathway activity (MPA) methods constitute a new paradigm that allows recoding poorly informative genomic measurements into cell activity quantitative values and relate them to phenotypes. Here we provide a review on the MPA methods available and explain their contribution to systems medicine approaches for addressing challenges in the diagnostic and treatment of complex diseases.

}, keywords = {Algorithms, Humans, Postmortem Changes, Signal Transduction, Systems biology, Transcriptome}, issn = {1477-4054}, doi = {10.1093/bib/bby040}, author = {Amadoz, Alicia and Hidalgo, Marta R and Cubuk, Cankut and Carbonell-Caballero, Jos{\'e} and Dopazo, Joaquin} } @article {422, title = {Differential metabolic activity and discovery of therapeutic targets using summarized metabolic pathway models.}, journal = {NPJ Syst Biol Appl}, volume = {5}, year = {2019}, month = {2019}, pages = {7}, abstract = {

In spite of the increasing availability of genomic and transcriptomic data, there is still a gap between the detection of perturbations in gene expression and the understanding of their contribution to the molecular mechanisms that ultimately account for the phenotype studied. Alterations in the metabolism are behind the initiation and progression of many diseases, including cancer. The wealth of available knowledge on metabolic processes can therefore be used to derive mechanistic models that link gene expression perturbations to changes in metabolic activity that provide relevant clues on molecular mechanisms of disease and drug modes of action (MoA). In particular, pathway modules, which recapitulate the main aspects of metabolism, are especially suitable for this type of modeling. We present Metabolizer, a web-based application that offers an intuitive, easy-to-use interactive interface to analyze differences in pathway metabolic module activities that can also be used for class prediction and in silico prediction of knock-out (KO) effects. Moreover, Metabolizer can automatically predict the optimal KO intervention for restoring a diseased phenotype. We provide different types of validations of some of the predictions made by Metabolizer. Metabolizer is a web tool that allows understanding molecular mechanisms of disease or the MoA of drugs within the context of the metabolism by using gene expression measurements. In addition, this tool automatically suggests potential therapeutic targets for individualized therapeutic interventions.

}, keywords = {Computational Biology, Computer Simulation, Drug discovery, Gene Regulatory Networks, Humans, Internet, Metabolic Networks and Pathways, Models, Biological, Neoplasms, Phenotype, Software, Transcriptome}, issn = {2056-7189}, doi = {10.1038/s41540-019-0087-2}, author = {Cubuk, Cankut and Hidalgo, Marta R and Amadoz, Alicia and Rian, Kinza and Salavert, Francisco and Pujana, Miguel A and Mateo, Francesca and Herranz, Carmen and Carbonell-Caballero, Jos{\'e} and Dopazo, Joaquin} } @article {397, title = {The effects of death and post-mortem cold ischemia on human tissue transcriptomes.}, journal = {Nat Commun}, volume = {9}, year = {2018}, month = {2018 02 13}, pages = {490}, abstract = {

Post-mortem tissues samples are a key resource for investigating patterns of gene expression. However, the processes triggered by death and the post-mortem interval (PMI) can significantly alter physiologically normal RNA levels. We investigate the impact of PMI on gene expression using data from multiple tissues of post-mortem donors obtained from the GTEx project. We find that many genes change expression over relatively short PMIs in a tissue-specific manner, but this potentially confounding effect in a biological analysis can be minimized by taking into account appropriate covariates. By comparing ante- and post-mortem blood samples, we identify the cascade of transcriptional events triggered by death of the organism. These events do not appear to simply reflect stochastic variation resulting from mRNA degradation, but active and ongoing regulation of transcription. Finally, we develop a model to predict the time since death from the analysis of the transcriptome of a few readily accessible tissues.

}, keywords = {Blood, Cold Ischemia, Death, Female, gene expression, Humans, Models, Biological, Postmortem Changes, RNA, Messenger, Stochastic Processes, Transcriptome}, issn = {2041-1723}, doi = {10.1038/s41467-017-02772-x}, author = {Ferreira, Pedro G and Mu{\~n}oz-Aguirre, Manuel and Reverter, Ferran and S{\'a} Godinho, Caio P and Sousa, Abel and Amadoz, Alicia and Sodaei, Reza and Hidalgo, Marta R and Pervouchine, Dmitri and Carbonell-Caballero, Jos{\'e} and Nurtdinov, Ramil and Breschi, Alessandra and Amador, Raziel and Oliveira, Patr{\'\i}cia and Cubuk, Cankut and Curado, Jo{\~a}o and Aguet, Fran{\c c}ois and Oliveira, Carla and Dopazo, Joaquin and Sammeth, Michael and Ardlie, Kristin G and Guig{\'o}, Roderic} } @article {405, title = {Gene Expression Integration into Pathway Modules Reveals a Pan-Cancer Metabolic Landscape.}, journal = {Cancer Res}, volume = {78}, year = {2018}, month = {2018 11 01}, pages = {6059-6072}, abstract = {

Metabolic reprogramming plays an important role in cancer development and progression and is a well-established hallmark of cancer. Despite its inherent complexity, cellular metabolism can be decomposed into functional modules that represent fundamental metabolic processes. Here, we performed a pan-cancer study involving 9,428 samples from 25 cancer types to reveal metabolic modules whose individual or coordinated activity predict cancer type and outcome, in turn highlighting novel therapeutic opportunities. Integration of gene expression levels into metabolic modules suggests that the activity of specific modules differs between cancers and the corresponding tissues of origin. Some modules may cooperate, as indicated by the positive correlation of their activity across a range of tumors. The activity of many metabolic modules was significantly associated with prognosis at a stronger magnitude than any of their constituent genes. Thus, modules may be classified as tumor suppressors and oncomodules according to their potential impact on cancer progression. Using this modeling framework, we also propose novel potential therapeutic targets that constitute alternative ways of treating cancer by inhibiting their reprogrammed metabolism. Collectively, this study provides an extensive resource of predicted cancer metabolic profiles and dependencies. Combining gene expression with metabolic modules identifies molecular mechanisms of cancer undetected on an individual gene level and allows discovery of new potential therapeutic targets. .

}, keywords = {Cell Line, Tumor, Cluster Analysis, Disease Progression, Gene Expression Profiling, Gene Expression Regulation, Neoplastic, Gene Regulatory Networks, Humans, Kaplan-Meier Estimate, Metabolome, mutation, Neoplasms, Oncogenes, Phenotype, Prognosis, RNA, Small Interfering, Sequence Analysis, RNA, Transcriptome, Treatment Outcome}, issn = {1538-7445}, doi = {10.1158/0008-5472.CAN-17-2705}, author = {Cubuk, Cankut and Hidalgo, Marta R and Amadoz, Alicia and Pujana, Miguel A and Mateo, Francesca and Herranz, Carmen and Carbonell-Caballero, Jos{\'e} and Dopazo, Joaquin} } @article {404, title = {Models of cell signaling uncover molecular mechanisms of high-risk neuroblastoma and predict disease outcome.}, journal = {Biol Direct}, volume = {13}, year = {2018}, month = {2018 08 22}, pages = {16}, abstract = {

BACKGROUND: Despite the progress in neuroblastoma therapies the mortality of high-risk patients is still high (40-50\%) and the molecular basis of the disease remains poorly known. Recently, a mathematical model was used to demonstrate that the network regulating stress signaling by the c-Jun N-terminal kinase pathway played a crucial role in survival of patients with neuroblastoma irrespective of their MYCN amplification status. This demonstrates the enormous potential of computational models of biological modules for the discovery of underlying molecular mechanisms of diseases.

RESULTS: Since signaling is known to be highly relevant in cancer, we have used a computational model of the whole cell signaling network to understand the molecular determinants of bad prognostic in neuroblastoma. Our model produced a comprehensive view of the molecular mechanisms of neuroblastoma tumorigenesis and progression.

CONCLUSION: We have also shown how the activity of signaling circuits can be considered a reliable model-based prognostic biomarker.

REVIEWERS: This article was reviewed by Tim Beissbarth, Wenzhong Xiao and Joanna Polanska. For the full reviews, please go to the Reviewers{\textquoteright} comments section.

}, keywords = {Computational Biology, Gene Expression Regulation, Neoplastic, Humans, JNK Mitogen-Activated Protein Kinases, Models, Theoretical, Neuroblastoma, Signal Transduction}, issn = {1745-6150}, doi = {10.1186/s13062-018-0219-4}, author = {Hidalgo, Marta R and Amadoz, Alicia and Cubuk, Cankut and Carbonell-Caballero, Jos{\'e} and Dopazo, Joaquin} } @article {434, title = {High throughput estimation of functional cell activities reveals disease mechanisms and predicts relevant clinical outcomes.}, journal = {Oncotarget}, volume = {8}, year = {2017}, month = {2017 Jan 17}, pages = {5160-5178}, abstract = {

Understanding the aspects of the cell functionality that account for disease or drug action mechanisms is a main challenge for precision medicine. Here we propose a new method that models cell signaling using biological knowledge on signal transduction. The method recodes individual gene expression values (and/or gene mutations) into accurate measurements of changes in the activity of signaling circuits, which ultimately constitute high-throughput estimations of cell functionalities caused by gene activity within the pathway. Moreover, such estimations can be obtained either at cohort-level, in case/control comparisons, or personalized for individual patients. The accuracy of the method is demonstrated in an extensive analysis involving 5640 patients from 12 different cancer types. Circuit activity measurements not only have a high diagnostic value but also can be related to relevant disease outcomes such as survival, and can be used to assess therapeutic interventions.

}, keywords = {Computational Biology, gene expression, Gene Regulatory Networks, Humans, mutation, Neoplasms, Precision Medicine, Sequence Analysis, RNA, Signal Transduction}, issn = {1949-2553}, doi = {10.18632/oncotarget.14107}, author = {Hidalgo, Marta R and Cubuk, Cankut and Amadoz, Alicia and Salavert, Francisco and Carbonell-Caballero, Jos{\'e} and Dopazo, Joaquin} } @article {388, title = {Reference genome assessment from a population scale perspective: an accurate profile of variability and noise.}, journal = {Bioinformatics}, volume = {33}, year = {2017}, month = {2017 Nov 15}, pages = {3511-3517}, abstract = {

Motivation: Current plant and animal genomic studies are often based on newly assembled genomes that have not been properly consolidated. In this scenario, misassembled regions can easily lead to false-positive findings. Despite quality control scores are included within genotyping protocols, they are usually employed to evaluate individual sample quality rather than reference sequence reliability. We propose a statistical model that combines quality control scores across samples in order to detect incongruent patterns at every genomic region. Our model is inherently robust since common artifact signals are expected to be shared between independent samples over misassembled regions of the genome.

Results: The reliability of our protocol has been extensively tested through different experiments and organisms with accurate results, improving state-of-the-art methods. Our analysis demonstrates synergistic relations between quality control scores and allelic variability estimators, that improve the detection of misassembled regions, and is able to find strong artifact signals even within the human reference assembly. Furthermore, we demonstrated how our model can be trained to properly rank the confidence of a set of candidate variants obtained from new independent samples.

Availability and implementation: This tool is freely available at http://gitlab.com/carbonell/ces.

Contact: jcarbonell.cipf@gmail.com or joaquin.dopazo@juntadeandalucia.es.

Supplementary information: Supplementary data are available at Bioinformatics online.

}, keywords = {Animals, Genetic Variation, Genome, Genomics, Genotype, Humans, Models, Statistical, Quality Control, Reproducibility of Results, Software}, issn = {1367-4811}, doi = {10.1093/bioinformatics/btx482}, url = {https://academic.oup.com/bioinformatics/article-lookup/doi/10.1093/bioinformatics/btx482}, author = {Carbonell-Caballero, Jos{\'e} and Amadoz, Alicia and Alonso, Roberto and Hidalgo, Marta R and Cubuk, Cankut and Conesa, David and L{\'o}pez-Qu{\'\i}lez, Antonio and Dopazo, Joaquin} } @article {1184, title = {267 Spanish exomes reveal population-specific differences in disease-related genetic variation.}, journal = {Molecular biology and evolution}, year = {2016}, month = {2016 Jan 13}, abstract = {Recent results from large-scale genomic projects suggest that allele frequencies, which are highly relevant for medical purposes, differ considerably across different populations. The need for a detailed catalogue of local variability motivated the whole exome sequencing of 267 unrelated individuals, representative of the healthy Spanish population. Like in other studies, a considerable number of rare variants were found (almost one third of the described variants). There were also relevant differences in allelic frequencies in polymorphic variants, including about 10,000 polymorphisms private to the Spanish population. The allelic frequencies of variants conferring susceptibility to complex diseases (including cancer, schizophrenia, Alzheimer disease, type 2 diabetes and other pathologies) were overall similar to those of other populations. However, the trend is the opposite for variants linked to Mendelian and rare diseases (including several retinal degenerative dystrophies and cardiomyopathies) that show marked frequency differences between populations. Interestingly, a correspondence between differences in allelic frequencies and disease prevalence was found, highlighting the relevance of frequency differences in disease risk. These differences are also observed in variants that disrupt known drug binding sites, suggesting an important role for local variability in population-specific drug resistances or adverse effects. We have made the Spanish population variant server web page that contains population frequency information for the complete list of 170,888 variant positions we found publicly available (http://spv.babelomics.org/), We show that it if fundamental to determine population-specific variant frequencies in order to distinguish real disease associations from population-specific polymorphisms.}, keywords = {disease, NGS, polymorphisms, Population genomics, prioritization, SNP}, issn = {1537-1719}, doi = {10.1093/molbev/msw005}, url = {https://mbe.oxfordjournals.org/content/early/2016/02/17/molbev.msw005.full}, author = {Joaqu{\'\i}n Dopazo and Amadoz, Alicia and Bleda, Marta and Garc{\'\i}a-Alonso, Luz and Alem{\'a}n, Alejandro and Garcia-Garcia, Francisco and Rodriguez, Juan A and Daub, Josephine T and Muntan{\'e}, Gerard and Antonio Rueda and Vela-Boza, Alicia and L{\'o}pez-Domingo, Francisco J and Florido, Javier P and Arce, Pablo and Ruiz-Ferrer, Macarena and M{\'e}ndez-Vidal, Cristina and Arnold, Todd E and Spleiss, Olivia and Alvarez-Tejado, Miguel and Navarro, Arcadi and Bhattacharya, Shomi S and Borrego, Salud and Santoyo-L{\'o}pez, Javier and Anti{\v n}olo, Guillermo} } @article {1203, title = {Actionable pathways: interactive discovery of therapeutic targets using signaling pathway models.}, journal = {Nucleic acids research}, year = {2016}, month = {2016 May 2}, abstract = {The discovery of actionable targets is crucial for targeted therapies and is also a constituent part of the drug discovery process. The success of an intervention over a target depends critically on its contribution, within the complex network of gene interactions, to the cellular processes responsible for disease progression or therapeutic response. Here we present PathAct, a web server that predicts the effect that interventions over genes (inhibitions or activations that simulate knock-outs, drug treatments or over-expressions) can have over signal transmission within signaling pathways and, ultimately, over the cell functionalities triggered by them. PathAct implements an advanced graphical interface that provides a unique interactive working environment in which the suitability of potentially actionable genes, that could eventually become drug targets for personalized or individualized therapies, can be easily tested. The PathAct tool can be found at: http://pathact.babelomics.org.}, keywords = {actionable genes, Disease mechanism, drug action mechanism, Drug discovery, pathway analysis, personalized medicine, signalling, therapeutic targets}, issn = {1362-4962}, doi = {10.1093/nar/gkw369}, url = {http://nar.oxfordjournals.org/content/early/2016/05/02/nar.gkw369.full}, author = {Salavert, Francisco and Hidago, Marta R and Amadoz, Alicia and Cubuk, Cankut and Medina, Ignacio and Crespo, Daniel and Carbonell-Caballero, Jos{\'e} and Joaqu{\'\i}n Dopazo} } @article {1129, title = {Babelomics 5.0: functional interpretation for new generations of genomic data.}, journal = {Nucleic acids research}, volume = {43}, number = {W1}, year = {2015}, month = {2015 Apr 20}, pages = {W117-W121}, abstract = {Babelomics has been running for more than one decade offering a user-friendly interface for the functional analysis of gene expression and genomic data. Here we present its fifth release, which includes support for Next Generation Sequencing data including gene expression (RNA-seq), exome or genome resequencing. Babelomics has simplified its interface, being now more intuitive. Improved visualization options, such as a genome viewer as well as an interactive network viewer, have been implemented. New technical enhancements at both, client and server sides, makes the user experience faster and more dynamic. Babelomics offers user-friendly access to a full range of methods that cover: (i) primary data analysis, (ii) a variety of tests for different experimental designs and (iii) different enrichment and network analysis algorithms for the interpretation of the results of such tests in the proper functional context. In addition to the public server, local copies of Babelomics can be downloaded and installed. Babelomics is freely available at: http://www.babelomics.org.}, keywords = {babelomics, data integration, gene set analysis, interactome, network analysis, NGS, RNA-seq, Systems biology, transcriptomics}, issn = {1362-4962}, doi = {10.1093/nar/gkv384}, url = {http://nar.oxfordjournals.org/content/43/W1/W117}, author = {Alonso, Roberto and Salavert, Francisco and Garcia-Garcia, Francisco and Carbonell-Caballero, Jos{\'e} and Bleda, Marta and Garc{\'\i}a-Alonso, Luz and Sanchis-Juan, Alba and Perez-Gil, Daniel and Marin-Garcia, Pablo and S{\'a}nchez, Rub{\'e}n and Cubuk, Cankut and Hidalgo, Marta R and Amadoz, Alicia and Hernansaiz-Ballesteros, Rosa D and Alem{\'a}n, Alejandro and T{\'a}rraga, Joaqu{\'\i}n and Montaner, David and Medina, Ignacio and Dopazo, Joaquin} } @article {474, title = {Using activation status of signaling pathways as mechanism-based biomarkers to predict drug sensitivity.}, journal = {Sci Rep}, volume = {5}, year = {2015}, month = {2015 Dec 18}, pages = {18494}, abstract = {

Many complex traits, as drug response, are associated with changes in biological pathways rather than being caused by single gene alterations. Here, a predictive framework is presented in which gene expression data are recoded into activity statuses of signal transduction circuits (sub-pathways within signaling pathways that connect receptor proteins to final effector proteins that trigger cell actions). Such activity values are used as features by a prediction algorithm which can efficiently predict a continuous variable such as the IC50 value. The main advantage of this prediction method is that the features selected by the predictor, the signaling circuits, are themselves rich-informative, mechanism-based biomarkers which provide insight into or drug molecular mechanisms of action (MoA).

}, keywords = {Algorithms, Antineoplastic Agents, biomarkers, Cell Line, Tumor, Cell Survival, gene expression, Humans, Lethal Dose 50, Neoplasms, Phosphorylation, Proteins, Signal Transduction}, issn = {2045-2322}, doi = {10.1038/srep18494}, author = {Amadoz, Alicia and Sebasti{\'a}n-Leon, Patricia and Vidal, Enrique and Salavert, Francisco and Dopazo, Joaquin} } @article {565, title = {Understanding disease mechanisms with models of signaling pathway activities}, journal = {BMC systems biology}, volume = {8}, year = {2014}, month = {10}, pages = {121}, doi = {10.1186/s12918-014-0121-3}, author = {Sebasti{\'a}n-Leon, Patricia and Vidal, Enrique and Minguez, Pablo and Conesa, Ana and Tarazona, Sonia and Amadoz, Alicia and Armero, Carmen and Salavert Torres, Francisco and Vidal-Puig, Antonio and Montaner, David and Dopazo, Joaquin} } @article {1093, title = {Understanding disease mechanisms with models of signaling pathway activities.}, journal = {BMC systems biology}, volume = {8}, year = {2014}, month = {2014 Oct 25}, pages = {121}, abstract = {BackgroundUnderstanding the aspects of the cell functionality that account for disease or drug action mechanisms is one of the main challenges in the analysis of genomic data and is on the basis of the future implementation of precision medicine.ResultsHere we propose a simple probabilistic model in which signaling pathways are separated into elementary sub-pathways or signal transmission circuits (which ultimately trigger cell functions) and then transforms gene expression measurements into probabilities of activation of such signal transmission circuits. Using this model, differential activation of such circuits between biological conditions can be estimated. Thus, circuit activation statuses can be interpreted as biomarkers that discriminate among the compared conditions. This type of mechanism-based biomarkers accounts for cell functional activities and can easily be associated to disease or drug action mechanisms. The accuracy of the proposed model is demonstrated with simulations and real datasets.ConclusionsThe proposed model provides detailed information that enables the interpretation disease mechanisms as a consequence of the complex combinations of altered gene expression values. Moreover, it offers a framework for suggesting possible ways of therapeutic intervention in a pathologically perturbed system.}, keywords = {Disease mechanism, pathway, signalling, Systems biology}, issn = {1752-0509}, doi = {10.1186/s12918-014-0121-3}, url = {http://www.biomedcentral.com/1752-0509/8/121/abstract}, author = {Sebasti{\'a}n-Leon, Patricia and Vidal, Enrique and Minguez, Pablo and Ana Conesa and Sonia Tarazona and Amadoz, Alicia and Armero, Carmen and Salavert, Francisco and Vidal-Puig, Antonio and Montaner, David and Joaqu{\'\i}n Dopazo} } @article {512, title = {Discovering the hidden sub-network component in a ranked list of genes or proteins derived from genomic experiments.}, journal = {Nucleic Acids Res}, volume = {40}, year = {2012}, month = {2012 Nov 01}, pages = {e158}, abstract = {

Genomic experiments (e.g. differential gene expression, single-nucleotide polymorphism association) typically produce ranked list of genes. We present a simple but powerful approach which uses protein-protein interaction data to detect sub-networks within such ranked lists of genes or proteins. We performed an exhaustive study of network parameters that allowed us concluding that the average number of components and the average number of nodes per component are the parameters that best discriminate between real and random networks. A novel aspect that increases the efficiency of this strategy in finding sub-networks is that, in addition to direct connections, also connections mediated by intermediate nodes are considered to build up the sub-networks. The possibility of using of such intermediate nodes makes this approach more robust to noise. It also overcomes some limitations intrinsic to experimental designs based on differential expression, in which some nodes are invariant across conditions. The proposed approach can also be used for candidate disease-gene prioritization. Here, we demonstrate the usefulness of the approach by means of several case examples that include a differential expression analysis in Fanconi Anemia, a genome-wide association study of bipolar disorder and a genome-scale study of essentiality in cancer genes. An efficient and easy-to-use web interface (available at http://www.babelomics.org) based on HTML5 technologies is also provided to run the algorithm and represent the network.

}, keywords = {Bipolar Disorder, Fanconi Anemia, Gene Regulatory Networks, Genes, Neoplasm, Genome-Wide Association Study, Genomics, Humans, Protein Interaction Mapping}, issn = {1362-4962}, doi = {10.1093/nar/gks699}, author = {Garc{\'\i}a-Alonso, Luz and Alonso, Roberto and Vidal, Enrique and Amadoz, Alicia and De Maria, Alejandro and Minguez, Pablo and Medina, Ignacio and Dopazo, Joaquin} }