Ung, Choong Yong

Abstract Glioblastoma is the most common and lethal primary brain tumor in adults, with a median survival less than 15 months despite current standard of care. Oncolytic immunovirotherapies (OVs) have shown promise in early clinical trials. OVs are engineered to specifically target cancer cells, trigger autophagic programmed cell death, and release tumor antigens in the context of T-cell activating cytokines. Cellular communication network factor 1 (CCN1) is an extracellular matrix protein expressed in 70% of glioblastomas that has been shown to drastically reduce OV efficacy, particularly herpes simplex virus type 1 (HSV-1). Our objectives were to decode protein-protein interaction (PPI) networks activated by CCN1, identify critical nodes in resistance networks, and use this analysis to design OVs that overcome CCN1-mediated resistance. We used NetDecoder to elucidate phenotype-specific PPI subnetworks in LN229 human glioblastoma cells with tetracycline-inducible CCN1 expression. Publicly available microarray data on CCN1-induced and control samples established differentially expressed genes, which served as sources in our PPI network to derive prioritized context-specific subnetworks. We found that 11 source genes collaborate via 39 deep nodes to confer OV resistance to CCN1-expressing glioblastomas. Of these, a router (IKBKE) and a sink (YBX1) have been previously implicated in glioblastoma pathogenesis. We conclude that autophagy regulator HSP90AA1, a critical node, may be targeted to improve OV efficacy in CCN1-expressing glioblastomas. Furthermore, WebGestalt overrepresentation analysis of subnetwork nodes suggests that oncolytic adenovirus may be more effective than HSV-1 as a synthetic biology chassis.

Abstract Background: Large scale genetic perturbations screens have been instrumental for many key discoveries. In the last decade sequence specific genetic perturbations created RNAi and CRISPR (PMID: 25446316). Recently, specific genome targeting approaches using single guide RNAs (sgRNAs) for the CRISPR-Cas9 system have been repurposed to induce sequence-specific repression or activation of gene expression at a genome-wide scale. Different perturbation approaches have different effects and encompass a trade-off between specificity and duration in establishing a measurable perturbation which can result in different biological readouts. However, bench-side strategies are time consuming and expensive. The method described here aims to build a systems-based computational platform to assess pharmacological and biological perturbations at genome-wide network level.Methods: Our network tool NetDecoder (Nucleic Acids Res. 2016; 44(10):e100) that models genome-wide information flow will be used to build context specific protein-protein interaction networks. These networks will be subsequently subjected to in silico perturbations and their impacts on network properties will be assessed. AVANA and Demeter datasets will be used during method validation procedures. Results: Using publicly available expression datasets we generated context specific Breast Cancer networks and individually manipulated nodes and edges of context-specific networks associated with pharmacological FDA approved targets. Our perturbation networks captured cellular rewiring by engulfing new gene players, alterations of the network size and its connectivity. Conclusions: Our in-silico genome wide perturbation platform provides a flexible approach to fine-tune single and tandem manipulations, allowing for inexpensively monitoring a range of pharmacological interventions using large omics data. Citation Format: Cristina Correia, Choong Yong Ung, Cheng Zhang, Scott H. Kaufmann, Hu Li. A genome wide pharmacological perturbation platform [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 1226.

Abstract Current understandings of individual disease etiology and therapeutics are limited despite great need. To fill the gap, we propose a novel computational pipeline which collects potent disease gene cooperative pathways to envision individualized disease etiology and therapies. Our algorithm constructs individualized disease modules de novo which enable us to elucidate the importance of mutated genes in specific patients and to understand the synthetic penetrance of these genes across patients. We reveal that importance of notorious cancer drivers TP53 and PIK3CA fluctuate widely across breast cancers and peak in tumors with distinct numbers of mutations, and that rarely mutated genes such as XPO1 and PLEKHA1 have high disease module importance in specific individuals. Furthermore, individualized module disruption enables us to devise customized singular and combinatorial target therapies which were highly varied across patients demonstrating the need for precision therapeutics pipelines. As the first analysis of de novo individualized disease modules, we illustrate the power of individualized disease modules for precision medicine by providing deep novel insights on the activity of diseased genes in individuals. Citation Format: Taylor Weiskittel, Cristina Correia, Choong Ung, Cheng Zhang, Hu Li. De novo individualized disease modules reveal the synthetic penetrance of genes and inform personalized treatment regimens [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 490.

Together, single-cell technologies and systems biology have been used to investigate previously unanswerable questions in biomedicine with unparalleled detail. Despite these advances, gaps in analytical capacity remain. Machine learning, which has revolutionized biomedical imaging analysis, drug discovery, and systems biology, is an ideal strategy to fill these gaps in single-cell studies. Machine learning additionally has proven to be remarkably synergistic with single-cell data because it remedies unique challenges while capitalizing on the positive aspects of single-cell data. In this review, we describe how systems-biology algorithms have layered machine learning with biological components to provide systems level analyses of single-cell omics data, thus elucidating complex biological mechanisms. Accordingly, we highlight the trifecta of single-cell, systems-biology, and machine-learning approaches and illustrate how this trifecta can significantly contribute to five key areas of scientific research: cell trajectory and identity, individualized medicine, pharmacology, spatial omics, and multi-omics. Given its success to date, the systems-biology, single-cell omics, and machine-learning trifecta has proven to be a potent combination that will further advance biomedical research.

Abstract Glioblastoma multiforme (GBM) is a particularly aggressive primary brain tumor that is relatively common in adults. The disease comes with a dismal prognosis, often with expected survival just over a year when treated with surgery, chemotherapy, and radiation. A promising addition to our GBM treatment arsenal are oncolytic viruses (OVs), a class of immunotherapy now in clinical trials. The most effective OVs home in on tumor cells, killing them and releasing tumor-associated antigens that elicit antitumor immunity. This effectiveness, however, varies depending on the cell state—a property that may be influenced by the composition of the extracellular matrix (ECM) of GBM. Cellular communication network factor 1 (CCN1) is found in the ECM of the majority of GBMs and is predictive of resistance to OVs, particularly those derived from herpes simplex virus type 1 (HSV-1). In this study, our aim is to understand how extracellular CCN1 alters the GBM intracellular state in a manner that confers resistance to OVs. We began by using our NetDecoder platform to build CCN1-specific prioritized protein-protein interaction (PPI) networks from published LN229 human GBM transcriptomic datasets. From this network, we identified 12 PPI edges that dominate the network flow in CCN1-high LN229 cells when compared to CCN1-low controls. The 21 unique nodes that comprise these binary PPI edges were extracted and subjected to downstream biological process overrepresentation analyses using Gene Ontology. We found that the cytokine-mediated signaling pathway (FDR = 4.2 x 10-11; enrichment ratio = 15.8) and the cellular response to cytokine stimulus (FDR = 3.1 x 10-9; enrichment ratio = 10.9) were highly enriched when mapped to these nodes. We conclude that CCN1, acting as a paracrine signal in the GBM tumor microenvironment, primes the cell to inhibit HSV-1 and other candidate OV replication prior to any virus interactions, thus contributing to immunovirotherapy resistance. While CCN1 and its cell surface integrin binding partners are ostensible drug targets for overcoming this resistance, pharmacologic interventions that target CCN1-specific high flow edges in our prioritized PPI network may be more effective. Citation Format: Dileep D. Monie, Cristina Correia, Cheng Zhang, Choong Ung, Richard G. Vile, Hu Li. Innate immunity and cytokine signaling pathways drive resistance to HSV-1 oncolytic immunovirotherapy in glioblastoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 1679.

Abstract One of the greatest barriers to curative treatment of neuroblastoma is its frequent metastatic outgrowth prior to diagnosis, especially in cases driven by amplification of the MYCN oncogene. However, only a limited number of regulatory proteins that contribute to this complex MYCN-mediated process have been elucidated. Here we show that the growth arrest-specific 7 (GAS7) gene, located at chromosome band 17p13.1, is preferentially deleted in high-risk MYCN-driven neuroblastoma. GAS7 expression was also suppressed in MYCN-amplified neuroblastoma lacking 17p deletion. GAS7 deficiency led to accelerated metastasis in both zebrafish and mammalian models of neuroblastoma with overexpression or amplification of MYCN. Analysis of expression profiles and the ultrastructure of zebrafish neuroblastoma tumors with MYCN overexpression identified that GAS7 deficiency led to (i) downregulation of genes involved in cell–cell interaction, (ii) loss of contact among tumor cells as critical determinants of accelerated metastasis, and (iii) increased levels of MYCN protein. These results provide the first genetic evidence that GAS7 depletion is a critical early step in the cascade of events culminating in neuroblastoma metastasis in the context of MYCN overexpression. Significance: Heterozygous deletion or MYCN-mediated repression of GAS7 in neuroblastoma releases an important brake on tumor cell dispersion and migration to distant sites, providing a novel mechanism underlying tumor metastasis in MYCN-driven neuroblastoma. See related commentary by Menard, p. 2815

Abstract With the emergence of genome editing technologies and synthetic biology, it is now possible to engineer genetic circuits driving a cell's phenotypic response to a stressor. However, capturing a continuous response, rather than simply a binary ‘on’ or ‘off’ response, remains a bioengineering challenge. No tools currently exist to identify gene candidates responsible for predetermining and fine-tuning cell response phenotypes. To address this gap, we devised a novel Regulostat Inferelator (RSI) algorithm to decipher intrinsic molecular devices or networks that predetermine cellular phenotypic responses. The RSI algorithm is designed to extract gene expression patterns from basal transcriptomic data in order to identify ‘regulostat’ constituent gene pairs, which exhibit rheostat-like mode-of-cooperation capable of fine-tuning cellular response. Our proof-of-concept study provides computational evidence for the existence of regulostats and that these networks predetermine cellular response prior to exposure to a stressor or drug. In addition, our work, for the first time, provides evidence of context-specific, drug–regulostat interactions in predetermining drug response phenotypes in cancer cells. Given RSI-inferred regulostat networks offer insights for prioritizing gene candidates capable of rendering a resistant phenotype sensitive to a given drug, we envision that this tool will be of great value in bioengineering and medicine.

Pediatric cancers rarely exhibit recurrent mutational events when compared to most adult cancers. This poses a challenge in understanding how cancers initiate, progress, and metastasize in early childhood. Also, due to limited detected driver mutations, it is difficult to benchmark key genes for drug development. In this review, we use neuroblastoma, a pediatric solid tumor of neural crest origin, as a paradigm for exploring “big data” applications in pediatric oncology. Computational strategies derived from big data science–network- and machine learning-based modeling and drug repositioning—hold the promise of shedding new light on the molecular mechanisms driving neuroblastoma pathogenesis and identifying potential therapeutics to combat this devastating disease. These strategies integrate robust data input, from genomic and transcriptomic studies, clinical data, and in vivo and in vitro experimental models specific to neuroblastoma and other types of cancers that closely mimic its biological characteristics. We discuss contexts in which “big data” and computational approaches, especially network-based modeling, may advance neuroblastoma research, describe currently available data and resources, and propose future models of strategic data collection and analyses for neuroblastoma and other related diseases.