%0 Journal Article %J Int J Mol Sci %D 2023 %T Crosstalk between Metabolite Production and Signaling Activity in Breast Cancer. %A Cubuk, Cankut %A Loucera, Carlos %A Peña-Chilet, Maria %A Dopazo, Joaquin %X

The reprogramming of metabolism is a recognized cancer hallmark. It is well known that different signaling pathways regulate and orchestrate this reprogramming that contributes to cancer initiation and development. However, recent evidence is accumulating, suggesting that several metabolites could play a relevant role in regulating signaling pathways. To assess the potential role of metabolites in the regulation of signaling pathways, both metabolic and signaling pathway activities of Breast invasive Carcinoma (BRCA) have been modeled using mechanistic models. Gaussian Processes, powerful machine learning methods, were used in combination with SHapley Additive exPlanations (SHAP), a recent methodology that conveys causality, to obtain potential causal relationships between the production of metabolites and the regulation of signaling pathways. A total of 317 metabolites were found to have a strong impact on signaling circuits. The results presented here point to the existence of a complex crosstalk between signaling and metabolic pathways more complex than previously was thought.

%B Int J Mol Sci %V 24 %8 2023 Apr 18 %G eng %N 8 %R 10.3390/ijms24087450 %0 Journal Article %J Cell Rep %D 2023 %T Defective extracellular matrix remodeling in brown adipose tissue is associated with fibro-inflammation and reduced diet-induced thermogenesis. %A Pellegrinelli, Vanessa %A Figueroa-Juárez, Elizabeth %A Samuelson, Isabella %A U-Din, Mueez %A Rodriguez-Fdez, Sonia %A Virtue, Samuel %A Leggat, Jennifer %A Cubuk, Cankut %A Peirce, Vivian J %A Niemi, Tarja %A Campbell, Mark %A Rodriguez-Cuenca, Sergio %A Dopazo, Joaquin %A Carobbio, Stefania %A Virtanen, Kirsi A %A Vidal-Puig, Antonio %X

The relevance of extracellular matrix (ECM) remodeling is reported in white adipose tissue (AT) and obesity-related dysfunctions, but little is known about the importance of ECM remodeling in brown AT (BAT) function. Here, we show that a time course of high-fat diet (HFD) feeding progressively impairs diet-induced thermogenesis concomitantly with the development of fibro-inflammation in BAT. Higher markers of fibro-inflammation are associated with lower cold-induced BAT activity in humans. Similarly, when mice are housed at thermoneutrality, inactivated BAT features fibro-inflammation. We validate the pathophysiological relevance of BAT ECM remodeling in response to temperature challenges and HFD using a model of a primary defect in the collagen turnover mediated by partial ablation of the Pepd prolidase. Pepd-heterozygous mice display exacerbated dysfunction and BAT fibro-inflammation at thermoneutrality and in HFD. Our findings show the relevance of ECM remodeling in BAT activation and provide a mechanism for BAT dysfunction in obesity.

%B Cell Rep %V 42 %P 112640 %8 2023 Jun 13 %G eng %N 6 %R 10.1016/j.celrep.2023.112640 %0 Journal Article %J Computational and Structural Biotechnology Journal %D 2021 %T Genome-scale mechanistic modeling of signaling pathways made easy: A bioconductor/cytoscape/web server framework for the analysis of omic data %A Rian, Kinza %A Hidalgo, Marta R. %A Cubuk, Cankut %A Falco, Matias M. %A Loucera, Carlos %A Esteban-Medina, Marina %A Alamo-Alvarez, Inmaculada %A Peña-Chilet, Maria %A Dopazo, Joaquin %B Computational and Structural Biotechnology Journal %V 19 %P 2968 - 2978 %8 Jan-01-2021 %G eng %U https://linkinghub.elsevier.com/retrieve/pii/S2001037021002038 %! Computational and Structural Biotechnology Journal %R 10.1016/j.csbj.2021.05.022 %0 Journal Article %J BioData Min %D 2021 %T Mechanistic modeling of the SARS-CoV-2 disease map. %A Rian, Kinza %A Esteban-Medina, Marina %A Hidalgo, Marta R %A Cubuk, Cankut %A Falco, Matias M %A Loucera, Carlos %A Gunyel, Devrim %A Ostaszewski, Marek %A Peña-Chilet, Maria %A Dopazo, Joaquin %X

Here we present a web interface that implements a comprehensive mechanistic model of the SARS-CoV-2 disease map. In this framework, the detailed activity of the human signaling circuits related to the viral infection, covering from the entry and replication mechanisms to the downstream consequences as inflammation and antigenic response, can be inferred from gene expression experiments. Moreover, the effect of potential interventions, such as knock-downs, or drug effects (currently the system models the effect of more than 8000 DrugBank drugs) can be studied. This freely available tool not only provides an unprecedentedly detailed view of the mechanisms of viral invasion and the consequences in the cell but has also the potential of becoming an invaluable asset in the search for efficient antiviral treatments.

%B BioData Min %V 14 %P 5 %8 2021 Jan 21 %G eng %N 1 %1 https://www.ncbi.nlm.nih.gov/pubmed/33478554?dopt=Abstract %R 10.1186/s13040-021-00234-1 %0 Journal Article %J Mol Med %D 2021 %T Taxonomic variations in the gut microbiome of gout patients with and without tophi might have a functional impact on urate metabolism. %A Méndez-Salazar, Eder Orlando %A Vázquez-Mellado, Janitzia %A Casimiro-Soriguer, Carlos S %A Dopazo, Joaquin %A Cubuk, Cankut %A Zamudio-Cuevas, Yessica %A Francisco-Balderas, Adriana %A Martínez-Flores, Karina %A Fernández-Torres, Javier %A Lozada-Pérez, Carlos %A Pineda, Carlos %A Sánchez-González, Austreberto %A Silveira, Luis H %A Burguete-García, Ana I %A Orbe-Orihuela, Citlalli %A Lagunas-Martínez, Alfredo %A Vazquez-Gomez, Alonso %A López-Reyes, Alberto %A Palacios-González, Berenice %A Martínez-Nava, Gabriela Angélica %K Biodiversity %K Computational Biology %K Dysbiosis %K Gastrointestinal Microbiome %K Gout %K Humans %K Metagenome %K metagenomics %K Protein Interaction Mapping %K Protein Interaction Maps %K Uric Acid %X

OBJECTIVE: To evaluate the taxonomic composition of the gut microbiome in gout patients with and without tophi formation, and predict bacterial functions that might have an impact on urate metabolism.

METHODS: Hypervariable V3-V4 regions of the bacterial 16S rRNA gene from fecal samples of gout patients with and without tophi (n = 33 and n = 25, respectively) were sequenced and compared to fecal samples from 53 healthy controls. We explored predictive functional profiles using bioinformatics in order to identify differences in taxonomy and metabolic pathways.

RESULTS: We identified a microbiome characterized by the lowest richness and a higher abundance of Phascolarctobacterium, Bacteroides, Akkermansia, and Ruminococcus_gnavus_group genera in patients with gout without tophi when compared to controls. The Proteobacteria phylum and the Escherichia-Shigella genus were more abundant in patients with tophaceous gout than in controls. Fold change analysis detected nine genera enriched in healthy controls compared to gout groups (Bifidobacterium, Butyricicoccus, Oscillobacter, Ruminococcaceae_UCG_010, Lachnospiraceae_ND2007_group, Haemophilus, Ruminococcus_1, Clostridium_sensu_stricto_1, and Ruminococcaceae_UGC_013). We found that the core microbiota of both gout groups shared Bacteroides caccae, Bacteroides stercoris ATCC 43183, and Bacteroides coprocola DSM 17136. These bacteria might perform functions linked to one-carbon metabolism, nucleotide binding, amino acid biosynthesis, and purine biosynthesis. Finally, we observed differences in key bacterial enzymes involved in urate synthesis, degradation, and elimination.

CONCLUSION: Our findings revealed that taxonomic variations in the gut microbiome of gout patients with and without tophi might have a functional impact on urate metabolism.

%B Mol Med %V 27 %P 50 %8 2021 05 24 %G eng %N 1 %1 https://www.ncbi.nlm.nih.gov/pubmed/34030623?dopt=Abstract %R 10.1186/s10020-021-00311-5 %0 Journal Article %J Cells %D 2020 %T Mechanistic Models of Signaling Pathways Reveal the Drug Action Mechanisms behind Gender-Specific Gene Expression for Cancer Treatments. %A Cubuk, Cankut %A Can, Fatma E %A Peña-Chilet, Maria %A Dopazo, Joaquin %K Female %K Gene Expression Regulation, Neoplastic %K Humans %K Male %K Neoplasms %K Signal Transduction %X

Despite the existence of differences in gene expression across numerous genes between males and females having been known for a long time, these have been mostly ignored in many studies, including drug development and its therapeutic use. In fact, the consequences of such differences over the disease mechanisms or the drug action mechanisms are completely unknown. Here we applied mechanistic mathematical models of signaling activity to reveal the ultimate functional consequences that gender-specific gene expression activities have over cell functionality and fate. Moreover, we also used the mechanistic modeling framework to simulate the drug interventions and unravel how drug action mechanisms are affected by gender-specific differential gene expression. Interestingly, some cancers have many biological processes significantly affected by these gender-specific differences (e.g., bladder or head and neck carcinomas), while others (e.g., glioblastoma or rectum cancer) are almost insensitive to them. We found that many of these gender-specific differences affect cancer-specific pathways or in physiological signaling pathways, also involved in cancer origin and development. Finally, mechanistic models have the potential to be used for finding alternative therapeutic interventions on the pathways targeted by the drug, which lead to similar results compensating the downstream consequences of gender-specific differences in gene expression.

%B Cells %V 9 %8 2020 06 29 %G eng %N 7 %1 https://www.ncbi.nlm.nih.gov/pubmed/32610626?dopt=Abstract %R 10.3390/cells9071579 %0 Journal Article %J Brief Bioinform %D 2019 %T A comparison of mechanistic signaling pathway activity analysis methods. %A Amadoz, Alicia %A Hidalgo, Marta R %A Cubuk, Cankut %A Carbonell-Caballero, José %A Dopazo, Joaquin %K Algorithms %K Humans %K Postmortem Changes %K Signal Transduction %K Systems biology %K Transcriptome %X

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.

%B Brief Bioinform %V 20 %P 1655-1668 %8 2019 09 27 %G eng %N 5 %1 https://www.ncbi.nlm.nih.gov/pubmed/29868818?dopt=Abstract %R 10.1093/bib/bby040 %0 Journal Article %J NPJ Syst Biol Appl %D 2019 %T Differential metabolic activity and discovery of therapeutic targets using summarized metabolic pathway models. %A Cubuk, Cankut %A Hidalgo, Marta R %A Amadoz, Alicia %A Rian, Kinza %A Salavert, Francisco %A Pujana, Miguel A %A Mateo, Francesca %A Herranz, Carmen %A Carbonell-Caballero, José %A Dopazo, Joaquin %K Computational Biology %K Computer Simulation %K Drug discovery %K Gene Regulatory Networks %K Humans %K Internet %K Metabolic Networks and Pathways %K Models, Biological %K Neoplasms %K Phenotype %K Software %K Transcriptome %X

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.

%B NPJ Syst Biol Appl %V 5 %P 7 %8 2019 %G eng %1 https://www.ncbi.nlm.nih.gov/pubmed/30854222?dopt=Abstract %R 10.1038/s41540-019-0087-2 %0 Journal Article %J Nat Commun %D 2018 %T The effects of death and post-mortem cold ischemia on human tissue transcriptomes. %A Ferreira, Pedro G %A Muñoz-Aguirre, Manuel %A Reverter, Ferran %A Sá Godinho, Caio P %A Sousa, Abel %A Amadoz, Alicia %A Sodaei, Reza %A Hidalgo, Marta R %A Pervouchine, Dmitri %A Carbonell-Caballero, José %A Nurtdinov, Ramil %A Breschi, Alessandra %A Amador, Raziel %A Oliveira, Patrícia %A Cubuk, Cankut %A Curado, João %A Aguet, François %A Oliveira, Carla %A Dopazo, Joaquin %A Sammeth, Michael %A Ardlie, Kristin G %A Guigó, Roderic %K Blood %K Cold Ischemia %K Death %K Female %K gene expression %K Humans %K Models, Biological %K Postmortem Changes %K RNA, Messenger %K Stochastic Processes %K Transcriptome %X

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.

%B Nat Commun %V 9 %P 490 %8 2018 02 13 %G eng %N 1 %1 https://www.ncbi.nlm.nih.gov/pubmed/29440659?dopt=Abstract %R 10.1038/s41467-017-02772-x %0 Journal Article %J Cancer Res %D 2018 %T Gene Expression Integration into Pathway Modules Reveals a Pan-Cancer Metabolic Landscape. %A Cubuk, Cankut %A Hidalgo, Marta R %A Amadoz, Alicia %A Pujana, Miguel A %A Mateo, Francesca %A Herranz, Carmen %A Carbonell-Caballero, José %A Dopazo, Joaquin %K Cell Line, Tumor %K Cluster Analysis %K Disease Progression %K Gene Expression Profiling %K Gene Expression Regulation, Neoplastic %K Gene Regulatory Networks %K Humans %K Kaplan-Meier Estimate %K Metabolome %K mutation %K Neoplasms %K Oncogenes %K Phenotype %K Prognosis %K RNA, Small Interfering %K Sequence Analysis, RNA %K Transcriptome %K Treatment Outcome %X

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

%B Cancer Res %V 78 %P 6059-6072 %8 2018 11 01 %G eng %N 21 %1 https://www.ncbi.nlm.nih.gov/pubmed/30135189?dopt=Abstract %R 10.1158/0008-5472.CAN-17-2705 %0 Journal Article %J Biol Direct %D 2018 %T Models of cell signaling uncover molecular mechanisms of high-risk neuroblastoma and predict disease outcome. %A Hidalgo, Marta R %A Amadoz, Alicia %A Cubuk, Cankut %A Carbonell-Caballero, José %A Dopazo, Joaquin %K Computational Biology %K Gene Expression Regulation, Neoplastic %K Humans %K JNK Mitogen-Activated Protein Kinases %K Models, Theoretical %K Neuroblastoma %K Signal Transduction %X

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' comments section.

%B Biol Direct %V 13 %P 16 %8 2018 08 22 %G eng %N 1 %1 https://www.ncbi.nlm.nih.gov/pubmed/30134948?dopt=Abstract %R 10.1186/s13062-018-0219-4 %0 Journal Article %J Oncotarget %D 2017 %T High throughput estimation of functional cell activities reveals disease mechanisms and predicts relevant clinical outcomes. %A Hidalgo, Marta R %A Cubuk, Cankut %A Amadoz, Alicia %A Salavert, Francisco %A Carbonell-Caballero, José %A Dopazo, Joaquin %K Computational Biology %K gene expression %K Gene Regulatory Networks %K Humans %K mutation %K Neoplasms %K Precision Medicine %K Sequence Analysis, RNA %K Signal Transduction %X

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.

%B Oncotarget %V 8 %P 5160-5178 %8 2017 Jan 17 %G eng %N 3 %1 https://www.ncbi.nlm.nih.gov/pubmed/28042959?dopt=Abstract %R 10.18632/oncotarget.14107 %0 Journal Article %J Bioinformatics %D 2017 %T Reference genome assessment from a population scale perspective: an accurate profile of variability and noise. %A Carbonell-Caballero, José %A Amadoz, Alicia %A Alonso, Roberto %A Hidalgo, Marta R %A Cubuk, Cankut %A Conesa, David %A López-Quílez, Antonio %A Dopazo, Joaquin %K Animals %K Genetic Variation %K Genome %K Genomics %K Genotype %K Humans %K Models, Statistical %K Quality Control %K Reproducibility of Results %K Software %X

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.

%B Bioinformatics %V 33 %P 3511-3517 %8 2017 Nov 15 %G eng %U https://academic.oup.com/bioinformatics/article-lookup/doi/10.1093/bioinformatics/btx482 %N 22 %1 https://www.ncbi.nlm.nih.gov/pubmed/28961772?dopt=Abstract %R 10.1093/bioinformatics/btx482 %0 Journal Article %J Nucleic acids research %D 2016 %T Actionable pathways: interactive discovery of therapeutic targets using signaling pathway models. %A Salavert, Francisco %A Hidago, Marta R %A Amadoz, Alicia %A Cubuk, Cankut %A Medina, Ignacio %A Crespo, Daniel %A Carbonell-Caballero, José %A Joaquín Dopazo %K actionable genes %K Disease mechanism %K drug action mechanism %K Drug discovery %K pathway analysis %K personalized medicine %K signalling %K therapeutic targets %X 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. %B Nucleic acids research %8 2016 May 2 %G eng %U http://nar.oxfordjournals.org/content/early/2016/05/02/nar.gkw369.full %R 10.1093/nar/gkw369 %0 Journal Article %J Stress (Amsterdam, Netherlands) %D 2016 %T Chronic subordination stress selectively downregulates the insulin signaling pathway in liver and skeletal muscle but not in adipose tissue of male mice. %A Sanghez, Valentina %A Cubuk, Cankut %A Sebastián-Leon, Patricia %A Carobbio, Stefania %A Dopazo, Joaquin %A Vidal-Puig, Antonio %A Bartolomucci, Alessandro %K Adipose tissue %K insulin %K IRS1 %K IRS2 %K metabolic syndrome %K obesity %K pathway analysis %X Chronic stress has been associated with obesity, glucose intolerance, and insulin resistance. We developed a model of chronic psychosocial stress (CPS) in which subordinate mice are vulnerable to obesity and the metabolic-like syndrome while dominant mice exhibit a healthy metabolic phenotype. Here we tested the hypothesis that the metabolic difference between subordinate and dominant mice is associated with changes in functional pathways relevant for insulin sensitivity, glucose and lipid homeostasis. Male mice were exposed to CPS for four weeks and fed either a standard diet or a high-fat diet (HFD). We first measured, by real-time PCR candidate genes, in the liver, skeletal muscle, and the perigonadal white adipose tissue (pWAT). Subsequently, we used a probabilistic analysis approach to analyze different ways in which signals can be transmitted across the pathways in each tissue. Results showed that subordinate mice displayed a drastic downregulation of the insulin pathway in liver and muscle, indicative of insulin resistance, already on standard diet. Conversely, pWAT showed molecular changes suggestive of facilitated fat deposition in an otherwise insulin-sensitive tissue. The molecular changes in subordinate mice fed a standard diet were greater compared to HFD-fed controls. Finally, dominant mice maintained a substantially normal metabolic and molecular phenotype even when fed a HFD. Overall, our data demonstrate that subordination stress is a potent stimulus for the downregulation of the insulin signaling pathway in liver and muscle and a major risk factor for the development of obesity, insulin resistance, and type 2 diabetes mellitus. %B Stress (Amsterdam, Netherlands) %P 1-11 %8 2016 Mar 7 %G eng %U http://www.tandfonline.com/doi/abs/10.3109/10253890.2016.1151491?journalCode=ists20 %R 10.3109/10253890.2016.1151491 %0 Journal Article %J Mol Metab %D 2016 %T Stress-induced activation of brown adipose tissue prevents obesity in conditions of low adaptive thermogenesis. %A Razzoli, Maria %A Frontini, Andrea %A Gurney, Allison %A Mondini, Eleonora %A Cubuk, Cankut %A Katz, Liora S %A Cero, Cheryl %A Bolan, Patrick J %A Dopazo, Joaquin %A Vidal-Puig, Antonio %A Cinti, Saverio %A Bartolomucci, Alessandro %X

BACKGROUND: Stress-associated conditions such as psychoemotional reactivity and depression have been paradoxically linked to either weight gain or weight loss. This bi-directional effect of stress is not understood at the functional level. Here we tested the hypothesis that pre-stress level of adaptive thermogenesis and brown adipose tissue (BAT) functions explain the vulnerability or resilience to stress-induced obesity.

METHODS: We used wt and triple β1,β2,β3-Adrenergic Receptors knockout (β-less) mice exposed to a model of chronic subordination stress (CSS) at either room temperature (22 °C) or murine thermoneutrality (30 °C). A combined behavioral, physiological, molecular, and immunohistochemical analysis was conducted to determine stress-induced modulation of energy balance and BAT structure and function. Immortalized brown adipocytes were used for in vitro assays.

RESULTS: Departing from our initial observation that βARs are dispensable for cold-induced BAT browning, we demonstrated that under physiological conditions promoting low adaptive thermogenesis and BAT activity (e.g. thermoneutrality or genetic deletion of the βARs), exposure to CSS acted as a stimulus for BAT activation and thermogenesis, resulting in resistance to diet-induced obesity despite the presence of hyperphagia. Conversely, in wt mice acclimatized to room temperature, and therefore characterized by sustained BAT function, exposure to CSS increased vulnerability to obesity. Exposure to CSS enhanced the sympathetic innervation of BAT in wt acclimatized to thermoneutrality and in β-less mice. Despite increased sympathetic innervation suggesting adrenergic-mediated browning, norepinephrine did not promote browning in βARs knockout brown adipocytes, which led us to identify an alternative sympathetic/brown adipocytes purinergic pathway in the BAT. This pathway is downregulated under conditions of low adaptive thermogenesis requirements, is induced by stress, and elicits activation of UCP1 in wt and β-less brown adipocytes. Importantly, this purinergic pathway is conserved in human BAT.

CONCLUSION: Our findings demonstrate that thermogenesis and BAT function are determinant of the resilience or vulnerability to stress-induced obesity. Our data support a model in which adrenergic and purinergic pathways exert complementary/synergistic functions in BAT, thus suggesting an alternative to βARs agonists for the activation of human BAT.

%B Mol Metab %V 5 %P 19-33 %8 2016 Jan %G eng %N 1 %1 https://www.ncbi.nlm.nih.gov/pubmed/26844204?dopt=Abstract %R 10.1016/j.molmet.2015.10.005 %0 Journal Article %J Nucleic acids research %D 2015 %T Babelomics 5.0: functional interpretation for new generations of genomic data. %A Alonso, Roberto %A Salavert, Francisco %A Garcia-Garcia, Francisco %A Carbonell-Caballero, José %A Bleda, Marta %A García-Alonso, Luz %A Sanchis-Juan, Alba %A Perez-Gil, Daniel %A Marin-Garcia, Pablo %A Sánchez, Rubén %A Cubuk, Cankut %A Hidalgo, Marta R %A Amadoz, Alicia %A Hernansaiz-Ballesteros, Rosa D %A Alemán, Alejandro %A Tárraga, Joaquín %A Montaner, David %A Medina, Ignacio %A Dopazo, Joaquin %K babelomics %K data integration %K gene set analysis %K interactome %K network analysis %K NGS %K RNA-seq %K Systems biology %K transcriptomics %X 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. %B Nucleic acids research %V 43 %P W117-W121 %8 2015 Apr 20 %G eng %U http://nar.oxfordjournals.org/content/43/W1/W117 %R 10.1093/nar/gkv384