@article {810, title = {drexml: A command line tool and Python package for drug repurposing.}, journal = {Comput Struct Biotechnol J}, volume = {23}, year = {2024}, month = {2024 Dec}, pages = {1129-1143}, abstract = {

We introduce drexml, a command line tool and Python package for rational data-driven drug repurposing. The package employs machine learning and mechanistic signal transduction modeling to identify drug targets capable of regulating a particular disease. In addition, it employs explainability tools to contextualize potential drug targets within the functional landscape of the disease. The methodology is validated in Fanconi Anemia and Familial Melanoma, two distinct rare diseases where there is a pressing need for solutions. In the Fanconi Anemia case, the model successfully predicts previously validated repurposed drugs, while in the Familial Melanoma case, it identifies a promising set of drugs for further investigation.

}, issn = {2001-0370}, doi = {10.1016/j.csbj.2024.02.027}, author = {Esteban-Medina, Marina and de la Oliva Roque, V{\'\i}ctor Manuel and Herr{\'a}iz-Gil, Sara and Pe{\~n}a-Chilet, Maria and Dopazo, Joaquin and Loucera, Carlos} } @article {805, title = {Drug-target identification in COVID-19 disease mechanisms using computational systems biology approaches.}, journal = {Front Immunol}, volume = {14}, year = {2024}, month = {2023}, pages = {1282859}, abstract = {

INTRODUCTION: The COVID-19 Disease Map project is a large-scale community effort uniting 277 scientists from 130 Institutions around the globe. We use high-quality, mechanistic content describing SARS-CoV-2-host interactions and develop interoperable bioinformatic pipelines for novel target identification and drug repurposing.

METHODS: Extensive community work allowed an impressive step forward in building interfaces between Systems Biology tools and platforms. Our framework can link biomolecules from omics data analysis and computational modelling to dysregulated pathways in a cell-, tissue- or patient-specific manner. Drug repurposing using text mining and AI-assisted analysis identified potential drugs, chemicals and microRNAs that could target the identified key factors.

RESULTS: Results revealed drugs already tested for anti-COVID-19 efficacy, providing a mechanistic context for their mode of action, and drugs already in clinical trials for treating other diseases, never tested against COVID-19.

DISCUSSION: The key advance is that the proposed framework is versatile and expandable, offering a significant upgrade in the arsenal for virus-host interactions and other complex pathologies.

}, keywords = {Computer Simulation, COVID-19, drug repositioning, Humans, SARS-CoV-2, Systems biology}, issn = {1664-3224}, doi = {10.3389/fimmu.2023.1282859}, author = {Niarakis, Anna and Ostaszewski, Marek and Mazein, Alexander and Kuperstein, Inna and Kutmon, Martina and Gillespie, Marc E and Funahashi, Akira and Acencio, Marcio Luis and Hemedan, Ahmed and Aichem, Michael and Klein, Karsten and Czauderna, Tobias and Burtscher, Felicia and Yamada, Takahiro G and Hiki, Yusuke and Hiroi, Noriko F and Hu, Finterly and Pham, Nhung and Ehrhart, Friederike and Willighagen, Egon L and Valdeolivas, Alberto and Dugourd, Aur{\'e}lien and Messina, Francesco and Esteban-Medina, Marina and Pe{\~n}a-Chilet, Maria and Rian, Kinza and Soliman, Sylvain and Aghamiri, Sara Sadat and Puniya, Bhanwar Lal and Naldi, Aur{\'e}lien and Helikar, Tom{\'a}{\v s} and Singh, Vidisha and Fern{\'a}ndez, Marco Fari{\~n}as and Bermudez, Viviam and Tsirvouli, Eirini and Montagud, Arnau and No{\"e}l, Vincent and Ponce-de-Leon, Miguel and Maier, Dieter and Bauch, Angela and Gyori, Benjamin M and Bachman, John A and Luna, Augustin and Pi{\~n}ero, Janet and Furlong, Laura I and Balaur, Irina and Rougny, Adrien and Jarosz, Yohan and Overall, Rupert W and Phair, Robert and Perfetto, Livia and Matthews, Lisa and Rex, Devasahayam Arokia Balaya and Orlic-Milacic, Marija and Gomez, Luis Cristobal Monraz and De Meulder, Bertrand and Ravel, Jean Marie and Jassal, Bijay and Satagopam, Venkata and Wu, Guanming and Golebiewski, Martin and Gawron, Piotr and Calzone, Laurence and Beckmann, Jacques S and Evelo, Chris T and D{\textquoteright}Eustachio, Peter and Schreiber, Falk and Saez-Rodriguez, Julio and Dopazo, Joaquin and Kuiper, Martin and Valencia, Alfonso and Wolkenhauer, Olaf and Kitano, Hiroaki and Barillot, Emmanuel and Auffray, Charles and Balling, Rudi and Schneider, Reinhard} } @article {796, title = {Defective extracellular matrix remodeling in brown adipose tissue is associated with fibro-inflammation and reduced diet-induced thermogenesis.}, journal = {Cell Rep}, volume = {42}, year = {2023}, month = {2023 Jun 13}, pages = {112640}, abstract = {

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.

}, issn = {2211-1247}, doi = {10.1016/j.celrep.2023.112640}, author = {Pellegrinelli, Vanessa and Figueroa-Ju{\'a}rez, Elizabeth and Samuelson, Isabella and U-Din, Mueez and Rodriguez-Fdez, Sonia and Virtue, Samuel and Leggat, Jennifer and Cubuk, Cankut and Peirce, Vivian J and Niemi, Tarja and Campbell, Mark and Rodriguez-Cuenca, Sergio and Dopazo, Joaquin and Carobbio, Stefania and Virtanen, Kirsi A and Vidal-Puig, Antonio} } @article {766, title = {Detection of High Level of Co-Infection and the Emergence of Novel SARS CoV-2 Delta-Omicron and Omicron-Omicron Recombinants in the Epidemiological Surveillance of Andalusia.}, journal = {Int J Mol Sci}, volume = {24}, year = {2023}, month = {2023 Jan 26}, abstract = {

Recombination is an evolutionary strategy to quickly acquire new viral properties inherited from the parental lineages. The systematic survey of the SARS-CoV-2 genome sequences of the Andalusian genomic surveillance strategy has allowed the detection of an unexpectedly high number of co-infections, which constitute the ideal scenario for the emergence of new recombinants. Whole genome sequence of SARS-CoV-2 has been carried out as part of the genomic surveillance programme. Sample sources included the main hospitals in the Andalusia region. In addition to the increase of co-infections and known recombinants, three novel SARS-CoV-2 delta-omicron and omicron-omicron recombinant variants with two break points have been detected. Our observations document an epidemiological scenario in which co-infection and recombination are detected more frequently. Finally, we describe a family case in which co-infection is followed by the detection of a recombinant made from the two co-infecting variants. This increased number of recombinants raises the risk of emergence of recombinant variants with increased transmissibility and pathogenicity.

}, issn = {1422-0067}, doi = {10.3390/ijms24032419}, author = {Perez-Florido, Javier and Casimiro-Soriguer, Carlos S and Ortuno, Francisco and Fernandez-Rueda, Jose L and Aguado, Andrea and Lara, Mar{\'\i}a and Riazzo, Cristina and Rodriguez-Iglesias, Manuel A and Camacho-Martinez, Pedro and Merino-Diaz, Laura and Pupo-Ledo, Inmaculada and de Salazar, Adolfo and Vi{\~n}uela, Laura and Fuentes, Ana and Chueca, Natalia and Garc{\'\i}a, Federico and Dopazo, Joaquin and Lepe, Jose A} } @article {751, title = {Discovering potential interactions between rare diseases and COVID-19 by combining mechanistic models of viral infection with statistical modeling.}, journal = {Hum Mol Genet}, year = {2022}, month = {2022 Jan 12}, abstract = {

Recent studies have demonstrated a relevant role of the host genetics in the COVID-19 prognosis. Most of the 7000 rare diseases described to date have a genetic component, typically highly penetrant. However, this vast spectrum of genetic variability remains yet unexplored with respect to possible interactions with COVID-19. Here, a mathematical mechanistic model of the COVID-19 molecular disease mechanism has been used to detect potential interactions between rare disease genes and the COVID-19 infection process and downstream consequences. Out of the 2518 disease genes analyzed, causative of 3854 rare diseases, a total of 254 genes have a direct effect on the COVID-19 molecular disease mechanism and 207 have an indirect effect revealed by a significant strong correlation. This remarkable potential of interaction occurs for more than 300 rare diseases. Mechanistic modeling of COVID-19 disease map has allowed a holistic systematic analysis of the potential interactions between the loss of function in known rare disease genes and the pathological consequences of COVID-19 infection. The results identify links between disease genes and COVID-19 hallmarks and demonstrate the usefulness of the proposed approach for future preventive measures in some rare diseases.

}, issn = {1460-2083}, doi = {10.1093/hmg/ddac007}, author = {L{\'o}pez-S{\'a}nchez, Macarena and Loucera, Carlos and Pe{\~n}a-Chilet, Maria and Dopazo, Joaquin} } @article {743, title = {Deciphering Genomic Heterogeneity and the Internal Composition of Tumour Activities through a Hierarchical Factorisation Model}, journal = {Mathematics}, volume = {9}, year = {2021}, month = {Jan-11-2021}, pages = {2833}, doi = {10.3390/math9212833}, url = {https://www.mdpi.com/2227-7390/9/21/2833https://www.mdpi.com/2227-7390/9/21/2833/pdf}, author = {Carbonell-Caballero, Jos{\'e} and L{\'o}pez-Qu{\'\i}lez, Antonio and Conesa, David and Dopazo, Joaquin} } @article {720, title = {A DNA damage repair gene-associated signature predicts responses of patients with advanced soft-tissue sarcoma to treatment with trabectedin.}, journal = {Mol Oncol}, volume = {15}, year = {2021}, month = {2021 12}, pages = {3691-3705}, abstract = {

Predictive biomarkers of trabectedin represent an unmet need in advanced soft-tissue sarcomas (STS). DNA damage repair (DDR) genes, involved in homologous recombination or nucleotide excision repair, had been previously described as biomarkers of trabectedin resistance or sensitivity, respectively. The majority of these studies only focused on specific factors (ERCC1, ERCC5, and BRCA1) and did not evaluate several other DDR-related genes that could have a relevant role for trabectedin efficacy. In this retrospective translational study, 118 genes involved in DDR were evaluated to determine, by transcriptomics, a predictive gene signature of trabectedin efficacy. A six-gene predictive signature of trabectedin efficacy was built in a series of 139 tumor samples from patients with advanced STS. Patients in the high-risk gene signature group showed a significantly worse progression-free survival compared with patients in the low-risk group (2.1 vs 6.0 months, respectively). Differential gene expression analysis defined new potential predictive biomarkers of trabectedin sensitivity (PARP3 and CCNH) or resistance (DNAJB11 and PARP1). Our study identified a new gene signature that significantly predicts patients with higher probability to respond to treatment with trabectedin. Targeting some genes of this signature emerges as a potential strategy to enhance trabectedin efficacy.

}, issn = {1878-0261}, doi = {10.1002/1878-0261.12996}, author = {Moura, David S and Pe{\~n}a-Chilet, Maria and Cordero Varela, Juan Antonio and Alvarez-Alegret, Ramiro and Agra-Pujol, Carolina and Izquierdo, Francisco and Ramos, Rafael and Ortega-Medina, Luis and Martin-Davila, Francisco and Castilla-Ramirez, Carolina and Hernandez-Leon, Carmen Nieves and Romagosa, Cleofe and Vaz Salgado, Maria Angeles and Lavernia, Javier and Bagu{\'e}, Silvia and Mayodormo-Aranda, Empar and Vicioso, Luis and Hern{\'a}ndez Barcel{\'o}, Jose Emilio and Rubio-Casadevall, Jordi and de Juan, Ana and Fia{\~n}o-Valverde, Maria Concepcion and Hindi, Nadia and Lopez-Alvarez, Maria and Lacerenza, Serena and Dopazo, Joaquin and Gutierrez, Antonio and Alvarez, Rosa and Valverde, Claudia and Martinez-Trufero, Javier and Martin-Broto, Javier} } @article {713, title = {Drug repurposing for COVID-19 using machine learning and mechanistic models of signal transduction circuits related to SARS-CoV-2 infection.}, journal = {Signal Transduct Target Ther}, volume = {5}, year = {2020}, month = {2020 12 11}, pages = {290}, keywords = {Computational Chemistry, COVID-19, drug repositioning, Humans, Machine Learning, Molecular Docking Simulation, Molecular Targeted Therapy, Proteins, SARS-CoV-2, Signal Transduction}, issn = {2059-3635}, doi = {10.1038/s41392-020-00417-y}, author = {Loucera, Carlos and Esteban-Medina, Marina and Rian, Kinza and Falco, Matias M and Dopazo, Joaquin and Pe{\~n}a-Chilet, Maria} } @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 {436, title = {Dysfunctional mitochondrial fission impairs cell reprogramming.}, journal = {Cell Cycle}, volume = {15}, year = {2016}, month = {2016 Dec}, pages = {3240-3250}, abstract = {

We have recently shown that mitochondrial fission is induced early in reprogramming in a Drp1-dependent manner; however, the identity of the factors controlling Drp1 recruitment to mitochondria was unexplored. To investigate this, we used a panel of RNAi targeting factors involved in the regulation of mitochondrial dynamics and we observed that MiD51, Gdap1 and, to a lesser extent, Mff were found to play key roles in this process. Cells derived from Gdap1-null mice were used to further explore the role of this factor in cell reprogramming. Microarray data revealed a prominent down-regulation of cell cycle pathways in Gdap1-null cells early in reprogramming and cell cycle profiling uncovered a G2/M growth arrest in Gdap1-null cells undergoing reprogramming. High-Content analysis showed that this growth arrest was DNA damage-independent. We propose that lack of efficient mitochondrial fission impairs cell reprogramming by interfering with cell cycle progression in a DNA damage-independent manner.

}, keywords = {Animals, Cell Cycle Checkpoints, Cellular Reprogramming, DNA Damage, G2 Phase, Gene Knockdown Techniques, Mice, Mitochondrial Dynamics, Mitosis, Nerve Tissue Proteins, Pluripotent Stem Cells, Transcription Factors}, issn = {1551-4005}, doi = {10.1080/15384101.2016.1241930}, author = {Prieto, Javier and Le{\'o}n, Marian and Ponsoda, Xavier and Garcia-Garcia, Francisco and Bort, Roque and Serna, Eva and Barneo-Mu{\~n}oz, Manuela and Palau, Francesc and Dopazo, Joaquin and L{\'o}pez-Garc{\'\i}a, Carlos and Torres, Josema} } @article {486, title = {Deciphering intrafamilial phenotypic variability by exome sequencing in a Bardet-Biedl family.}, journal = {Mol Genet Genomic Med}, volume = {2}, year = {2014}, month = {2014 Mar}, pages = {124-33}, abstract = {

Bardet-Biedl syndrome (BBS) is a model ciliopathy characterized by a wide range of clinical variability. The heterogeneity of this condition is reflected in the number of underlying gene defects and the epistatic interactions between the proteins encoded. BBS is generally inherited in an autosomal recessive trait. However, in some families, mutations across different loci interact to modulate the expressivity of the phenotype. In order to investigate the magnitude of epistasis in one BBS family with remarkable intrafamilial phenotypic variability, we designed an exome sequencing-based approach using SOLID 5500xl platform. This strategy allowed the reliable detection of the primary causal mutations in our family consisting of two novel compound heterozygous mutations in McKusick-Kaufman syndrome (MKKS) gene (p.D90G and p.V396F). Additionally, exome sequencing enabled the detection of one novel heterozygous NPHP4 variant which is predicted to activate a cryptic acceptor splice site and is only present in the most severely affected patient. Here, we provide an exome sequencing analysis of a BBS family and show the potential utility of this tool, in combination with network analysis, to detect disease-causing mutations and second-site modifiers. Our data demonstrate how next-generation sequencing (NGS) can facilitate the dissection of epistatic phenomena, and shed light on the genetic basis of phenotypic variability.

}, issn = {2324-9269}, doi = {10.1002/mgg3.50}, author = {Gonz{\'a}lez-del Pozo, Mar{\'\i}a and M{\'e}ndez-Vidal, Cristina and Santoyo-L{\'o}pez, Javier and Vela-Boza, Alicia and Bravo-Gil, Nereida and Rueda, Antonio and Garc{\'\i}a-Alonso, Luz and V{\'a}zquez-Marouschek, Carmen and Dopazo, Joaquin and Borrego, Salud and Anti{\v n}olo, Guillermo} } @article {503, title = {Defining the genomic signature of totipotency and pluripotency during early human development.}, journal = {PLoS One}, volume = {8}, year = {2013}, month = {2013}, pages = {e62135}, abstract = {

The genetic mechanisms governing human pre-implantation embryo development and the in vitro counterparts, human embryonic stem cells (hESCs), still remain incomplete. Previous global genome studies demonstrated that totipotent blastomeres from day-3 human embryos and pluripotent inner cell masses (ICMs) from blastocysts, display unique and differing transcriptomes. Nevertheless, comparative gene expression analysis has revealed that no significant differences exist between hESCs derived from blastomeres versus those obtained from ICMs, suggesting that pluripotent hESCs involve a new developmental progression. To understand early human stages evolution, we developed an undifferentiation network signature (UNS) and applied it to a differential gene expression profile between single blastomeres from day-3 embryos, ICMs and hESCs. This allowed us to establish a unique signature composed of highly interconnected genes characteristic of totipotency (61 genes), in vivo pluripotency (20 genes), and in vitro pluripotency (107 genes), and which are also proprietary according to functional analysis. This systems biology approach has led to an improved understanding of the molecular and signaling processes governing human pre-implantation embryo development, as well as enabling us to comprehend how hESCs might adapt to in vitro culture conditions.

}, keywords = {Blastocyst Inner Cell Mass, Blastomeres, Cell Differentiation, Embryonic Development, Embryonic Stem Cells, Gene Expression Profiling, Gene Regulatory Networks, Genome, Human, Humans, Molecular Sequence Annotation, Pluripotent Stem Cells, Totipotent Stem Cells}, issn = {1932-6203}, doi = {10.1371/journal.pone.0062135}, author = {Galan, Amparo and Diaz-Gimeno, Patricia and Poo, Maria Eugenia and Valbuena, Diana and Sanchez, Eva and Ruiz, Veronica and Dopazo, Joaquin and Montaner, David and Conesa, Ana and Simon, Carlos} } @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} } @article {513, title = {Diversification of the expanded teleost-specific toll-like receptor family in Atlantic cod, Gadus morhua.}, journal = {BMC Evol Biol}, volume = {12}, year = {2012}, month = {2012 Dec 29}, pages = {256}, abstract = {

BACKGROUND: Toll-like receptors (Tlrs) are major molecular pattern recognition receptors of the innate immune system. Atlantic cod (Gadus morhua) is the first vertebrate known to have lost most of the mammalian Tlr orthologues, particularly all bacterial recognising and other cell surface Tlrs. On the other hand, its genome encodes a unique repertoire of teleost-specific Tlrs. The aim of this study was to investigate if these duplicate Tlrs have been retained through adaptive evolution to compensate for the lack of other cell surface Tlrs in the cod genome.

RESULTS: In this study, one tlr21, 12 tlr22 and two tlr23 genes representing the teleost-specific Tlr family have been cloned and characterised in cod. Phylogenetic analysis grouped all tlr22 genes under a single clade, indicating that the multiple cod paralogues have arisen through lineage-specific duplications. All tlrs examined were transcribed in immune-related tissues as well as in stomach, gut and gonads of adult cod and were differentially expressed during early development. These tlrs were also differentially regulated following immune challenge by immersion with Vibrio anguillarum, indicating their role in the immune response. An increase in water temperature from 4 to 12{\textdegree}C was associated with a 5.5-fold down-regulation of tlr22d transcript levels in spleen. Maximum likelihood analysis with different evolution models revealed that tlr22 genes are under positive selection. A total of 24 codons were found to be positively selected, of which 19 are in the ligand binding region of ectodomain.

CONCLUSION: Positive selection pressure coupled with experimental evidence of differential expression strongly support the hypothesis that teleost-specific tlr paralogues in cod are undergoing neofunctionalisation and can recognise bacterial pathogen-associated molecular patterns to compensate for the lack of other cell surface Tlrs.

}, keywords = {Amino Acid Sequence, Animals, Binding Sites, Evolution, Molecular, Fish Diseases, Fish Proteins, Gadus morhua, Gene Expression Profiling, Genetic Variation, Gills, Head Kidney, Host-Pathogen Interactions, Models, Molecular, Molecular Sequence Data, Multigene Family, Phylogeny, Protein Structure, Tertiary, Reverse Transcriptase Polymerase Chain Reaction, Selection, Genetic, Sequence Analysis, DNA, Sequence Homology, Amino Acid, Temperature, Toll-Like Receptors, Vibrio}, issn = {1471-2148}, doi = {10.1186/1471-2148-12-256}, author = {Sundaram, Arvind Y M and Kiron, Viswanath and Dopazo, Joaquin and Fernandes, Jorge M O} } @article {529, title = {Differential expression in RNA-seq: a matter of depth.}, journal = {Genome Res}, volume = {21}, year = {2011}, month = {2011 Dec}, pages = {2213-23}, abstract = {

Next-generation sequencing (NGS) technologies are revolutionizing genome research, and in particular, their application to transcriptomics (RNA-seq) is increasingly being used for gene expression profiling as a replacement for microarrays. However, the properties of RNA-seq data have not been yet fully established, and additional research is needed for understanding how these data respond to differential expression analysis. In this work, we set out to gain insights into the characteristics of RNA-seq data analysis by studying an important parameter of this technology: the sequencing depth. We have analyzed how sequencing depth affects the detection of transcripts and their identification as differentially expressed, looking at aspects such as transcript biotype, length, expression level, and fold-change. We have evaluated different algorithms available for the analysis of RNA-seq and proposed a novel approach--NOISeq--that differs from existing methods in that it is data-adaptive and nonparametric. Our results reveal that most existing methodologies suffer from a strong dependency on sequencing depth for their differential expression calls and that this results in a considerable number of false positives that increases as the number of reads grows. In contrast, our proposed method models the noise distribution from the actual data, can therefore better adapt to the size of the data set, and is more effective in controlling the rate of false discoveries. This work discusses the true potential of RNA-seq for studying regulation at low expression ranges, the noise within RNA-seq data, and the issue of replication.

}, keywords = {Algorithms, Expressed Sequence Tags, Gene Expression Profiling, Gene Expression Regulation, Humans, Models, Genetic, Oligonucleotide Array Sequence Analysis}, issn = {1549-5469}, doi = {10.1101/gr.124321.111}, author = {Tarazona, Sonia and Garc{\'\i}a-Alcalde, Fernando and Dopazo, Joaquin and Ferrer, Alberto and Conesa, Ana} } @article {21266330, title = {Differential Lipid Partitioning Between Adipocytes and Tissue Macrophages Modulates Macrophage Lipotoxicity and M2/M1 Polarization in Obese Mice.}, journal = {Diabetes}, volume = {60}, number = {3}, year = {2011}, month = {2011 Jan 24}, pages = {797-809}, abstract = {

OBJECTIVE Obesity-associated insulin resistance is characterized by a state of chronic, low-grade inflammation that is associated with the accumulation of M1 proinflammatory macrophages in adipose tissue. Although different evidence explains the mechanisms linking the expansion of adipose tissue and adipose tissue macrophage (ATM) polarization, in the current study we investigated the concept of lipid-induced toxicity as the pathogenic link that could explain the trigger of this response. RESEARCH DESIGN AND METHODS We addressed this question using isolated ATMs and adipocytes from genetic and diet-induced murine models of obesity. Through transcriptomic and lipidomic analysis, we created a model integrating transcript and lipid species networks simultaneously occurring in adipocytes and ATMs and their reversibility by thiazolidinedione treatment. RESULTS We show that polarization of ATMs is associated with lipid accumulation and the consequent formation of foam cell-like cells in adipose tissue. Our study reveals that early stages of adipose tissue expansion are characterized by M2-polarized ATMs and that progressive lipid accumulation within ATMs heralds the M1 polarization, a macrophage phenotype associated with severe obesity and insulin resistance. Furthermore, rosiglitazone treatment, which promotes redistribution of lipids toward adipocytes and extends the M2 ATM polarization state, prevents the lipid alterations associated with M1 ATM polarization. CONCLUSIONS Our data indicate that the M1 ATM polarization in obesity might be a macrophage-specific manifestation of a more general lipotoxic pathogenic mechanism. This indicates that strategies to optimize fat deposition and repartitioning toward adipocytes might improve insulin sensitivity by preventing ATM lipotoxicity and M1 polarization.

}, author = {Prieur, Xavier and Mok, Crystal Y L and Velagapudi, Vidya R and N{\'u}{\~n}ez, Vanessa and Fuentes, Luc{\'\i}a and Montaner, David and Ishikawa, Ko and Camacho, Alberto and Barbarroja, Nuria and O{\textquoteright}Rahilly, Stephen and Sethi, Jaswinder and Dopazo, Joaquin and Oresic, Matej and Ricote, Mercedes and Vidal-Puig, Antonio} } @article {547, title = {DNA methylation epigenotypes in breast cancer molecular subtypes.}, journal = {Breast Cancer Res}, volume = {12}, year = {2010}, month = {2010}, pages = {R77}, abstract = {

INTRODUCTION: Identification of gene expression based breast cancer subtypes is considered as a critical means of prognostication. Genetic mutations along with epigenetic alterations contribute to gene expression changes occurring in breast cancer. So far, these epigenetic contributions to sporadic breast cancer subtypes have not been well characterized, and there is only a limited understanding of the epigenetic mechanisms affected in those particular breast cancer subtypes. The present study was undertaken to dissect the breast cancer methylome and deliver specific epigenotypes associated with particular breast cancer subtypes.

METHODS: Using a microarray approach we analyzed DNA methylation in regulatory regions of 806 cancer related genes in 28 breast cancer paired samples. We subsequently performed substantial technical and biological validation by Pyrosequencing, investigating the top qualifying 19 CpG regions in independent cohorts encompassing 47 basal-like, 44 ERBB2+ overexpressing, 48 luminal A and 48 luminal B paired breast cancer/adjacent tissues. Using all-subset selection method, we identified the most subtype predictive methylation profiles in multivariable logistic regression analysis.

RESULTS: The approach efficiently recognized 15 individual CpG loci differentially methylated in breast cancer tumor subtypes. We further identify novel subtype specific epigenotypes which clearly demonstrate the differences in the methylation profiles of basal-like and human epidermal growth factor 2 (HER2)-overexpressing tumors.

CONCLUSIONS: Our results provide evidence that well defined DNA methylation profiles enables breast cancer subtype prediction and support the utilization of this biomarker for prognostication and therapeutic stratification of patients with breast cancer.

}, keywords = {Aged, Breast Neoplasms, CpG Islands, DNA Methylation, Epigenesis, Genetic, Female, Gene Expression Profiling, Genes, p53, Genotype, Humans, Ki-67 Antigen, Middle Aged, mutation, Neoplasm Grading, Oligonucleotide Array Sequence Analysis, Receptor, ErbB-2, Tumor Suppressor Protein p53}, issn = {1465-542X}, doi = {10.1186/bcr2721}, author = {Bediaga, Naiara G and Acha-Sagredo, Amelia and Guerra, Isabel and Viguri, Amparo and Albaina, Carmen and Ruiz Diaz, Irune and Rezola, Ricardo and Alberdi, Maria Jesus and Dopazo, Joaquin and Montaner, David and Renobales, Mertxe and Fernandez, Agustin F and Field, John K and Fraga, Mario F and Liloglou, Triantafillos and de Pancorbo, Marian M} } @article {590, title = {Direct functional assessment of the composite phenotype through multivariate projection strategies.}, journal = {Genomics}, volume = {92}, year = {2008}, month = {2008 Dec}, pages = {373-83}, abstract = {

We present a novel approach for the analysis of transcriptomics data that integrates functional annotation of gene sets with expression values in a multivariate fashion, and directly assesses the relation of functional features to a multivariate space of response phenotypical variables. Multivariate projection methods are used to obtain new correlated variables for a set of genes that share a given function. These new functional variables are then related to the response variables of interest. The analysis of the principal directions of the multivariate regression allows for the identification of gene function features correlated with the phenotype. Two different transcriptomics studies are used to illustrate the statistical and interpretative aspects of the methodology. We demonstrate the superiority of the proposed method over equivalent approaches.

}, keywords = {Breast Neoplasms, Computational Biology, Databases, Genetic, Female, Gene Expression Profiling, Humans, Mathematical Computing, Multivariate Analysis, Phenotype}, issn = {1089-8646}, doi = {10.1016/j.ygeno.2008.05.015}, author = {Conesa, Ana and Bro, Rasmus and Garcia-Garcia, Francisco and Prats, Jos{\'e} Manuel and G{\"o}tz, Stefan and Kjeldahl, Karin and Montaner, David and Dopazo, Joaquin} } @article {603, title = {DBAli tools: mining the protein structure space.}, journal = {Nucleic Acids Res}, volume = {35}, year = {2007}, month = {2007 Jul}, pages = {W393-7}, abstract = {

The DBAli tools use a comprehensive set of structural alignments in the DBAli database to leverage the structural information deposited in the Protein Data Bank (PDB). These tools include (i) the DBAlit program that allows users to input the 3D coordinates of a protein structure for comparison by MAMMOTH against all chains in the PDB; (ii) the AnnoLite and AnnoLyze programs that annotate a target structure based on its stored relationships to other structures; (iii) the ModClus program that clusters structures by sequence and structure similarities; (iv) the ModDom program that identifies domains as recurrent structural fragments and (v) an implementation of the COMPARER method in the SALIGN command in MODELLER that creates a multiple structure alignment for a set of related protein structures. Thus, the DBAli tools, which are freely accessible via the World Wide Web at http://salilab.org/DBAli/, allow users to mine the protein structure space by establishing relationships between protein structures and their functions.

}, keywords = {Algorithms, Amino Acid Sequence, Computational Biology, Data Interpretation, Statistical, Databases, Protein, Internet, Molecular Sequence Data, Protein Conformation, Proteins, Pseudomonas aeruginosa, Sequence Alignment, Sequence Analysis, Protein, Sequence Homology, Amino Acid, Software, Structure-Activity Relationship}, issn = {1362-4962}, doi = {10.1093/nar/gkm236}, author = {Marti-Renom, Marc A and Pieper, Ursula and Madhusudhan, M S and Rossi, Andrea and Eswar, Narayanan and Davis, Fred P and Al-Shahrour, F{\'a}tima and Dopazo, Joaquin and Sali, Andrej} }