1
|
Hossain I, Fanfani V, Fischer J, Quackenbush J, Burkholz R. Biologically informed NeuralODEs for genome-wide regulatory dynamics. Genome Biol 2024; 25:127. [PMID: 38773638 DOI: 10.1186/s13059-024-03264-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2023] [Accepted: 04/30/2024] [Indexed: 05/24/2024] Open
Abstract
BACKGROUND Gene regulatory network (GRN) models that are formulated as ordinary differential equations (ODEs) can accurately explain temporal gene expression patterns and promise to yield new insights into important cellular processes, disease progression, and intervention design. Learning such gene regulatory ODEs is challenging, since we want to predict the evolution of gene expression in a way that accurately encodes the underlying GRN governing the dynamics and the nonlinear functional relationships between genes. Most widely used ODE estimation methods either impose too many parametric restrictions or are not guided by meaningful biological insights, both of which impede either scalability, explainability, or both. RESULTS We developed PHOENIX, a modeling framework based on neural ordinary differential equations (NeuralODEs) and Hill-Langmuir kinetics, that overcomes limitations of other methods by flexibly incorporating prior domain knowledge and biological constraints to promote sparse, biologically interpretable representations of GRN ODEs. We tested the accuracy of PHOENIX in a series of in silico experiments, benchmarking it against several currently used tools. We demonstrated PHOENIX's flexibility by modeling regulation of oscillating expression profiles obtained from synchronized yeast cells. We also assessed the scalability of PHOENIX by modeling genome-scale GRNs for breast cancer samples ordered in pseudotime and for B cells treated with Rituximab. CONCLUSIONS PHOENIX uses a combination of user-defined prior knowledge and functional forms from systems biology to encode biological "first principles" as soft constraints on the GRN allowing us to predict subsequent gene expression patterns in a biologically explainable manner.
Collapse
Affiliation(s)
| | - Viola Fanfani
- Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | - Jonas Fischer
- Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | | | - Rebekka Burkholz
- CISPA Helmholtz Center for Information Security, Saarbrücken, Germany
| |
Collapse
|
2
|
Stone K, Platig J, Quackenbush J, Fagny M. Complex Traits Heritability is Highly Clustered in the eQTL Bipartite Network. bioRxiv 2024:2024.02.27.582063. [PMID: 38464142 PMCID: PMC10925220 DOI: 10.1101/2024.02.27.582063] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/12/2024]
Abstract
Single Nucleotide Polymorphisms (SNPs) associated with traits typically explain a small part of the trait genetic heritability-with the remainder thought to be distributed throughout the genome. Such SNPs are likely to alter expression levels of biologically relevant genes. Expression Quantitative Trait Locus (eQTL) networks analysis has helped to functionally characterize such variants. We systematically analyze the distribution of SNP heritability for ten traits across 29 tissue-specific eQTL networks. We find that heritability is clustered in a small number or tissue-specific, functionally relevant SNP-gene modules and that the greatest occurs in local "hubs" that are both the cornerstone of the network's modules and tissue-specific regulatory elements. The network structure could thus both amplify the genotype-phenotype connection and buffer the deleterious effect of the genetic variations on other traits. Together, these results define a conceptual framework for understanding complex trait architecture and identifying key mutations carrying most of the heritability.
Collapse
Affiliation(s)
- Katherine Stone
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, United States
- Department of Data Science and Center for Cancer Computational Biology, Dana-Farber Cancer Institute, Boston, Massachusetts, USA
| | - John Platig
- Center for Public Health Genomics, University of Virginia, Charlottesville, Virginia, USA
- Department of Public Health Sciences, University of Virginia, Charlottesville, Virginia, USA
- Department of Biomedical Engineering, University of Virginia, Charlottesville, Virginia, USA
| | - John Quackenbush
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, United States
- Department of Data Science and Center for Cancer Computational Biology, Dana-Farber Cancer Institute, Boston, Massachusetts, USA
- Channing Division of Network Medicine, Brigham and Women's Hospital, Boston, Massachusetts, United States
| | - Maud Fagny
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, United States
- Department of Data Science and Center for Cancer Computational Biology, Dana-Farber Cancer Institute, Boston, Massachusetts, USA
- Université Paris-Saclay, INRAE, CNRS, AgroParisTech, Genetique Quantitative et Evolution - Le Moulon, Gif-sur-Yvette 91190 France
| |
Collapse
|
3
|
Hossain I, Fanfani V, Fischer J, Quackenbush J, Burkholz R. Biologically informed NeuralODEs for genome-wide regulatory dynamics. bioRxiv 2024:2023.02.24.529835. [PMID: 36909563 PMCID: PMC10002636 DOI: 10.1101/2023.02.24.529835] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
Modeling dynamics of gene regulatory networks using ordinary differential equations (ODEs) allow a deeper understanding of disease progression and response to therapy, thus aiding in intervention optimization. Although there exist methods to infer regulatory ODEs, these are generally limited to small networks, rely on dimensional reduction, or impose non-biological parametric restrictions - all impeding scalability and explainability. PHOENIX is a neural ODE framework incorporating prior domain knowledge as soft constraints to infer sparse, biologically interpretable dynamics. Extensive experiments - on simulated and real data - demonstrate PHOENIX's unique ability to learn key regulatory dynamics while scaling to the whole genome.
Collapse
|
4
|
Saha E, Fanfani V, Mandros P, Ben-Guebila M, Fischer J, Hoff-Shutta K, Glass K, DeMeo DL, Lopes-Ramos C, Quackenbush J. Bayesian Optimized sample-specific Networks Obtained By Omics data (BONOBO). bioRxiv 2023:2023.11.16.567119. [PMID: 38014256 PMCID: PMC10680741 DOI: 10.1101/2023.11.16.567119] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/29/2023]
Abstract
Gene regulatory networks (GRNs) are effective tools for inferring complex interactions between molecules that regulate biological processes and hence can provide insights into drivers of biological systems. Inferring co-expression networks is a critical element of GRN inference as the correlation between expression patterns may indicate that genes are coregulated by common factors. However, methods that estimate co-expression networks generally derive an aggregate network representing the mean regulatory properties of the population and so fail to fully capture population heterogeneity. To address these concerns, we introduce BONOBO (Bayesian Optimized Networks Obtained By assimilating Omics data), a scalable Bayesian model for deriving individual sample-specific co-expression networks by recognizing variations in molecular interactions across individuals. For every sample, BONOBO assumes a Gaussian distribution on the log-transformed centered gene expression and a conjugate prior distribution on the sample-specific co-expression matrix constructed from all other samples in the data. Combining the sample-specific gene expression with the prior distribution, BONOBO yields a closed-form solution for the posterior distribution of the sample-specific co-expression matrices, thus making the method extremely scalable. We demonstrate the utility of BONOBO in several contexts, including analyzing gene regulation in yeast transcription factor knockout studies, prognostic significance of miRNA-mRNA interaction in human breast cancer subtypes, and sex differences in gene regulation within human thyroid tissue. We find that BONOBO outperforms other sample-specific co-expression network inference methods and provides insight into individual differences in the drivers of biological processes.
Collapse
Affiliation(s)
- Enakshi Saha
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA
| | - Viola Fanfani
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA
| | - Panagiotis Mandros
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA
| | - Marouen Ben-Guebila
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA
| | - Jonas Fischer
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA
| | - Katherine Hoff-Shutta
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA
- Channing Division of Network Medicine, Brigham and Women’s Hospital, Boston, MA, USA
| | - Kimberly Glass
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA
- Channing Division of Network Medicine, Brigham and Women’s Hospital, Boston, MA, USA
- Department of Medicine, Harvard Medical School, Boston, MA, USA
| | - Dawn Lisa DeMeo
- Channing Division of Network Medicine, Brigham and Women’s Hospital, Boston, MA, USA
- Department of Medicine, Harvard Medical School, Boston, MA, USA
| | - Camila Lopes-Ramos
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA
- Channing Division of Network Medicine, Brigham and Women’s Hospital, Boston, MA, USA
- Department of Medicine, Harvard Medical School, Boston, MA, USA
| | - John Quackenbush
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA
- Channing Division of Network Medicine, Brigham and Women’s Hospital, Boston, MA, USA
- Department of Data Science, Dana-Farber Cancer Institute, Boston, MA, USA
| |
Collapse
|
5
|
Saha E, Guebila MB, Fanfani V, Fischer J, Shutta KH, Mandros P, DeMeo DL, Quackenbush J, Lopes-Ramos CM. Gene regulatory Networks Reveal Sex Difference in Lung Adenocarcinoma. bioRxiv 2023:2023.09.22.559001. [PMID: 37790409 PMCID: PMC10543009 DOI: 10.1101/2023.09.22.559001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
Abstract
Lung adenocarcinoma (LUAD) has been observed to have significant sex differences in incidence, prognosis, and response to therapy. However, the molecular mechanisms responsible for these disparities have not been investigated extensively. Sample-specific gene regulatory network methods were used to analyze RNA sequencing data from non-cancerous human lung samples from The Genotype Tissue Expression Project (GTEx) and lung adenocarcinoma primary tumor samples from The Cancer Genome Atlas (TCGA); results were validated on independent data. We observe that genes associated with key biological pathways including cell proliferation, immune response and drug metabolism are differentially regulated between males and females in both healthy lung tissue, as well as in tumor, and that these regulatory differences are further perturbed by tobacco smoking. We also uncovered significant sex bias in transcription factor targeting patterns of clinically actionable oncogenes and tumor suppressor genes, including AKT2 and KRAS. Using differentially regulated genes between healthy and tumor samples in conjunction with a drug repurposing tool, we identified several small-molecule drugs that might have sex-biased efficacy as cancer therapeutics and further validated this observation using an independent cell line database. These findings underscore the importance of including sex as a biological variable and considering gene regulatory processes in developing strategies for disease prevention and management.
Collapse
Affiliation(s)
- Enakshi Saha
- Department of Biostatistics, Harvard T. H. Chan School of Public Health, Boston, MA 02115, USA
| | - Marouen Ben Guebila
- Department of Biostatistics, Harvard T. H. Chan School of Public Health, Boston, MA 02115, USA
| | - Viola Fanfani
- Department of Biostatistics, Harvard T. H. Chan School of Public Health, Boston, MA 02115, USA
| | - Jonas Fischer
- Department of Biostatistics, Harvard T. H. Chan School of Public Health, Boston, MA 02115, USA
| | - Katherine H Shutta
- Department of Biostatistics, Harvard T. H. Chan School of Public Health, Boston, MA 02115, USA
- Channing Division of Network Medicine, Brigham and Women's Hospital, Boston, MA, USA 02115
| | - Panagiotis Mandros
- Department of Biostatistics, Harvard T. H. Chan School of Public Health, Boston, MA 02115, USA
| | - Dawn L DeMeo
- Channing Division of Network Medicine, Brigham and Women's Hospital, Boston, MA, USA 02115
- Department of Medicine, Harvard Medical School, Boston, MA 02115, USA
| | - John Quackenbush
- Department of Biostatistics, Harvard T. H. Chan School of Public Health, Boston, MA 02115, USA
- Channing Division of Network Medicine, Brigham and Women's Hospital, Boston, MA, USA 02115
- Department of Data Science, Dana-Farber Cancer Institute, Boston, MA 02115, USA
| | - Camila M Lopes-Ramos
- Department of Biostatistics, Harvard T. H. Chan School of Public Health, Boston, MA 02115, USA
- Channing Division of Network Medicine, Brigham and Women's Hospital, Boston, MA, USA 02115
- Department of Medicine, Harvard Medical School, Boston, MA 02115, USA
| |
Collapse
|
6
|
Hossain I, Fanfani V, Quackenbush J, Burkholz R. Biologically informed NeuralODEs for genome-wide regulatory dynamics. Res Sq 2023:rs.3.rs-2675584. [PMID: 36993392 PMCID: PMC10055646 DOI: 10.21203/rs.3.rs-2675584/v1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Indexed: 06/19/2023]
Abstract
Models that are formulated as ordinary differential equations (ODEs) can accurately explain temporal gene expression patterns and promise to yield new insights into important cellular processes, disease progression, and intervention design. Learning such ODEs is challenging, since we want to predict the evolution of gene expression in a way that accurately encodes the causal gene-regulatory network (GRN) governing the dynamics and the nonlinear functional relationships between genes. Most widely used ODE estimation methods either impose too many parametric restrictions or are not guided by meaningful biological insights, both of which impedes scalability and/or explainability. To overcome these limitations, we developed PHOENIX, a modeling framework based on neural ordinary differential equations (NeuralODEs) and Hill-Langmuir kinetics, that can flexibly incorporate prior domain knowledge and biological constraints to promote sparse, biologically interpretable representations of ODEs. We test accuracy of PHOENIX in a series of in silico experiments benchmarking it against several currently used tools for ODE estimation. We also demonstrate PHOENIX's flexibility by studying oscillating expression data from synchronized yeast cells and assess its scalability by modelling genome-scale breast cancer expression for samples ordered in pseudotime. Finally, we show how the combination of user-defined prior knowledge and functional forms from systems biology allows PHOENIX to encode key properties of the underlying GRN, and subsequently predict expression patterns in a biologically explainable way.
Collapse
Affiliation(s)
- Intekhab Hossain
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | - Viola Fanfani
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | - John Quackenbush
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | - Rebekka Burkholz
- Helmholtz Center for Information Security (CISPA), Saarbrücken, Germany
| |
Collapse
|
7
|
Ben Guebila M, Wang T, Lopes-Ramos CM, Fanfani V, Weighill D, Burkholz R, Schlauch D, Paulson JN, Altenbuchinger M, Shutta KH, Sonawane AR, Lim J, Calderer G, van IJzendoorn DGP, Morgan D, Marin A, Chen CY, Song Q, Saha E, DeMeo DL, Padi M, Platig J, Kuijjer ML, Glass K, Quackenbush J. The Network Zoo: a multilingual package for the inference and analysis of gene regulatory networks. Genome Biol 2023; 24:45. [PMID: 36894939 PMCID: PMC9999668 DOI: 10.1186/s13059-023-02877-1] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2022] [Accepted: 02/15/2023] [Indexed: 03/11/2023] Open
Abstract
Inference and analysis of gene regulatory networks (GRNs) require software that integrates multi-omic data from various sources. The Network Zoo (netZoo; netzoo.github.io) is a collection of open-source methods to infer GRNs, conduct differential network analyses, estimate community structure, and explore the transitions between biological states. The netZoo builds on our ongoing development of network methods, harmonizing the implementations in various computing languages and between methods to allow better integration of these tools into analytical pipelines. We demonstrate the utility using multi-omic data from the Cancer Cell Line Encyclopedia. We will continue to expand the netZoo to incorporate additional methods.
Collapse
Affiliation(s)
- Marouen Ben Guebila
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | - Tian Wang
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
- Present Address: Biology Department, Boston College, Chestnut Hill, MA, USA
| | - Camila M Lopes-Ramos
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
- Channing Division of Network Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA
| | - Viola Fanfani
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | - Des Weighill
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
- Present Address: Lineberger Comprehensive Cancer Center, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA
| | - Rebekka Burkholz
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
- Present Address: CISPA Helmholtz Center for Information Security, Saarbrücken, Germany
| | - Daniel Schlauch
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
- Present Address: Genospace, LLC, Boston, MA, USA
| | - Joseph N Paulson
- Department of Biochemistry and Molecular Biology, Pennsylvania State University College of Medicine, Hershey, PA, USA
| | - Michael Altenbuchinger
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
- Present Address: Department of Medical Bioinformatics, University Medical Center Göttingen, Göttingen, Germany
| | - Katherine H Shutta
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
- Channing Division of Network Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA
| | - Abhijeet R Sonawane
- Channing Division of Network Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA
- Present Address: Center for Interdisciplinary Cardiovascular Sciences, Division of Cardiovascular Medicine, Department of Medicine, Brigham and Women's Hospital, Boston, MA, USA
| | - James Lim
- Department of Molecular and Cellular Biology, University of Arizona, Tucson, AZ, USA
- Present Address: Monoceros Biosystems, LLC, San Diego, CA, USA
| | - Genis Calderer
- Center for Molecular Medicine Norway, Nordic EMBL Partnership, University of Oslo, Oslo, Norway
| | - David G P van IJzendoorn
- Department of Pathology, Leiden University Medical Center, Leiden, The Netherlands
- Present Address: Department of Pathology, Stanford University School of Medicine, Palo Alto, CA, USA
| | - Daniel Morgan
- Channing Division of Network Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA
- Present Address: School of Biomedical Sciences, Hong Kong University, Pokfulam, Hong Kong
| | | | - Cho-Yi Chen
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
- Dana-Farber Cancer Institute, Boston, MA, USA
- Present Address: Institute of Biomedical Informatics, National Yang Ming Chiao Tung University, Taipei, 112, Taiwan
| | - Qi Song
- Channing Division of Network Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA
- Present Address: Computational Biology Department, Carnegie Mellon University, Pittsburgh, PA, USA
| | - Enakshi Saha
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | - Dawn L DeMeo
- Channing Division of Network Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA
| | - Megha Padi
- Department of Molecular and Cellular Biology, University of Arizona, Tucson, AZ, USA
| | - John Platig
- Channing Division of Network Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA
| | - Marieke L Kuijjer
- Center for Molecular Medicine Norway, Nordic EMBL Partnership, University of Oslo, Oslo, Norway
- Department of Pathology, Leiden University Medical Center, Leiden, The Netherlands
- Leiden Center for Computational Oncology, Leiden University, Leiden, The Netherlands
| | - Kimberly Glass
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
- Channing Division of Network Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA
| | - John Quackenbush
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA.
- Channing Division of Network Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA.
- Dana-Farber Cancer Institute, Boston, MA, USA.
| |
Collapse
|
8
|
Shutta KH, Weighill D, Burkholz R, Guebila M, DeMeo DL, Zacharias HU, Quackenbush J, Altenbuchinger M. DRAGON: Determining Regulatory Associations using Graphical models on multi-Omic Networks. Nucleic Acids Res 2022; 51:e15. [PMID: 36533448 PMCID: PMC9943674 DOI: 10.1093/nar/gkac1157] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2022] [Revised: 11/08/2022] [Accepted: 11/23/2022] [Indexed: 12/23/2022] Open
Abstract
The increasing quantity of multi-omic data, such as methylomic and transcriptomic profiles collected on the same specimen or even on the same cell, provides a unique opportunity to explore the complex interactions that define cell phenotype and govern cellular responses to perturbations. We propose a network approach based on Gaussian Graphical Models (GGMs) that facilitates the joint analysis of paired omics data. This method, called DRAGON (Determining Regulatory Associations using Graphical models on multi-Omic Networks), calibrates its parameters to achieve an optimal trade-off between the network's complexity and estimation accuracy, while explicitly accounting for the characteristics of each of the assessed omics 'layers.' In simulation studies, we show that DRAGON adapts to edge density and feature size differences between omics layers, improving model inference and edge recovery compared to state-of-the-art methods. We further demonstrate in an analysis of joint transcriptome - methylome data from TCGA breast cancer specimens that DRAGON can identify key molecular mechanisms such as gene regulation via promoter methylation. In particular, we identify Transcription Factor AP-2 Beta (TFAP2B) as a potential multi-omic biomarker for basal-type breast cancer. DRAGON is available as open-source code in Python through the Network Zoo package (netZooPy v0.8; netzoo.github.io).
Collapse
Affiliation(s)
| | | | - Rebekka Burkholz
- CISPA Helmholtz Center for Information Security, Saarbrücken, Germany
| | - Marouen Ben Guebila
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | - Dawn L DeMeo
- Channing Division of Network Medicine, Brigham and Women’s Hospital, and Department of Medicine, Harvard Medical School, Boston, MA, USA
| | - Helena U Zacharias
- Department of Internal Medicine I, University Medical Center Schleswig-Holstein, Campus Kiel, Kiel, Germany,Institute of Clinical Molecular Biology, Kiel University and University Medical Center Schleswig-Holstein, Campus Kiel, Kiel, Germany,Peter L. Reichertz Institute for Medical Informatics of TU Braunschweig and Hannover Medical School, Hannover, Germany
| | | | - Michael Altenbuchinger
- To whom correspondence should be addressed. Tel: +49 551 39 61788; Fax: +49 551 39 61783;
| |
Collapse
|
9
|
Yao S, Campbell PT, Ugai T, Gierach G, Abubakar M, Adalsteinsson V, Almeida J, Brennan P, Chanock S, Golub T, Hanash S, Harris C, Hathaway CA, Kelsey K, Landi MT, Mahmood F, Newton C, Quackenbush J, Rodig S, Schultz N, Tearney G, Tworoger SS, Wang M, Zhang X, Garcia-Closas M, Rebbeck TR, Ambrosone CB, Ogino S. Proceedings of the fifth international Molecular Pathological Epidemiology (MPE) meeting. Cancer Causes Control 2022; 33:1107-1120. [PMID: 35759080 PMCID: PMC9244289 DOI: 10.1007/s10552-022-01594-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2021] [Accepted: 05/20/2022] [Indexed: 01/19/2023]
Abstract
Cancer heterogeneities hold the key to a deeper understanding of cancer etiology and progression and the discovery of more precise cancer therapy. Modern pathological and molecular technologies offer a powerful set of tools to profile tumor heterogeneities at multiple levels in large patient populations, from DNA to RNA, protein and epigenetics, and from tumor tissues to tumor microenvironment and liquid biopsy. When coupled with well-validated epidemiologic methodology and well-characterized epidemiologic resources, the rich tumor pathological and molecular tumor information provide new research opportunities at an unprecedented breadth and depth. This is the research space where Molecular Pathological Epidemiology (MPE) emerged over a decade ago and has been thriving since then. As a truly multidisciplinary field, MPE embraces collaborations from diverse fields including epidemiology, pathology, immunology, genetics, biostatistics, bioinformatics, and data science. Since first convened in 2013, the International MPE Meeting series has grown into a dynamic and dedicated platform for experts from these disciplines to communicate novel findings, discuss new research opportunities and challenges, build professional networks, and educate the next-generation scientists. Herein, we share the proceedings of the Fifth International MPE meeting, held virtually online, on May 24 and 25, 2021. The meeting consisted of 21 presentations organized into the three main themes, which were recent integrative MPE studies, novel cancer profiling technologies, and new statistical and data science approaches. Looking forward to the near future, the meeting attendees anticipated continuous expansion and fruition of MPE research in many research fronts, particularly immune-epidemiology, mutational signatures, liquid biopsy, and health disparities.
Collapse
Affiliation(s)
- Song Yao
- Department of Cancer Prevention and Control, Roswell Park Comprehensive Cancer Center, Elm & Carlton Streets, Buffalo, NY, 14263, USA.
| | - Peter T Campbell
- Department of Epidemiology and Population Health, Albert Einstein College of Medicine, Bronx, NY, USA
| | - Tomotaka Ugai
- Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | - Gretchen Gierach
- Division of Cancer Epidemiology and Genetics, National Cancer Institute, Bethesda, MD, USA
| | - Mustapha Abubakar
- Division of Cancer Epidemiology and Genetics, National Cancer Institute, Bethesda, MD, USA
| | | | - Jonas Almeida
- Division of Cancer Epidemiology and Genetics, National Cancer Institute, Bethesda, MD, USA
| | - Paul Brennan
- International Agency for Research On Cancer (IARC/WHO), Genomic Epidemiology Branch, Lyon, France
| | - Stephen Chanock
- Division of Cancer Epidemiology and Genetics, National Cancer Institute, Bethesda, MD, USA
| | - Todd Golub
- Broad Institute of MIT and Harvard, Boston, MA, USA
- Department of Pediatric Oncology, Dana-Farber Cancer Institute, Boston, MA, USA
| | - Samir Hanash
- Department of Clinical Cancer Prevention, MD Anderson Cancer Institute, Houston, TX, USA
| | - Curtis Harris
- Center for Cancer Research, National Cancer Institute, Bethesda, MD, USA
| | - Cassandra A Hathaway
- Department of Cancer Epidemiology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, USA
| | - Karl Kelsey
- Department of Epidemiology, Brown School of Public Health, Brown University, Providence, RI, USA
| | - Maria Teresa Landi
- Division of Cancer Epidemiology and Genetics, National Cancer Institute, Bethesda, MD, USA
| | - Faisal Mahmood
- Department of Pathology, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA
| | - Christina Newton
- Department of Population Science, American Cancer Society, Atlanta, GA, USA
| | - John Quackenbush
- Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital, Harvard Medical School, Boston, MA, USA
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | - Scott Rodig
- Department of Pathology, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA
| | - Nikolaus Schultz
- Department of Epidemiology and Biostatistics, Memorial Sloan Kettering Cancer Center, New York, NY, USA
| | - Guillermo Tearney
- Department of Pathology and Wellman Center for Photomedicine, Massachusetts General Hospital, Boston, MA, USA
| | - Shelley S Tworoger
- Department of Cancer Epidemiology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, USA
| | - Molin Wang
- Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, MA, USA
- Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital, Harvard Medical School, Boston, MA, USA
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | - Xuehong Zhang
- Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital, Harvard Medical School, Boston, MA, USA
- Department of Nutrition, Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | | | - Timothy R Rebbeck
- Zhu Family Center for Global Cancer Prevention, Harvard T.H. Chan School of Public Health and Dana-Farber Cancer Institute, Boston, MA, USA
| | - Christine B Ambrosone
- Department of Cancer Prevention and Control, Roswell Park Comprehensive Cancer Center, Elm & Carlton Streets, Buffalo, NY, 14263, USA
| | - Shuji Ogino
- Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, MA, USA.
- Broad Institute of MIT and Harvard, Boston, MA, USA.
- Department of Pathology, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA.
| |
Collapse
|
10
|
Moll M, Hobbs BD, Menon A, Ghosh AJ, Putman RK, Hino T, Hata A, Silverman EK, Quackenbush J, Castaldi PJ, Hersh CP, McGeachie MJ, Sin DD, Tal-Singer R, Nishino M, Hatabu H, Hunninghake GM, Cho MH. Blood gene expression risk profiles and interstitial lung abnormalities: COPDGene and ECLIPSE cohort studies. Respir Res 2022; 23:157. [PMID: 35715807 PMCID: PMC9204872 DOI: 10.1186/s12931-022-02077-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2022] [Accepted: 06/03/2022] [Indexed: 12/05/2022] Open
Abstract
BACKGROUND Interstitial lung abnormalities (ILA) are radiologic findings that may progress to idiopathic pulmonary fibrosis (IPF). Blood gene expression profiles can predict IPF mortality, but whether these same genes associate with ILA and ILA outcomes is unknown. This study evaluated if a previously described blood gene expression profile associated with IPF mortality is associated with ILA and all-cause mortality. METHODS In COPDGene and ECLIPSE study participants with visual scoring of ILA and gene expression data, we evaluated the association of a previously described IPF mortality score with ILA and mortality. We also trained a new ILA score, derived using genes from the IPF score, in a subset of COPDGene. We tested the association with ILA and mortality on the remainder of COPDGene and ECLIPSE. RESULTS In 1469 COPDGene (training n = 734; testing n = 735) and 571 ECLIPSE participants, the IPF score was not associated with ILA or mortality. However, an ILA score derived from IPF score genes was associated with ILA (meta-analysis of test datasets OR 1.4 [95% CI: 1.2-1.6]) and mortality (HR 1.25 [95% CI: 1.12-1.41]). Six of the 11 genes in the ILA score had discordant directions of effects compared to the IPF score. The ILA score partially mediated the effects of age on mortality (11.8% proportion mediated). CONCLUSIONS An ILA gene expression score, derived from IPF mortality-associated genes, identified genes with concordant and discordant effects on IPF mortality and ILA. These results suggest shared, and unique biologic processes, amongst those with ILA, IPF, aging, and death.
Collapse
Affiliation(s)
- Matthew Moll
- Channing Division for Network Medicine, Brigham and Women's Hospital, Boston, MA, 02115, USA
- Division of Pulmonary and Critical Care Medicine, Brigham and Women's Hospital, Boston, MA, 02115, USA
- Harvard Medical School, Boston, MA, 02115, USA
| | - Brian D Hobbs
- Channing Division for Network Medicine, Brigham and Women's Hospital, Boston, MA, 02115, USA
- Division of Pulmonary and Critical Care Medicine, Brigham and Women's Hospital, Boston, MA, 02115, USA
- Harvard Medical School, Boston, MA, 02115, USA
| | - Aravind Menon
- Division of Pulmonary and Critical Care Medicine, Brigham and Women's Hospital, Boston, MA, 02115, USA
- Harvard Medical School, Boston, MA, 02115, USA
| | - Auyon J Ghosh
- Channing Division for Network Medicine, Brigham and Women's Hospital, Boston, MA, 02115, USA
- Division of Pulmonary and Critical Care Medicine, Brigham and Women's Hospital, Boston, MA, 02115, USA
- Harvard Medical School, Boston, MA, 02115, USA
| | - Rachel K Putman
- Division of Pulmonary and Critical Care Medicine, Brigham and Women's Hospital, Boston, MA, 02115, USA
- Harvard Medical School, Boston, MA, 02115, USA
| | - Takuya Hino
- Harvard Medical School, Boston, MA, 02115, USA
- Department of Radiology, Center for Pulmonary Functional Imaging, Brigham and Women's Hospital, Boston, MA, 02115, USA
| | - Akinori Hata
- Harvard Medical School, Boston, MA, 02115, USA
- Department of Radiology, Center for Pulmonary Functional Imaging, Brigham and Women's Hospital, Boston, MA, 02115, USA
| | - Edwin K Silverman
- Channing Division for Network Medicine, Brigham and Women's Hospital, Boston, MA, 02115, USA
- Division of Pulmonary and Critical Care Medicine, Brigham and Women's Hospital, Boston, MA, 02115, USA
- Harvard Medical School, Boston, MA, 02115, USA
| | - John Quackenbush
- Channing Division for Network Medicine, Brigham and Women's Hospital, Boston, MA, 02115, USA
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, 02115, USA
| | - Peter J Castaldi
- Channing Division for Network Medicine, Brigham and Women's Hospital, Boston, MA, 02115, USA
- Harvard Medical School, Boston, MA, 02115, USA
- Division of General Internal Medicine and Primary Care, Department of Medicine, Brigham and Women's Hospital, Boston, MA, 02115, Canada
| | - Craig P Hersh
- Channing Division for Network Medicine, Brigham and Women's Hospital, Boston, MA, 02115, USA
- Division of Pulmonary and Critical Care Medicine, Brigham and Women's Hospital, Boston, MA, 02115, USA
- Harvard Medical School, Boston, MA, 02115, USA
| | - Michael J McGeachie
- Channing Division for Network Medicine, Brigham and Women's Hospital, Boston, MA, 02115, USA
| | - Don D Sin
- Centre for Heart Lung Innovation, St. Paul's Hospital, and Department of Medicine (Respiratory Division), University of British Columbia, Vancouver, BC, Canada
| | | | - Mizuki Nishino
- Harvard Medical School, Boston, MA, 02115, USA
- Department of Radiology, Center for Pulmonary Functional Imaging, Brigham and Women's Hospital, Boston, MA, 02115, USA
| | - Hiroto Hatabu
- Harvard Medical School, Boston, MA, 02115, USA
- Department of Radiology, Center for Pulmonary Functional Imaging, Brigham and Women's Hospital, Boston, MA, 02115, USA
| | - Gary M Hunninghake
- Division of Pulmonary and Critical Care Medicine, Brigham and Women's Hospital, Boston, MA, 02115, USA
- Harvard Medical School, Boston, MA, 02115, USA
| | - Michael H Cho
- Channing Division for Network Medicine, Brigham and Women's Hospital, Boston, MA, 02115, USA.
- Division of Pulmonary and Critical Care Medicine, Brigham and Women's Hospital, Boston, MA, 02115, USA.
- Harvard Medical School, Boston, MA, 02115, USA.
| |
Collapse
|
11
|
Gaynor SM, Fagny M, Lin X, Platig J, Quackenbush J. Connectivity in eQTL networks dictates reproducibility and genomic properties. Cell Rep Methods 2022; 2:100218. [PMID: 35637906 PMCID: PMC9142682 DOI: 10.1016/j.crmeth.2022.100218] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/22/2021] [Revised: 02/08/2022] [Accepted: 04/25/2022] [Indexed: 01/11/2023]
Abstract
Expression quantitative trait locus (eQTL) analysis associates SNPs with gene expression; these relationships can be represented as a bipartite network with association strength as "edge weights" between SNPs and genes. However, most eQTL networks use binary edge weights based on thresholded FDR estimates: definitions that influence reproducibility and downstream analyses. We constructed twenty-nine tissue-specific eQTL networks using GTEx data and evaluated a comprehensive set of network specifications based on false discovery rates, test statistics, and p values, focusing on the degree centrality-a metric of an SNP or gene node's potential network influence. We found a thresholded Benjamini-Hochberg q value weighted by the Z-statistic balances metric reproducibility and computational efficiency. Our estimated gene degrees positively correlate with gene degrees in gene regulatory networks, demonstrating that these networks are complementary in understanding regulation. Gene degrees also correlate with genetic diversity, and heritability analyses show that highly connected nodes are enriched for tissue-relevant traits.
Collapse
Affiliation(s)
- Sheila M. Gaynor
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA 02115, USA
- Department of Biostatistics and Computational Biology and Center for Cancer Computational Biology, Dana-Farber Cancer Institute, Boston, MA 02115, USA
| | - Maud Fagny
- Department of Biostatistics and Computational Biology and Center for Cancer Computational Biology, Dana-Farber Cancer Institute, Boston, MA 02115, USA
- Université Paris-Saclay, INRAE, CNRS, AgroParisTech, GQE - Le Moulon, 91190 Gif-sur-Yvette, France
| | - Xihong Lin
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA 02115, USA
- Department of Statistics, Harvard University, Cambridge, MA 02138, USA
| | - John Platig
- Department of Biostatistics and Computational Biology and Center for Cancer Computational Biology, Dana-Farber Cancer Institute, Boston, MA 02115, USA
- Channing Division of Network Medicine, Brigham and Women’s Hospital, Boston, MA 02115, USA
- Department of Medicine, Harvard Medical School, Boston, MA 02115, USA
| | - John Quackenbush
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA 02115, USA
- Department of Biostatistics and Computational Biology and Center for Cancer Computational Biology, Dana-Farber Cancer Institute, Boston, MA 02115, USA
- Channing Division of Network Medicine, Brigham and Women’s Hospital, Boston, MA 02115, USA
| |
Collapse
|
12
|
Weighill D, Ben Guebila M, Glass K, Quackenbush J, Platig J. Predicting genotype-specific gene regulatory networks. Genome Res 2022; 32:524-533. [PMID: 35193937 PMCID: PMC8896459 DOI: 10.1101/gr.275107.120] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2021] [Accepted: 01/11/2022] [Indexed: 11/25/2022]
Abstract
Understanding how each person's unique genotype influences their individual patterns of gene regulation has the potential to improve our understanding of human health and development, and to refine genotype-specific disease risk assessments and treatments. However, the effects of genetic variants are not typically considered when constructing gene regulatory networks, despite the fact that many disease-associated genetic variants are thought to have regulatory effects, including the disruption of transcription factor (TF) binding. We developed EGRET (Estimating the Genetic Regulatory Effect on TFs), which infers a genotype-specific gene regulatory network for each individual in a study population. EGRET begins by constructing a genotype-informed TF-gene prior network derived using TF motif predictions, expression quantitative trait locus (eQTL) data, individual genotypes, and the predicted effects of genetic variants on TF binding. It then uses a technique known as message passing to integrate this prior network with gene expression and TF protein–protein interaction data to produce a refined, genotype-specific regulatory network. We used EGRET to infer gene regulatory networks for two blood-derived cell lines and identified genotype-associated, cell line–specific regulatory differences that we subsequently validated using allele-specific expression, chromatin accessibility QTLs, and differential ChIP-seq TF binding. We also inferred EGRET networks for three cell types from each of 119 individuals and identified cell type–specific regulatory differences associated with diseases related to those cell types. EGRET is, to our knowledge, the first method that infers networks reflective of individual genetic variation in a way that provides insight into the genetic regulatory associations driving complex phenotypes.
Collapse
Affiliation(s)
- Deborah Weighill
- Harvard T.H. Chan School of Public Health, Boston, Massachusetts 02115, USA
| | | | - Kimberly Glass
- Harvard T.H. Chan School of Public Health, Boston, Massachusetts 02115, USA.,Channing Division of Network Medicine, Brigham and Women's Hospital, Boston, Massachusetts 02115, USA.,Harvard Medical School, Boston, Massachusetts 02115, USA
| | - John Quackenbush
- Harvard T.H. Chan School of Public Health, Boston, Massachusetts 02115, USA.,Channing Division of Network Medicine, Brigham and Women's Hospital, Boston, Massachusetts 02115, USA
| | - John Platig
- Channing Division of Network Medicine, Brigham and Women's Hospital, Boston, Massachusetts 02115, USA.,Harvard Medical School, Boston, Massachusetts 02115, USA
| |
Collapse
|
13
|
Guebila MB, Morgan DC, Glass K, Kuijjer ML, DeMeo DL, Quackenbush J. gpuZoo: Cost-effective estimation of gene regulatory networks using the Graphics Processing Unit. NAR Genom Bioinform 2022; 4:lqac002. [PMID: 35156023 PMCID: PMC8826808 DOI: 10.1093/nargab/lqac002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2021] [Revised: 12/28/2021] [Accepted: 02/02/2022] [Indexed: 11/14/2022] Open
Abstract
Gene regulatory network inference allows for the modeling of genome-scale regulatory processes that are altered during development, in disease, and in response to perturbations. Our group has developed a collection of tools to model various regulatory processes, including transcriptional (PANDA, SPIDER) and post-transcriptional (PUMA) gene regulation, as well as gene regulation in individual samples (LIONESS). These methods work by postulating a network structure and then optimizing that structure to be consistent with multiple lines of biological evidence through repeated operations on data matrices. Although our methods are widely used, the corresponding computational complexity, and the associated costs and run times, do limit some applications. To improve the cost/time performance of these algorithms, we developed gpuZoo which implements GPU-accelerated calculations, dramatically improving the performance of these algorithms. The runtime of the gpuZoo implementation in MATLAB and Python is up to 61 times faster and 28 times less expensive than multi-core CPU implementation of the same methods. gpuZoo is available in MATLAB through the netZooM package https://github.com/netZoo/netZooM and in Python through the netZooPy package https://github.com/netZoo/netZooPy.
Collapse
Affiliation(s)
- Marouen Ben Guebila
- Department of Biostatistics, Harvard School of Public Health, Boston, MA, USA
| | - Daniel C Morgan
- Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA
| | - Kimberly Glass
- Department of Biostatistics, Harvard School of Public Health, Boston, MA, USA
| | - Marieke L Kuijjer
- Centre for Molecular Medicine Norway (NCMM), Nordic EMBL Partnership, University of Oslo, Oslo, Norway
| | - Dawn L DeMeo
- Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA
| | - John Quackenbush
- Department of Biostatistics, Harvard School of Public Health, Boston, MA, USA
| |
Collapse
|
14
|
Ben Guebila M, Lopes-Ramos CM, Weighill D, Sonawane A, Burkholz R, Shamsaei B, Platig J, Glass K, Kuijjer M, Quackenbush J. GRAND: a database of gene regulatory network models across human conditions. Nucleic Acids Res 2022; 50:D610-D621. [PMID: 34508353 PMCID: PMC8728257 DOI: 10.1093/nar/gkab778] [Citation(s) in RCA: 23] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2021] [Revised: 08/17/2021] [Accepted: 09/08/2021] [Indexed: 12/14/2022] Open
Abstract
Gene regulation plays a fundamental role in shaping tissue identity, function, and response to perturbation. Regulatory processes are controlled by complex networks of interacting elements, including transcription factors, miRNAs and their target genes. The structure of these networks helps to determine phenotypes and can ultimately influence the development of disease or response to therapy. We developed GRAND (https://grand.networkmedicine.org) as a database for computationally-inferred, context-specific gene regulatory network models that can be compared between biological states, or used to predict which drugs produce changes in regulatory network structure. The database includes 12 468 genome-scale networks covering 36 human tissues, 28 cancers, 1378 unperturbed cell lines, as well as 173 013 TF and gene targeting scores for 2858 small molecule-induced cell line perturbation paired with phenotypic information. GRAND allows the networks to be queried using phenotypic information and visualized using a variety of interactive tools. In addition, it includes a web application that matches disease states to potentially therapeutic small molecule drugs using regulatory network properties.
Collapse
Affiliation(s)
- Marouen Ben Guebila
- Department of Biostatistics, Harvard School of Public Health, Boston, MA, USA
| | | | - Deborah Weighill
- Department of Biostatistics, Harvard School of Public Health, Boston, MA, USA
| | - Abhijeet Rajendra Sonawane
- Center for Interdisciplinary Cardiovascular Sciences, Cardiovascular Medicine, Brigham and Women's Hospital, Harvard Medical School, Boston, MA02115, USA
| | - Rebekka Burkholz
- Department of Biostatistics, Harvard School of Public Health, Boston, MA, USA
| | - Behrouz Shamsaei
- Division of Biostatistics and Bioinformatics, Department of Environmental and Public Health Sciences, University of Cincinnati College of Medicine, Cincinnati, OH, USA
| | - John Platig
- Channing Division of Network Medicine, Department of Medicine, Harvard Medical School and Brigham and Women’s Hospital, Boston, MA, USA
| | - Kimberly Glass
- Department of Biostatistics, Harvard School of Public Health, Boston, MA, USA
- Channing Division of Network Medicine, Department of Medicine, Harvard Medical School and Brigham and Women’s Hospital, Boston, MA, USA
| | - Marieke L Kuijjer
- Center for Molecular Medicine Norway, Faculty of Medicine, University of Oslo, Oslo, Norway
- Leiden University Medical Center, Leiden, The Netherlands
| | - John Quackenbush
- Department of Biostatistics, Harvard School of Public Health, Boston, MA, USA
- Channing Division of Network Medicine, Department of Medicine, Harvard Medical School and Brigham and Women’s Hospital, Boston, MA, USA
| |
Collapse
|
15
|
Lopes-Ramos CM, Belova T, Brunner TH, Ben Guebila M, Osorio D, Quackenbush J, Kuijjer ML. Regulatory Network of PD1 Signaling Is Associated with Prognosis in Glioblastoma Multiforme. Cancer Res 2021; 81:5401-5412. [PMID: 34493595 PMCID: PMC8563450 DOI: 10.1158/0008-5472.can-21-0730] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2021] [Revised: 07/20/2021] [Accepted: 09/02/2021] [Indexed: 01/07/2023]
Abstract
Glioblastoma is an aggressive cancer of the brain and spine. While analysis of glioblastoma 'omics data has somewhat improved our understanding of the disease, it has not led to direct improvement in patient survival. Cancer survival is often characterized by differences in gene expression, but the mechanisms that drive these differences are generally unknown. We therefore set out to model the regulatory mechanisms associated with glioblastoma survival. We inferred individual patient gene regulatory networks using data from two different expression platforms from The Cancer Genome Atlas. We performed comparative network analysis between patients with long- and short-term survival. Seven pathways were identified as associated with survival, all of them involved in immune signaling; differential regulation of PD1 signaling was validated to correspond with outcome in an independent dataset from the German Glioma Network. In this pathway, transcriptional repression of genes for which treatment options are available was lost in short-term survivors; this was independent of mutational burden and only weakly associated with T-cell infiltration. Collectively, these results provide a new way to stratify patients with glioblastoma that uses network features as biomarkers to predict survival. They also identify new potential therapeutic interventions, underscoring the value of analyzing gene regulatory networks in individual patients with cancer. SIGNIFICANCE: Genome-wide network modeling of individual glioblastomas identifies dysregulation of PD1 signaling in patients with poor prognosis, indicating this approach can be used to understand how gene regulation influences cancer progression.
Collapse
Affiliation(s)
- Camila M. Lopes-Ramos
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, Massachusetts
| | - Tatiana Belova
- Centre for Molecular Medicine Norway, University of Oslo, Oslo, Norway
| | | | - Marouen Ben Guebila
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, Massachusetts
| | - Daniel Osorio
- Centre for Molecular Medicine Norway, University of Oslo, Oslo, Norway
| | - John Quackenbush
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, Massachusetts.,Department of Cancer Biology, Dana-Farber Cancer Institute, Boston, Massachusetts.,Channing Division of Network Medicine, Harvard Medical School, Boston, Massachusetts
| | - Marieke L. Kuijjer
- Centre for Molecular Medicine Norway, University of Oslo, Oslo, Norway.,Department of Pathology, Leiden University Medical Center, Leiden, the Netherlands.,Corresponding Author: Marieke L. Kuijjer, Centre for Molecular Medicine Norway, University of Oslo, Guastadalléen 21, Oslo 0318, Norway. Phone: 47-22840528; E-mail:
| |
Collapse
|
16
|
Neri de Souza Reis V, Tahira AC, Daguano Gastaldi V, Mari P, Portolese J, Feio dos Santos AC, Lisboa B, Mari J, Caetano SC, Brunoni D, Bordini D, Silvestre de Paula C, Vêncio RZN, Quackenbush J, Brentani H. Environmental Influences Measured by Epigenetic Clock and Vulnerability Components at Birth Impact Clinical ASD Heterogeneity. Genes (Basel) 2021; 12:genes12091433. [PMID: 34573415 PMCID: PMC8467464 DOI: 10.3390/genes12091433] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2021] [Revised: 09/03/2021] [Accepted: 09/08/2021] [Indexed: 11/16/2022] Open
Abstract
Although Autism Spectrum Disorders (ASD) is recognized as being heavily influenced by genetic factors, the role of epigenetic and environmental factors is still being established. This study aimed to identify ASD vulnerability components based on familial history and intrauterine environmental stress exposure, explore possible vulnerability subgroups, access DNA methylation age acceleration (AA) as a proxy of stress exposure during life, and evaluate the association of ASD vulnerability components and AA to phenotypic severity measures. Principal Component Analysis (PCA) was used to search the vulnerability components from 67 mothers of autistic children. We found that PC1 had a higher correlation with psychosocial stress (maternal stress, maternal education, and social class), and PC2 had a higher correlation with biological factors (psychiatric family history and gestational complications). Comparing the methylome between above and below PC1 average subgroups we found 11,879 statistically significant differentially methylated probes (DMPs, p < 0.05). DMPs CpG sites were enriched in variably methylated regions (VMRs), most showing environmental and genetic influences. Hypermethylated probes presented higher rates in different regulatory regions associated with functional SNPs, indicating that the subgroups may have different affected regulatory regions and their liability to disease explained by common variations. Vulnerability components score moderated by epigenetic clock AA was associated with Vineland Total score (p = 0.0036, adjR2 = 0.31), suggesting risk factors with stress burden can influence ASD phenotype.
Collapse
Affiliation(s)
- Viviane Neri de Souza Reis
- Departamento & Instituto de Psiquiatria, Faculdade de Medicina FMUSP, Universidade de São Paulo, São Paulo 05403-903, SP, Brazil; (V.N.d.S.R.); (A.C.T.); (V.D.G.); (P.M.); (J.P.); (A.C.F.d.S.); (B.L.)
| | - Ana Carolina Tahira
- Departamento & Instituto de Psiquiatria, Faculdade de Medicina FMUSP, Universidade de São Paulo, São Paulo 05403-903, SP, Brazil; (V.N.d.S.R.); (A.C.T.); (V.D.G.); (P.M.); (J.P.); (A.C.F.d.S.); (B.L.)
- Instituto Butantan, São Paulo 05503-900, SP, Brazil
| | - Vinícius Daguano Gastaldi
- Departamento & Instituto de Psiquiatria, Faculdade de Medicina FMUSP, Universidade de São Paulo, São Paulo 05403-903, SP, Brazil; (V.N.d.S.R.); (A.C.T.); (V.D.G.); (P.M.); (J.P.); (A.C.F.d.S.); (B.L.)
| | - Paula Mari
- Departamento & Instituto de Psiquiatria, Faculdade de Medicina FMUSP, Universidade de São Paulo, São Paulo 05403-903, SP, Brazil; (V.N.d.S.R.); (A.C.T.); (V.D.G.); (P.M.); (J.P.); (A.C.F.d.S.); (B.L.)
| | - Joana Portolese
- Departamento & Instituto de Psiquiatria, Faculdade de Medicina FMUSP, Universidade de São Paulo, São Paulo 05403-903, SP, Brazil; (V.N.d.S.R.); (A.C.T.); (V.D.G.); (P.M.); (J.P.); (A.C.F.d.S.); (B.L.)
| | - Ana Cecilia Feio dos Santos
- Departamento & Instituto de Psiquiatria, Faculdade de Medicina FMUSP, Universidade de São Paulo, São Paulo 05403-903, SP, Brazil; (V.N.d.S.R.); (A.C.T.); (V.D.G.); (P.M.); (J.P.); (A.C.F.d.S.); (B.L.)
- Laboratório de Pesquisas Básicas em Malária—Entomologia, Seção de Parasitologia—Instituto Evandro Chagas/SVS/MS, Ananindeua 66093-020, PA, Brazil
| | - Bianca Lisboa
- Departamento & Instituto de Psiquiatria, Faculdade de Medicina FMUSP, Universidade de São Paulo, São Paulo 05403-903, SP, Brazil; (V.N.d.S.R.); (A.C.T.); (V.D.G.); (P.M.); (J.P.); (A.C.F.d.S.); (B.L.)
| | - Jair Mari
- Departamento de Psiquiatria, Universidade Federal de São Paulo (UNIFESP), São Paulo 04023-062, SP, Brazil; (J.M.); (S.C.C.); (D.B.); (C.S.d.P.)
| | - Sheila C. Caetano
- Departamento de Psiquiatria, Universidade Federal de São Paulo (UNIFESP), São Paulo 04023-062, SP, Brazil; (J.M.); (S.C.C.); (D.B.); (C.S.d.P.)
| | - Décio Brunoni
- Centro de Ciências Biológicas e da Saúde, Universidade Presbiteriana Mackenzie (UPM), São Paulo 01302-907, SP, Brazil;
| | - Daniela Bordini
- Departamento de Psiquiatria, Universidade Federal de São Paulo (UNIFESP), São Paulo 04023-062, SP, Brazil; (J.M.); (S.C.C.); (D.B.); (C.S.d.P.)
| | - Cristiane Silvestre de Paula
- Departamento de Psiquiatria, Universidade Federal de São Paulo (UNIFESP), São Paulo 04023-062, SP, Brazil; (J.M.); (S.C.C.); (D.B.); (C.S.d.P.)
- Centro de Ciências Biológicas e da Saúde, Universidade Presbiteriana Mackenzie (UPM), São Paulo 01302-907, SP, Brazil;
| | - Ricardo Z. N. Vêncio
- Departamento de Computação e Matemática FFCLRP-USP, Universidade de São Paulo, Ribeirão Preto 14040-901, SP, Brazil;
| | - John Quackenbush
- Center for Cancer Computational Biology, Department of Biostatistics and Computational Biology, Dana-Farber Cancer Institute, Boston, MA 02115, USA; or
- Department of Biostatistics, Harvard T. H. Chan School of Public Health, Boston, MA 02115, USA
| | - Helena Brentani
- Departamento & Instituto de Psiquiatria, Faculdade de Medicina FMUSP, Universidade de São Paulo, São Paulo 05403-903, SP, Brazil; (V.N.d.S.R.); (A.C.T.); (V.D.G.); (P.M.); (J.P.); (A.C.F.d.S.); (B.L.)
- Correspondence: ; Tel.: +55-(11)-99-931-4349
| |
Collapse
|
17
|
Snyder JM, Huang RY, Bai H, Rao VR, Cornes S, Barnholtz-Sloan JS, Gutman D, Fasano R, Van Meir EG, Brat D, Eschbacher J, Quackenbush J, Wen PY, Lee JW. Analysis of morphological characteristics of IDH-mutant/wildtype brain tumors using whole-lesion phenotype analysis. Neurooncol Adv 2021; 3:vdab088. [PMID: 34409295 PMCID: PMC8367280 DOI: 10.1093/noajnl/vdab088] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
Background Although IDH-mutant tumors aggregate to the frontotemporal regions, the clustering pattern of IDH-wildtype tumors is less clear. As voxel-based lesion-symptom mapping (VLSM) has several limitations for solid lesion mapping, a new technique, whole-lesion phenotype analysis (WLPA), is developed. We utilize WLPA to assess spatial clustering of tumors with IDH mutation from The Cancer Genome Atlas and The Cancer Imaging Archive. Methods The degree of tumor clustering segmented from T1 weighted images is measured to every other tumor by a function of lesion similarity to each other via the Hausdorff distance. Each tumor is ranked according to the degree to which its neighboring tumors show identical phenotypes, and through a permutation technique, significant tumors are determined. VLSM was applied through a previously described method. Results A total of 244 patients of mixed-grade gliomas (WHO II-IV) are analyzed, of which 150 were IDH-wildtype and 139 were glioblastomas. VLSM identifies frontal lobe regions that are more likely associated with the presence of IDH mutation but no regions where IDH-wildtype was more likely to be present. WLPA identifies both IDH-mutant and -wildtype tumors exhibit statistically significant spatial clustering. Conclusion WLPA may provide additional statistical power when compared with VLSM without making several potentially erroneous assumptions. WLPA identifies tumors most likely to exhibit particular phenotypes, rather than producing anatomical maps, and may be used in conjunction with VLSM to understand the relationship between tumor morphology and biologically relevant tumor phenotypes.
Collapse
Affiliation(s)
- James M Snyder
- Departments of Neurosurgery and Neurology, Henry Ford Health System, Detroit, Michigan, USA
| | - Raymond Y Huang
- Department of Radiology, Brigham and Women's Hospital, Harvard Medical School, Boston, Massachusetts, USA
| | - Harrison Bai
- Department of Diagnostic Imaging, Rhode Island Hospital and Warren Alpert Medical School of Brown University, Providence, Rhode Island, USA
| | - Vikram R Rao
- Department of Neurology and Weill Institute for Neurosciences, University of California, San Francisco, San Francisco, California, USA
| | - Susannah Cornes
- Department of Neurology and Weill Institute for Neurosciences, University of California, San Francisco, San Francisco, California, USA
| | - Jill S Barnholtz-Sloan
- Department of Population and Quantitative Health Sciences, School of Medicine Case Western Reserve University and University Hospitals of Cleveland, Cleveland, Ohio, USA
| | - David Gutman
- Winship Cancer Institute, Emory University, Atlanta, Georgia, USA
| | - Rebecca Fasano
- Department of Neurology, Emory University, Atlanta, Georgia, USA
| | - Erwin G Van Meir
- O'Neal Comprehensive Cancer Center, University of Alabama at Birmingham (UAB), Birmingham, Alabama, USA
| | - Daniel Brat
- Department of Pathology, Feinberg School of Medicine, Northwestern University, Chicago, Illinois, USA
| | | | - John Quackenbush
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Center for Cancer Computational Biology, Dana Farber Cancer Institute, Boston, Massachusetts, USA
| | - Patrick Y Wen
- Center for Neuro-Oncology, Dana-Farber Cancer Institute, Boston, Massachusetts, USA
| | - Jong Woo Lee
- Department of Neurology, Brigham and Women's Hospital, Harvard Medical School, Boston, Massachusetts, USA
| |
Collapse
|
18
|
Paulson JN, Williams BL, Hehnly C, Mishra N, Sinnar SA, Zhang L, Ssentongo P, Mbabazi-Kabachelor E, Wijetunge DSS, von Bredow B, Mulondo R, Kiwanuka J, Bajunirwe F, Bazira J, Bebell LM, Burgoine K, Couto-Rodriguez M, Ericson JE, Erickson T, Ferrari M, Gladstone M, Guo C, Haran M, Hornig M, Isaacs AM, Kaaya BN, Kangere SM, Kulkarni AV, Kumbakumba E, Li X, Limbrick DD, Magombe J, Morton SU, Mugamba J, Ng J, Olupot-Olupot P, Onen J, Peterson MR, Roy F, Sheldon K, Townsend R, Weeks AD, Whalen AJ, Quackenbush J, Ssenyonga P, Galperin MY, Almeida M, Atkins H, Warf BC, Lipkin WI, Broach JR, Schiff SJ. Paenibacillus infection with frequent viral coinfection contributes to postinfectious hydrocephalus in Ugandan infants. Sci Transl Med 2021; 12:12/563/eaba0565. [PMID: 32998967 DOI: 10.1126/scitranslmed.aba0565] [Citation(s) in RCA: 27] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2019] [Accepted: 05/06/2020] [Indexed: 12/14/2022]
Abstract
Postinfectious hydrocephalus (PIH), which often follows neonatal sepsis, is the most common cause of pediatric hydrocephalus worldwide, yet the microbial pathogens underlying this disease remain to be elucidated. Characterization of the microbial agents causing PIH would enable a shift from surgical palliation of cerebrospinal fluid (CSF) accumulation to prevention of the disease. Here, we examined blood and CSF samples collected from 100 consecutive infant cases of PIH and control cases comprising infants with non-postinfectious hydrocephalus in Uganda. Genomic sequencing of samples was undertaken to test for bacterial, fungal, and parasitic DNA; DNA and RNA sequencing was used to identify viruses; and bacterial culture recovery was used to identify potential causative organisms. We found that infection with the bacterium Paenibacillus, together with frequent cytomegalovirus (CMV) coinfection, was associated with PIH in our infant cohort. Assembly of the genome of a facultative anaerobic bacterial isolate recovered from cultures of CSF samples from PIH cases identified a strain of Paenibacillus thiaminolyticus This strain, designated Mbale, was lethal when injected into mice in contrast to the benign reference Paenibacillus strain. These findings show that an unbiased pan-microbial approach enabled characterization of Paenibacillus in CSF samples from PIH cases, and point toward a pathway of more optimal treatment and prevention for PIH and other proximate neonatal infections.
Collapse
Affiliation(s)
- Joseph N Paulson
- Department of Biostatistics, Product Development, Genentech Inc., South San Francisco, CA 94080, USA
| | - Brent L Williams
- Center for Infection and Immunity, Mailman School of Public Health, Columbia University, New York, NY 10032, USA.,Department of Epidemiology, Columbia University Mailman School of Public Health, New York, NY 10032, USA
| | - Christine Hehnly
- Institute for Personalized Medicine, Department of Biochemistry and Molecular Biology, Pennsylvania State University College of Medicine, Hershey, PA 17033, USA
| | - Nischay Mishra
- Center for Infection and Immunity, Mailman School of Public Health, Columbia University, New York, NY 10032, USA.,Department of Epidemiology, Columbia University Mailman School of Public Health, New York, NY 10032, USA
| | - Shamim A Sinnar
- Center for Neural Engineering, Pennsylvania State University, University Park, PA 16802, USA.,Department of Medicine, Pennsylvania State University College of Medicine, Hershey, PA 17033, USA
| | - Lijun Zhang
- Institute for Personalized Medicine, Department of Biochemistry and Molecular Biology, Pennsylvania State University College of Medicine, Hershey, PA 17033, USA
| | - Paddy Ssentongo
- Center for Neural Engineering, Pennsylvania State University, University Park, PA 16802, USA.,Department of Engineering Science and Mechanics, Pennsylvania State University, University Park, PA 16802, USA.,Department of Public Health Sciences, Pennsylvania State University College of Medicine, Hershey, PA 17033, USA
| | | | - Dona S S Wijetunge
- Department of Pathology, Pennsylvania State University College of Medicine, Hershey, PA 17033, USA
| | - Benjamin von Bredow
- Department of Pathology, Pennsylvania State University College of Medicine, Hershey, PA 17033, USA
| | - Ronnie Mulondo
- CURE Children's Hospital of Uganda, Plot 97-105, Bugwere Road, P.O. Box 903 Mbale, Uganda
| | - Julius Kiwanuka
- Department of Pediatrics, Mbarara University of Science and Technology, P.O. Box 1410 Mbarara, Uganda
| | - Francis Bajunirwe
- Department of Epidemiology, Mbarara University of Science and Technology, P.O. Box 1410, Mbarara, Uganda
| | - Joel Bazira
- Department of Microbiology, Mbarara University of Science and Technology, P.O. Box 1410 Mbarara, Uganda
| | - Lisa M Bebell
- Division of Infectious Disease, Massachusetts Genereal Hospital, Harvard Medical School, 55 Fruit St, GRJ-504, Boston, MA 02114, USA
| | - Kathy Burgoine
- Neonatal Unit, Department of Paediatrics and Child Health, Mbale Regional Referral Hospital, Plot 29-33 Pallisa Road, P.O. Box 1966, Mbale, Uganda.,Mbale Clinical Research Institute, Mbale Regional Referral Hospital, Plot 29-33 Pallisa Road, P.O. Box 1966 Mbale, Uganda.,University of Liverpool, Liverpool, L69 3BX, UK
| | - Mara Couto-Rodriguez
- Center for Infection and Immunity, Mailman School of Public Health, Columbia University, New York, NY 10032, USA.,Biotia, 100 6th avenue, New York, NY 10013, USA
| | - Jessica E Ericson
- Division of Pediatric Infectious Disease, Pennsylvania State University College of Medicine, Hershey, PA 17033, USA
| | - Tim Erickson
- CURE Children's Hospital of Uganda, Plot 97-105, Bugwere Road, P.O. Box 903 Mbale, Uganda
| | - Matthew Ferrari
- Center for Infectious Disease Dynamics, Pennsylvania State University, University Park, PA 16802, USA.,Department of Biology, Pennsylvania State University, University Park, PA 16802, USA.,Department of Statistics, Pennsylvania State University, University Park, PA 16802, USA
| | - Melissa Gladstone
- Institute for Translational Medicine, University of Liverpool, Liverpool, L12 2AP, UK
| | - Cheng Guo
- Center for Infection and Immunity, Mailman School of Public Health, Columbia University, New York, NY 10032, USA
| | - Murali Haran
- Department of Statistics, Pennsylvania State University, University Park, PA 16802, USA
| | - Mady Hornig
- Department of Epidemiology, Columbia University Mailman School of Public Health, New York, NY 10032, USA
| | - Albert M Isaacs
- Department of Neuroscience, Washington University School of Medicine, St. Louis, MO 63130, USA
| | - Brian Nsubuga Kaaya
- CURE Children's Hospital of Uganda, Plot 97-105, Bugwere Road, P.O. Box 903 Mbale, Uganda
| | - Sheila M Kangere
- CURE Children's Hospital of Uganda, Plot 97-105, Bugwere Road, P.O. Box 903 Mbale, Uganda
| | - Abhaya V Kulkarni
- Division of Neurosurgery, Hospital for Sick Children, University of Toronto, Toronto, Ontario, M5G 1X8, Canada
| | - Elias Kumbakumba
- Department of Pediatrics, Mbarara University of Science and Technology, P.O. Box 1410 Mbarara, Uganda
| | - Xiaoxiao Li
- Institute for Translational Medicine, University of Liverpool, Liverpool, L12 2AP, UK
| | - David D Limbrick
- Department of Neurological Surgery, Washington University School of Medicine, St. Louis, MO 63130, USA
| | - Joshua Magombe
- CURE Children's Hospital of Uganda, Plot 97-105, Bugwere Road, P.O. Box 903 Mbale, Uganda
| | - Sarah U Morton
- Division of Newborn Medicine, Boston Children's Hospital and Department of Pediatrics, Harvard Medical School, Boston MA 02115, USA
| | - John Mugamba
- CURE Children's Hospital of Uganda, Plot 97-105, Bugwere Road, P.O. Box 903 Mbale, Uganda
| | - James Ng
- Center for Infection and Immunity, Mailman School of Public Health, Columbia University, New York, NY 10032, USA
| | - Peter Olupot-Olupot
- Mbale Clinical Research Institute, Mbale Regional Referral Hospital, Plot 29-33 Pallisa Road, P.O. Box 1966 Mbale, Uganda.,Busitema University, Mbale Campus, Plot 29-33 Pallisa Road, P.O. Box 1966, Mbale, Uganda
| | - Justin Onen
- CURE Children's Hospital of Uganda, Plot 97-105, Bugwere Road, P.O. Box 903 Mbale, Uganda
| | - Mallory R Peterson
- Center for Neural Engineering, Pennsylvania State University, University Park, PA 16802, USA.,Department of Engineering Science and Mechanics, Pennsylvania State University, University Park, PA 16802, USA
| | - Farrah Roy
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA 02115, USA
| | - Kathryn Sheldon
- Institute for Personalized Medicine, Department of Biochemistry and Molecular Biology, Pennsylvania State University College of Medicine, Hershey, PA 17033, USA
| | - Reid Townsend
- Department of Medicine, Washington University School of Medicine , St. Louis, MO 63130, USA
| | - Andrew D Weeks
- Sanyu Research Unit, Liverpool Women's Hospital, University of Liverpool, Liverpool L8 7SS, UK
| | - Andrew J Whalen
- Department of Mechanical Engineering, Pennsylvania State University, University Park, PA 16802, USA
| | - John Quackenbush
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA 02115, USA
| | - Peter Ssenyonga
- CURE Children's Hospital of Uganda, Plot 97-105, Bugwere Road, P.O. Box 903 Mbale, Uganda
| | - Michael Y Galperin
- National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, MD 20894, USA
| | - Mathieu Almeida
- Université Paris-Saclay, INRAE, MGP, Jouy-en-Josas, 78350, France
| | - Hannah Atkins
- Department of Comparative Medicine, Pennsylvania State University College of Medicine, Hershey, PA 17033, USA
| | - Benjamin C Warf
- Department of Neurosurgery, Boston Children's Hospital, Harvard Medical School, Boston, MA 02115, USA
| | - W Ian Lipkin
- Center for Infection and Immunity, Mailman School of Public Health, Columbia University, New York, NY 10032, USA.,Department of Epidemiology, Columbia University Mailman School of Public Health, New York, NY 10032, USA
| | - James R Broach
- Institute for Personalized Medicine, Department of Biochemistry and Molecular Biology, Pennsylvania State University College of Medicine, Hershey, PA 17033, USA
| | - Steven J Schiff
- Center for Neural Engineering, Pennsylvania State University, University Park, PA 16802, USA. .,Department of Engineering Science and Mechanics, Pennsylvania State University, University Park, PA 16802, USA.,Center for Infectious Disease Dynamics, Pennsylvania State University, University Park, PA 16802, USA.,Department of Neurosurgery, Pennsylvania State University College of Medicine, Hershey, PA 17033, USA.,Department of Physics, Pennsylvania State University, University Park, PA 16802, USA
| |
Collapse
|
19
|
Burkholz R, Quackenbush J, Bojar D. Using graph convolutional neural networks to learn a representation for glycans. Cell Rep 2021; 35:109251. [PMID: 34133929 PMCID: PMC9208909 DOI: 10.1016/j.celrep.2021.109251] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2021] [Revised: 05/05/2021] [Accepted: 05/24/2021] [Indexed: 02/06/2023] Open
Abstract
As the only nonlinear and the most diverse biological sequence, glycans offer substantial challenges for computational biology. These carbohydrates participate in nearly all biological processes—from protein folding to viral cell entry—yet are still not well understood. There are few computational methods to link glycan sequences to functions, and they do not fully leverage all available information about glycans. SweetNet is a graph convolutional neural network that uses graph representation learning to facilitate a computational understanding of glycobiology. SweetNet explicitly incorporates the nonlinear nature of glycans and establishes a framework to map any glycan sequence to a representation. We show that SweetNet outperforms other computational methods in predicting glycan properties on all reported tasks. More importantly, we show that glycan representations, learned by SweetNet, are predictive of organismal phenotypic and environmental properties. Finally, we use glycan-focused machine learning to predict viral glycan binding, which can be used to discover viral receptors. Burkholz et al. develop an analysis platform for glycans, using graph convolutional neural networks, that considers the branched nature of these carbohydrates. They demonstrate that glycan-focused machine learning can be employed for various purposes, such as to cluster species according to their glycomic similarity or to identify viral receptors.
Collapse
Affiliation(s)
- Rebekka Burkholz
- Department of Biostatistics, Harvard School of Public Health, Boston, MA, USA
| | - John Quackenbush
- Department of Biostatistics, Harvard School of Public Health, Boston, MA, USA; Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA; Department of Biostatistics and Computational Biology, Dana-Farber Cancer Institute, Boston, MA, USA
| | - Daniel Bojar
- Department of Chemistry and Molecular Biology, University of Gothenburg, Gothenburg, Sweden; Wallenberg Centre for Molecular and Translational Medicine, University of Gothenburg, Gothenburg, Sweden.
| |
Collapse
|
20
|
Lopes-Ramos CM, Chen CY, Kuijjer ML, Paulson JN, Sonawane AR, Fagny M, Platig J, Glass K, Quackenbush J, DeMeo DL. Sex Differences in Gene Expression and Regulatory Networks across 29 Human Tissues. Cell Rep 2021; 31:107795. [PMID: 32579922 DOI: 10.1016/j.celrep.2020.107795] [Citation(s) in RCA: 148] [Impact Index Per Article: 49.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2016] [Revised: 04/01/2020] [Accepted: 05/29/2020] [Indexed: 11/25/2022] Open
Abstract
Sex differences manifest in many diseases and may drive sex-specific therapeutic responses. To understand the molecular basis of sex differences, we evaluated sex-biased gene regulation by constructing sample-specific gene regulatory networks in 29 human healthy tissues using 8,279 whole-genome expression profiles from the Genotype-Tissue Expression (GTEx) project. We find sex-biased regulatory network structures in each tissue. Even though most transcription factors (TFs) are not differentially expressed between males and females, many have sex-biased regulatory targeting patterns. In each tissue, genes that are differentially targeted by TFs between the sexes are enriched for tissue-related functions and diseases. In brain tissue, for example, genes associated with Parkinson's disease and Alzheimer's disease are targeted by different sets of TFs in each sex. Our systems-based analysis identifies a repertoire of TFs that play important roles in sex-specific architecture of gene regulatory networks, and it underlines sex-specific regulatory processes in both health and disease.
Collapse
Affiliation(s)
| | - Cho-Yi Chen
- Institute of Biomedical Informatics, National Yang-Ming University, Taipei, Taiwan
| | - Marieke L Kuijjer
- Centre for Molecular Medicine Norway (NCMM), Nordic EMBL Partnership, University of Oslo, Oslo, Norway
| | - Joseph N Paulson
- Department of Biostatistics, Product Development, Genentech Inc., San Francisco, CA, USA
| | - Abhijeet R Sonawane
- Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA
| | - Maud Fagny
- Genetique Quantitative et Evolution-Le Moulon, Universite Paris-Saclay, Institut National de Recherche pour l'Agriculture, l'Alimentation et l'Environnement, Centre National de la Recherche Scientifique, AgroParisTech, Gif-sur-Yvette, France
| | - John Platig
- Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA
| | - Kimberly Glass
- Department of Biostatistics, Harvard School of Public Health, Boston, MA, USA; Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA
| | - John Quackenbush
- Department of Biostatistics, Harvard School of Public Health, Boston, MA, USA; Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA; Department of Biostatistics and Computational Biology, Dana-Farber Cancer Institute, Boston, MA, USA.
| | - Dawn L DeMeo
- Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA; Division of Pulmonary and Critical Care Medicine, Brigham and Women's Hospital, Boston, MA, USA.
| |
Collapse
|
21
|
Abstract
Profiling of whole transcriptomes has become a cornerstone of molecular biology and an invaluable tool for the characterization of clinical phenotypes and the identification of disease subtypes. Analyses of these data are becoming ever more sophisticated as we move beyond simple comparisons to consider networks of higher-order interactions and associations. Gene regulatory networks (GRNs) model the regulatory relationships of transcription factors and genes and have allowed the identification of differentially regulated processes in disease systems. In this perspective, we discuss gene targeting scores, which measure changes in inferred regulatory network interactions, and their use in identifying disease-relevant processes. In addition, we present an example analysis for pancreatic ductal adenocarcinoma (PDAC), demonstrating the power of gene targeting scores to identify differential processes between complex phenotypes, processes that would have been missed by only performing differential expression analysis. This example demonstrates that gene targeting scores are an invaluable addition to gene expression analysis in the characterization of diseases and other complex phenotypes.
Collapse
Affiliation(s)
- Deborah Weighill
- Department of Biostatistics, Harvard T. H. Chan School of Public Health, Harvard University, Boston, MA, United States
| | - Marouen Ben Guebila
- Department of Biostatistics, Harvard T. H. Chan School of Public Health, Harvard University, Boston, MA, United States
| | - Kimberly Glass
- Department of Biostatistics, Harvard T. H. Chan School of Public Health, Harvard University, Boston, MA, United States
- Channing Division of Network Medicine, Brigham and Women’s Hospital, Boston, MA, United States
- Harvard Medical School, Harvard University, Boston, MA, United States
| | - John Platig
- Channing Division of Network Medicine, Brigham and Women’s Hospital, Boston, MA, United States
- Harvard Medical School, Harvard University, Boston, MA, United States
| | - Jen Jen Yeh
- Departments of Surgery and Pharmacology, Lineberger Comprehensive Cancer Center, University of North Carolina at Chapel Hill, Chapel Hill, NC, United States
| | - John Quackenbush
- Department of Biostatistics, Harvard T. H. Chan School of Public Health, Harvard University, Boston, MA, United States
- Channing Division of Network Medicine, Brigham and Women’s Hospital, Boston, MA, United States
| |
Collapse
|
22
|
Kuijjer ML, Fagny M, Marin A, Quackenbush J, Glass K. PUMA: PANDA Using MicroRNA Associations. Bioinformatics 2021; 36:4765-4773. [PMID: 32860050 PMCID: PMC7750953 DOI: 10.1093/bioinformatics/btaa571] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2020] [Revised: 05/19/2020] [Accepted: 06/10/2020] [Indexed: 12/27/2022] Open
Abstract
Motivation Conventional methods to analyze genomic data do not make use of the interplay between multiple factors, such as between microRNAs (miRNAs) and the messenger RNA (mRNA) transcripts they regulate, and thereby often fail to identify the cellular processes that are unique to specific tissues. We developed PUMA (PANDA Using MicroRNA Associations), a computational tool that uses message passing to integrate a prior network of miRNA target predictions with target gene co-expression information to model genome-wide gene regulation by miRNAs. We applied PUMA to 38 tissues from the Genotype-Tissue Expression project, integrating RNA-Seq data with two different miRNA target predictions priors, built on predictions from TargetScan and miRanda, respectively. We found that while target predictions obtained from these two different resources are considerably different, PUMA captures similar tissue-specific miRNA–target regulatory interactions in the different network models. Furthermore, the tissue-specific functions of miRNAs we identified based on regulatory profiles (available at: https://kuijjer.shinyapps.io/puma_gtex/) are highly similar between networks modeled on the two target prediction resources. This indicates that PUMA consistently captures important tissue-specific miRNA regulatory processes. In addition, using PUMA we identified miRNAs regulating important tissue-specific processes that, when mutated, may result in disease development in the same tissue. Availability and implementation PUMA is available in C++, MATLAB and Python on GitHub (https://github.com/kuijjerlab and https://netzoo.github.io/). Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Marieke L Kuijjer
- Centre for Molecular Medicine Norway, University of Oslo, Oslo 0318, Norway
| | - Maud Fagny
- UMR7206 Eco-Anthropologie, Muséum National d'Histoire Naturelle, Centre National de la Recherche Scientifique, Université de Paris, Paris 75016, France
| | - Alessandro Marin
- Centre for Computing in Science Education, Department of Physics, University of Oslo, Oslo 0316, Norway
| | - John Quackenbush
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA 02115, USA.,Department of Data Science, Dana-Farber Cancer Institute, Boston, MA 02215, USA.,Channing Division of Network Medicine, Harvard Medical School, Boston, MA 02115, USA
| | - Kimberly Glass
- Channing Division of Network Medicine, Harvard Medical School, Boston, MA 02115, USA
| |
Collapse
|
23
|
Weighill D, Guebila MB, Lopes-Ramos C, Glass K, Quackenbush J, Platig J, Burkholz R. Gene regulatory network inference as relaxed graph matching. Proc AAAI Conf Artif Intell 2021; 35:10263-10272. [PMID: 34707916 PMCID: PMC8546743] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Subscribe] [Scholar Register] [Indexed: 05/17/2023]
Abstract
Bipartite network inference is a ubiquitous problem across disciplines. One important example in the field molecular biology is gene regulatory network inference. Gene regulatory networks are an instrumental tool aiding in the discovery of the molecular mechanisms driving diverse diseases, including cancer. However, only noisy observations of the projections of these regulatory networks are typically assayed. In an effort to better estimate regulatory networks from their noisy projections, we formulate a non-convex but analytically tractable optimization problem called OTTER. This problem can be interpreted as relaxed graph matching between the two projections of the bipartite network. OTTER's solutions can be derived explicitly and inspire a spectral algorithm, for which we provide network recovery guarantees. We also provide an alternative approach based on gradient descent that is more robust to noise compared to the spectral algorithm. Interestingly, this gradient descent approach resembles the message passing equations of an established gene regulatory network inference method, PANDA. Using three cancer-related data sets, we show that OTTER outperforms state-of-the-art inference methods in predicting transcription factor binding to gene regulatory regions. To encourage new graph matching applications to this problem, we have made all networks and validation data publicly available.
Collapse
Affiliation(s)
- Deborah Weighill
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA 02115
| | - Marouen Ben Guebila
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA 02115
| | - Camila Lopes-Ramos
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA 02115
| | - Kimberly Glass
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA 02115
- Channing Division of Network Medicine, Brigham and Women's Hospital
- Harvard Medical School, Boston, MA 02115
| | - John Quackenbush
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA 02115
- Channing Division of Network Medicine, Brigham and Women's Hospital
- Harvard Medical School, Boston, MA 02115
| | - John Platig
- Channing Division of Network Medicine, Brigham and Women's Hospital
- Harvard Medical School, Boston, MA 02115
| | - Rebekka Burkholz
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA 02115
| |
Collapse
|
24
|
Young AT, Carette X, Helmel M, Steen H, Husson RN, Quackenbush J, Platig J. Multi-omic regulatory networks capture downstream effects of kinase inhibition in Mycobacterium tuberculosis. NPJ Syst Biol Appl 2021; 7:8. [PMID: 33514755 PMCID: PMC7846781 DOI: 10.1038/s41540-020-00164-4] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2020] [Accepted: 12/07/2020] [Indexed: 11/30/2022] Open
Abstract
The ability of Mycobacterium tuberculosis (Mtb) to adapt to diverse stresses in its host environment is crucial for pathogenesis. Two essential Mtb serine/threonine protein kinases, PknA and PknB, regulate cell growth in response to environmental stimuli, but little is known about their downstream effects. By combining RNA-Seq data, following treatment with either an inhibitor of both PknA and PknB or an inactive control, with publicly available ChIP-Seq and protein–protein interaction data for transcription factors, we show that the Mtb transcription factor (TF) regulatory network propagates the effects of kinase inhibition and leads to widespread changes in regulatory programs involved in cell wall integrity, stress response, and energy production, among others. We also observe that changes in TF regulatory activity correlate with kinase-specific phosphorylation of those TFs. In addition to characterizing the downstream regulatory effects of PknA/PknB inhibition, this demonstrates the need for regulatory network approaches that can incorporate signal-driven transcription factor modifications.
Collapse
Affiliation(s)
- Albert T Young
- School of Medicine, University of California, San Francisco, USA
| | - Xavier Carette
- Division of Infectious Diseases, Boston Children's Hospital, Boston, USA.,Harvard Medical School, Boston, USA
| | - Michaela Helmel
- Harvard Medical School, Boston, USA.,Department of Pathology, Boston Children's Hospital, Boston, USA
| | - Hanno Steen
- Division of Infectious Diseases, Boston Children's Hospital, Boston, USA.,Harvard Medical School, Boston, USA.,Department of Pathology, Boston Children's Hospital, Boston, USA
| | - Robert N Husson
- Division of Infectious Diseases, Boston Children's Hospital, Boston, USA.,Harvard Medical School, Boston, USA
| | - John Quackenbush
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, USA.,Channing Division of Network Medicine, Brigham and Women's Hospital, Boston, USA
| | - John Platig
- Harvard Medical School, Boston, USA. .,Channing Division of Network Medicine, Brigham and Women's Hospital, Boston, USA.
| |
Collapse
|
25
|
Abstract
Despite their known importance in clinical medicine, differences based on sex and gender are among the least studied factors affecting cancer susceptibility, progression, survival, and therapeutic response. In particular, the molecular mechanisms driving sex differences are poorly understood and so most approaches to precision medicine use mutational or other genomic data to assign therapy without considering how the sex of the individual might influence therapeutic efficacy. The mandate by the National Institutes of Health that research studies include sex as a biological variable has begun to expand our understanding on its importance. Sex differences in cancer may arise due to a combination of environmental, genetic, and epigenetic factors, as well as differences in gene regulation, and expression. Extensive sex differences occur genome-wide, and ultimately influence cancer biology and outcomes. In this review, we summarize the current state of knowledge about sex-specific genetic and genome-wide influences in cancer, describe how differences in response to environmental exposures and genetic and epigenetic alterations alter the trajectory of the disease, and provide insights into the importance of integrative analyses in understanding the interplay of sex and genomics in cancer. In particular, we will explore some of the emerging analytical approaches, such as the use of network methods, that are providing a deeper understanding of the drivers of differences based on sex and gender. Better understanding these complex factors and their interactions will improve cancer prevention, treatment, and outcomes for all individuals.
Collapse
Affiliation(s)
- Camila M Lopes-Ramos
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, United States
| | - John Quackenbush
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, United States.,Department of Data Science, Dana-Farber Cancer Institute, Boston, MA, United States.,Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, United States
| | - Dawn L DeMeo
- Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, United States.,Division of Pulmonary and Critical Care Medicine, Brigham and Women's Hospital, Boston, MA, United States
| |
Collapse
|
26
|
Lam HC, Cloonan SM, Bhashyam AR, Haspel JA, Singh A, Sathirapongsasuti JF, Cervo M, Yao H, Chung AL, Mizumura K, An CH, Shan B, Franks JM, Haley KJ, Owen CA, Tesfaigzi Y, Washko GR, Quackenbush J, Silverman EK, Rahman I, Kim HP, Mahmood A, Biswal SS, Ryter SW, Choi AM. Histone deacetylase 6-mediated selective autophagy regulates COPD-associated cilia dysfunction. J Clin Invest 2020; 130:6189. [PMID: 33136096 DOI: 10.1172/jci143863] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022] Open
|
27
|
Haibe-Kains B, Adam GA, Hosny A, Khodakarami F, Waldron L, Wang B, McIntosh C, Goldenberg A, Kundaje A, Greene CS, Broderick T, Hoffman MM, Leek JT, Korthauer K, Huber W, Brazma A, Pineau J, Tibshirani R, Hastie T, Ioannidis JPA, Quackenbush J, Aerts HJWL. Transparency and reproducibility in artificial intelligence. Nature 2020; 586:E14-E16. [PMID: 33057217 PMCID: PMC8144864 DOI: 10.1038/s41586-020-2766-y] [Citation(s) in RCA: 140] [Impact Index Per Article: 35.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2020] [Accepted: 08/10/2020] [Indexed: 01/15/2023]
Abstract
Breakthroughs in artificial intelligence (AI) hold enormous potential as it can automate complex tasks and go even beyond human performance. In their study, McKinney et al. showed the high potential of AI for breast cancer screening. However, the lack of methods’ details and algorithm code undermines its scientific value. Here, we identify obstacles hindering transparent and reproducible AI research as faced by McKinney et al., and provide solutions to these obstacles with implications for the broader field.
Collapse
Affiliation(s)
- Benjamin Haibe-Kains
- Princess Margaret Cancer Centre, University Health Network, Toronto, Ontario, Canada.
- Department of Medical Biophysics, University of Toronto, Toronto, Ontario, Canada.
- Department of Computer Science, University of Toronto, Toronto, Ontario, Canada.
- Ontario Institute for Cancer Research, Toronto, Ontario, Canada.
- Vector Institute for Artificial Intelligence, Toronto, Ontario, Canada.
| | - George Alexandru Adam
- Department of Computer Science, University of Toronto, Toronto, Ontario, Canada
- Vector Institute for Artificial Intelligence, Toronto, Ontario, Canada
| | - Ahmed Hosny
- Artificial Intelligence in Medicine (AIM) Program, Brigham and Women's Hospital, Harvard Medical School, Boston, MA, USA
- Radiation Oncology and Radiology, Dana-Farber Cancer Institute, Brigham and Women's Hospital, Harvard Medical School, Boston, MA, USA
| | - Farnoosh Khodakarami
- Princess Margaret Cancer Centre, University Health Network, Toronto, Ontario, Canada
- Department of Medical Biophysics, University of Toronto, Toronto, Ontario, Canada
| | - Levi Waldron
- Department of Epidemiology and Biostatistics and Institute for Implementation Science in Population Health, CUNY Graduate School of Public Health and Health Policy, New York, NY, USA
| | - Bo Wang
- Department of Medical Biophysics, University of Toronto, Toronto, Ontario, Canada
- Department of Computer Science, University of Toronto, Toronto, Ontario, Canada
- Vector Institute for Artificial Intelligence, Toronto, Ontario, Canada
- Peter Munk Cardiac Centre, University Health Network, Toronto, Ontario, Canada
- Department of Laboratory Medicine and Pathobiology, University of Toronto, Ontario, Canada
| | - Chris McIntosh
- Department of Medical Biophysics, University of Toronto, Toronto, Ontario, Canada
- Vector Institute for Artificial Intelligence, Toronto, Ontario, Canada
- Peter Munk Cardiac Centre, University Health Network, Toronto, Ontario, Canada
| | - Anna Goldenberg
- Department of Computer Science, University of Toronto, Toronto, Ontario, Canada
- Vector Institute for Artificial Intelligence, Toronto, Ontario, Canada
- SickKids Research Institute, Toronto, Ontario, Canada
- Child and Brain Development Program, CIFAR, Toronto, Ontario, Canada
| | - Anshul Kundaje
- Department of Genetics, Stanford University School of Medicine, Stanford, CA, USA
- Department of Computer Science, Stanford University, Stanford, CA, USA
| | - Casey S Greene
- Dept. of Systems Pharmacology and Translational Therapeutics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
- Childhood Cancer Data Lab, Alex's Lemonade Stand Foundation, Philadelphia, PA, USA
| | - Tamara Broderick
- Department of Electrical Engineering and Computer Science, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Michael M Hoffman
- Princess Margaret Cancer Centre, University Health Network, Toronto, Ontario, Canada
- Department of Medical Biophysics, University of Toronto, Toronto, Ontario, Canada
- Department of Computer Science, University of Toronto, Toronto, Ontario, Canada
- Vector Institute for Artificial Intelligence, Toronto, Ontario, Canada
| | - Jeffrey T Leek
- Department of Biostatistics, Johns Hopkins Bloomberg School of Public Health, Baltimore, MD, USA
| | - Keegan Korthauer
- Department of Statistics, University of British Columbia, Vancouver, British Columbia, Canada
- BC Children's Hospital Research Institute, Vancouver, British Columbia, Canada
| | - Wolfgang Huber
- European Molecular Biology Laboratory, Genome Biology Unit, Heidelberg, Germany
| | - Alvis Brazma
- European Molecular Biology Laboratory, European Bioinformatics Institute, EMBL-EBI, Hinxton, UK
| | - Joelle Pineau
- McGill University, Montreal, Quebec, Canada
- Montreal Institute for Learning Algorithms, Quebec, Canada
| | - Robert Tibshirani
- Department of Statistics, Stanford University School of Humanities and Sciences, Stanford, CA, USA
- Department of Biomedical Data Science, Stanford University School of Medicine, Stanford, CA, USA
| | - Trevor Hastie
- Department of Statistics, Stanford University School of Humanities and Sciences, Stanford, CA, USA
- Department of Biomedical Data Science, Stanford University School of Medicine, Stanford, CA, USA
| | - John P A Ioannidis
- Department of Statistics, Stanford University School of Humanities and Sciences, Stanford, CA, USA
- Department of Biomedical Data Science, Stanford University School of Medicine, Stanford, CA, USA
- Department of Medicine, Stanford University School of Medicine, Stanford, CA, USA
- Meta-Research Innovation Center at Stanford (METRICS), Stanford, CA, USA
- Department of Epidemiology and Population Health, Stanford University School of Medicine, Stanford, CA, USA
| | - John Quackenbush
- Department of Biostatistics, Harvard T.H Chan School of Public Health, Boston, MA, USA
- Channing Division of Network Medicine, Brigham and Women's Hospital, Boston, MA, USA
- Department of Data Science, Dana-Farber Cancer Institute, Boston, MA, USA
| | - Hugo J W L Aerts
- Artificial Intelligence in Medicine (AIM) Program, Brigham and Women's Hospital, Harvard Medical School, Boston, MA, USA
- Radiation Oncology and Radiology, Dana-Farber Cancer Institute, Brigham and Women's Hospital, Harvard Medical School, Boston, MA, USA
- Radiology and Nuclear Medicine, Maastricht University, Maastricht, The Netherlands
- Cardiovascular Imaging Research Center, Massachusetts General Hospital, Harvard Medical School, Boston, MA, USA
| |
Collapse
|
28
|
Lopes-Ramos CM, Kuijjer M, Glass K, DeMeo D, Quackenbush J. Abstract 6569: Regulatory networks of liver carcinoma reveal sex specific patterns of gene regulation. Cancer Res 2020. [DOI: 10.1158/1538-7445.am2020-6569] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Abstract
Despite pronounced sex differences in cancer incidence, severity, and response to treatment, most current approaches to clinical management, as well as therapeutics development and selection, are sex-independent. For most cancer types, including liver cancer, males have a higher risk of developing the disease and a lower survival rate than women. However, the molecular features that drive these sex differences are poorly understood. We inferred patient-specific regulatory networks of liver hepatocellular carcinoma using data from TCGA. By comparing the female and male networks, we found marked sex differences in transcriptional regulatory processes relevant to disease development, progression, and response to therapy. We found that oncogenes have significantly higher regulatory targeting in males, while tumor suppressor genes have significantly higher targeting in females. Many “hallmark” cancer pathways, including the HEDGEHOG, WNT, and TGF-β signaling pathways, were significantly more highly targeted in males, while drug metabolism and immune related pathways were enriched for genes highly targeted in females. We also evaluated sex-biased somatic mutation patterns using mutation and copy number alteration data in TCGA. By summarizing mutations found in genes into pathway mutation scores, we found sex-biased mutation profiles for many pathways, providing additional support for biological sex differences associated with WNT, NOTCH, TGF-β, and HEDGEHOG signaling pathways. Our analysis uncovered patterns of gene regulation that differentiate male and female liver cancer and may be associated with sex differences in prognosis and treatment response. These findings provide insight into the mechanisms that drive clinically observed sex differences and underscore the importance of considering sex as a factor influencing disease etiology and in developing and prescribing therapies. Our network approach can provide insights into why some therapies have differential effect in males and females, and suggest new ways to optimize drug response in each sex.
Citation Format: Camila M. Lopes-Ramos, Marieke Kuijjer, Kimberly Glass, Dawn DeMeo, John Quackenbush. Regulatory networks of liver carcinoma reveal sex specific patterns of gene regulation [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 6569.
Collapse
|
29
|
Morrow JD, Make B, Regan E, Han M, Hersh CP, Tal-Singer R, Quackenbush J, Choi AMK, Silverman EK, DeMeo DL. DNA Methylation Is Predictive of Mortality in Current and Former Smokers. Am J Respir Crit Care Med 2020; 201:1099-1109. [PMID: 31995399 DOI: 10.1164/rccm.201902-0439oc] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
Rationale: Smoking results in at least a decade lower life expectancy. Mortality among current smokers is two to three times as high as never smokers. DNA methylation is an epigenetic modification of the human genome that has been associated with both cigarette smoking and mortality.Objectives: We sought to identify DNA methylation marks in blood that are predictive of mortality in a subset of the COPDGene (Genetic Epidemiology of COPD) study, representing 101 deaths among 667 current and former smokers.Methods: We assayed genome-wide DNA methylation in non-Hispanic white smokers with and without chronic obstructive pulmonary disease (COPD) using blood samples from the COPDGene enrollment visit. We tested whether DNA methylation was associated with mortality in models adjusted for COPD status, age, sex, current smoking status, and pack-years of cigarette smoking. Replication was performed in a subset of 231 individuals from the ECLIPSE (Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints) study.Measurements and Main Results: We identified seven CpG sites associated with mortality (false discovery rate < 20%) that replicated in the ECLIPSE cohort (P < 0.05). None of these marks were associated with longitudinal lung function decline in survivors, smoking history, or current smoking status. However, differential methylation of two replicated PIK3CD (phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit delta) sites were associated with lung function at enrollment (P < 0.05). We also observed associations between DNA methylation and gene expression for the PIK3CD sites.Conclusions: This study is the first to identify variable DNA methylation associated with all-cause mortality in smokers with and without COPD. Evaluating predictive epigenomic marks of smokers in peripheral blood may allow for targeted risk stratification and aid in delivery of future tailored therapeutic interventions.
Collapse
Affiliation(s)
| | - Barry Make
- National Jewish Health, Denver, Colorado
| | | | - MeiLan Han
- Division of Pulmonary and Critical Care Medicine, University of Michigan, Ann Arbor, Michigan
| | - Craig P Hersh
- Channing Division of Network Medicine and.,Division of Pulmonary and Critical Care Medicine, Brigham and Women's Hospital, Boston, Massachusetts
| | | | - John Quackenbush
- Department of Biostatistics and Computational Biology, Dana-Farber Cancer Institute, Boston, Massachusetts; and
| | - Augustine M K Choi
- Department of Medicine, NewYork-Presbyterian/Weill Cornell Medical Center, New York, New York
| | - Edwin K Silverman
- Channing Division of Network Medicine and.,Division of Pulmonary and Critical Care Medicine, Brigham and Women's Hospital, Boston, Massachusetts
| | - Dawn L DeMeo
- Channing Division of Network Medicine and.,Division of Pulmonary and Critical Care Medicine, Brigham and Women's Hospital, Boston, Massachusetts
| |
Collapse
|
30
|
Gaynor SM, Sun R, Lin X, Quackenbush J. Identification of differentially expressed gene sets using the Generalized Berk-Jones statistic. Bioinformatics 2020; 35:4568-4576. [PMID: 31062858 DOI: 10.1093/bioinformatics/btz277] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2018] [Revised: 02/14/2019] [Accepted: 04/23/2019] [Indexed: 01/09/2023] Open
Abstract
MOTIVATION Cancer genomics studies frequently aim to identify genes that are differentially expressed between clinically distinct patient subgroups, generally by testing single genes one at a time. However, the results of any individual transcriptomic study are often not fully reproducible. A particular challenge impeding statistical analysis is the difficulty of distinguishing between differential expression comprising part of the genomic disease etiology and that induced by downstream effects. More robust analytical approaches that are well-powered to detect potentially causative genes, are less prone to discovering spurious associations, and can deliver reproducible findings across different studies are needed. RESULTS We propose a set-based procedure for testing of differential expression and show that this set-based approach can produce more robust results by aggregating information across multiple, correlated genomic markers. Specifically, we adapt the Generalized Berk-Jones statistic to test for the transcription factors that may contribute to the progression of estrogen receptor positive breast cancer. We demonstrate the ability of our method to produce reproducible findings by applying the same analysis to 21 publicly available datasets, producing a similar list of significant transcription factors across most studies. Our Generalized Berk-Jones approach produces results that show improved consistency over three set-based testing algorithms: Generalized Higher Criticism, Gene Set Analysis and Gene Set Enrichment Analysis. AVAILABILITY AND IMPLEMENTATION Data are in the MetaGxBreast R package. Code is available at github.com/ryanrsun/gaynor_sun_GBJ_breast_cancer. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Sheila M Gaynor
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Harvard University, Boston, MA 02115, USA.,Department of Biostatistics and Computational Biology and Center for Cancer Computational Biology, Dana-Farber Cancer Institute, Boston, MA 02215, USA
| | - Ryan Sun
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Harvard University, Boston, MA 02115, USA
| | - Xihong Lin
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Harvard University, Boston, MA 02115, USA
| | - John Quackenbush
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Harvard University, Boston, MA 02115, USA.,Department of Biostatistics and Computational Biology and Center for Cancer Computational Biology, Dana-Farber Cancer Institute, Boston, MA 02215, USA
| |
Collapse
|
31
|
Chauhan K, Nadkarni GN, Fleming F, McCullough J, He CJ, Quackenbush J, Murphy B, Donovan MJ, Coca SG, Bonventre JV. Initial Validation of a Machine Learning-Derived Prognostic Test (KidneyIntelX) Integrating Biomarkers and Electronic Health Record Data To Predict Longitudinal Kidney Outcomes. ACTA ACUST UNITED AC 2020; 1:731-739. [DOI: 10.34067/kid.0002252020] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2020] [Accepted: 06/25/2020] [Indexed: 11/27/2022]
Abstract
BackgroundIndividuals with type 2 diabetes (T2D) or the apolipoprotein L1 high-risk (APOL1-HR) genotypes are at increased risk of rapid kidney function decline (RKFD) and kidney failure. We hypothesized that a prognostic test using machine learning integrating blood biomarkers and longitudinal electronic health record (EHR) data would improve risk stratification.MethodsWe selected two cohorts from the Mount Sinai BioMe Biobank: T2D (n=871) and African ancestry with APOL1-HR (n=498). We measured plasma tumor necrosis factor receptors (TNFR) 1 and 2 and kidney injury molecule-1 (KIM-1) and used random forest algorithms to integrate biomarker and EHR data to generate a risk score for a composite outcome: RKFD (eGFR decline of ≥5 ml/min per year), or 40% sustained eGFR decline, or kidney failure. We compared performance to a validated clinical model and applied thresholds to assess the utility of the prognostic test (KidneyIntelX) to accurately stratify patients into risk categories.ResultsOverall, 23% of those with T2D and 18% of those with APOL1-HR experienced the composite kidney end point over a median follow-up of 4.6 and 5.9 years, respectively. The area under the receiver operator characteristic curve (AUC) of KidneyIntelX was 0.77 (95% CI, 0.75 to 0.79) in T2D, and 0.80 (95% CI, 0.77 to 0.83) in APOL1-HR, outperforming the clinical models (AUC, 0.66 [95% CI, 0.65 to 0.67] and 0.72 [95% CI, 0.71 to 0.73], respectively; P<0.001). The positive predictive values for KidneyIntelX were 62% and 62% versus 46% and 39% for the clinical models (P<0.01) in high-risk (top 15%) stratum for T2D and APOL1-HR, respectively. The negative predictive values for KidneyIntelX were 92% in T2D and 96% for APOL1-HR versus 85% and 93% for the clinical model, respectively (P=0.76 and 0.93, respectively), in low-risk stratum (bottom 50%).ConclusionsIn patients with T2D or APOL1-HR, a prognostic test (KidneyIntelX) integrating biomarker levels with longitudinal EHR data significantly improved prediction of a composite kidney end point of RKFD, 40% decline in eGFR, or kidney failure over validated clinical models.
Collapse
|
32
|
Altenbuchinger M, Weihs A, Quackenbush J, Grabe HJ, Zacharias HU. Gaussian and Mixed Graphical Models as (multi-)omics data analysis tools. Biochim Biophys Acta Gene Regul Mech 2020; 1863:194418. [PMID: 31639475 PMCID: PMC7166149 DOI: 10.1016/j.bbagrm.2019.194418] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/31/2019] [Revised: 08/21/2019] [Accepted: 08/21/2019] [Indexed: 11/30/2022]
Abstract
Gaussian Graphical Models (GGMs) are tools to infer dependencies between biological variables. Popular applications are the reconstruction of gene, protein, and metabolite association networks. GGMs are an exploratory research tool that can be useful to discover interesting relations between genes (functional clusters) or to identify therapeutically interesting genes, but do not necessarily infer a network in the mechanistic sense. Although GGMs are well investigated from a theoretical and applied perspective, important extensions are not well known within the biological community. GGMs assume, for instance, multivariate normal distributed data. If this assumption is violated Mixed Graphical Models (MGMs) can be the better choice. In this review, we provide the theoretical foundations of GGMs, present extensions such as MGMs or multi-class GGMs, and illustrate how those methods can provide insight in biological mechanisms. We summarize several applications and present user-friendly estimation software. This article is part of a Special Issue entitled: Transcriptional Profiles and Regulatory Gene Networks edited by Dr. Dr. Federico Manuel Giorgi and Dr. Shaun Mahony.
Collapse
Affiliation(s)
- Michael Altenbuchinger
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, MA Boston, 02115, USA.
| | - Antoine Weihs
- Department of Psychiatry and Psychotherapy, University Medicine Greifswald, 17475 Greifswald, Germany
| | - John Quackenbush
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, MA Boston, 02115, USA; Channing Division of Network Medicine, Brigham and Women's Hospital, Boston, MA 02115, USA; Department of Medicine, Harvard Medical School, Boston, MA 02115, USA
| | - Hans Jörgen Grabe
- Department of Psychiatry and Psychotherapy, University Medicine Greifswald, 17475 Greifswald, Germany; German Center for Neurodegenerative Diseases DZNE, Site Rostock/Greifswald, 17475 Greifswald, Germany
| | - Helena U Zacharias
- Department of Psychiatry and Psychotherapy, University Medicine Greifswald, 17475 Greifswald, Germany.
| |
Collapse
|
33
|
Silverman EK, Schmidt HHHW, Anastasiadou E, Altucci L, Angelini M, Badimon L, Balligand JL, Benincasa G, Capasso G, Conte F, Di Costanzo A, Farina L, Fiscon G, Gatto L, Gentili M, Loscalzo J, Marchese C, Napoli C, Paci P, Petti M, Quackenbush J, Tieri P, Viggiano D, Vilahur G, Glass K, Baumbach J. Molecular networks in Network Medicine: Development and applications. Wiley Interdiscip Rev Syst Biol Med 2020; 12:e1489. [PMID: 32307915 DOI: 10.1002/wsbm.1489] [Citation(s) in RCA: 84] [Impact Index Per Article: 21.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/23/2019] [Revised: 02/29/2020] [Accepted: 03/20/2020] [Indexed: 12/14/2022]
Abstract
Network Medicine applies network science approaches to investigate disease pathogenesis. Many different analytical methods have been used to infer relevant molecular networks, including protein-protein interaction networks, correlation-based networks, gene regulatory networks, and Bayesian networks. Network Medicine applies these integrated approaches to Omics Big Data (including genetics, epigenetics, transcriptomics, metabolomics, and proteomics) using computational biology tools and, thereby, has the potential to provide improvements in the diagnosis, prognosis, and treatment of complex diseases. We discuss briefly the types of molecular data that are used in molecular network analyses, survey the analytical methods for inferring molecular networks, and review efforts to validate and visualize molecular networks. Successful applications of molecular network analysis have been reported in pulmonary arterial hypertension, coronary heart disease, diabetes mellitus, chronic lung diseases, and drug development. Important knowledge gaps in Network Medicine include incompleteness of the molecular interactome, challenges in identifying key genes within genetic association regions, and limited applications to human diseases. This article is categorized under: Models of Systems Properties and Processes > Mechanistic Models Translational, Genomic, and Systems Medicine > Translational Medicine Analytical and Computational Methods > Analytical Methods Analytical and Computational Methods > Computational Methods.
Collapse
Affiliation(s)
- Edwin K Silverman
- Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, Massachusetts, USA
| | - Harald H H W Schmidt
- Department of Pharmacology and Personalized Medicine, School of Mental Health and Neuroscience, Faculty of Health, Medicine and Life Science, Maastricht University, Maastricht, The Netherlands
| | - Eleni Anastasiadou
- Department of Experimental Medicine, Sapienza University of Rome, Rome, Italy
| | - Lucia Altucci
- Department of Precision Medicine, University of Campania 'Luigi Vanvitelli', Naples, Italy
| | - Marco Angelini
- Department of Computer, Control and Management Engineering, Sapienza University of Rome, Rome, Italy
| | - Lina Badimon
- Cardiovascular Program-ICCC, IR-Hospital de la Santa Creu i Sant Pau, CiberCV, IIB-Sant Pau, Autonomous University of Barcelona, Barcelona, Spain
| | - Jean-Luc Balligand
- Pole of Pharmacology and Therapeutics (FATH), Institute for Clinical and Experimental Research (IREC), UCLouvain, Brussels, Belgium
| | - Giuditta Benincasa
- Department of Advanced Clinical and Surgical Sciences, University of Campania "Luigi Vanvitelli", Naples, Italy
| | - Giovambattista Capasso
- Department of Translational Medical Sciences, University of Campania "L. Vanvitelli", Naples, Italy.,BIOGEM, Ariano Irpino, Italy
| | - Federica Conte
- Institute for Systems Analysis and Computer Science "Antonio Ruberti", National Research Council, Rome, Italy
| | - Antonella Di Costanzo
- Department of Precision Medicine, University of Campania 'Luigi Vanvitelli', Naples, Italy
| | - Lorenzo Farina
- Department of Computer, Control and Management Engineering, Sapienza University of Rome, Rome, Italy
| | - Giulia Fiscon
- Institute for Systems Analysis and Computer Science "Antonio Ruberti", National Research Council, Rome, Italy
| | - Laurent Gatto
- de Duve Institute, Brussels, Belgium.,Institute for Experimental and Clinical Research (IREC), UCLouvain, Brussels, Belgium
| | - Michele Gentili
- Department of Computer, Control and Management Engineering, Sapienza University of Rome, Rome, Italy
| | - Joseph Loscalzo
- Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, Massachusetts, USA.,Division of Cardiovascular Medicine, Brigham and Women's Hospital, Boston, Massachusetts, USA
| | - Cinzia Marchese
- Department of Experimental Medicine, Sapienza University of Rome, Rome, Italy
| | - Claudio Napoli
- Department of Advanced Clinical and Surgical Sciences, University of Campania "Luigi Vanvitelli", Naples, Italy
| | - Paola Paci
- Department of Computer, Control and Management Engineering, Sapienza University of Rome, Rome, Italy
| | - Manuela Petti
- Department of Computer, Control and Management Engineering, Sapienza University of Rome, Rome, Italy
| | - John Quackenbush
- Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, Massachusetts, USA.,Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, USA
| | - Paolo Tieri
- CNR National Research Council of Italy, IAC Institute for Applied Computing, Rome, Italy
| | - Davide Viggiano
- BIOGEM, Ariano Irpino, Italy.,Department of Medicine and Health Sciences, University of Molise, Campobasso, Italy
| | - Gemma Vilahur
- Cardiovascular Program-ICCC, IR-Hospital de la Santa Creu i Sant Pau, CiberCV, IIB-Sant Pau, Autonomous University of Barcelona, Barcelona, Spain
| | - Kimberly Glass
- Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, Massachusetts, USA.,Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, USA
| | - Jan Baumbach
- Department of Experimental Bioinformatics, TUM School of Life Sciences Weihenstephan, Technical University of Munich, Maximus-von-Imhof-Forum 3, Freising, Germany.,Institute of Mathematics and Computer Science, University of Southern Denmark, Odense, Denmark
| |
Collapse
|
34
|
Qiang J, Ding W, Kuijjer M, Quackenbush J, Chen P. Clustering Sparse Data With Feature Correlation With Application to Discover Subtypes in Cancer. IEEE Access 2020; 8:67775-67789. [PMID: 36329870 PMCID: PMC9629797 DOI: 10.1109/access.2020.2982569] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
In this paper, given data with high-dimensional features, we study this problem of how to calculate the similarity between two samples by considering feature interaction network, where a feature interaction network represents the relationship between features. This is different from some traditional methods, those of which learn similarities based on a sample network that represents the relationship between samples. Therefore, we propose a novel network-based similarity metric for computing the similarity between samples, which incorporates the knowledge of feature interaction network, in order to overcome the data sparseness problem. Our similarity metric uses a new Feature Alignment Similarity measure, which does not directly compute the similarities among samples, but projects each sample into a feature interaction network and measures the similarities between two samples using the similarities between the vertices of the samples in the network. As such, when two samples do not share any common features, they are likely to have higher similarity values when their features share the similar network regions. For ensuring that the metric is useful in a real-world application, we apply our metric to discover subtypes in tumor mutational data by incorporating the information of the gene interaction network. Our experimental results from using synthetic data and real-world tumor mutational data show that our approach outperforms the top competitors in cancer subtype discovery. Furthermore, our approach can identify cancer subtypes that cannot be detected by other clustering algorithms in real cancer data.
Collapse
Affiliation(s)
- Jipeng Qiang
- Department of Computer Science, Yangzhou University, Yangzhou 225127, China
- Department of Computer Science, University of Massachusetts Boston, Boston, MA 02125, USA
| | - Wei Ding
- Department of Computer Science, University of Massachusetts Boston, Boston, MA 02125, USA
| | - Marieke Kuijjer
- Centre for Molecular Medicine Norway, University of Oslo Faculty of Medicine, 0318 Oslo, Norway
| | - John Quackenbush
- Department of Biostatistics, Harvard T. H. Chan School of Public Health, Boston, MA 02115, USA
| | - Ping Chen
- Department of Computer Science, University of Massachusetts Boston, Boston, MA 02125, USA
| |
Collapse
|
35
|
Lietz CE, Garbutt C, Barry WT, Deshpande V, Chen YL, Lozano-Calderon SA, Wang Y, Lawney B, Ebb D, Cote GM, Duan Z, Hornicek FJ, Choy E, Petur Nielsen G, Haibe-Kains B, Quackenbush J, Spentzos D. MicroRNA-mRNA networks define translatable molecular outcome phenotypes in osteosarcoma. Sci Rep 2020; 10:4409. [PMID: 32157112 PMCID: PMC7064533 DOI: 10.1038/s41598-020-61236-3] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2019] [Accepted: 02/03/2020] [Indexed: 12/30/2022] Open
Abstract
There is a lack of well validated prognostic biomarkers in osteosarcoma, a rare, recalcitrant disease for which treatment standards have not changed in over 20 years. We performed microRNA sequencing in 74 frozen osteosarcoma biopsy samples, constituting the largest single center translationally analyzed osteosarcoma cohort to date, and we separately analyzed a multi-omic dataset from a large NCI supported national cooperative group cohort. We validated the prognostic value of candidate microRNA signatures and contextualized them in relevant transcriptomic and epigenomic networks. Our results reveal the existence of molecularly defined phenotypes associated with outcome independent of clinicopathologic features. Through machine learning based integrative pharmacogenomic analysis, the microRNA biomarkers identify novel therapeutics for stratified application in osteosarcoma. The previously unrecognized osteosarcoma subtypes with distinct clinical courses and response to therapy could be translatable for discerning patients appropriate for more intensified, less intensified, or alternate therapeutic regimens.
Collapse
Affiliation(s)
- Christopher E Lietz
- Department of Orthopaedic Surgery, Massachusetts General Hospital, Harvard Medical School, Boston, MA, United States
| | - Cassandra Garbutt
- Department of Orthopaedic Surgery, Massachusetts General Hospital, Harvard Medical School, Boston, MA, United States
- Illumina, Inc., San Diego, United States
| | - William T Barry
- Department of Biostatistics and Computational Biology, Dana Farber Cancer Institute, Harvard Medical School, Boston, MA, United States
| | - Vikram Deshpande
- Department of Pathology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, United States
| | - Yen-Lin Chen
- Department of Radiation Oncology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, United States
| | - Santiago A Lozano-Calderon
- Department of Orthopaedic Surgery, Massachusetts General Hospital, Harvard Medical School, Boston, MA, United States
| | - Yaoyu Wang
- Department of Biostatistics and Computational Biology, Dana Farber Cancer Institute, Harvard Medical School, Boston, MA, United States
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, United States
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, United States
| | - Brian Lawney
- Department of Biostatistics and Computational Biology, Dana Farber Cancer Institute, Harvard Medical School, Boston, MA, United States
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, United States
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, United States
| | - David Ebb
- Pediatric Hematology-Oncology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, United States
| | - Gregory M Cote
- Department of Hematology/Oncology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, United States
| | - Zhenfeng Duan
- Department of Orthopaedic Surgery, UCLA, Los Angeles, CA, United States
| | | | - Edwin Choy
- Department of Hematology/Oncology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, United States
| | - G Petur Nielsen
- Department of Pathology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, United States
| | - Benjamin Haibe-Kains
- Department of Medical Biophysics, Princess Margaret Cancer Centre, University of Toronto, Toronto, Canada
| | - John Quackenbush
- Department of Biostatistics and Computational Biology, Dana Farber Cancer Institute, Harvard Medical School, Boston, MA, United States
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, United States
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, United States
| | - Dimitrios Spentzos
- Department of Orthopaedic Surgery, Massachusetts General Hospital, Harvard Medical School, Boston, MA, United States.
| |
Collapse
|
36
|
Schwede M, Waldron L, Mok SC, Wei W, Basunia A, Merritt MA, Mitsiades CS, Parmigiani G, Harrington DP, Quackenbush J, Birrer MJ, Culhane AC. The Impact of Stroma Admixture on Molecular Subtypes and Prognostic Gene Signatures in Serous Ovarian Cancer. Cancer Epidemiol Biomarkers Prev 2019; 29:509-519. [PMID: 31871106 DOI: 10.1158/1055-9965.epi-18-1359] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2018] [Revised: 04/26/2019] [Accepted: 12/06/2019] [Indexed: 12/18/2022] Open
Abstract
BACKGROUND Recent efforts to improve outcomes for high-grade serous ovarian cancer, a leading cause of cancer death in women, have focused on identifying molecular subtypes and prognostic gene signatures, but existing subtypes have poor cross-study robustness. We tested the contribution of cell admixture in published ovarian cancer molecular subtypes and prognostic gene signatures. METHODS Gene signatures of tumor and stroma were developed using paired microdissected tissue from two independent studies. Stromal genes were investigated in two molecular subtype classifications and 61 published gene signatures. Prognostic performance of gene signatures of stromal admixture was evaluated in 2,527 ovarian tumors (16 studies). Computational simulations of increasing stromal cell proportion were performed by mixing gene-expression profiles of paired microdissected ovarian tumor and stroma. RESULTS Recently described ovarian cancer molecular subtypes are strongly associated with the cell admixture. Tumors were classified as different molecular subtypes in simulations where the percentage of stromal cells increased. Stromal gene expression in bulk tumors was associated with overall survival (hazard ratio, 1.17; 95% confidence interval, 1.11-1.23), and in one data set, increased stroma was associated with anatomic sampling location. Five published prognostic gene signatures were no longer prognostic in a multivariate model that adjusted for stromal content. CONCLUSIONS Cell admixture affects the interpretation and reproduction of ovarian cancer molecular subtypes and gene signatures derived from bulk tissue. Elucidating the role of stroma in the tumor microenvironment and in prognosis is important. IMPACT Single-cell analyses may be required to refine the molecular subtypes of high-grade serous ovarian cancer.
Collapse
Affiliation(s)
- Matthew Schwede
- Department of Data Sciences, Dana-Farber Cancer Institute, Boston, Massachusetts
| | - Levi Waldron
- Biostatistics, CUNY Graduate School of Public Health and Health Policy, New York, New York
| | - Samuel C Mok
- Department of Gynecologic Oncology and Reproductive Medicine, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Wei Wei
- Pfizer, Andover, Massachusetts
| | - Azfar Basunia
- Department of Data Sciences, Dana-Farber Cancer Institute, Boston, Massachusetts.,Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, Massachusetts
| | | | | | - Giovanni Parmigiani
- Department of Data Sciences, Dana-Farber Cancer Institute, Boston, Massachusetts.,Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, Massachusetts
| | - David P Harrington
- Department of Data Sciences, Dana-Farber Cancer Institute, Boston, Massachusetts.,Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, Massachusetts
| | - John Quackenbush
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, Massachusetts
| | - Michael J Birrer
- Division of Hematology-Oncology, University of Alabama at Birmingham, Birmingham, Alabama.
| | - Aedín C Culhane
- Department of Data Sciences, Dana-Farber Cancer Institute, Boston, Massachusetts. .,Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, Massachusetts
| |
Collapse
|
37
|
Fagny M, Platig J, Kuijjer ML, Lin X, Quackenbush J. Nongenic cancer-risk SNPs affect oncogenes, tumour-suppressor genes, and immune function. Br J Cancer 2019; 122:569-577. [PMID: 31806877 PMCID: PMC7028992 DOI: 10.1038/s41416-019-0614-3] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2019] [Revised: 09/23/2019] [Accepted: 10/07/2019] [Indexed: 12/31/2022] Open
Abstract
Background Genome-wide association studies (GWASes) have identified many noncoding germline single-nucleotide polymorphisms (SNPs) that are associated with an increased risk of developing cancer. However, how these SNPs affect cancer risk is still largely unknown. Methods We used a systems biology approach to analyse the regulatory role of cancer-risk SNPs in thirteen tissues. By using data from the Genotype-Tissue Expression (GTEx) project, we performed an expression quantitative trait locus (eQTL) analysis. We represented both significant cis- and trans-eQTLs as edges in tissue-specific eQTL bipartite networks. Results Each tissue-specific eQTL network is organised into communities that group sets of SNPs and functionally related genes. When mapping cancer-risk SNPs to these networks, we find that in each tissue, these SNPs are significantly overrepresented in communities enriched for immune response processes, as well as tissue-specific functions. Moreover, cancer-risk SNPs are more likely to be ‘cores’ of their communities, influencing the expression of many genes within the same biological processes. Finally, cancer-risk SNPs preferentially target oncogenes and tumour-suppressor genes, suggesting that they may alter the expression of these key cancer genes. Conclusions This approach provides a new way of understanding genetic effects on cancer risk and provides a biological context for interpreting the results of GWAS cancer studies.
Collapse
Affiliation(s)
- Maud Fagny
- Genetique Quantitative et Evolution-Le Moulon, Institut National de la Recherche agronomique, Université Paris-Sud, Centre National de la Recherche Scientifique, AgroParisTech, Université Paris-Saclay, Paris, France
| | - John Platig
- Channing Division of Network Medicine, Brigham and Women's Hospital, Boston, MA, USA.,Harvard Medical School, Boston, MA, USA
| | - Marieke Lydia Kuijjer
- Department of Biostatistics and Computational Biology, Dana-Farber Cancer Institute, Boston, MA, USA.,Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA.,Centre for Molecular Medicine Norway, University of Oslo, Oslo, Norway
| | - Xihong Lin
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | - John Quackenbush
- Channing Division of Network Medicine, Brigham and Women's Hospital, Boston, MA, USA. .,Department of Biostatistics and Computational Biology, Dana-Farber Cancer Institute, Boston, MA, USA. .,Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA. .,Department of Cancer Biology, Dana-Farber Cancer Institute, Boston, MA, USA.
| |
Collapse
|
38
|
St Hilaire MA, Kristal BS, Rahman SA, Sullivan JP, Quackenbush J, Duffy JF, Barger LK, Gooley JJ, Czeisler CA, Lockley SW. Using a Single Daytime Performance Test to Identify Most Individuals at High-Risk for Performance Impairment during Extended Wake. Sci Rep 2019; 9:16681. [PMID: 31723161 PMCID: PMC6853981 DOI: 10.1038/s41598-019-52930-y] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2018] [Accepted: 10/25/2019] [Indexed: 12/20/2022] Open
Abstract
We explored the predictive value of a neurobehavioral performance assessment under rested baseline conditions (evaluated at 8 hours awake following 8 hours of sleep) on neurobehavioral response to moderate sleep loss (evaluated at 20 hours awake two days later) in 151 healthy young participants (18-30 years). We defined each participant's response-to-sleep-loss phenotype based on the number of attentional failures on a 10-min visual psychomotor vigilance task taken at 20 hours awake (resilient: less than 6 attentional failures, n = 26 participants; non-resilient: 6 or more attentional failures, n = 125 participants). We observed that 97% of rested participants with 2 or more attentional failures (n = 73 of 151) and 100% of rested participants with 3 or more attentional failures (n = 57 of 151) were non-resilient after moderate sleep loss. Our approach can accurately identify a significant proportion of individuals who are at high risk for neurobehavioral performance impairment from staying up late with a single neurobehavioral performance assessment conducted during rested conditions. Additional methods are needed to predict the future performance of individuals who are not identified as high risk during baseline.
Collapse
Affiliation(s)
- Melissa A St Hilaire
- Division of Sleep and Circadian Disorders, Brigham & Women's Hospital, 221 Longwood Avenue, Boston, MA, 02115, USA.
- Division of Sleep Medicine, Harvard Medical School, 221 Longwood Avenue, Boston, MA, 02115, USA.
| | - Bruce S Kristal
- Division of Sleep and Circadian Disorders, Brigham & Women's Hospital, 221 Longwood Avenue, Boston, MA, 02115, USA
- Division of Sleep Medicine, Harvard Medical School, 221 Longwood Avenue, Boston, MA, 02115, USA
| | - Shadab A Rahman
- Division of Sleep and Circadian Disorders, Brigham & Women's Hospital, 221 Longwood Avenue, Boston, MA, 02115, USA
- Division of Sleep Medicine, Harvard Medical School, 221 Longwood Avenue, Boston, MA, 02115, USA
| | - Jason P Sullivan
- Division of Sleep and Circadian Disorders, Brigham & Women's Hospital, 221 Longwood Avenue, Boston, MA, 02115, USA
| | - John Quackenbush
- Biostatistics and Computational Biology, Dana-Farber Cancer Institute, 450 Brookline Avenue, Boston, MA, 02215, USA
| | - Jeanne F Duffy
- Division of Sleep and Circadian Disorders, Brigham & Women's Hospital, 221 Longwood Avenue, Boston, MA, 02115, USA
- Division of Sleep Medicine, Harvard Medical School, 221 Longwood Avenue, Boston, MA, 02115, USA
| | - Laura K Barger
- Division of Sleep and Circadian Disorders, Brigham & Women's Hospital, 221 Longwood Avenue, Boston, MA, 02115, USA
- Division of Sleep Medicine, Harvard Medical School, 221 Longwood Avenue, Boston, MA, 02115, USA
| | - Joshua J Gooley
- Programme in Neuroscience and Behavioural Disorders, Duke-National University of Singapore Medical School, 8 College Road, Singapore, 169857, Singapore
| | - Charles A Czeisler
- Division of Sleep and Circadian Disorders, Brigham & Women's Hospital, 221 Longwood Avenue, Boston, MA, 02115, USA
- Division of Sleep Medicine, Harvard Medical School, 221 Longwood Avenue, Boston, MA, 02115, USA
| | - Steven W Lockley
- Division of Sleep and Circadian Disorders, Brigham & Women's Hospital, 221 Longwood Avenue, Boston, MA, 02115, USA
- Division of Sleep Medicine, Harvard Medical School, 221 Longwood Avenue, Boston, MA, 02115, USA
| |
Collapse
|
39
|
Abstract
BACKGROUND In biomedical research, network inference algorithms are typically used to infer complex association patterns between biological entities, such as between genes or proteins, using data from a population. This resulting aggregate network, in essence, averages over the networks of those individuals in the population. LIONESS (Linear Interpolation to Obtain Network Estimates for Single Samples) is a method that can be used together with a network inference algorithm to extract networks for individual samples in a population. The method's key characteristic is that, by modeling networks for individual samples in a data set, it can capture network heterogeneity in a population. LIONESS was originally made available as a function within the PANDA (Passing Attributes between Networks for Data Assimilation) regulatory network reconstruction framework. However, the LIONESS algorithm is generalizable and can be used to model single sample networks based on a wide range of network inference algorithms. RESULTS In this software article, we describe lionessR, an R implementation of LIONESS that can be applied to any network inference method in R that outputs a complete, weighted adjacency matrix. As an example, we provide a vignette of an application of lionessR to model single sample networks based on correlated gene expression in a bone cancer dataset. We show how the tool can be used to identify differential patterns of correlation between two groups of patients. CONCLUSIONS We developed lionessR, an open source R package to model single sample networks. We show how lionessR can be used to inform us on potential precision medicine applications in cancer. The lionessR package is a user-friendly tool to perform such analyses. The package, which includes a vignette describing the application, is freely available at: https://github.com/kuijjerlab/lionessR and at: http://bioconductor.org/packages/lionessR .
Collapse
Affiliation(s)
- Marieke L Kuijjer
- Centre for Molecular Medicine Norway (NCMM), Nordic EMBL Partnership, University of Oslo, Gaustadalléen 21, Oslo, 0318, Norway.
| | - Ping-Han Hsieh
- Centre for Molecular Medicine Norway (NCMM), Nordic EMBL Partnership, University of Oslo, Gaustadalléen 21, Oslo, 0318, Norway
| | - John Quackenbush
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, 677 Huntington Ave, Boston, 02215, USA.,Channing Division of Network Medicine, Brigham and Women's Hospital, 181 Longwood Avenue, Boston, 02215, USA.,Department of Medicine, Harvard Medical School, 25 Shattuck Street, Boston, 02215, USA
| | - Kimberly Glass
- Channing Division of Network Medicine, Brigham and Women's Hospital, 181 Longwood Avenue, Boston, 02215, USA.,Department of Medicine, Harvard Medical School, 25 Shattuck Street, Boston, 02215, USA
| |
Collapse
|
40
|
Campbell PT, Ambrosone CB, Nishihara R, Aerts HJWL, Bondy M, Chatterjee N, Garcia-Closas M, Giannakis M, Golden JA, Heng YJ, Kip NS, Koshiol J, Liu XS, Lopes-Ramos CM, Mucci LA, Nowak JA, Phipps AI, Quackenbush J, Schoen RE, Sholl LM, Tamimi RM, Wang M, Weijenberg MP, Wu CJ, Wu K, Yao S, Yu KH, Zhang X, Rebbeck TR, Ogino S. Proceedings of the fourth international molecular pathological epidemiology (MPE) meeting. Cancer Causes Control 2019; 30:799-811. [PMID: 31069578 PMCID: PMC6614001 DOI: 10.1007/s10552-019-01177-z] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2018] [Accepted: 04/27/2019] [Indexed: 02/06/2023]
Abstract
An important premise of epidemiology is that individuals with the same disease share similar underlying etiologies and clinical outcomes. In the past few decades, our knowledge of disease pathogenesis has improved, and disease classification systems have evolved to the point where no complex disease processes are considered homogenous. As a result, pathology and epidemiology have been integrated into the single, unified field of molecular pathological epidemiology (MPE). Advancing integrative molecular and population-level health sciences and addressing the unique research challenges specific to the field of MPE necessitates assembling experts in diverse fields, including epidemiology, pathology, biostatistics, computational biology, bioinformatics, genomics, immunology, and nutritional and environmental sciences. Integrating these seemingly divergent fields can lead to a greater understanding of pathogenic processes. The International MPE Meeting Series fosters discussion that addresses the specific research questions and challenges in this emerging field. The purpose of the meeting series is to: discuss novel methods to integrate pathology and epidemiology; discuss studies that provide pathogenic insights into population impact; and educate next-generation scientists. Herein, we share the proceedings of the Fourth International MPE Meeting, held in Boston, MA, USA, on 30 May-1 June, 2018. Major themes of this meeting included 'integrated genetic and molecular pathologic epidemiology', 'immunology-MPE', and 'novel disease phenotyping'. The key priority areas for future research identified by meeting attendees included integration of tumor immunology and cancer disparities into epidemiologic studies, further collaboration between computational and population-level scientists to gain new insight on exposure-disease associations, and future pooling projects of studies with comparable data.
Collapse
Affiliation(s)
- Peter T Campbell
- Behavioral and Epidemiology Research Group, American Cancer Society, 250 Williams Street NW, Atlanta, GA, 30303, USA.
| | - Christine B Ambrosone
- Department of Cancer Prevention and Control, Roswell Park Comprehensive Cancer Center, Buffalo, NY, USA
| | - Reiko Nishihara
- Program in MPE Molecular Pathological Epidemiology, Department of Pathology, Brigham and Women's Hospital and Harvard Medical School, 450 Brookline Ave, Room SM1036, Boston, MA, 02215, USA
- Department of Nutrition, Harvard T.H. Chan School of Public Health, Boston, MA, USA
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | - Hugo J W L Aerts
- Departments of Radiation Oncology and Radiology, Dana-Farber Cancer Institute, Brigham and Women's Hospital, Harvard Medical School, Boston, MA, USA
| | - Melissa Bondy
- Cancer Prevention and Population Sciences, Dan L. Duncan Comprehensive Cancer Center, Epidemiology and Population Sciences, Department of Medicine, Baylor College of Medicine, Houston, TX, USA
| | - Nilanjan Chatterjee
- Department of Biostatistics, Bloomberg School of Public Health, Department of Oncology, School of Medicine, Johns Hopkins University, Baltimore, MD, USA
| | - Montserrat Garcia-Closas
- Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Bethesda, MD, USA
| | - Marios Giannakis
- Department of Medical Oncology, Dana-Farber Cancer Institute and Harvard Medical School, Boston, MA, USA
- Broad Institute of Harvard & MIT, Cambridge, MA, USA
| | - Jeffrey A Golden
- Department of Pathology, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA
| | - Yujing J Heng
- Department of Pathology, Beth Israel Deaconess Medical Center, Harvard Medical School, Boston, MA, USA
| | - N Sertac Kip
- Sema4, Mount Sinai Icahn School of Medicine, Genetics & Genomic Sciences and Pathology, Branford, CT, USA
| | - Jill Koshiol
- Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Bethesda, MD, USA
| | - X Shirley Liu
- Department of Biostatistics and Computational Biology, Dana-Farber Cancer Institute and Harvard T.H. Chan School of Public Health, Boston, MA, 02215, USA
| | - Camila M Lopes-Ramos
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | - Lorelei A Mucci
- Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | - Jonathan A Nowak
- Department of Pathology, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA
| | - Amanda I Phipps
- Department of Epidemiology, University of Washington, Seattle, WA, USA
| | - John Quackenbush
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | - Robert E Schoen
- Departments of Medicine and Epidemiology, The University of Pittsburgh, Pittsburgh, PA, USA
| | - Lynette M Sholl
- Department of Pathology, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA
| | - Rulla M Tamimi
- Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, MA, USA
- Channing Division of Network Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA
| | - Molin Wang
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
- Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, MA, USA
- Channing Division of Network Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA
| | - Matty P Weijenberg
- Department of Epidemiology, GROW-School for Oncology and Developmental Biology, Maastricht University, Maastricht, The Netherlands
| | - Catherine J Wu
- Department of Medical Oncology, Dana-Farber Cancer Institute and Harvard Medical School, Boston, MA, USA
- Broad Institute of Harvard & MIT, Cambridge, MA, USA
| | - Kana Wu
- Department of Nutrition, Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | - Song Yao
- Department of Cancer Prevention and Control, Roswell Park Comprehensive Cancer Center, Buffalo, NY, USA
| | - Kun-Hsing Yu
- Department of Pathology, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA
- Department of Biomedical Informatics, Harvard Medical School, Boston, MA, USA
| | - Xuehong Zhang
- Channing Division of Network Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA
| | - Timothy R Rebbeck
- Department of Medical Oncology, Dana-Farber Cancer Institute and Harvard Medical School, Boston, MA, USA
- Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | - Shuji Ogino
- Program in MPE Molecular Pathological Epidemiology, Department of Pathology, Brigham and Women's Hospital and Harvard Medical School, 450 Brookline Ave, Room SM1036, Boston, MA, 02215, USA.
- Broad Institute of Harvard & MIT, Cambridge, MA, USA.
- Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, MA, USA.
- Department of Oncologic Pathology, Dana-Farber Cancer Institute, Boston, MA, USA.
| |
Collapse
|
41
|
Hicks SC, Okrah K, Paulson JN, Quackenbush J, Irizarry RA, Bravo HC. Smooth quantile normalization. Biostatistics 2019; 19:185-198. [PMID: 29036413 DOI: 10.1093/biostatistics/kxx028] [Citation(s) in RCA: 56] [Impact Index Per Article: 11.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2016] [Accepted: 05/07/2017] [Indexed: 11/14/2022] Open
Abstract
Between-sample normalization is a critical step in genomic data analysis to remove systematic bias and unwanted technical variation in high-throughput data. Global normalization methods are based on the assumption that observed variability in global properties is due to technical reasons and are unrelated to the biology of interest. For example, some methods correct for differences in sequencing read counts by scaling features to have similar median values across samples, but these fail to reduce other forms of unwanted technical variation. Methods such as quantile normalization transform the statistical distributions across samples to be the same and assume global differences in the distribution are induced by only technical variation. However, it remains unclear how to proceed with normalization if these assumptions are violated, for example, if there are global differences in the statistical distributions between biological conditions or groups, and external information, such as negative or control features, is not available. Here, we introduce a generalization of quantile normalization, referred to as smooth quantile normalization (qsmooth), which is based on the assumption that the statistical distribution of each sample should be the same (or have the same distributional shape) within biological groups or conditions, but allowing that they may differ between groups. We illustrate the advantages of our method on several high-throughput datasets with global differences in distributions corresponding to different biological conditions. We also perform a Monte Carlo simulation study to illustrate the bias-variance tradeoff and root mean squared error of qsmooth compared to other global normalization methods. A software implementation is available from https://github.com/stephaniehicks/qsmooth.
Collapse
Affiliation(s)
- Stephanie C Hicks
- Department of Biostatistics and Computational Biology, Dana-Farber Cancer Institute, 450 Brookline Ave, Boston, MA 02215, USA and Department of Biostatistics, Harvard T.H. Chan School of Public Health, 677 Huntington Ave, Boston, MA 02115, USA
| | - Kwame Okrah
- Genetech, Product Development Biostatistics, 1 DNA Way, South San Francisco, CA 94080, USA
| | - Joseph N Paulson
- Department of Biostatistics and Computational Biology, Dana-Farber Cancer Institute, 450 Brookline Ave, Boston, MA 02215, USA and Department of Biostatistics, Harvard T.H. Chan School of Public Health, 677 Huntington Ave, Boston, MA 02115, USA
| | - John Quackenbush
- Department of Biostatistics and Computational Biology, Dana-Farber Cancer Institute, 450 Brookline Ave, Boston, MA 02215, USA and Department of Biostatistics, Harvard T.H. Chan School of Public Health, 677 Huntington Ave, Boston, MA 02115, USA
| | - Rafael A Irizarry
- Department of Biostatistics and Computational Biology, Dana-Farber Cancer Institute, 450 Brookline Ave, Boston, MA 02215, USA and Department of Biostatistics, Harvard T.H. Chan School of Public Health, 677 Huntington Ave, Boston, MA 02115, USA
| | - Héctor Corrada Bravo
- Department of Computer Science, University of Maryland, College Park, USA and Center for Bioinformatics and Computational Biology, Institute of Advanced Computer Studies, University of Maryland, 8314 Paint Branch Dr., College Park, MD 20742, College Park, USA
| |
Collapse
|
42
|
Kuijjer ML, Tung MG, Yuan G, Quackenbush J, Glass K. Estimating Sample-Specific Regulatory Networks. iScience 2019; 14:226-240. [PMID: 30981959 PMCID: PMC6463816 DOI: 10.1016/j.isci.2019.03.021] [Citation(s) in RCA: 78] [Impact Index Per Article: 15.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2018] [Revised: 01/30/2019] [Accepted: 03/21/2019] [Indexed: 10/28/2022] Open
Abstract
Biological systems are driven by intricate interactions among molecules. Many methods have been developed that draw on large numbers of expression samples to infer connections between genes (or their products). The result is an aggregate network representing a single estimate for the likelihood of each interaction, or "edge," in the network. Although informative, aggregate models fail to capture population heterogeneity. Here we propose a method to reverse engineer sample-specific networks from aggregate networks. We demonstrate our approach in several contexts, including simulated, yeast microarray, and human lymphoblastoid cell line RNA sequencing data. We use these sample-specific networks to study changes in network topology across time and to characterize shifts in gene regulation that were not apparent in the expression data. We believe that generating sample-specific networks will greatly facilitate the application of network methods to large, complex, and heterogeneous multi-omic datasets, supporting the emerging field of precision network medicine.
Collapse
Affiliation(s)
- Marieke Lydia Kuijjer
- Centre for Molecular Medicine Norway (NCMM), Nordic EMBL Partnership, University of Oslo, 0318 Oslo, Norway
| | - Matthew George Tung
- Department of Anesthesiology, Critical Care, and Pain Medicine, Massachusetts General Hospital, Boston, MA 02114, USA
| | - GuoCheng Yuan
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA 02115, USA; Department of Biostatistics and Computational Biology, Dana-Farber Cancer Institute, Boston, MA 02215, USA
| | - John Quackenbush
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA 02115, USA; Channing Division of Network Medicine, Brigham and Women's Hospital, Boston, MA 02115, USA; Department of Medicine, Harvard Medical School, Boston, MA 02115, USA
| | - Kimberly Glass
- Channing Division of Network Medicine, Brigham and Women's Hospital, Boston, MA 02115, USA; Department of Medicine, Harvard Medical School, Boston, MA 02115, USA.
| |
Collapse
|
43
|
Lopes-Ramos CM, Kuijjer ML, Ogino S, Fuchs CS, DeMeo DL, Glass K, Quackenbush J. Gene Regulatory Network Analysis Identifies Sex-Linked Differences in Colon Cancer Drug Metabolism. Cancer Res 2018; 78:5538-5547. [PMID: 30275053 PMCID: PMC6169995 DOI: 10.1158/0008-5472.can-18-0454] [Citation(s) in RCA: 58] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2018] [Revised: 06/04/2018] [Accepted: 07/20/2018] [Indexed: 12/12/2022]
Abstract
Understanding sex differences in colon cancer is essential to advance disease prevention, diagnosis, and treatment. Males have a higher risk of developing colon cancer and a lower survival rate than women. However, the molecular features that drive these sex differences are poorly understood. In this study, we use both transcript-based and gene regulatory network methods to analyze RNA-seq data from The Cancer Genome Atlas for 445 patients with colon cancer. We compared gene expression between tumors in men and women and observed significant sex differences in sex chromosome genes only. We then inferred patient-specific gene regulatory networks and found significant regulatory differences between males and females, with drug and xenobiotics metabolism via cytochrome P450 pathways more strongly targeted in females. This finding was validated in a dataset of 1,193 patients from five independent studies. While targeting, the drug metabolism pathway did not change overall survival for males treated with adjuvant chemotherapy, females with greater targeting showed an increase in 10-year overall survival probability, 89% [95% confidence interval (CI), 78-100] survival compared with 61% (95% CI, 45-82) for women with lower targeting, respectively (P = 0.034). Our network analysis uncovers patterns of transcriptional regulation that differentiate male and female colon cancer and identifies differences in regulatory processes involving the drug metabolism pathway associated with survival in women who receive adjuvant chemotherapy. This approach can be used to investigate the molecular features that drive sex differences in other cancers and complex diseases.Significance: A network-based approach reveals that sex-specific patterns of gene targeting by transcriptional regulators are associated with survival outcome in colon cancer. This approach can be used to understand how sex influences progression and response to therapies in other cancers. Cancer Res; 78(19); 5538-47. ©2018 AACR.
Collapse
Affiliation(s)
- Camila M Lopes-Ramos
- Department of Biostatistics and Computational Biology, Dana-Farber Cancer Institute, Boston, Massachusetts
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, Massachusetts
| | - Marieke L Kuijjer
- Department of Biostatistics and Computational Biology, Dana-Farber Cancer Institute, Boston, Massachusetts
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, Massachusetts
| | - Shuji Ogino
- Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, Massachusetts
- Program in MPE Molecular Pathological Epidemiology, Department of Pathology, Brigham and Women's Hospital and Harvard Medical School, Boston, Massachusetts
- Department of Oncologic Pathology, Dana-Farber Cancer Institute, Boston, Massachusetts
- Broad Institute of MIT and Harvard, Cambridge, Massachusetts
| | - Charles S Fuchs
- Yale Cancer Center, New Haven, Connecticut
- Department of Medicine, Yale School of Medicine, New Haven, Connecticut
- Smilow Cancer Hospital, New Haven, Connecticut
| | - Dawn L DeMeo
- Channing Division of Network Medicine, Brigham and Women's Hospital, and Harvard Medical School, Boston, Massachusetts
- Division of Pulmonary and Critical Care Medicine, Brigham and Women's Hospital, Boston, Massachusetts
| | - Kimberly Glass
- Channing Division of Network Medicine, Brigham and Women's Hospital, and Harvard Medical School, Boston, Massachusetts
| | - John Quackenbush
- Department of Biostatistics and Computational Biology, Dana-Farber Cancer Institute, Boston, Massachusetts.
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, Massachusetts
- Department of Cancer Biology, Dana-Farber Cancer Institute, Boston, Massachusetts
| |
Collapse
|
44
|
Parmar C, Barry JD, Hosny A, Quackenbush J, Aerts HJWL. Data Analysis Strategies in Medical Imaging. Clin Cancer Res 2018; 24:3492-3499. [PMID: 29581134 PMCID: PMC6082690 DOI: 10.1158/1078-0432.ccr-18-0385] [Citation(s) in RCA: 93] [Impact Index Per Article: 15.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2018] [Revised: 02/26/2018] [Accepted: 03/22/2018] [Indexed: 12/27/2022]
Abstract
Radiographic imaging continues to be one of the most effective and clinically useful tools within oncology. Sophistication of artificial intelligence has allowed for detailed quantification of radiographic characteristics of tissues using predefined engineered algorithms or deep learning methods. Precedents in radiology as well as a wealth of research studies hint at the clinical relevance of these characteristics. However, critical challenges are associated with the analysis of medical imaging data. Although some of these challenges are specific to the imaging field, many others like reproducibility and batch effects are generic and have already been addressed in other quantitative fields such as genomics. Here, we identify these pitfalls and provide recommendations for analysis strategies of medical imaging data, including data normalization, development of robust models, and rigorous statistical analyses. Adhering to these recommendations will not only improve analysis quality but also enhance precision medicine by allowing better integration of imaging data with other biomedical data sources. Clin Cancer Res; 24(15); 3492-9. ©2018 AACR.
Collapse
Affiliation(s)
- Chintan Parmar
- Department of Radiation Oncology, Dana-Farber Cancer Institute, Brigham and Women's Hospital, Harvard Medical School, Boston, Massachusetts
| | - Joseph D Barry
- Department of Biostatistics & Computational Biology, Dana-Farber Cancer Institute, Boston, Massachusetts
| | - Ahmed Hosny
- Department of Radiation Oncology, Dana-Farber Cancer Institute, Brigham and Women's Hospital, Harvard Medical School, Boston, Massachusetts
| | - John Quackenbush
- Department of Biostatistics & Computational Biology, Dana-Farber Cancer Institute, Boston, Massachusetts
- Department of Cancer Biology, Dana-Farber Cancer Institute, Boston, Massachusetts
| | - Hugo J W L Aerts
- Department of Radiation Oncology, Dana-Farber Cancer Institute, Brigham and Women's Hospital, Harvard Medical School, Boston, Massachusetts.
- Department of Radiology, Dana-Farber Cancer Institute, Brigham and Women's Hospital, Harvard Medical School, Boston, Massachusetts
| |
Collapse
|
45
|
Abstract
Artificial intelligence (AI) algorithms, particularly deep learning, have demonstrated remarkable progress in image-recognition tasks. Methods ranging from convolutional neural networks to variational autoencoders have found myriad applications in the medical image analysis field, propelling it forward at a rapid pace. Historically, in radiology practice, trained physicians visually assessed medical images for the detection, characterization and monitoring of diseases. AI methods excel at automatically recognizing complex patterns in imaging data and providing quantitative, rather than qualitative, assessments of radiographic characteristics. In this Opinion article, we establish a general understanding of AI methods, particularly those pertaining to image-based tasks. We explore how these methods could impact multiple facets of radiology, with a general focus on applications in oncology, and demonstrate ways in which these methods are advancing the field. Finally, we discuss the challenges facing clinical implementation and provide our perspective on how the domain could be advanced.
Collapse
Affiliation(s)
- Ahmed Hosny
- Department of Radiation Oncology, Dana-Farber Cancer Institute, Harvard Medical School, Boston, MA, USA
| | - Chintan Parmar
- Department of Radiation Oncology, Dana-Farber Cancer Institute, Harvard Medical School, Boston, MA, USA
| | - John Quackenbush
- Department of Biostatistics & Computational Biology, Dana-Farber Cancer Institute, Boston, MA, USA
- Department of Cancer Biology, Dana-Farber Cancer Institute, Boston, MA, USA
| | - Lawrence H Schwartz
- Department of Radiology, Columbia University College of Physicians and Surgeons, New York, NY, USA
- Department of Radiology, New York Presbyterian Hospital, New York, NY, USA
| | - Hugo J W L Aerts
- Department of Radiation Oncology, Dana-Farber Cancer Institute, Harvard Medical School, Boston, MA, USA.
- Department of Radiology, Dana-Farber Cancer Institute, Brigham and Women's Hospital, Harvard Medical School, Boston, MA, USA.
| |
Collapse
|
46
|
Abstract
Artificial intelligence (AI) algorithms, particularly deep learning, have demonstrated remarkable progress in image-recognition tasks. Methods ranging from convolutional neural networks to variational autoencoders have found myriad applications in the medical image analysis field, propelling it forward at a rapid pace. Historically, in radiology practice, trained physicians visually assessed medical images for the detection, characterization and monitoring of diseases. AI methods excel at automatically recognizing complex patterns in imaging data and providing quantitative, rather than qualitative, assessments of radiographic characteristics. In this Opinion article, we establish a general understanding of AI methods, particularly those pertaining to image-based tasks. We explore how these methods could impact multiple facets of radiology, with a general focus on applications in oncology, and demonstrate ways in which these methods are advancing the field. Finally, we discuss the challenges facing clinical implementation and provide our perspective on how the domain could be advanced.
Collapse
Affiliation(s)
- Ahmed Hosny
- Department of Radiation Oncology, Dana-Farber Cancer Institute, Harvard Medical School, Boston, MA, USA
| | - Chintan Parmar
- Department of Radiation Oncology, Dana-Farber Cancer Institute, Harvard Medical School, Boston, MA, USA
| | - John Quackenbush
- Department of Biostatistics & Computational Biology, Dana-Farber Cancer Institute, Boston, MA, USA
- Department of Cancer Biology, Dana-Farber Cancer Institute, Boston, MA, USA
| | - Lawrence H Schwartz
- Department of Radiology, Columbia University College of Physicians and Surgeons, New York, NY, USA
- Department of Radiology, New York Presbyterian Hospital, New York, NY, USA
| | - Hugo J W L Aerts
- Department of Radiation Oncology, Dana-Farber Cancer Institute, Harvard Medical School, Boston, MA, USA.
- Department of Radiology, Dana-Farber Cancer Institute, Brigham and Women's Hospital, Harvard Medical School, Boston, MA, USA.
| |
Collapse
|
47
|
Barry JD, Fagny M, Paulson JN, Aerts HJWL, Platig J, Quackenbush J. Histopathological Image QTL Discovery of Immune Infiltration Variants. iScience 2018; 5:80-89. [PMID: 30240647 PMCID: PMC6123851 DOI: 10.1016/j.isci.2018.07.001] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2018] [Revised: 05/30/2018] [Accepted: 07/03/2018] [Indexed: 12/20/2022] Open
Abstract
Genotype-to-phenotype association studies typically use macroscopic physiological measurements or molecular readouts as quantitative traits. There are comparatively few suitable quantitative traits available between cell and tissue length scales, a limitation that hinders our ability to identify variants affecting phenotype at many clinically informative levels. Here we show that quantitative image features, automatically extracted from histopathological imaging data, can be used for image quantitative trait loci (iQTLs) mapping and variant discovery. Using thyroid pathology images, clinical metadata, and genomics data from the Genotype-Tissue Expression (GTEx) project, we establish and validate a quantitative imaging biomarker for immune cell infiltration. A total of 100,215 variants were selected for iQTL profiling and tested for genotype-phenotype associations with our quantitative imaging biomarker. Significant associations were found in HDAC9 and TXNDC5. We validated the TXNDC5 association using GTEx cis-expression QTL data and an independent hypothyroidism dataset from the Electronic Medical Records and Genomics network.
Collapse
Affiliation(s)
- Joseph D Barry
- Center for Cancer Computational Biology and Department of Biostatistics and Computational Biology, Dana-Farber Cancer Institute, 450 Brookline Avenue, Boston, MA 02215, USA; Department of Biostatistics, Harvard T.H. Chan School of Public Health, 667 Huntington Avenue, Boston, MA 02115, USA.
| | - Maud Fagny
- Center for Cancer Computational Biology and Department of Biostatistics and Computational Biology, Dana-Farber Cancer Institute, 450 Brookline Avenue, Boston, MA 02215, USA; Department of Biostatistics, Harvard T.H. Chan School of Public Health, 667 Huntington Avenue, Boston, MA 02115, USA
| | - Joseph N Paulson
- Center for Cancer Computational Biology and Department of Biostatistics and Computational Biology, Dana-Farber Cancer Institute, 450 Brookline Avenue, Boston, MA 02215, USA; Department of Biostatistics, Harvard T.H. Chan School of Public Health, 667 Huntington Avenue, Boston, MA 02115, USA
| | - Hugo J W L Aerts
- Department of Radiology, Dana-Farber Cancer Institute, 450 Brookline Avenue, Boston, MA 02215, USA
| | - John Platig
- Center for Cancer Computational Biology and Department of Biostatistics and Computational Biology, Dana-Farber Cancer Institute, 450 Brookline Avenue, Boston, MA 02215, USA; Department of Biostatistics, Harvard T.H. Chan School of Public Health, 667 Huntington Avenue, Boston, MA 02115, USA
| | - John Quackenbush
- Center for Cancer Computational Biology and Department of Biostatistics and Computational Biology, Dana-Farber Cancer Institute, 450 Brookline Avenue, Boston, MA 02215, USA; Department of Biostatistics, Harvard T.H. Chan School of Public Health, 667 Huntington Avenue, Boston, MA 02115, USA
| |
Collapse
|
48
|
Morrow JD, Glass K, Cho MH, Hersh CP, Pinto-Plata V, Celli B, Marchetti N, Criner G, Bueno R, Washko G, Choi AMK, Quackenbush J, Silverman EK, DeMeo DL. Human Lung DNA Methylation Quantitative Trait Loci Colocalize with Chronic Obstructive Pulmonary Disease Genome-Wide Association Loci. Am J Respir Crit Care Med 2018; 197:1275-1284. [PMID: 29313708 PMCID: PMC5955059 DOI: 10.1164/rccm.201707-1434oc] [Citation(s) in RCA: 47] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2017] [Accepted: 01/03/2018] [Indexed: 12/23/2022] Open
Abstract
RATIONALE As the third leading cause of death in the United States, the impact of chronic obstructive pulmonary disease (COPD) makes identification of its molecular mechanisms of great importance. Genome-wide association studies (GWASs) have identified multiple genomic regions associated with COPD. However, genetic variation only explains a small fraction of the susceptibility to COPD, and sub-genome-wide significant loci may play a role in pathogenesis. OBJECTIVES Regulatory annotation with epigenetic evidence may give priority for further investigation, particularly for GWAS associations in noncoding regions. We performed integrative genomics analyses using DNA methylation profiling and genome-wide SNP genotyping from lung tissue samples from 90 subjects with COPD and 36 control subjects. METHODS We performed methylation quantitative trait loci (mQTL) analyses, testing for SNPs associated with percent DNA methylation and assessed the colocalization of these results with previous COPD GWAS findings using Bayesian methods in the R package coloc to highlight potential regulatory features of the loci. MEASUREMENTS AND MAIN RESULTS We identified 942,068 unique SNPs and 33,996 unique CpG sites among the significant (5% false discovery rate) cis-mQTL results. The genome-wide significant and subthreshold (P < 10-4) GWAS SNPs were enriched in the significant mQTL SNPs (hypergeometric test P < 0.00001). We observed enrichment for sites located in CpG shores and shelves, but not CpG islands. Using Bayesian colocalization, we identified loci in regions near KCNK3, EEFSEC, PIK3CD, DCDC2C, TCERG1L, FRMD4B, and IL27. CONCLUSIONS Colocalization of mQTL and GWAS loci provides regulatory characterization of significant and subthreshold GWAS findings, supporting a role for genetic control of methylation in COPD pathogenesis.
Collapse
Affiliation(s)
| | | | - Michael H. Cho
- Channing Division of Network Medicine
- Division of Pulmonary and Critical Care Medicine, and
| | - Craig P. Hersh
- Channing Division of Network Medicine
- Division of Pulmonary and Critical Care Medicine, and
| | | | | | - Nathaniel Marchetti
- Division of Pulmonary and Critical Care Medicine, Temple University, Philadelphia, Pennsylvania
| | - Gerard Criner
- Division of Pulmonary and Critical Care Medicine, Temple University, Philadelphia, Pennsylvania
| | - Raphael Bueno
- Division of Thoracic Surgery, Brigham and Women’s Hospital, Boston, Massachusetts
| | - George Washko
- Division of Pulmonary and Critical Care Medicine, and
| | - Augustine M. K. Choi
- Department of Medicine, New York Presbyterian/Weill Cornell Medical Center, New York, New York; and
| | - John Quackenbush
- Department of Biostatistics and Computational Biology, Dana-Farber Cancer Institute, Boston, Massachusetts
| | - Edwin K. Silverman
- Channing Division of Network Medicine
- Division of Pulmonary and Critical Care Medicine, and
| | - Dawn L. DeMeo
- Channing Division of Network Medicine
- Division of Pulmonary and Critical Care Medicine, and
| |
Collapse
|
49
|
Abstract
SummaryTo review the current state of the art in computational methods for the analysis of DNA microarray data.The review considers methods of microarray data collection, transformation and representation, comparisons and predictions of gene expression from the data, their mechanistic analysis, related systems biology, and the application of clustering techniques.Functional genomics approaches have greatly increased the rate at which data on biological systems is generated, leading to corresponding challenges in analyzing the data through advanced computational techniques . The paper compares and contrasts the application of computational clustering for discovery, comparison, and prediction of gene expression classes, together with their evaluation and relation to mechanistic analyses of biological systems.Methods for assaying gene expression levels by DNA microarray experiments produce considerably more data than other techniques, and require a wide variety of computational techniques for identifying patterns of expression that may be biologically significant. These will have to be verified and validated by comparison to results from other methods, integrated with other systems data, and provide the feedback for further experimentation for testing mechanistic or other biological hypotheses.
Collapse
|
50
|
Domenyuk V, Gatalica Z, Santhanam R, Wei X, Stark A, Kennedy P, Toussaint B, Levenberg S, Wang R, Xiao N, Greil R, Rinnerthaler G, Gampenrieder S, Heimberger AB, Berry DJ, Barker A, Demetri GD, Quackenbush J, Marshall JL, Poste G, Vacirca JL, Vidal GA, Schwartzberg LS, Halbert DD, Voss A, Miglarese MR, Famulok M, Mayer G, Spetzler D. Abstract P2-09-09: Polyligand profiling differentiates cancer patients according to their benefit of treatment. Cancer Res 2018. [DOI: 10.1158/1538-7445.sabcs17-p2-09-09] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Abstract
Introduction: Deconvolution of multi-nodal perturbations in cancer network architecture demands highly multiplexed profiling assays. We demonstrate the value of polyligand profiling of tumor systems states using libraries of single stranded oligodeoxynucleotides (ssODN) to distinguish between tumor tissue from breast cancer patients who did or did not derive benefit from treatment regimens containing trastuzumab.
Methods: This study included cases from women with invasive breast cancer who received chemotherapy+ trastuzumab (C+T) or trastuzumab monotherapy with available retrospective data on the time to next treatment (TTNT). A library of 2x1012 unique ssODN was exposed to FFPE tissues from patients who benefited (B) or not (NB) from trastuzumab-based regimens in several rounds of positive and negative selection. Two enriched libraries were screened on independent set of 42 B and 19 NB cases using a modified IHC protocol for detection of bound ssODNs. Poly-Ligand Profiles (PLP) were scored by a blinded pathologist. Two libraries, EL-NB and EL-B, showed significant p-values between groups of responders and non-responders. A Cox-PH model was fitted using either tumors' HER2 status or PLP test results as the independent variable. Median survival time was calculated from the Kaplan-Meier estimate. A separate group of 63 cases with TTNT data from chemotherapy without trastuzumab was used as a control to distinguish prognostic from predictive performance.
Results: The PLP scores of EL-NB and EL-B were assessed by receiver operating characteristic (ROC) curves and resulted in a combined AUC value of 0.81. EL-NB and EL-B were able to effectively classify B and NB patients with either HER2-negative/equivocal (AUC = 0.73) or HER2-positive cancers (AUC = 0.84). In contrast, HER2 status alone yielded an AUC value of 0.47. The combined PLP scores for the independent set of 63 patients treated with C excluding trastuzumab resulted in an AUC value of 0.53, indicating that the assay was predictive and not simply prognostic. Kaplan-Meier curves analysis shows that PLP+ cases have 429 days median TTNT, while PLP- cases have 129 days (HR = 0.38, log-rank p = 0.001). Analysis based on HER2 status showed no significant difference in TTNT between patients that were HER2+ (280 days) or HER2-negative/equivocal (336 days, HR = 1.27, log-rank p =0.45).
Summary: Performance of the PLP assay in differentiating patients who did or did not benefit from trastuzumab therapy outperforms the standard IHC assay for HER2 status. These results represent a promising step towards the development of a CDx to identify the 50-70% of HER2+ patients who will not benefit from trastuzumab. In addition, PLP also has the potential to identify the HER2-negative/equivocal patients who may benefit from trastuzumab-containing regimens.
Citation Format: Domenyuk V, Gatalica Z, Santhanam R, Wei X, Stark A, Kennedy P, Toussaint B, Levenberg S, Wang R, Xiao N, Greil R, Rinnerthaler G, Gampenrieder S, Heimberger AB, Berry DJ, Barker A, Demetri GD, Quackenbush J, Marshall JL, Poste G, Vacirca JL, Vidal GA, Schwartzberg LS, Halbert DD, Voss A, Miglarese MR, Famulok M, Mayer G, Spetzler D. Polyligand profiling differentiates cancer patients according to their benefit of treatment [abstract]. In: Proceedings of the 2017 San Antonio Breast Cancer Symposium; 2017 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2018;78(4 Suppl):Abstract nr P2-09-09.
Collapse
Affiliation(s)
- V Domenyuk
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| | - Z Gatalica
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| | - R Santhanam
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| | - X Wei
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| | - A Stark
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| | - P Kennedy
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| | - B Toussaint
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| | - S Levenberg
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| | - R Wang
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| | - N Xiao
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| | - R Greil
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| | - G Rinnerthaler
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| | - S Gampenrieder
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| | - AB Heimberger
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| | - DJ Berry
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| | - A Barker
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| | - GD Demetri
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| | - J Quackenbush
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| | - JL Marshall
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| | - G Poste
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| | - JL Vacirca
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| | - GA Vidal
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| | - LS Schwartzberg
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| | - DD Halbert
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| | - A Voss
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| | - MR Miglarese
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| | - M Famulok
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| | - G Mayer
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| | - D Spetzler
- Caris Life Sciences, Phoenix, AZ; Paracelsus Medical University Salzburg, Austria and Salzburg Cancer Research Institute, and Cancer Cluster Salzburg, Salzburg, Austria; University of Texas MD Anderson Cancer Center, Houston, TX; Complex Adaptive Systems Initiative, Arizona State University, Scottsdale, AZ; Dana-Farber Cancer Institute and Ludwig Center at Harvard Medical School, Boston, MA; Dana-Farber Cancer Institute, Boston, Boston, MA; Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, Washington, DC; North Shore Hematology Oncology Associates Cancer Center, New York, NY; University of Tennessee Health Science Center, Memphis, TN; LIMES Program Unit Chemical Biology & Medicinal Chemistry, c/o Kekulé Institute for Organic Chemistry and Biochemistry, University of Bonn, Bonn, Germany; Chemical Biology Max-Planck-Fellowship Group, Center of Advanced European Studies and Research (CAESAR, Bonn, Germany; Center of Aptamer Research and Development, University of Bonn, Bonn, Germany
| |
Collapse
|