Integrated network analyses for functional genomic studies in cancer

The online Tabloid Proteome: an annotated database of protein associations

A complete knowledge of the proteome can only be attained by determining the associations between proteins, along with the nature of these associations (e.g. physical contact in protein–protein interactions, participation in complex formation or different roles in the same pathway). Despite extensive efforts in elucidating direct protein interactions, our knowledge on the complete spectrum of protein associations remains limited. We therefore developed a new approach that detects protein associations from identifications obtained after re-processing of large-scale, public mass spectrometry-based proteomics data. Our approach infers protein association based on the co-occurrence of proteins across many different proteomics experiments, and provides information that is almost completely complementary to traditional direct protein interaction studies. We here present a web interface to query and explore the associations derived from this method, called the online Tabloid Proteome. The online Tabloid Proteome also integrates biological knowledge from several existing resources to annotate our derived protein associations. The online Tabloid Proteome is freely available through a user-friendly web interface, which provides intuitive navigation and data exploration options for the user at http://iomics.ugent.be/tabloidproteome.

Functional Genomics Approach Identifies Novel Signaling Regulators of TGFα Ectodomain Shedding

Ectodomain shedding of cell-surface precursor proteins by metalloproteases generates important cellular signaling molecules. Of importance for disease is the release of ligands that activate the EGFR, such as TGFα, which is mostly carried out by ADAM17 [a member of the A-disintegrin and metalloprotease (ADAM) domain family]. EGFR ligand shedding has been linked to many diseases, in particular cancer development, growth and metastasis, as well as resistance to cancer therapeutics. Excessive EGFR ligand release can outcompete therapeutic EGFR inhibition or the inhibition of other growth factor pathways by providing bypass signaling via EGFR activation. Drugging metalloproteases directly have failed clinically because it indiscriminately affected shedding of numerous substrates. It is therefore essential to identify regulators for EGFR ligand cleavage. Here, integration of a functional shRNA genomic screen, computational network analysis, and dedicated validation tests succeeded in identifying several key signaling pathways as novel regulators of TGFα shedding in cancer cells. Most notably, a cluster of genes with NFκB pathway regulatory functions was found to strongly influence TGFα release, albeit independent of their NFκB regulatory functions. Inflammatory regulators thus also govern cancer cell growth-promoting ectodomain cleavage, lending mechanistic understanding to the well-known connection between inflammation and cancer.Implications: Using genomic screens and network analysis, this study defines targets that regulate ectodomain shedding and suggests new treatment opportunities for EGFR-driven cancers. Mol Cancer Res; 16(1); 147-61. ©2017 AACR.

Next generation of network medicine: interdisciplinary signaling approaches

In the last decade, network approaches have transformed our understanding of biological systems. Network analyses and visualizations have allowed us to identify essential molecules and modules in biological systems, and improved our understanding of how changes in cellular processes can lead to complex diseases, such as cancer, infectious and neurodegenerative diseases. "Network medicine" involves unbiased large-scale network-based analyses of diverse data describing interactions between genes, diseases, phenotypes, drug targets, drug transport, drug side-effects, disease trajectories and more. In terms of drug discovery, network medicine exploits our understanding of the network connectivity and signaling system dynamics to help identify optimal, often novel, drug targets. Contrary to initial expectations, however, network approaches have not yet delivered a revolution in molecular medicine. In this review, we propose that a key reason for the limited impact, so far, of network medicine is a lack of quantitative multi-disciplinary studies involving scientists from different backgrounds. To support this argument, we present existing approaches from structural biology, 'omics' technologies (e.g., genomics, proteomics, lipidomics) and computational modeling that point towards how multi-disciplinary efforts allow for important new insights. We also highlight some breakthrough studies as examples of the potential of these approaches, and suggest ways to make greater use of the power of interdisciplinarity. This review reflects discussions held at an interdisciplinary signaling workshop which facilitated knowledge exchange from experts from several different fields, including in silico modelers, computational biologists, biochemists, geneticists, molecular and cell biologists as well as cancer biologists and pharmacologists.

Multitype Network-Guided Target Controllability in Phenotypically Characterized Osteosarcoma: Role of Tumor Microenvironment

This study highlights the relevance of network-guided controllability analysis as a precision oncology tool. Target controllability through networks is potentially relevant to cancer research for the identification of therapeutic targets. With reference to a recent study on multiple phenotypes from 22 osteosarcoma (OS) cell lines characterized both in vitro and in vivo, we found that a variety of critical proteins in OS regulation circuits were in part phenotype specific and in part shared. To generalize our inference approach and match cancer phenotypic heterogeneity, we employed multitype networks and identified targets in correspondence with protein sub-complexes. Therefore, we established the relevance for diagnostic and therapeutic purposes of inspecting interactive targets, namely those enriched by significant connectivity patterns in protein sub-complexes. Emerging targets appeared with reference to the OS microenvironment, and relatively to small leucine-rich proteoglycan members and D-type cyclins, among other collagen, laminin, and keratin proteins. These described were evidences shared across all phenotypes; instead, specific evidences were provided by critical proteins including IGFBP7 and PDGFRA in the invasive phenotype, and FGFR3 and THBS1 in the colony forming phenotype.

Time-course gene profiling and networks in demethylated retinoblastoma cell line

Retinoblastoma, a very aggressive cancer of the developing retina, initiatiates by the biallelic loss of RB1 gene, and progresses very quickly following RB1 inactivation. While its genome is stable, multiple pathways are deregulated, also epigenetically. After reviewing the main findings in relation with recently validated markers, we propose an integrative bioinformatics approach to include in the previous group new markers obtained from the analysis of a single cell line subject to epigenetic treatment. In particular, differentially expressed genes are identified from time course microarray experiments on the WERI-RB1 cell line treated with 5-Aza-2′-deoxycytidine (decitabine; DAC). By inducing demethylation of CpG island in promoter genes that are involved in biological processes, for instance apoptosis, we performed the following main integrative analysis steps: i) Gene expression profiling at 48h, 72h and 96h after DAC treatment; ii) Time differential gene co-expression networks and iii) Context-driven marker association (transcriptional factor regulated protein networks, master regulatory paths). The observed DAC-driven temporal profiles and regulatory connectivity patterns are obtained by the application of computational tools, with support from curated literature. It is worth emphasizing the capacity of networks to reconcile multi-type evidences, thus generating testable hypotheses made available by systems scale predictive inference power. Despite our small experimental setting, we propose through such integrations valuable impacts of epigenetic treatment in terms of gene expression measurements, and then validate evidenced apoptotic effects.

Pathway-based network modeling finds hidden genes in shRNA screen for regulators of acute lymphoblastic leukemia

Data integration stands to improve interpretation of RNAi screens which, as a result of off-target effects, typically yield numerous gene hits of which only a few validate. These off-target effects can result from seed matches to unintended gene targets (reagent-based) or cellular pathways, which can compensate for gene perturbations (biology-based). We focus on the biology-based effects and use network modeling tools to discover pathways de novo around RNAi hits. By looking at hits in a functional context, we can uncover novel biology not identified from any individual 'omics measurement. We leverage multiple 'omic measurements using the Simultaneous Analysis of Multiple Networks (SAMNet) computational framework to model a genome scale shRNA screen investigating Acute Lymphoblastic Leukemia (ALL) progression in vivo. Our network model is enriched for cellular processes associated with hematopoietic differentiation and homeostasis even though none of the individual 'omic sets showed this enrichment. The model identifies genes associated with the TGF-beta pathway and predicts a role in ALL progression for many genes without this functional annotation. We further experimentally validate the hidden genes - Wwp1, a ubiquitin ligase, and Hgs, a multi-vesicular body associated protein - for their role in ALL progression. Our ALL pathway model includes genes with roles in multiple types of leukemia and roles in hematological development. We identify a tumor suppressor role for Wwp1 in ALL progression. This work demonstrates that network integration approaches can compensate for off-target effects, and that these methods can uncover novel biology retroactively on existing screening data. We anticipate that this framework will be valuable to multiple functional genomic technologies - siRNA, shRNA, and CRISPR - generally, and will improve the utility of functional genomic studies.

Identification of potential glucocorticoid receptor therapeutic targets in multiple myeloma

Glucocorticoids (GC) are a cornerstone of combination therapies for multiple myeloma. However, patients ultimately develop resistance to GCs frequently based on decreased glucocorticoid receptor (GR) expression. An understanding of the direct targets of GC actions, which induce cell death, is expected to culminate in potential therapeutic strategies for inducing cell death by regulating downstream targets in the absence of a functional GR. The specific goal of our research is to identify primary GR targets that contribute to GC-induced cell death, with the ultimate goal of developing novel therapeutics around these targets that can be used to overcome resistance to GCs in the absence of GR. Using the MM.1S glucocorticoid-sensitive human myeloma cell line, we began with the broad platform of gene expression profiling to identify glucocorticoid-regulated genes further refined by combination treatment with phosphatidylinositol-3’-kinase inhibition (PI3Ki). To further refine the search to distinguish direct and indirect targets of GR that respond to the combination GC and PI3Ki treatment of MM.1S cells, we integrated 1) gene expression profiles of combination GC treatment with PI3Ki, which induces synergistic cell death; 2) negative correlation between genes inhibited by combination treatment in MM.1S cells and genes over-expressed in myeloma patients to establish clinical relevance and 3) GR chromatin immunoprecipitation with massively parallel sequencing (ChIP-Seq) in myeloma cells to identify global chromatin binding for the glucocorticoid receptor (GR). Using established bioinformatics platforms, we have integrated these data sets to identify a subset of candidate genes that may form the basis for a comprehensive picture of glucocorticoid actions in multiple myeloma. As a proof of principle, we have verified two targets, namely RRM2 and BCL2L1, as primary functional targets of GR involved in GC-induced cell death.

Systematic evaluation of underlying defects in DNA repair as an approach to case-only assessment of familial prostate cancer

Risk assessment for prostate cancer is challenging due to its genetic heterogeneity. In this study, our goal was to develop an operational framework to select and evaluate gene variants that may contribute to familial prostate cancer risk. Drawing on orthogonal sources, we developed a candidate list of genes relevant to prostate cancer, then analyzed germline exomes from 12 case-only prostate cancer patients from high-risk families to identify patterns of protein-damaging gene variants. We described an average of 5 potentially disruptive variants in each individual and annotated them in the context of public databases representing human variation. Novel damaging variants were found in several genes of relevance to prostate cancer. Almost all patients had variants associated with defects in DNA damage response. Many also had variants linked to androgen signaling. Treatment of primary T-lymphocytes from these prostate cancer patients versus controls with DNA damaging agents showed elevated levels of the DNA double strand break (DSB) marker γH2AX (p < 0.05), supporting the idea of an underlying defect in DNA repair. This work suggests the value of focusing on underlying defects in DNA damage in familial prostate cancer risk assessment and demonstrates an operational framework for exome sequencing in case-only prostate cancer genetic evaluation.

An argument for mechanism-based statistical inference in cancer

Cancer is perhaps the prototypical systems disease, and as such has been the focus of extensive study in quantitative systems biology. However, translating these programs into personalized clinical care remains elusive and incomplete. In this perspective, we argue that realizing this agenda—in particular, predicting disease phenotypes, progression and treatment response for individuals—requires going well beyond standard computational and bioinformatics tools and algorithms. It entails designing global mathematical models over network-scale configurations of genomic states and molecular concentrations, and learning the model parameters from limited available samples of high-dimensional and integrative omics data. As such, any plausible design should accommodate: biological mechanism, necessary for both feasible learning and interpretable decision making; stochasticity, to deal with uncertainty and observed variation at many scales; and a capacity for statistical inference at the patient level. This program, which requires a close, sustained collaboration between mathematicians and biologists, is illustrated in several contexts, including learning biomarkers, metabolism, cell signaling, network inference and tumorigenesis.

In vivo shRNA screens in solid tumors

Loss-of-function (LOF) experiments targeting multiple genes during tumorigenesis can be implemented using pooled shRNA libraries. RNAi screens in animal models rely on the use of multiple shRNAs to simultaneously disrupt gene function, as well as to serve as barcodes for cell fate outcomes during tumorigenesis. Here we provide a protocol for performing RNAi screens in orthotopic mouse tumor models, referring to glioma and lung adenocarcinoma as specific examples. The protocol aims to provide guidelines for applying RNAi to a diverse spectrum of solid tumors and to highlight crucial considerations when designing and performing these studies. It covers shRNA library assembly and packaging into lentiviral particles, and transduction into tumor-initiating cells (TICs), followed by in vivo transplantation, tumor DNA recovery, sequencing and analysis. Depending on the target genes and tumor model, tumor suppressors and oncogenes can be identified or biological pathways can be dissected in 6-9 weeks.

Approaches to in vitro tissue regeneration with application for human disease modeling and drug development

Reliable in vitro human disease models that capture the complexity of in vivo tissue behaviors are crucial to gain mechanistic insights into human disease and enable the development of treatments that are effective across broad patient populations. The integration of stem cell technologies, tissue engineering, emerging biomaterials strategies and microfabrication processes, as well as computational and systems biology approaches, is enabling new tools to generate reliable in vitro systems to study the molecular basis of human disease and facilitate drug development. In this review, we discuss these recently developed tools and emphasize opportunities and challenges involved in combining these technologies toward regenerative science.

Tumor suppressing properties of rodent parvovirus NS1 proteins and their derivatives

Cancer chemotherapy with monospecific agents is often hampered by the rapid development of tumor resistance to the drug used. Therefore, combination treatments aiming at several different targets are sought. Viral regulatory proteins, modified or not, appear ideal for this purpose because of their multimodal killing action against neoplastically transformed cells. The large nonstructural protein NS1 of rodent parvoviruses is an excellent candidate for an anticancer agent, shown to interfere specifically with cancer cell growth and survival. The present review describes the structure, functions, and regulation of the multifunctional protein NS1, its specific interference with cell processes and cell protein activities, and what is known so far about the mechanisms underlying NS1 interference with cancer growth. It further outlines prospects for the development of new, multimodal cancer toxins and their potential applications.