1
|
Chen CK. Inference of gene networks using gene expression data with applications. Heliyon 2024; 10:e26065. [PMID: 38449656 PMCID: PMC10915353 DOI: 10.1016/j.heliyon.2024.e26065] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2024] [Revised: 02/06/2024] [Accepted: 02/07/2024] [Indexed: 03/08/2024] Open
Abstract
Gene networks (GNs) use graphs to represent the interaction relationships between genes. Large-scale GNs are often sparse and contain hub genes that interact with many other genes. In this paper, we propose a novel method called NetARD, which utilizes Automatic Relevance Determination (ARD) to estimate partial correlations, to infer GNs with the hub genes from gene expression data. We test NetARD on simulated GNs and in silico GNs, and it outperforms existing methods. In our high-throughput gene expression data analysis, we integrate the NetARD into a method called GN Co-expression Extension (GNCE). This approach infers the GNs of co-expressed genes, with genes from a predefined GN serving as hub genes. We validate this approach by extending the core GN of transcription factor genes of E. coli using microarray data. In an application example, we identify biological process (BP) Gene Ontology (GO) terms that are significantly involved in cancer progression. This task is accomplished by analyzing the GN inferred through GNCE using the core GN associated with the colorectal cancer pathway and RNA-seq data.
Collapse
Affiliation(s)
- Chi-Kan Chen
- Department of Applied Mathematics, National Chung Hsing University, 145 Xingda Rd., South Dist., Taichung City, 40227, Taiwan
| |
Collapse
|
2
|
Pasquier C, Guerlais V, Pallez D, Rapetti-Mauss R, Soriani O. A network embedding approach to identify active modules in biological interaction networks. Life Sci Alliance 2023; 6:e202201550. [PMID: 37339804 PMCID: PMC10282331 DOI: 10.26508/lsa.202201550] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2022] [Revised: 06/06/2023] [Accepted: 06/06/2023] [Indexed: 06/22/2023] Open
Abstract
The identification of condition-specific gene sets from transcriptomic experiments is important to reveal regulatory and signaling mechanisms associated with a given cellular response. Statistical methods of differential expression analysis, designed to assess individual gene variations, have trouble highlighting modules of small varying genes whose interaction is essential to characterize phenotypic changes. To identify these highly informative gene modules, several methods have been proposed in recent years, but they have many limitations that make them of little use to biologists. Here, we propose an efficient method for identifying these active modules that operates on a data embedding combining gene expressions and interaction data. Applications carried out on real datasets show that our method can identify new groups of genes of high interest corresponding to functions not revealed by traditional approaches. Software is available at https://github.com/claudepasquier/amine.
Collapse
Affiliation(s)
- Claude Pasquier
- Laboratoire d'Informatique, Signaux et Systèmes de Sophia-Antipolis, I3S - UMR7271 - UNS CNRS, Les Algorithmes - bât. Euclide B, Sophia Antipolis, France
| | - Vincent Guerlais
- Laboratoire d'Informatique, Signaux et Systèmes de Sophia-Antipolis, I3S - UMR7271 - UNS CNRS, Les Algorithmes - bât. Euclide B, Sophia Antipolis, France
| | - Denis Pallez
- Laboratoire d'Informatique, Signaux et Systèmes de Sophia-Antipolis, I3S - UMR7271 - UNS CNRS, Les Algorithmes - bât. Euclide B, Sophia Antipolis, France
| | - Raphaël Rapetti-Mauss
- iBV - Institut de Biologie Valrose, Université Nice Sophia Antipolis, Faculté des Sciences, Parc Valrose, Nice cedex 2, France
| | - Olivier Soriani
- iBV - Institut de Biologie Valrose, Université Nice Sophia Antipolis, Faculté des Sciences, Parc Valrose, Nice cedex 2, France
| |
Collapse
|
3
|
Corria-Osorio J, Carmona SJ, Stefanidis E, Andreatta M, Ortiz-Miranda Y, Muller T, Rota IA, Crespo I, Seijo B, Castro W, Jimenez-Luna C, Scarpellino L, Ronet C, Spill A, Lanitis E, Romero P, Luther SA, Irving M, Coukos G. Orthogonal cytokine engineering enables novel synthetic effector states escaping canonical exhaustion in tumor-rejecting CD8 + T cells. Nat Immunol 2023; 24:869-883. [PMID: 37081150 PMCID: PMC10154250 DOI: 10.1038/s41590-023-01477-2] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2022] [Accepted: 03/01/2023] [Indexed: 04/22/2023]
Abstract
To date, no immunotherapy approaches have managed to fully overcome T-cell exhaustion, which remains a mandatory fate for chronically activated effector cells and a major therapeutic challenge. Understanding how to reprogram CD8+ tumor-infiltrating lymphocytes away from exhausted effector states remains an elusive goal. Our work provides evidence that orthogonal gene engineering of T cells to secrete an interleukin (IL)-2 variant binding the IL-2Rβγ receptor and the alarmin IL-33 reprogrammed adoptively transferred T cells to acquire a novel, synthetic effector state, which deviated from canonical exhaustion and displayed superior effector functions. These cells successfully overcame homeostatic barriers in the host and led-in the absence of lymphodepletion or exogenous cytokine support-to high levels of engraftment and tumor regression. Our work unlocks a new opportunity of rationally engineering synthetic CD8+ T-cell states endowed with the ability to avoid exhaustion and control advanced solid tumors.
Collapse
Affiliation(s)
- Jesus Corria-Osorio
- Ludwig Institute for Cancer Research, Lausanne Branch, University of Lausanne; and Department of Oncology, Lausanne University Hospital, Epalinges, Switzerland.
- AGORA Cancer Research Center, Lausanne, Switzerland.
| | - Santiago J Carmona
- Ludwig Institute for Cancer Research, Lausanne Branch, University of Lausanne; and Department of Oncology, Lausanne University Hospital, Epalinges, Switzerland
- AGORA Cancer Research Center, Lausanne, Switzerland
| | - Evangelos Stefanidis
- Ludwig Institute for Cancer Research, Lausanne Branch, University of Lausanne; and Department of Oncology, Lausanne University Hospital, Epalinges, Switzerland
| | - Massimo Andreatta
- Ludwig Institute for Cancer Research, Lausanne Branch, University of Lausanne; and Department of Oncology, Lausanne University Hospital, Epalinges, Switzerland
- AGORA Cancer Research Center, Lausanne, Switzerland
| | - Yaquelin Ortiz-Miranda
- Ludwig Institute for Cancer Research, Lausanne Branch, University of Lausanne; and Department of Oncology, Lausanne University Hospital, Epalinges, Switzerland
- AGORA Cancer Research Center, Lausanne, Switzerland
| | - Tania Muller
- Ludwig Institute for Cancer Research, Lausanne Branch, University of Lausanne; and Department of Oncology, Lausanne University Hospital, Epalinges, Switzerland
- AGORA Cancer Research Center, Lausanne, Switzerland
| | - Ioanna A Rota
- Ludwig Institute for Cancer Research, Lausanne Branch, University of Lausanne; and Department of Oncology, Lausanne University Hospital, Epalinges, Switzerland
- AGORA Cancer Research Center, Lausanne, Switzerland
| | - Isaac Crespo
- Ludwig Institute for Cancer Research, Lausanne Branch, University of Lausanne; and Department of Oncology, Lausanne University Hospital, Epalinges, Switzerland
- AGORA Cancer Research Center, Lausanne, Switzerland
| | - Bili Seijo
- Ludwig Institute for Cancer Research, Lausanne Branch, University of Lausanne; and Department of Oncology, Lausanne University Hospital, Epalinges, Switzerland
- AGORA Cancer Research Center, Lausanne, Switzerland
| | - Wilson Castro
- Ludwig Institute for Cancer Research, Lausanne Branch, University of Lausanne; and Department of Oncology, Lausanne University Hospital, Epalinges, Switzerland
- AGORA Cancer Research Center, Lausanne, Switzerland
| | - Cristina Jimenez-Luna
- Ludwig Institute for Cancer Research, Lausanne Branch, University of Lausanne; and Department of Oncology, Lausanne University Hospital, Epalinges, Switzerland
| | | | - Catherine Ronet
- Ludwig Institute for Cancer Research, Lausanne Branch, University of Lausanne; and Department of Oncology, Lausanne University Hospital, Epalinges, Switzerland
- AGORA Cancer Research Center, Lausanne, Switzerland
| | - Aodrenn Spill
- Ludwig Institute for Cancer Research, Lausanne Branch, University of Lausanne; and Department of Oncology, Lausanne University Hospital, Epalinges, Switzerland
- AGORA Cancer Research Center, Lausanne, Switzerland
| | - Evripidis Lanitis
- Ludwig Institute for Cancer Research, Lausanne Branch, University of Lausanne; and Department of Oncology, Lausanne University Hospital, Epalinges, Switzerland
| | - Pedro Romero
- Ludwig Institute for Cancer Research, Lausanne Branch, University of Lausanne; and Department of Oncology, Lausanne University Hospital, Epalinges, Switzerland
- AGORA Cancer Research Center, Lausanne, Switzerland
| | - Sanjiv A Luther
- Department of Immunobiology, University of Lausanne, Epalinges, Switzerland
| | - Melita Irving
- Ludwig Institute for Cancer Research, Lausanne Branch, University of Lausanne; and Department of Oncology, Lausanne University Hospital, Epalinges, Switzerland
| | - George Coukos
- Ludwig Institute for Cancer Research, Lausanne Branch, University of Lausanne; and Department of Oncology, Lausanne University Hospital, Epalinges, Switzerland.
- AGORA Cancer Research Center, Lausanne, Switzerland.
| |
Collapse
|
4
|
Zhang L, Fan S, Vera J, Lai X. A network medicine approach for identifying diagnostic and prognostic biomarkers and exploring drug repurposing in human cancer. Comput Struct Biotechnol J 2022; 21:34-45. [PMID: 36514340 PMCID: PMC9732137 DOI: 10.1016/j.csbj.2022.11.037] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2022] [Revised: 11/18/2022] [Accepted: 11/18/2022] [Indexed: 12/03/2022] Open
Abstract
Cancer is a heterogeneous disease mainly driven by abnormal gene perturbations in regulatory networks. Therefore, it is appealing to identify the common and specific perturbed genes from multiple cancer networks. We developed an integrative network medicine approach to identify novel biomarkers and investigate drug repurposing across cancer types. We used a network-based method to prioritize genes in cancer-specific networks reconstructed using human transcriptome and interactome data. The prioritized genes show extensive perturbation and strong regulatory interaction with other highly perturbed genes, suggesting their vital contribution to tumorigenesis and tumor progression, and are therefore regarded as cancer genes. The cancer genes detected show remarkable performances in discriminating tumors from normal tissues and predicting survival times of cancer patients. Finally, we developed a network proximity approach to systematically screen drugs and identified dozens of candidates with repurposable potential in several cancer types. Taken together, we demonstrated the power of the network medicine approach to identify novel biomarkers and repurposable drugs in multiple cancer types. We have also made the data and code freely accessible to ensure reproducibility and reusability of the developed computational workflow.
Collapse
Affiliation(s)
- Le Zhang
- College of Computer Science, Sichuan University, Chengdu, China
| | - Shiwei Fan
- College of Computer Science, Sichuan University, Chengdu, China
| | - Julio Vera
- Laboratory of Systems Tumor Immunology, Department of Dermatology, Universitätsklinikum Erlangen and Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen, Germany,Deutsches Zentrum Immuntherapie, Erlangen, Germany,Comprehensive Cancer Center Erlangen, Erlangen, Germany
| | - Xin Lai
- Laboratory of Systems Tumor Immunology, Department of Dermatology, Universitätsklinikum Erlangen and Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen, Germany,Deutsches Zentrum Immuntherapie, Erlangen, Germany,Comprehensive Cancer Center Erlangen, Erlangen, Germany,BioMediTech, Faculty of Medicine and Health Technology, Tampere University, Tampere, Finland,Corresponding author at: Universitätsklinikum Erlangen, Erlangen, Germany; Tampere University, Tampere, Finland.
| |
Collapse
|
5
|
Sahoo A, Pechmann S. Functional network motifs defined through integration of protein-protein and genetic interactions. PeerJ 2022; 10:e13016. [PMID: 35223214 PMCID: PMC8877332 DOI: 10.7717/peerj.13016] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2021] [Accepted: 02/06/2022] [Indexed: 01/11/2023] Open
Abstract
Cells are enticingly complex systems. The identification of feedback regulation is critically important for understanding this complexity. Network motifs defined as small graphlets that occur more frequently than expected by chance have revolutionized our understanding of feedback circuits in cellular networks. However, with their definition solely based on statistical over-representation, network motifs often lack biological context, which limits their usefulness. Here, we define functional network motifs (FNMs) through the systematic integration of genetic interaction data that directly inform on functional relationships between genes and encoded proteins. Occurring two orders of magnitude less frequently than conventional network motifs, we found FNMs significantly enriched in genes known to be functionally related. Moreover, our comprehensive analyses of FNMs in yeast showed that they are powerful at capturing both known and putative novel regulatory interactions, thus suggesting a promising strategy towards the systematic identification of feedback regulation in biological networks. Many FNMs appeared as excellent candidates for the prioritization of follow-up biochemical characterization, which is a recurring bottleneck in the targeting of complex diseases. More generally, our work highlights a fruitful avenue for integrating and harnessing genomic network data.
Collapse
Affiliation(s)
- Amruta Sahoo
- Département de Biochimie, Université de Montréal, Montréal, QC, Canada
| | | |
Collapse
|
6
|
Chen CK. Inference of genetic regulatory networks with regulatory hubs using vector autoregressions and automatic relevance determination with model selections. Stat Appl Genet Mol Biol 2021; 20:121-143. [PMID: 34963205 DOI: 10.1515/sagmb-2020-0054] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2020] [Accepted: 11/15/2021] [Indexed: 12/11/2022]
Abstract
The inference of genetic regulatory networks (GRNs) reveals how genes interact with each other. A few genes can regulate many genes as targets to control cell functions. We present new methods based on the order-1 vector autoregression (VAR1) for inferring GRNs from gene expression time series. The methods use the automatic relevance determination (ARD) to incorporate the regulatory hub structure into the estimation of VAR1 in a Bayesian framework. Several sparse approximation schemes are applied to the estimated regression weights or VAR1 model to generate the sparse weighted adjacency matrices representing the inferred GRNs. We apply the proposed and several widespread reference methods to infer GRNs with up to 100 genes using simulated, DREAM4 in silico and experimental E. coli gene expression time series. We show that the proposed methods are efficient on simulated hub GRNs and scale-free GRNs using short time series simulated by VAR1s and outperform reference methods on small-scale DREAM4 in silico GRNs and E. coli GRNs. They can utilize the known major regulatory hubs to improve the performance on larger DREAM4 in silico GRNs and E. coli GRNs. The impact of nonlinear time series data on the performance of proposed methods is discussed.
Collapse
Affiliation(s)
- Chi-Kan Chen
- Department of Applied Mathematics, National Chung Hsing University, 145 Xingda Rd., South District, Taichung City, Taiwan, ROC
| |
Collapse
|
7
|
Kurz D, Sánchez CS, Axenie C. Data-Driven Discovery of Mathematical and Physical Relations in Oncology Data Using Human-Understandable Machine Learning. Front Artif Intell 2021; 4:713690. [PMID: 34901835 PMCID: PMC8655230 DOI: 10.3389/frai.2021.713690] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2021] [Accepted: 10/08/2021] [Indexed: 11/19/2022] Open
Abstract
For decades, researchers have used the concepts of rate of change and differential equations to model and forecast neoplastic processes. This expressive mathematical apparatus brought significant insights in oncology by describing the unregulated proliferation and host interactions of cancer cells, as well as their response to treatments. Now, these theories have been given a new life and found new applications. With the advent of routine cancer genome sequencing and the resulting abundance of data, oncology now builds an "arsenal" of new modeling and analysis tools. Models describing the governing physical laws of tumor-host-drug interactions can be now challenged with biological data to make predictions about cancer progression. Our study joins the efforts of the mathematical and computational oncology community by introducing a novel machine learning system for data-driven discovery of mathematical and physical relations in oncology. The system utilizes computational mechanisms such as competition, cooperation, and adaptation in neural networks to simultaneously learn the statistics and the governing relations between multiple clinical data covariates. Targeting an easy adoption in clinical oncology, the solutions of our system reveal human-understandable properties and features hidden in the data. As our experiments demonstrate, our system can describe nonlinear conservation laws in cancer kinetics and growth curves, symmetries in tumor's phenotypic staging transitions, the preoperative spatial tumor distribution, and up to the nonlinear intracellular and extracellular pharmacokinetics of neoadjuvant therapies. The primary goal of our work is to enhance or improve the mechanistic understanding of cancer dynamics by exploiting heterogeneous clinical data. We demonstrate through multiple instantiations that our system is extracting an accurate human-understandable representation of the underlying dynamics of physical interactions central to typical oncology problems. Our results and evaluation demonstrate that, using simple-yet powerful-computational mechanisms, such a machine learning system can support clinical decision-making. To this end, our system is a representative tool of the field of mathematical and computational oncology and offers a bridge between the data, the modeler, the data scientist, and the practicing clinician.
Collapse
Affiliation(s)
- Daria Kurz
- Interdisziplinäres Brustzentrum, Helios Klinikum München West, Akademisches Lehrkrankenhaus der Ludwig-Maximilians Universität München, Munich, Germany
| | - Carlos Salort Sánchez
- TUM School of Life Sciences Weihenstephan, Technical University of Munich, Freising, Germany
| | - Cristian Axenie
- Audi Konfuzius-Institut Ingolstadt Laboratory, Technische Hochschule Ingolstadt, Ingolstadt, Germany
| |
Collapse
|
8
|
Abstract
This perspective article gathers the latest developments in mathematical and computational oncology tools that exploit network approaches for the mathematical modelling, analysis, and simulation of cancer development and therapy design. It instigates the community to explore new paths and synergies under the umbrella of the Special Issue “Networks in Cancer: From Symmetry Breaking to Targeted Therapy”. The focus of the perspective is to demonstrate how networks can model the physics, analyse the interactions, and predict the evolution of the multiple processes behind tumour-host encounters across multiple scales. From agent-based modelling and mechano-biology to machine learning and predictive modelling, the perspective motivates a methodology well suited to mathematical and computational oncology and suggests approaches that mark a viable path towards adoption in the clinic.
Collapse
|
9
|
Caldu-Primo JL, Verduzco-Martínez JA, Alvarez-Buylla ER, Davila-Velderrain J. In vivo and in vitro human gene essentiality estimations capture contrasting functional constraints. NAR Genom Bioinform 2021; 3:lqab063. [PMID: 34268495 PMCID: PMC8276763 DOI: 10.1093/nargab/lqab063] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2021] [Revised: 06/18/2021] [Accepted: 07/07/2021] [Indexed: 11/28/2022] Open
Abstract
Gene essentiality estimation is a popular empirical approach to link genotypes to phenotypes. In humans, essentiality is estimated based on loss-of-function (LoF) mutation intolerance, either from population exome sequencing (in vivo) data or CRISPR-based in vitro perturbation experiments. Both approaches identify genes presumed to have detrimental consequences on the organism upon mutation. Are these genes constrained by having key cellular/organismal roles? Do in vivo and in vitro estimations equally recover these constraints? Insights into these questions have important implications in generalizing observations from cell models and interpreting disease risk genes. To empirically address these questions, we integrate genome-scale datasets and compare structural, functional and evolutionary features of essential genes versus genes with extremely high mutational tolerance. We found that essentiality estimates do recover functional constraints. However, the organismal or cellular context of estimation leads to functionally contrasting properties underlying the constraint. Our results suggest that depletion of LoF mutations in human populations effectively captures organismal-level functional constraints not experimentally accessible through CRISPR-based screens. Finally, we identify a set of genes (OrgEssential), which are mutationally intolerant in vivo but highly tolerant in vitro. These genes drive observed functional constraint differences and have an unexpected preference for nervous system expression.
Collapse
Affiliation(s)
- Jose Luis Caldu-Primo
- Instituto de Ecología, Universidad Nacional Autónoma de México, Cd. Universitaria, CDMX., 04510, México
| | - Jorge Armando Verduzco-Martínez
- Departamento de Biología Celular y Genética, Facultad de Ciencias Biológicas, Universidad Autónoma de Nuevo León, San Nicolás de los Garza, Nuevo León, 66400, México
| | - Elena R Alvarez-Buylla
- Instituto de Ecología, Universidad Nacional Autónoma de México, Cd. Universitaria, CDMX., 04510, México
| | | |
Collapse
|
10
|
Lai X, Dreyer FS, Cantone M, Eberhardt M, Gerer KF, Jaitly T, Uebe S, Lischer C, Ekici A, Wittmann J, Jäck HM, Schaft N, Dörrie J, Vera J. Network- and systems-based re-engineering of dendritic cells with non-coding RNAs for cancer immunotherapy. Theranostics 2021; 11:1412-1428. [PMID: 33391542 PMCID: PMC7738891 DOI: 10.7150/thno.53092] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2020] [Accepted: 10/15/2020] [Indexed: 12/12/2022] Open
Abstract
Dendritic cells (DCs) are professional antigen-presenting cells that induce and regulate adaptive immunity by presenting antigens to T cells. Due to their coordinative role in adaptive immune responses, DCs have been used as cell-based therapeutic vaccination against cancer. The capacity of DCs to induce a therapeutic immune response can be enhanced by re-wiring of cellular signalling pathways with microRNAs (miRNAs). Methods: Since the activation and maturation of DCs is controlled by an interconnected signalling network, we deploy an approach that combines RNA sequencing data and systems biology methods to delineate miRNA-based strategies that enhance DC-elicited immune responses. Results: Through RNA sequencing of IKKβ-matured DCs that are currently being tested in a clinical trial on therapeutic anti-cancer vaccination, we identified 44 differentially expressed miRNAs. According to a network analysis, most of these miRNAs regulate targets that are linked to immune pathways, such as cytokine and interleukin signalling. We employed a network topology-oriented scoring model to rank the miRNAs, analysed their impact on immunogenic potency of DCs, and identified dozens of promising miRNA candidates, with miR-15a and miR-16 as the top ones. The results of our analysis are presented in a database that constitutes a tool to identify DC-relevant miRNA-gene interactions with therapeutic potential (https://www.synmirapy.net/dc-optimization). Conclusions: Our approach enables the systematic analysis and identification of functional miRNA-gene interactions that can be experimentally tested for improving DC immunogenic potency.
Collapse
Affiliation(s)
- Xin Lai
- Laboratory of Systems Tumor Immunology, Department of Dermatology, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU) and Universitätsklinikum Erlangen, Erlangen, Germany
- Deutsches Zentrum Immuntherapie (DZI), Erlangen, Germany
- Comprehensive Cancer Center (CCC) Erlangen, Erlangen, Germany
| | - Florian S. Dreyer
- Laboratory of Systems Tumor Immunology, Department of Dermatology, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU) and Universitätsklinikum Erlangen, Erlangen, Germany
- Deutsches Zentrum Immuntherapie (DZI), Erlangen, Germany
- Comprehensive Cancer Center (CCC) Erlangen, Erlangen, Germany
| | - Martina Cantone
- Laboratory of Systems Tumor Immunology, Department of Dermatology, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU) and Universitätsklinikum Erlangen, Erlangen, Germany
- Deutsches Zentrum Immuntherapie (DZI), Erlangen, Germany
- Comprehensive Cancer Center (CCC) Erlangen, Erlangen, Germany
| | - Martin Eberhardt
- Laboratory of Systems Tumor Immunology, Department of Dermatology, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU) and Universitätsklinikum Erlangen, Erlangen, Germany
- Deutsches Zentrum Immuntherapie (DZI), Erlangen, Germany
- Comprehensive Cancer Center (CCC) Erlangen, Erlangen, Germany
| | - Kerstin F. Gerer
- RNA Group, Department of Dermatology, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU) and Universitätsklinikum Erlangen, Erlangen, Germany
- Deutsches Zentrum Immuntherapie (DZI), Erlangen, Germany
- Comprehensive Cancer Center (CCC) Erlangen, Erlangen, Germany
| | - Tanushree Jaitly
- Laboratory of Systems Tumor Immunology, Department of Dermatology, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU) and Universitätsklinikum Erlangen, Erlangen, Germany
- Deutsches Zentrum Immuntherapie (DZI), Erlangen, Germany
- Comprehensive Cancer Center (CCC) Erlangen, Erlangen, Germany
| | - Steffen Uebe
- Department of Human Genetics, Universitätsklinikum Erlangen, Erlangen, Germany
| | - Christopher Lischer
- Laboratory of Systems Tumor Immunology, Department of Dermatology, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU) and Universitätsklinikum Erlangen, Erlangen, Germany
- Deutsches Zentrum Immuntherapie (DZI), Erlangen, Germany
- Comprehensive Cancer Center (CCC) Erlangen, Erlangen, Germany
| | - Arif Ekici
- Department of Human Genetics, Universitätsklinikum Erlangen, Erlangen, Germany
| | - Jürgen Wittmann
- Division of Molecular Immunology, Department of Medicine 3, Universitätsklinikum Erlangen, Erlangen, Germany
| | - Hans-Martin Jäck
- Division of Molecular Immunology, Department of Medicine 3, Universitätsklinikum Erlangen, Erlangen, Germany
| | - Niels Schaft
- RNA Group, Department of Dermatology, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU) and Universitätsklinikum Erlangen, Erlangen, Germany
- Deutsches Zentrum Immuntherapie (DZI), Erlangen, Germany
- Comprehensive Cancer Center (CCC) Erlangen, Erlangen, Germany
| | - Jan Dörrie
- RNA Group, Department of Dermatology, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU) and Universitätsklinikum Erlangen, Erlangen, Germany
- Deutsches Zentrum Immuntherapie (DZI), Erlangen, Germany
- Comprehensive Cancer Center (CCC) Erlangen, Erlangen, Germany
| | - Julio Vera
- Laboratory of Systems Tumor Immunology, Department of Dermatology, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU) and Universitätsklinikum Erlangen, Erlangen, Germany
- Deutsches Zentrum Immuntherapie (DZI), Erlangen, Germany
- Comprehensive Cancer Center (CCC) Erlangen, Erlangen, Germany
| |
Collapse
|
11
|
Rubanova N, Pinna G, Kropp J, Campalans A, Radicella JP, Polesskaya A, Harel-Bellan A, Morozova N. MasterPATH: network analysis of functional genomics screening data. BMC Genomics 2020; 21:632. [PMID: 32928103 PMCID: PMC7491077 DOI: 10.1186/s12864-020-07047-2] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2019] [Accepted: 09/01/2020] [Indexed: 12/18/2022] Open
Abstract
Background Functional genomics employs several experimental approaches to investigate gene functions. High-throughput techniques, such as loss-of-function screening and transcriptome profiling, allow to identify lists of genes potentially involved in biological processes of interest (so called hit list). Several computational methods exist to analyze and interpret such lists, the most widespread of which aim either at investigating of significantly enriched biological processes, or at extracting significantly represented subnetworks. Results Here we propose a novel network analysis method and corresponding computational software that employs the shortest path approach and centrality measure to discover members of molecular pathways leading to the studied phenotype, based on functional genomics screening data. The method works on integrated interactomes that consist of both directed and undirected networks – HIPPIE, SIGNOR, SignaLink, TFactS, KEGG, TransmiR, miRTarBase. The method finds nodes and short simple paths with significant high centrality in subnetworks induced by the hit genes and by so-called final implementers – the genes that are involved in molecular events responsible for final phenotypic realization of the biological processes of interest. We present the application of the method to the data from miRNA loss-of-function screen and transcriptome profiling of terminal human muscle differentiation process and to the gene loss-of-function screen exploring the genes that regulates human oxidative DNA damage recognition. The analysis highlighted the possible role of several known myogenesis regulatory miRNAs (miR-1, miR-125b, miR-216a) and their targets (AR, NR3C1, ARRB1, ITSN1, VAV3, TDGF1), as well as linked two major regulatory molecules of skeletal myogenesis, MYOD and SMAD3, to their previously known muscle-related targets (TGFB1, CDC42, CTCF) and also to a number of proteins such as C-KIT that have not been previously studied in the context of muscle differentiation. The analysis also showed the role of the interaction between H3 and SETDB1 proteins for oxidative DNA damage recognition. Conclusion The current work provides a systematic methodology to discover members of molecular pathways in integrated networks using functional genomics screening data. It also offers a valuable instrument to explain the appearance of a set of genes, previously not associated with the process of interest, in the hit list of each particular functional genomics screening.
Collapse
Affiliation(s)
- Natalia Rubanova
- Institut des Hautes Etudes Scientifiques, Le Bois-Marie 35 rte de Chartres, 91440, Bures-sur-Yvette, France. .,Université Paris Diderot, Paris, France. .,Skolkovo Institute of Science and Technology, Skolkovo, Russia.
| | - Guillaume Pinna
- Institute for Integrative Biology of the Cell (I2BC), CEA, CNRS, Univ. Paris-Sud, Université Paris-Saclay, 91198, Gif-sur-Yvette cedex, France
| | - Jeremie Kropp
- Institut des Hautes Etudes Scientifiques, Le Bois-Marie 35 rte de Chartres, 91440, Bures-sur-Yvette, France
| | - Anna Campalans
- Institute of Molecular and Cellular Radiobiology, Institut François Jacob, CEA, F-92265, Fontenay-aux-Roses, France.,INSERM, U967, bâtiment 56 PC 103 18 route du Panorama, BP6 92265, Fontenay-aux-Roses Cedex, France.,Université Paris Sud, U967, bâtiment 56 PC 103 18 route du Panorama, BP6 92265, Fontenay-aux-Roses Cedex, France
| | - Juan Pablo Radicella
- Institute of Molecular and Cellular Radiobiology, Institut François Jacob, CEA, F-92265, Fontenay-aux-Roses, France.,INSERM, U967, bâtiment 56 PC 103 18 route du Panorama, BP6 92265, Fontenay-aux-Roses Cedex, France.,Université Paris Sud, U967, bâtiment 56 PC 103 18 route du Panorama, BP6 92265, Fontenay-aux-Roses Cedex, France
| | - Anna Polesskaya
- Ecole Polytechnique, Université Paris-Saclay, CNRS UMR 7654, Laboratoire de Biochimie, Ecole Polytechnique, 91128, Palaiseau, France
| | - Annick Harel-Bellan
- Institut des Hautes Etudes Scientifiques, Le Bois-Marie 35 rte de Chartres, 91440, Bures-sur-Yvette, France
| | - Nadya Morozova
- Institut des Hautes Etudes Scientifiques, Le Bois-Marie 35 rte de Chartres, 91440, Bures-sur-Yvette, France.,Institute for Integrative Biology of the Cell (I2BC), CEA, CNRS, Univ. Paris-Sud, Université Paris-Saclay, 91198, Gif-sur-Yvette cedex, France
| |
Collapse
|
12
|
Zhang B, Wang Y, Li H, Feng L, Li W, Cheng S. Identification of Prognostic Biomarkers for Multiple Solid Tumors Using a Human Villi Development Model. Front Cell Dev Biol 2020; 8:492. [PMID: 32656211 PMCID: PMC7325693 DOI: 10.3389/fcell.2020.00492] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2020] [Accepted: 05/25/2020] [Indexed: 01/03/2023] Open
Abstract
The processes of embryonic development that rely on epithelial-mesenchymal transition (EMT) for the implantation of trophoblast cells are co-opted by tumors, reflecting their inherent uncontrolled characteristics and leading to invasion and metastasis. Although tumorigenesis and embryogenesis have similar EMT characteristics, trophoblasts have been shown to exhibit "physiological metastasis" or be "pseudo-malignant," resulting in different outcomes. The gene co-expression network is the basis of embryonic development and tumorigenesis. We hypothesize that if the gene co-expression network in tumors is "off-track" from that in villi, it is more likely to develop into malignant tumors and have a worse prognosis, and we proposed the "off-track theory" for the first time. In this study, we examined gene co-expression networks in villi and multiple solid tumors. Through network functional enrichment analyses, we found that most tumors and villi exhibited a significantly enriched EMT, but the genes that performed this function were not identical. Then, we identified the "off-track genes" in the EMT-related gene interaction network using the "off-track theory," and through survival analysis, we discovered that the risk score of "off-track genes" was associated with poor survival of cancer patients. Our study indicated that villi development is a reliable and strictly regulated model that can illuminate the trajectory of human cancer development and that the gene co-expression networks in tumor development are "off-track" from those in villi. These "off-track genes" may have a substantial impact on tumor development and could reveal novel prognostic biomarkers.
Collapse
Affiliation(s)
- Botao Zhang
- Department of Neuro-oncology, Neurosurgery Center, Beijing Tiantan Hospital, Capital Medical University, Beijing, China.,State Key Laboratory of Molecular Oncology, Department of Etiology and Carcinogenesis, National Cancer Center/National Clinical Research Center for Cancer/Cancer Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China
| | - Yuanjing Wang
- Department of Gynecological Oncology, Beijing Obstetrics and Gynecology Hospital, Capital Medical University, Beijing, China
| | - Hongxia Li
- Department of Obstetrics and Gynecology, Beijing Shijitan Hospital, Capital Medical University, Beijing, China
| | - Lin Feng
- State Key Laboratory of Molecular Oncology, Department of Etiology and Carcinogenesis, National Cancer Center/National Clinical Research Center for Cancer/Cancer Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China
| | - Wenbin Li
- Department of Neuro-oncology, Neurosurgery Center, Beijing Tiantan Hospital, Capital Medical University, Beijing, China
| | - Shujun Cheng
- State Key Laboratory of Molecular Oncology, Department of Etiology and Carcinogenesis, National Cancer Center/National Clinical Research Center for Cancer/Cancer Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China
| |
Collapse
|
13
|
Werling DM, Pochareddy S, Choi J, An JY, Sheppard B, Peng M, Li Z, Dastmalchi C, Santpere G, Sousa AMM, Tebbenkamp ATN, Kaur N, Gulden FO, Breen MS, Liang L, Gilson MC, Zhao X, Dong S, Klei L, Cicek AE, Buxbaum JD, Adle-Biassette H, Thomas JL, Aldinger KA, O'Day DR, Glass IA, Zaitlen NA, Talkowski ME, Roeder K, State MW, Devlin B, Sanders SJ, Sestan N. Whole-Genome and RNA Sequencing Reveal Variation and Transcriptomic Coordination in the Developing Human Prefrontal Cortex. Cell Rep 2020; 31:107489. [PMID: 32268104 PMCID: PMC7295160 DOI: 10.1016/j.celrep.2020.03.053] [Citation(s) in RCA: 72] [Impact Index Per Article: 18.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2019] [Revised: 11/06/2019] [Accepted: 03/16/2020] [Indexed: 02/08/2023] Open
Abstract
Gene expression levels vary across developmental stage, cell type, and region in the brain. Genomic variants also contribute to the variation in expression, and some neuropsychiatric disorder loci may exert their effects through this mechanism. To investigate these relationships, we present BrainVar, a unique resource of paired whole-genome and bulk tissue RNA sequencing from the dorsolateral prefrontal cortex of 176 individuals across prenatal and postnatal development. Here we identify common variants that alter gene expression (expression quantitative trait loci [eQTLs]) constantly across development or predominantly during prenatal or postnatal stages. Both "constant" and "temporal-predominant" eQTLs are enriched for loci associated with neuropsychiatric traits and disorders and colocalize with specific variants. Expression levels of more than 12,000 genes rise or fall in a concerted late-fetal transition, with the transitional genes enriched for cell-type-specific genes and neuropsychiatric risk loci, underscoring the importance of cataloging developmental trajectories in understanding cortical physiology and pathology.
Collapse
Affiliation(s)
- Donna M Werling
- Department of Psychiatry, UCSF Weill Institute for Neurosciences, University of California, San Francisco, San Francisco, CA 94158, USA; Laboratory of Genetics, University of Wisconsin-Madison, Madison, WI 53706, USA
| | - Sirisha Pochareddy
- Department of Neuroscience and Kavli Institute for Neuroscience, Yale School of Medicine, New Haven, CT 06510, USA
| | - Jinmyung Choi
- Department of Neuroscience and Kavli Institute for Neuroscience, Yale School of Medicine, New Haven, CT 06510, USA
| | - Joon-Yong An
- Department of Psychiatry, UCSF Weill Institute for Neurosciences, University of California, San Francisco, San Francisco, CA 94158, USA; Department of Integrated Biomedical and Life Science, Korea University, Seoul 02841, Republic of Korea; School of Biosystem and Biomedical Science, College of Health Science, Korea University, Seoul 02841, Republic of Korea
| | - Brooke Sheppard
- Department of Psychiatry, UCSF Weill Institute for Neurosciences, University of California, San Francisco, San Francisco, CA 94158, USA
| | - Minshi Peng
- Department of Statistics and Data Science, Carnegie Mellon University, Pittsburgh, PA 15213, USA
| | - Zhen Li
- Department of Neuroscience and Kavli Institute for Neuroscience, Yale School of Medicine, New Haven, CT 06510, USA; Department of Neurosciences, University of California, San Diego, San Diego, CA 92093, USA
| | - Claudia Dastmalchi
- Department of Psychiatry, UCSF Weill Institute for Neurosciences, University of California, San Francisco, San Francisco, CA 94158, USA
| | - Gabriel Santpere
- Department of Neuroscience and Kavli Institute for Neuroscience, Yale School of Medicine, New Haven, CT 06510, USA; Neurogenomics Group, Research Programme on Biomedical Informatics, Hospital del Mar Medical Research Institute, Department of Experimental and Health Sciences, Universitat Pompeu Fabra, 08003 Barcelona, Catalonia, Spain
| | - André M M Sousa
- Department of Neuroscience and Kavli Institute for Neuroscience, Yale School of Medicine, New Haven, CT 06510, USA
| | - Andrew T N Tebbenkamp
- Department of Neuroscience and Kavli Institute for Neuroscience, Yale School of Medicine, New Haven, CT 06510, USA
| | - Navjot Kaur
- Department of Neuroscience and Kavli Institute for Neuroscience, Yale School of Medicine, New Haven, CT 06510, USA
| | - Forrest O Gulden
- Department of Neuroscience and Kavli Institute for Neuroscience, Yale School of Medicine, New Haven, CT 06510, USA
| | - Michael S Breen
- Seaver Autism Center for Research and Treatment, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA; Department of Psychiatry, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA; Mindich Child Health and Development Institute, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA; Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA
| | - Lindsay Liang
- Department of Psychiatry, UCSF Weill Institute for Neurosciences, University of California, San Francisco, San Francisco, CA 94158, USA
| | - Michael C Gilson
- Department of Psychiatry, UCSF Weill Institute for Neurosciences, University of California, San Francisco, San Francisco, CA 94158, USA
| | - Xuefang Zhao
- Center for Genomic Medicine and Department of Neurology, Massachusetts General Hospital, Boston, MA 02114, USA; Department of Neurology, Harvard Medical School, Boston, MA 02115, USA; Program in Medical and Population Genetics and Stanley Center for Psychiatric Research, Broad Institute, Cambridge, MA 02142, USA
| | - Shan Dong
- Department of Psychiatry, UCSF Weill Institute for Neurosciences, University of California, San Francisco, San Francisco, CA 94158, USA
| | - Lambertus Klei
- Department of Psychiatry, University of Pittsburgh School of Medicine, Pittsburgh, PA 15213, USA
| | - A Ercument Cicek
- Department of Computer Engineering, Bilkent University, Ankara 06800, Turkey; Computational Biology Department, Carnegie Mellon University, Pittsburgh, PA 15213, USA
| | - Joseph D Buxbaum
- Seaver Autism Center for Research and Treatment, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA; Department of Psychiatry, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA; Mindich Child Health and Development Institute, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA; Friedman Brain Institute, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA
| | - Homa Adle-Biassette
- Department of Pathology, Lariboisière Hospital, APHP, Biobank BB-0033-00064, and Université de Paris, 75006 Paris, France
| | - Jean-Leon Thomas
- Department of Neurology, Yale University School of Medicine, New Haven, CT 06511, USA; UMRS1127, Sorbonne Université, Institut du Cerveau et de la Moelle Épinière, 75013 Paris, France
| | - Kimberly A Aldinger
- Center for Integrative Brain Research, Seattle Children's Research Institute, Seattle, WA 98101, USA; Brotman Baty Institute for Precision Medicine, Seattle, WA 98195, USA
| | - Diana R O'Day
- Department of Pediatrics, University of Washington, Seattle, WA 98105, USA
| | - Ian A Glass
- Department of Pediatrics, University of Washington, Seattle, WA 98105, USA
| | - Noah A Zaitlen
- Department of Medicine, University of California, San Francisco, San Francisco, CA 94158, USA
| | - Michael E Talkowski
- Center for Genomic Medicine and Department of Neurology, Massachusetts General Hospital, Boston, MA 02114, USA; Department of Neurology, Harvard Medical School, Boston, MA 02115, USA; Program in Medical and Population Genetics and Stanley Center for Psychiatric Research, Broad Institute, Cambridge, MA 02142, USA
| | - Kathryn Roeder
- Department of Statistics and Data Science, Carnegie Mellon University, Pittsburgh, PA 15213, USA; Computational Biology Department, Carnegie Mellon University, Pittsburgh, PA 15213, USA
| | - Matthew W State
- Department of Psychiatry, UCSF Weill Institute for Neurosciences, University of California, San Francisco, San Francisco, CA 94158, USA; Institute for Human Genetics, University of California, San Francisco, San Francisco, CA 94158, USA
| | - Bernie Devlin
- Department of Psychiatry, University of Pittsburgh School of Medicine, Pittsburgh, PA 15213, USA
| | - Stephan J Sanders
- Department of Psychiatry, UCSF Weill Institute for Neurosciences, University of California, San Francisco, San Francisco, CA 94158, USA; Institute for Human Genetics, University of California, San Francisco, San Francisco, CA 94158, USA.
| | - Nenad Sestan
- Department of Neuroscience and Kavli Institute for Neuroscience, Yale School of Medicine, New Haven, CT 06510, USA; Department of Psychiatry, Yale University School of Medicine, New Haven, CT 06520, USA; Department of Genetics, Yale University School of Medicine, New Haven, CT 06520, USA; Department of Comparative Medicine, Program in Integrative Cell Signaling and Neurobiology of Metabolism, Yale School of Medicine, New Haven, CT 06510, USA; Program in Cellular Neuroscience, Neurodegeneration, and Repair and Yale Child Study Center, Yale School of Medicine, New Haven, CT 06510, USA.
| |
Collapse
|
14
|
Blatti C, Emad A, Berry MJ, Gatzke L, Epstein M, Lanier D, Rizal P, Ge J, Liao X, Sobh O, Lambert M, Post CS, Xiao J, Groves P, Epstein AT, Chen X, Srinivasan S, Lehnert E, Kalari KR, Wang L, Weinshilboum RM, Song JS, Jongeneel CV, Han J, Ravaioli U, Sobh N, Bushell CB, Sinha S. Knowledge-guided analysis of "omics" data using the KnowEnG cloud platform. PLoS Biol 2020; 18:e3000583. [PMID: 31971940 PMCID: PMC6977717 DOI: 10.1371/journal.pbio.3000583] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2019] [Accepted: 12/19/2019] [Indexed: 12/19/2022] Open
Abstract
We present Knowledge Engine for Genomics (KnowEnG), a free-to-use computational system for analysis of genomics data sets, designed to accelerate biomedical discovery. It includes tools for popular bioinformatics tasks such as gene prioritization, sample clustering, gene set analysis, and expression signature analysis. The system specializes in "knowledge-guided" data mining and machine learning algorithms, in which user-provided data are analyzed in light of prior information about genes, aggregated from numerous knowledge bases and encoded in a massive "Knowledge Network." KnowEnG adheres to "FAIR" principles (findable, accessible, interoperable, and reuseable): its tools are easily portable to diverse computing environments, run on the cloud for scalable and cost-effective execution, and are interoperable with other computing platforms. The analysis tools are made available through multiple access modes, including a web portal with specialized visualization modules. We demonstrate the KnowEnG system's potential value in democratization of advanced tools for the modern genomics era through several case studies that use its tools to recreate and expand upon the published analysis of cancer data sets.
Collapse
Affiliation(s)
- Charles Blatti
- Carl R. Woese Institute for Genomic Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Amin Emad
- Carl R. Woese Institute for Genomic Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
- Department of Electrical and Computer Engineering, McGill University, Montreal, Canada
| | - Matthew J. Berry
- National Center for Supercomputing Applications, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Lisa Gatzke
- National Center for Supercomputing Applications, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Milt Epstein
- Carl R. Woese Institute for Genomic Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Daniel Lanier
- Carl R. Woese Institute for Genomic Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Pramod Rizal
- National Center for Supercomputing Applications, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Jing Ge
- Carl R. Woese Institute for Genomic Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Xiaoxia Liao
- National Center for Supercomputing Applications, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Omar Sobh
- Carl R. Woese Institute for Genomic Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Mike Lambert
- National Center for Supercomputing Applications, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Corey S. Post
- Department of Computer Science, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Jinfeng Xiao
- Department of Computer Science, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Peter Groves
- National Center for Supercomputing Applications, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Aidan T. Epstein
- Carl R. Woese Institute for Genomic Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Xi Chen
- Carl R. Woese Institute for Genomic Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Subhashini Srinivasan
- Carl R. Woese Institute for Genomic Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Erik Lehnert
- Seven Bridges Genomics, Charlestown, Massachusetts, United States of America
| | - Krishna R. Kalari
- Department of Health Sciences Research, Mayo Clinic, Rochester, Minnesota, United States of America
| | - Liewei Wang
- Department of Molecular Pharmacology and Experimental Therapeutics, Mayo Clinic, Rochester, Minnesota, United States of America
| | - Richard M. Weinshilboum
- Department of Molecular Pharmacology and Experimental Therapeutics, Mayo Clinic, Rochester, Minnesota, United States of America
| | - Jun S. Song
- Carl R. Woese Institute for Genomic Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
- Cancer Center at Illinois, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - C. Victor Jongeneel
- Carl R. Woese Institute for Genomic Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Jiawei Han
- Department of Computer Science, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
- Cancer Center at Illinois, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Umberto Ravaioli
- Department of Electrical and Computer Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Nahil Sobh
- Carl R. Woese Institute for Genomic Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Colleen B. Bushell
- National Center for Supercomputing Applications, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Saurabh Sinha
- Carl R. Woese Institute for Genomic Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
- Department of Computer Science, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
- Cancer Center at Illinois, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
- * E-mail:
| |
Collapse
|
15
|
Thomas PD, Hill DP, Mi H, Osumi-Sutherland D, Van Auken K, Carbon S, Balhoff JP, Albou LP, Good B, Gaudet P, Lewis SE, Mungall CJ. Gene Ontology Causal Activity Modeling (GO-CAM) moves beyond GO annotations to structured descriptions of biological functions and systems. Nat Genet 2020; 51:1429-1433. [PMID: 31548717 DOI: 10.1038/s41588-019-0500-1] [Citation(s) in RCA: 68] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Paul D Thomas
- Department of Preventive Medicine, Division of Bioinformatics, University of Southern California, Los Angeles, CA, USA.
| | | | - Huaiyu Mi
- Department of Preventive Medicine, Division of Bioinformatics, University of Southern California, Los Angeles, CA, USA
| | - David Osumi-Sutherland
- European Bioinformatics Institute (EMBL-EBI), Wellcome Trust Genome Campus, Cambridge, UK
| | - Kimberly Van Auken
- Division of Biology and Biological Engineering, California Institute of Technology, Pasadena, CA, USA
| | - Seth Carbon
- Division of Environmental Genomics and Systems Biology, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - James P Balhoff
- Renaissance Computing Institute, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA
| | - Laurent-Philippe Albou
- Department of Preventive Medicine, Division of Bioinformatics, University of Southern California, Los Angeles, CA, USA
| | - Benjamin Good
- Division of Environmental Genomics and Systems Biology, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Pascale Gaudet
- SIB Swiss Institute of Bioinformatics, Geneva, Switzerland
| | - Suzanna E Lewis
- Division of Environmental Genomics and Systems Biology, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Christopher J Mungall
- Division of Environmental Genomics and Systems Biology, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| |
Collapse
|
16
|
Sulaimanov N, Kumar S, Burdet F, Ibberson M, Pagni M, Koeppl H. Inferring gene expression networks with hubs using a degree weighted Lasso approach. Bioinformatics 2019; 35:987-994. [PMID: 30165436 DOI: 10.1093/bioinformatics/bty716] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2017] [Revised: 06/08/2018] [Accepted: 08/25/2018] [Indexed: 11/14/2022] Open
Abstract
MOTIVATION Genome-scale gene networks contain regulatory genes called hubs that have many interaction partners. These genes usually play an essential role in gene regulation and cellular processes. Despite recent advancements in high-throughput technology, inferring gene networks with hub genes from high-dimensional data still remains a challenging problem. Novel statistical network inference methods are needed for efficient and accurate reconstruction of hub networks from high-dimensional data. RESULTS To address this challenge we propose DW-Lasso, a degree weighted Lasso (least absolute shrinkage and selection operator) method which infers gene networks with hubs efficiently under the low sample size setting. Our network reconstruction approach is formulated as a two stage procedure: first, the degree of networks is estimated iteratively, and second, the gene regulatory network is reconstructed using degree information. A useful property of the proposed method is that it naturally favors the accumulation of neighbors around hub genes and thereby helps in accurate modeling of the high-throughput data under the assumption that the underlying network exhibits hub structure. In a simulation study, we demonstrate good predictive performance of the proposed method in comparison to traditional Lasso type methods in inferring hub and scale-free graphs. We show the effectiveness of our method in an application to microarray data of Escherichia coli and RNA sequencing data of Kidney Clear Cell Carcinoma from The Cancer Genome Atlas datasets. AVAILABILITY AND IMPLEMENTATION Under the GNU General Public Licence at https://cran.r-project.org/package=DWLasso. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Nurgazy Sulaimanov
- Department of Electrical Engineering and Information Technology, Technische Universität Darmstadt, Darmstadt, Germany.,Department of Biology, Technische Universität Darmstadt, Darmstadt, Germany
| | - Sunil Kumar
- Department of Electrical Engineering and Information Technology, Technische Universität Darmstadt, Darmstadt, Germany.,Department of Biology, Technische Universität Darmstadt, Darmstadt, Germany
| | | | - Mark Ibberson
- Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Marco Pagni
- Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Heinz Koeppl
- Department of Electrical Engineering and Information Technology, Technische Universität Darmstadt, Darmstadt, Germany.,Department of Biology, Technische Universität Darmstadt, Darmstadt, Germany
| |
Collapse
|
17
|
Vodopivec M, Lah L, Narat M, Curk T. Metabolomic profiling of CHO fed-batch growth phases at 10, 100, and 1,000 L. Biotechnol Bioeng 2019; 116:2720-2729. [PMID: 31184374 DOI: 10.1002/bit.27087] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2019] [Revised: 05/17/2019] [Accepted: 06/03/2019] [Indexed: 01/08/2023]
Abstract
Established bioprocess monitoring is based on quick and reliable methods, including cell count and viability measurement, extracellular metabolite measurement, and the measurement of physicochemical qualities of the cultivation medium. These methods are sufficient for monitoring of process performance, but rarely give insight into the actual physiological states of the cell culture. However, understanding of the latter is essential for optimization of bioprocess development. Our study used LC-MS metabolomics as a tool for additional resolution of bioprocess monitoring and was designed at three bioreactors scales (10 L, 100 L, and 1,000 L) to gain insight into the basal metabolic states of the Chinese hamster ovary (CHO) cell culture during fed-batch. Metabolites characteristics of the four growth stages (early and late exponential phase, stationary phase, and the phase of decline) were identified by multivariate analysis. Enriched metabolic pathways were then established for each growth phase using the CHO metabolic network model. Biomass generation and nucleotide synthesis were enriched in early exponential phase, followed by increased protein production and imbalanced glutathione metabolism in late exponential phase. Glycolysis became downregulated in stationary phase and amino-acid metabolism increased. Phase of culture decline resulted in rise of oxidized glutathione and fatty acid concentrations. Intracellular metabolic profiles of the CHO fed-batch culture were also shown to be consistent with scale and thus demonstrate metabolomic profiling as an informative method to gain physiological insight into the cell culture states during bioprocess regardless of scale.
Collapse
Affiliation(s)
- Maja Vodopivec
- Bioprocess Development, Technical Development Biologics Mengeš, Novartis Technical Research & Development, Lek Pharmaceuticals d.d, Slovenia
| | - Ljerka Lah
- Bioprocess Development, Technical Development Biologics Mengeš, Novartis Technical Research & Development, Lek Pharmaceuticals d.d, Slovenia
| | - Mojca Narat
- Department of Animal Science, Biotechnical Faculty, University of Ljubljana, Ljubljana, Slovenia
| | - Tomaž Curk
- Bioinformatics Laboratory, Faculty of Computer and Information Science, University of Ljubljana, Ljubljana, Slovenija
| |
Collapse
|
18
|
Glaab E. Computational systems biology approaches for Parkinson's disease. Cell Tissue Res 2018; 373:91-109. [PMID: 29185073 PMCID: PMC6015628 DOI: 10.1007/s00441-017-2734-5] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2017] [Accepted: 11/06/2017] [Indexed: 12/26/2022]
Abstract
Parkinson's disease (PD) is a prime example of a complex and heterogeneous disorder, characterized by multifaceted and varied motor- and non-motor symptoms and different possible interplays of genetic and environmental risk factors. While investigations of individual PD-causing mutations and risk factors in isolation are providing important insights to improve our understanding of the molecular mechanisms behind PD, there is a growing consensus that a more complete understanding of these mechanisms will require an integrative modeling of multifactorial disease-associated perturbations in molecular networks. Identifying and interpreting the combinatorial effects of multiple PD-associated molecular changes may pave the way towards an earlier and reliable diagnosis and more effective therapeutic interventions. This review provides an overview of computational systems biology approaches developed in recent years to study multifactorial molecular alterations in complex disorders, with a focus on PD research applications. Strengths and weaknesses of different cellular pathway and network analyses, and multivariate machine learning techniques for investigating PD-related omics data are discussed, and strategies proposed to exploit the synergies of multiple biological knowledge and data sources. A final outlook provides an overview of specific challenges and possible next steps for translating systems biology findings in PD to new omics-based diagnostic tools and targeted, drug-based therapeutic approaches.
Collapse
Affiliation(s)
- Enrico Glaab
- Luxembourg Centre for Systems Biomedicine (LCSB), University of Luxembourg, 7 avenue des Hauts Fourneaux, L-4362, Esch-sur-Alzette, Luxembourg.
| |
Collapse
|
19
|
Li T, Kim A, Rosenbluh J, Horn H, Greenfeld L, An D, Zimmer A, Liberzon A, Bistline J, Natoli T, Li Y, Tsherniak A, Narayan R, Subramanian A, Liefeld T, Wong B, Thompson D, Calvo S, Carr S, Boehm J, Jaffe J, Mesirov J, Hacohen N, Regev A, Lage K. GeNets: a unified web platform for network-based genomic analyses. Nat Methods 2018; 15:543-546. [PMID: 29915188 PMCID: PMC6450090 DOI: 10.1038/s41592-018-0039-6] [Citation(s) in RCA: 43] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2017] [Accepted: 05/02/2018] [Indexed: 11/25/2022]
Abstract
Functional genomics networks are widely used to identify unexpected pathway relationships in large genomic datasets. However, it is challenging to quantitatively compare the signal-to-noise ratio of different networks, the biology they describe, and to identify the optimal network to interpret a particular genetic dataset. Via GeNets users can train a machine-learning model (Quack) to make such comparisons; and they can execute, store, and share analyses of genetic and RNA sequencing datasets.
Collapse
Affiliation(s)
- Taibo Li
- Department of Surgery, Massachusetts General Hospital, Boston, MA, USA.,Broad Institute of MIT and Harvard, Cambridge, MA, USA.,Department of Electrical Engineering & Computer Science, MIT, Cambridge, MA, USA.,Johns Hopkins School of Medicine, Baltimore, MD, USA
| | - April Kim
- Department of Surgery, Massachusetts General Hospital, Boston, MA, USA.,Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Joseph Rosenbluh
- Broad Institute of MIT and Harvard, Cambridge, MA, USA.,Department of Medical Oncology, Dana-Farber Cancer Institute, Harvard Medical School, Boston, MA, USA
| | - Heiko Horn
- Department of Surgery, Massachusetts General Hospital, Boston, MA, USA.,Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Liraz Greenfeld
- Department of Surgery, Massachusetts General Hospital, Boston, MA, USA.,Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - David An
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Andrew Zimmer
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | | | - Jon Bistline
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Ted Natoli
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Yang Li
- Broad Institute of MIT and Harvard, Cambridge, MA, USA.,Howard Hughes Medical Institute and Department of Molecular Biology, Massachusetts General Hospital, Boston, MA, USA.,Department of Statistics, Harvard University, Cambridge, MA, USA
| | | | - Rajiv Narayan
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | | | - Ted Liefeld
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Bang Wong
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Dawn Thompson
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Sarah Calvo
- Broad Institute of MIT and Harvard, Cambridge, MA, USA.,Howard Hughes Medical Institute and Department of Molecular Biology, Massachusetts General Hospital, Boston, MA, USA
| | - Steve Carr
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Jesse Boehm
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Jake Jaffe
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Jill Mesirov
- Broad Institute of MIT and Harvard, Cambridge, MA, USA.,Department of Medicine, University of California, San Diego, San Diego, CA, USA
| | - Nir Hacohen
- Broad Institute of MIT and Harvard, Cambridge, MA, USA.,Center for Immunology and Inflammatory Diseases, Massachusetts General Hospital, Boston, MA, USA.,Department of Surgery, Harvard Medical School, Boston, MA, USA
| | - Aviv Regev
- Broad Institute of MIT and Harvard, Cambridge, MA, USA.,Howard Hughes Medical Institute, Department of Biology, MIT, Cambridge, MA, USA
| | - Kasper Lage
- Department of Surgery, Massachusetts General Hospital, Boston, MA, USA. .,Broad Institute of MIT and Harvard, Cambridge, MA, USA. .,Department of Surgery, Harvard Medical School, Boston, MA, USA. .,Institute for Biological Psychiatry, Mental Health Center Sct. Hans, University of Copenhagen, Roskilde, Denmark.
| |
Collapse
|
20
|
Robinson S, Nevalainen J, Pinna G, Campalans A, Radicella JP, Guyon L. Incorporating interaction networks into the determination of functionally related hit genes in genomic experiments with Markov random fields. Bioinformatics 2018; 33:i170-i179. [PMID: 28881978 PMCID: PMC5870666 DOI: 10.1093/bioinformatics/btx244] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Abstract
Motivation Incorporating gene interaction data into the identification of ‘hit’ genes in genomic experiments is a well-established approach leveraging the ‘guilt by association’ assumption to obtain a network based hit list of functionally related genes. We aim to develop a method to allow for multivariate gene scores and multiple hit labels in order to extend the analysis of genomic screening data within such an approach. Results We propose a Markov random field-based method to achieve our aim and show that the particular advantages of our method compared with those currently used lead to new insights in previously analysed data as well as for our own motivating data. Our method additionally achieves the best performance in an independent simulation experiment. The real data applications we consider comprise of a survival analysis and differential expression experiment and a cell-based RNA interference functional screen. Availability and implementation We provide all of the data and code related to the results in the paper. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Sean Robinson
- CEA, BIG, Biologie à Grande Echelle, F-38054 Grenoble, France.,Université Grenoble-Alpes, F-38000 Grenoble, France.,INSERM, U1038, F-38054 Grenoble, France.,Department of Mathematics and Statistics, University of Turku, Turku, Finland.,Industrial Biotechnology, VTT Technical Research Centre of Finland, Turku, Finland
| | - Jaakko Nevalainen
- Department of Mathematics and Statistics, University of Turku, Turku, Finland.,School of Health Sciences, University of Tampere, Tampere, Finland
| | - Guillaume Pinna
- Plateforme ARN Interférence (PArI), DSV/ISVFJ/SBIGEM/UMR 9198 I2BC, CEA Saclay, Gif-sur-Yvette, France
| | - Anna Campalans
- Institute of Molecular and Cellular Radiobiology, CEA, Fontenay-aux-Roses, France.,INSERM, U967, Fontenay-aux-Roses, France.,Université Paris Diderot, U967, Fontenay-aux-Roses, France.,Université Paris Sud, U967, Fontenay-aux-Roses, France
| | - J Pablo Radicella
- Institute of Molecular and Cellular Radiobiology, CEA, Fontenay-aux-Roses, France.,INSERM, U967, Fontenay-aux-Roses, France.,Université Paris Diderot, U967, Fontenay-aux-Roses, France.,Université Paris Sud, U967, Fontenay-aux-Roses, France
| | - Laurent Guyon
- CEA, BIG, Biologie à Grande Echelle, F-38054 Grenoble, France.,Université Grenoble-Alpes, F-38000 Grenoble, France.,INSERM, U1038, F-38054 Grenoble, France
| |
Collapse
|
21
|
Agrawal M, Zitnik M, Leskovec J. Large-scale analysis of disease pathways in the human interactome. PACIFIC SYMPOSIUM ON BIOCOMPUTING. PACIFIC SYMPOSIUM ON BIOCOMPUTING 2018; 23:111-122. [PMID: 29218874 PMCID: PMC5731453] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Discovering disease pathways, which can be defined as sets of proteins associated with a given disease, is an important problem that has the potential to provide clinically actionable insights for disease diagnosis, prognosis, and treatment. Computational methods aid the discovery by relying on protein-protein interaction (PPI) networks. They start with a few known disease-associated proteins and aim to find the rest of the pathway by exploring the PPI network around the known disease proteins. However, the success of such methods has been limited, and failure cases have not been well understood. Here we study the PPI network structure of 519 disease pathways. We find that 90% of pathways do not correspond to single well-connected components in the PPI network. Instead, proteins associated with a single disease tend to form many separate connected components/regions in the network. We then evaluate state-of-the-art disease pathway discovery methods and show that their performance is especially poor on diseases with disconnected pathways. Thus, we conclude that network connectivity structure alone may not be sufficient for disease pathway discovery. However, we show that higher-order network structures, such as small subgraphs of the pathway, provide a promising direction for the development of new methods.
Collapse
Affiliation(s)
- Monica Agrawal
- Department of Computer Science, Stanford University, Stanford, CA, USA,
| | | | | |
Collapse
|
22
|
Pirkl M, Diekmann M, van der Wees M, Beerenwinkel N, Fröhlich H, Markowetz F. Inferring modulators of genetic interactions with epistatic nested effects models. PLoS Comput Biol 2017; 13:e1005496. [PMID: 28406896 PMCID: PMC5407847 DOI: 10.1371/journal.pcbi.1005496] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2016] [Revised: 04/27/2017] [Accepted: 04/03/2017] [Indexed: 12/27/2022] Open
Abstract
Maps of genetic interactions can dissect functional redundancies in cellular networks. Gene expression profiles as high-dimensional molecular readouts of combinatorial perturbations provide a detailed view of genetic interactions, but can be hard to interpret if different gene sets respond in different ways (called mixed epistasis). Here we test the hypothesis that mixed epistasis between a gene pair can be explained by the action of a third gene that modulates the interaction. We have extended the framework of Nested Effects Models (NEMs), a type of graphical model specifically tailored to analyze high-dimensional gene perturbation data, to incorporate logical functions that describe interactions between regulators on downstream genes and proteins. We benchmark our approach in the controlled setting of a simulation study and show high accuracy in inferring the correct model. In an application to data from deletion mutants of kinases and phosphatases in S. cerevisiae we show that epistatic NEMs can point to modulators of genetic interactions. Our approach is implemented in the R-package 'epiNEM' available from https://github.com/cbg-ethz/epiNEM and https://bioconductor.org/packages/epiNEM/.
Collapse
Affiliation(s)
- Martin Pirkl
- ETH Zurich, Department of Biosystems Science and Engineering, Basel, Switzerland
- SIB Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Madeline Diekmann
- ETH Zurich, Department of Biosystems Science and Engineering, Basel, Switzerland
- SIB Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | | | - Niko Beerenwinkel
- ETH Zurich, Department of Biosystems Science and Engineering, Basel, Switzerland
- SIB Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Holger Fröhlich
- Bonn-Aachen International Center for IT (B-IT), University of Bonn, Bonn, Germany
- UCB Biosciences GmbH, Monheim, Germany
| | - Florian Markowetz
- University of Cambridge, Cancer Research UK Cambridge Institute, Cambridge, United Kingdom
| |
Collapse
|
23
|
Kaushik A, Ali S, Gupta D. Altered Pathway Analyzer: A gene expression dataset analysis tool for identification and prioritization of differentially regulated and network rewired pathways. Sci Rep 2017; 7:40450. [PMID: 28084397 PMCID: PMC5233954 DOI: 10.1038/srep40450] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2016] [Accepted: 12/07/2016] [Indexed: 12/13/2022] Open
Abstract
Gene connection rewiring is an essential feature of gene network dynamics. Apart from its normal functional role, it may also lead to dysregulated functional states by disturbing pathway homeostasis. Very few computational tools measure rewiring within gene co-expression and its corresponding regulatory networks in order to identify and prioritize altered pathways which may or may not be differentially regulated. We have developed Altered Pathway Analyzer (APA), a microarray dataset analysis tool for identification and prioritization of altered pathways, including those which are differentially regulated by TFs, by quantifying rewired sub-network topology. Moreover, APA also helps in re-prioritization of APA shortlisted altered pathways enriched with context-specific genes. We performed APA analysis of simulated datasets and p53 status NCI-60 cell line microarray data to demonstrate potential of APA for identification of several case-specific altered pathways. APA analysis reveals several altered pathways not detected by other tools evaluated by us. APA analysis of unrelated prostate cancer datasets identifies sample-specific as well as conserved altered biological processes, mainly associated with lipid metabolism, cellular differentiation and proliferation. APA is designed as a cross platform tool which may be transparently customized to perform pathway analysis in different gene expression datasets. APA is freely available at http://bioinfo.icgeb.res.in/APA.
Collapse
Affiliation(s)
- Abhinav Kaushik
- Translational Bioinformatics Group, International Centre for Genetic Engineering and Biotechnology, New Delhi 110067, India
| | - Shakir Ali
- Department of Biochemistry, Jamia Hamdard, Deemed University, New Delhi 110062, India
| | - Dinesh Gupta
- Translational Bioinformatics Group, International Centre for Genetic Engineering and Biotechnology, New Delhi 110067, India
| |
Collapse
|
24
|
Loviglio MN, Beck CR, White JJ, Leleu M, Harel T, Guex N, Niknejad A, Bi W, Chen ES, Crespo I, Yan J, Charng WL, Gu S, Fang P, Coban-Akdemir Z, Shaw CA, Jhangiani SN, Muzny DM, Gibbs RA, Rougemont J, Xenarios I, Lupski JR, Reymond A. Identification of a RAI1-associated disease network through integration of exome sequencing, transcriptomics, and 3D genomics. Genome Med 2016; 8:105. [PMID: 27799067 PMCID: PMC5088687 DOI: 10.1186/s13073-016-0359-z] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2016] [Accepted: 09/16/2016] [Indexed: 02/13/2023] Open
Abstract
Background Smith-Magenis syndrome (SMS) is a developmental disability/multiple congenital anomaly disorder resulting from haploinsufficiency of RAI1. It is characterized by distinctive facial features, brachydactyly, sleep disturbances, and stereotypic behaviors. Methods We investigated a cohort of 15 individuals with a clinical suspicion of SMS who showed neither deletion in the SMS critical region nor damaging variants in RAI1 using whole exome sequencing. A combination of network analysis (co-expression and biomedical text mining), transcriptomics, and circularized chromatin conformation capture (4C-seq) was applied to verify whether modified genes are part of the same disease network as known SMS-causing genes. Results Potentially deleterious variants were identified in nine of these individuals using whole-exome sequencing. Eight of these changes affect KMT2D, ZEB2, MAP2K2, GLDC, CASK, MECP2, KDM5C, and POGZ, known to be associated with Kabuki syndrome 1, Mowat-Wilson syndrome, cardiofaciocutaneous syndrome, glycine encephalopathy, mental retardation and microcephaly with pontine and cerebellar hypoplasia, X-linked mental retardation 13, X-linked mental retardation Claes-Jensen type, and White-Sutton syndrome, respectively. The ninth individual carries a de novo variant in JAKMIP1, a regulator of neuronal translation that was recently found deleted in a patient with autism spectrum disorder. Analyses of co-expression and biomedical text mining suggest that these pathologies and SMS are part of the same disease network. Further support for this hypothesis was obtained from transcriptome profiling that showed that the expression levels of both Zeb2 and Map2k2 are perturbed in Rai1–/– mice. As an orthogonal approach to potentially contributory disease gene variants, we used chromatin conformation capture to reveal chromatin contacts between RAI1 and the loci flanking ZEB2 and GLDC, as well as between RAI1 and human orthologs of the genes that show perturbed expression in our Rai1–/– mouse model. Conclusions These holistic studies of RAI1 and its interactions allow insights into SMS and other disorders associated with intellectual disability and behavioral abnormalities. Our findings support a pan-genomic approach to the molecular diagnosis of a distinctive disorder. Electronic supplementary material The online version of this article (doi:10.1186/s13073-016-0359-z) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Maria Nicla Loviglio
- Center for Integrative Genomics, University of Lausanne, 1015, Lausanne, Switzerland
| | - Christine R Beck
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, 77030, USA
| | - Janson J White
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, 77030, USA
| | - Marion Leleu
- School of Life Sciences, EPFL (Ecole Polytechnique Fédérale de Lausanne), 1015, Lausanne, Switzerland.,Swiss Institute of Bioinformatics (SIB), 1015, Lausanne, Switzerland
| | - Tamar Harel
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, 77030, USA
| | - Nicolas Guex
- Swiss Institute of Bioinformatics (SIB), 1015, Lausanne, Switzerland
| | - Anne Niknejad
- Swiss Institute of Bioinformatics (SIB), 1015, Lausanne, Switzerland
| | - Weimin Bi
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, 77030, USA
| | - Edward S Chen
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, 77030, USA
| | - Isaac Crespo
- Swiss Institute of Bioinformatics (SIB), 1015, Lausanne, Switzerland
| | - Jiong Yan
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, 77030, USA.,Laboratory Medicine Program, UHN, Department of Laboratory Medicine and Pathobiology, University of Toronto, Toronto, ON, M5G 2C4, Canada
| | - Wu-Lin Charng
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, 77030, USA
| | - Shen Gu
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, 77030, USA
| | - Ping Fang
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, 77030, USA.,Present address: WuXiNextCODE, 101Main Street, Cambridge, MA, 02142, USA
| | - Zeynep Coban-Akdemir
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, 77030, USA
| | - Chad A Shaw
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, 77030, USA
| | - Shalini N Jhangiani
- Human Genome Sequencing Center, Baylor College of Medicine, Houston, TX, 77030, USA
| | - Donna M Muzny
- Human Genome Sequencing Center, Baylor College of Medicine, Houston, TX, 77030, USA
| | - Richard A Gibbs
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, 77030, USA.,Human Genome Sequencing Center, Baylor College of Medicine, Houston, TX, 77030, USA
| | - Jacques Rougemont
- School of Life Sciences, EPFL (Ecole Polytechnique Fédérale de Lausanne), 1015, Lausanne, Switzerland.,Swiss Institute of Bioinformatics (SIB), 1015, Lausanne, Switzerland
| | - Ioannis Xenarios
- Center for Integrative Genomics, University of Lausanne, 1015, Lausanne, Switzerland.,Swiss Institute of Bioinformatics (SIB), 1015, Lausanne, Switzerland
| | - James R Lupski
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, 77030, USA.,Human Genome Sequencing Center, Baylor College of Medicine, Houston, TX, 77030, USA.,Department of Pediatrics, Baylor College of Medicine, Houston, TX, 77030, USA.,Texas Children's Hospital, Houston, TX, 77030, USA
| | - Alexandre Reymond
- Center for Integrative Genomics, University of Lausanne, 1015, Lausanne, Switzerland.
| |
Collapse
|
25
|
Valentini G, Armano G, Frasca M, Lin J, Mesiti M, Re M. RANKS: a flexible tool for node label ranking and classification in biological networks. Bioinformatics 2016; 32:2872-4. [PMID: 27256314 DOI: 10.1093/bioinformatics/btw235] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2016] [Accepted: 04/22/2016] [Indexed: 11/12/2022] Open
Abstract
UNLABELLED RANKS is a flexible software package that can be easily applied to any bioinformatics task formalizable as ranking of nodes with respect to a property given as a label, such as automated protein function prediction, gene disease prioritization and drug repositioning. To this end RANKS provides an efficient and easy-to-use implementation of kernelized score functions, a semi-supervised algorithmic scheme embedding both local and global learning strategies for the analysis of biomolecular networks. To facilitate comparative assessment, baseline network-based methods, e.g. label propagation and random walk algorithms, have also been implemented. AVAILABILITY AND IMPLEMENTATION The package is available from CRAN: https://cran.r-project.org/ The package is written in R, except for the most computationally intensive functionalities which are implemented in C. CONTACT valentini@di.unimi.it SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Giorgio Valentini
- AnacletoLab, Department of Computer Science, University of Milan, 20135 Milan, Italy
| | - Giuliano Armano
- Department of Electric and Electronic Engineering, University of Cagliari, 09123 Cagliari, Italy
| | - Marco Frasca
- AnacletoLab, Department of Computer Science, University of Milan, 20135 Milan, Italy
| | - Jianyi Lin
- AnacletoLab, Department of Computer Science, University of Milan, 20135 Milan, Italy
| | - Marco Mesiti
- AnacletoLab, Department of Computer Science, University of Milan, 20135 Milan, Italy
| | - Matteo Re
- AnacletoLab, Department of Computer Science, University of Milan, 20135 Milan, Italy
| |
Collapse
|
26
|
Baryshnikova A. Systematic Functional Annotation and Visualization of Biological Networks. Cell Syst 2016; 2:412-21. [DOI: 10.1016/j.cels.2016.04.014] [Citation(s) in RCA: 58] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2015] [Revised: 02/09/2016] [Accepted: 04/18/2016] [Indexed: 10/21/2022]
|
27
|
Lyons J, Herring CA, Banerjee A, Simmons AJ, Lau KS. Multiscale analysis of the murine intestine for modeling human diseases. Integr Biol (Camb) 2016; 7:740-57. [PMID: 26040649 DOI: 10.1039/c5ib00030k] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
When functioning properly, the intestine is one of the key interfaces between the human body and its environment. It is responsible for extracting nutrients from our food and excreting our waste products. It provides an environment for a host of healthful microbes and serves as a first defense against pathogenic ones. These processes require tight homeostatic controls, which are provided by the interactions of a complex mix of epithelial, stromal, neural and immune cells, as well as the resident microflora. This homeostasis can be disrupted by invasive microbes, genetic lesions, and carcinogens, resulting in diseases such Clostridium difficile infection, inflammatory bowel disease (IBD) and cancer. Enormous strides have been made in understanding how this important organ functions in health and disease using everything from cell culture systems to animal models to human tissue samples. This has resulted in better therapies for all of these diseases, but there is still significant room for improvement. In the United States alone, 14,000 people per year die of C. difficile, up to 1.6 million people suffer from IBD, and more than 50,000 people die every year from colon cancer. Because these and other intestinal diseases arise from complex interactions between the different components of the gut ecosystem, we propose that systems approaches that address this complexity in an integrative manner may eventually lead to improved therapeutics that deliver lasting cures. This review will discuss the use of systems biology for studying intestinal diseases in vivo with particular emphasis on mouse models. Additionally, it will focus on established experimental techniques that have been used to drive this systems-level analysis, and emerging techniques that will push this field forward in the future.
Collapse
Affiliation(s)
- Jesse Lyons
- Department of Biological Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | | | | | | | | |
Collapse
|
28
|
Blatti C, Sinha S. Characterizing gene sets using discriminative random walks with restart on heterogeneous biological networks. Bioinformatics 2016; 32:2167-75. [PMID: 27153592 PMCID: PMC4937193 DOI: 10.1093/bioinformatics/btw151] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2015] [Accepted: 03/14/2016] [Indexed: 12/11/2022] Open
Abstract
Motivation: Analysis of co-expressed gene sets typically involves testing for enrichment of different annotations or ‘properties’ such as biological processes, pathways, transcription factor binding sites, etc., one property at a time. This common approach ignores any known relationships among the properties or the genes themselves. It is believed that known biological relationships among genes and their many properties may be exploited to more accurately reveal commonalities of a gene set. Previous work has sought to achieve this by building biological networks that combine multiple types of gene–gene or gene–property relationships, and performing network analysis to identify other genes and properties most relevant to a given gene set. Most existing network-based approaches for recognizing genes or annotations relevant to a given gene set collapse information about different properties to simplify (homogenize) the networks. Results: We present a network-based method for ranking genes or properties related to a given gene set. Such related genes or properties are identified from among the nodes of a large, heterogeneous network of biological information. Our method involves a random walk with restarts, performed on an initial network with multiple node and edge types that preserve more of the original, specific property information than current methods that operate on homogeneous networks. In this first stage of our algorithm, we find the properties that are the most relevant to the given gene set and extract a subnetwork of the original network, comprising only these relevant properties. We then re-rank genes by their similarity to the given gene set, based on a second random walk with restarts, performed on the above subnetwork. We demonstrate the effectiveness of this algorithm for ranking genes related to Drosophila embryonic development and aggressive responses in the brains of social animals. Availability and Implementation: DRaWR was implemented as an R package available at veda.cs.illinois.edu/DRaWR. Contact: blatti@illinois.edu Supplementary information:Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Charles Blatti
- Department of Computer Science, University of Illinois, Urbana, IL 61801, USA
| | - Saurabh Sinha
- Department of Computer Science, University of Illinois, Urbana, IL 61801, USA Institute of Genomic Biology, University of Illinois, Urbana, IL 61801, USA
| |
Collapse
|
29
|
Abstract
Motivation: Markov networks are undirected graphical models that are widely used to infer relations between genes from experimental data. Their state-of-the-art inference procedures assume the data arise from a Gaussian distribution. High-throughput omics data, such as that from next generation sequencing, often violates this assumption. Furthermore, when collected data arise from multiple related but otherwise nonidentical distributions, their underlying networks are likely to have common features. New principled statistical approaches are needed that can deal with different data distributions and jointly consider collections of datasets. Results: We present FuseNet, a Markov network formulation that infers networks from a collection of nonidentically distributed datasets. Our approach is computationally efficient and general: given any number of distributions from an exponential family, FuseNet represents model parameters through shared latent factors that define neighborhoods of network nodes. In a simulation study, we demonstrate good predictive performance of FuseNet in comparison to several popular graphical models. We show its effectiveness in an application to breast cancer RNA-sequencing and somatic mutation data, a novel application of graphical models. Fusion of datasets offers substantial gains relative to inference of separate networks for each dataset. Our results demonstrate that network inference methods for non-Gaussian data can help in accurate modeling of the data generated by emergent high-throughput technologies. Availability and implementation: Source code is at https://github.com/marinkaz/fusenet. Contact:blaz.zupan@fri.uni-lj.si Supplementary information:Supplementary information is available at Bioinformatics online.
Collapse
Affiliation(s)
- Marinka Žitnik
- Faculty of Computer and Information Science, University of Ljubljana, Ljubljana, Slovenia and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, USA
| | - Blaž Zupan
- Faculty of Computer and Information Science, University of Ljubljana, Ljubljana, Slovenia and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, USA Faculty of Computer and Information Science, University of Ljubljana, Ljubljana, Slovenia and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, USA
| |
Collapse
|
30
|
Cornish AJ, Filippis I, David A, Sternberg MJE. Exploring the cellular basis of human disease through a large-scale mapping of deleterious genes to cell types. Genome Med 2015; 7:95. [PMID: 26330083 PMCID: PMC4557825 DOI: 10.1186/s13073-015-0212-9] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2015] [Accepted: 07/31/2015] [Indexed: 01/09/2023] Open
Abstract
BACKGROUND Each cell type found within the human body performs a diverse and unique set of functions, the disruption of which can lead to disease. However, there currently exists no systematic mapping between cell types and the diseases they can cause. METHODS In this study, we integrate protein-protein interaction data with high-quality cell-type-specific gene expression data from the FANTOM5 project to build the largest collection of cell-type-specific interactomes created to date. We develop a novel method, called gene set compactness (GSC), that contrasts the relative positions of disease-associated genes across 73 cell-type-specific interactomes to map genes associated with 196 diseases to the cell types they affect. We conduct text-mining of the PubMed database to produce an independent resource of disease-associated cell types, which we use to validate our method. RESULTS The GSC method successfully identifies known disease-cell-type associations, as well as highlighting associations that warrant further study. This includes mast cells and multiple sclerosis, a cell population currently being targeted in a multiple sclerosis phase 2 clinical trial. Furthermore, we build a cell-type-based diseasome using the cell types identified as manifesting each disease, offering insight into diseases linked through etiology. CONCLUSIONS The data set produced in this study represents the first large-scale mapping of diseases to the cell types in which they are manifested and will therefore be useful in the study of disease systems. Overall, we demonstrate that our approach links disease-associated genes to the phenotypes they produce, a key goal within systems medicine.
Collapse
Affiliation(s)
- Alex J Cornish
- Department of Life Sciences, Imperial College London, Exhibition Road, London, SW7 2AZ, UK.
| | - Ioannis Filippis
- Department of Life Sciences, Imperial College London, Exhibition Road, London, SW7 2AZ, UK.
| | - Alessia David
- Department of Life Sciences, Imperial College London, Exhibition Road, London, SW7 2AZ, UK.
| | - Michael J E Sternberg
- Department of Life Sciences, Imperial College London, Exhibition Road, London, SW7 2AZ, UK.
| |
Collapse
|
31
|
Systems Analyses Reveal Shared and Diverse Attributes of Oct4 Regulation in Pluripotent Cells. Cell Syst 2015; 1:141-51. [PMID: 27135800 DOI: 10.1016/j.cels.2015.08.002] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2015] [Revised: 07/07/2015] [Accepted: 08/04/2015] [Indexed: 11/23/2022]
Abstract
We combine a genome-scale RNAi screen in mouse epiblast stem cells (EpiSCs) with genetic interaction, protein localization, and "protein-level dependency" studies-a systematic technique that uncovers post-transcriptional regulation-to delineate the network of factors that control the expression of Oct4, a key regulator of pluripotency. Our data signify that there are similarities, but also fundamental differences in Oct4 regulation in EpiSCs versus embryonic stem cells (ESCs). Through multiparametric data analyses, we predict that Tox4 is associating with the Paf1C complex, which maintains cell identity in both cell types, and validate that this protein-protein interaction exists in ESCs and EpiSCs. We also identify numerous knockdowns that increase Oct4 expression in EpiSCs, indicating that, in stark contrast to ESCs, Oct4 is under active repressive control in EpiSCs. These studies provide a framework for better understanding pluripotency and for dissecting the molecular events that govern the transition from the pre-implantation to the post-implantation state.
Collapse
|
32
|
Oh EY, Christensen SM, Ghanta S, Jeong JC, Bucur O, Glass B, Montaser-Kouhsari L, Knoblauch NW, Bertos N, Saleh SM, Haibe-Kains B, Park M, Beck AH. Extensive rewiring of epithelial-stromal co-expression networks in breast cancer. Genome Biol 2015; 16:128. [PMID: 26087699 PMCID: PMC4471934 DOI: 10.1186/s13059-015-0675-4] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2015] [Accepted: 05/13/2015] [Indexed: 02/06/2023] Open
Abstract
BACKGROUND Epithelial-stromal crosstalk plays a critical role in invasive breast cancer pathogenesis; however, little is known on a systems level about how epithelial-stromal interactions evolve during carcinogenesis. RESULTS We develop a framework for building genome-wide epithelial-stromal co-expression networks composed of pairwise co-expression relationships between mRNA levels of genes expressed in the epithelium and stroma across a population of patients. We apply this method to laser capture micro-dissection expression profiling datasets in the setting of breast carcinogenesis. Our analysis shows that epithelial-stromal co-expression networks undergo extensive rewiring during carcinogenesis, with the emergence of distinct network hubs in normal breast, and estrogen receptor-positive and estrogen receptor-negative invasive breast cancer, and the emergence of distinct patterns of functional network enrichment. In contrast to normal breast, the strongest epithelial-stromal co-expression relationships in invasive breast cancer mostly represent self-loops, in which the same gene is co-expressed in epithelial and stromal regions. We validate this observation using an independent laser capture micro-dissection dataset and confirm that self-loop interactions are significantly increased in cancer by performing computational image analysis of epithelial and stromal protein expression using images from the Human Protein Atlas. CONCLUSIONS Epithelial-stromal co-expression network analysis represents a new approach for systems-level analyses of spatially localized transcriptomic data. The analysis provides new biological insights into the rewiring of epithelial-stromal co-expression networks and the emergence of epithelial-stromal co-expression self-loops in breast cancer. The approach may facilitate the development of new diagnostics and therapeutics targeting epithelial-stromal interactions in cancer.
Collapse
Affiliation(s)
- Eun-Yeong Oh
- Cancer Research Institute, Beth Israel Deaconess Cancer Center, Boston, MA, 02215, USA. .,Department of Pathology, Beth Israel Deaconess Medical Center, Boston, MA, 02215, USA. .,Harvard Medical School, Boston, MA, 02215, USA.
| | - Stephen M Christensen
- Cancer Research Institute, Beth Israel Deaconess Cancer Center, Boston, MA, 02215, USA. .,Department of Pathology, Beth Israel Deaconess Medical Center, Boston, MA, 02215, USA. .,Harvard Medical School, Boston, MA, 02215, USA.
| | - Sindhu Ghanta
- Cancer Research Institute, Beth Israel Deaconess Cancer Center, Boston, MA, 02215, USA. .,Department of Pathology, Beth Israel Deaconess Medical Center, Boston, MA, 02215, USA. .,Harvard Medical School, Boston, MA, 02215, USA.
| | - Jong Cheol Jeong
- Cancer Research Institute, Beth Israel Deaconess Cancer Center, Boston, MA, 02215, USA. .,Department of Pathology, Beth Israel Deaconess Medical Center, Boston, MA, 02215, USA. .,Harvard Medical School, Boston, MA, 02215, USA.
| | - Octavian Bucur
- Cancer Research Institute, Beth Israel Deaconess Cancer Center, Boston, MA, 02215, USA. .,Department of Pathology, Beth Israel Deaconess Medical Center, Boston, MA, 02215, USA. .,Harvard Medical School, Boston, MA, 02215, USA.
| | - Benjamin Glass
- Cancer Research Institute, Beth Israel Deaconess Cancer Center, Boston, MA, 02215, USA. .,Department of Pathology, Beth Israel Deaconess Medical Center, Boston, MA, 02215, USA. .,Harvard Medical School, Boston, MA, 02215, USA.
| | - Laleh Montaser-Kouhsari
- Cancer Research Institute, Beth Israel Deaconess Cancer Center, Boston, MA, 02215, USA. .,Department of Pathology, Beth Israel Deaconess Medical Center, Boston, MA, 02215, USA. .,Harvard Medical School, Boston, MA, 02215, USA.
| | - Nicholas W Knoblauch
- Cancer Research Institute, Beth Israel Deaconess Cancer Center, Boston, MA, 02215, USA. .,Department of Pathology, Beth Israel Deaconess Medical Center, Boston, MA, 02215, USA. .,Harvard Medical School, Boston, MA, 02215, USA.
| | - Nicholas Bertos
- Goodman Cancer Research Centre, McGill University, Montreal, QC, Canada.
| | - Sadiq Mi Saleh
- Goodman Cancer Research Centre, McGill University, Montreal, QC, Canada.
| | - Benjamin Haibe-Kains
- Princess Margaret Cancer Centre, University Health Network, Toronto, ON, M5G 1L7, Canada. .,Department of Medical Biophysics, University of Toronto, Toronto, Ontario, M5G 1L7, Canada.
| | - Morag Park
- Goodman Cancer Research Centre, McGill University, Montreal, QC, Canada.
| | - Andrew H Beck
- Cancer Research Institute, Beth Israel Deaconess Cancer Center, Boston, MA, 02215, USA. .,Department of Pathology, Beth Israel Deaconess Medical Center, Boston, MA, 02215, USA. .,Harvard Medical School, Boston, MA, 02215, USA.
| |
Collapse
|
33
|
Di Lena P, Martelli PL, Fariselli P, Casadio R. NET-GE: a novel NETwork-based Gene Enrichment for detecting biological processes associated to Mendelian diseases. BMC Genomics 2015; 16 Suppl 8:S6. [PMID: 26110971 PMCID: PMC4480278 DOI: 10.1186/1471-2164-16-s8-s6] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Abstract
Background Enrichment analysis is a widely applied procedure for shedding light on the molecular mechanisms and functions at the basis of phenotypes, for enlarging the dataset of possibly related genes/proteins and for helping interpretation and prioritization of newly determined variations. Several standard and Network-based enrichment methods are available. Both approaches rely on the annotations that characterize the genes/proteins included in the input set; network based ones also include in different ways physical and functional relationships among different genes or proteins that can be extracted from the available biological networks of interactions. Results Here we describe a novel procedure based on the extraction from the STRING interactome of sub-networks connecting proteins that share the same Gene Ontology(GO) terms for Biological Process (BP). Enrichment analysis is performed by mapping the protein set to be analyzed on the sub-networks, and then by collecting the corresponding annotations. We test the ability of our enrichment method in finding annotation terms disregarded by other enrichment methods available. We benchmarked 244 sets of proteins associated to different Mendelian diseases, according to the OMIM web resource. In 143 cases (58%), the network-based procedure extracts GO terms neglected by the standard method, and in 86 cases (35%), some of the newly enriched GO terms are not included in the set of annotations characterizing the input proteins. We present in detail six cases where our network-based enrichment provides an insight into the biological basis of the diseases, outperforming other freely available network-based methods. Conclusions Considering a set of proteins in the context of their interaction network can help in better defining their functions. Our novel method exploits the information contained in the STRING database for building the minimal connecting network containing all the proteins annotated with the same GO term. The enrichment procedure is performed considering the GO-specific network modules and, when tested on the OMIM-derived benchmark sets, it is able to extract enrichment terms neglected by other methods. Our procedure is effective even when the size of the input protein set is small, requiring at least two input proteins.
Collapse
|