Search

Supplementary Table S1

Supplementary Figures and Supplementary Table captions

Purpose: Neuroblastoma is a pediatric malignancy of the sympathetic nervous system with diverse clinical behaviors. Genomic amplification of MYCN oncogene has been shown to drive neuroblastoma pathogenesis and correlate with aggressive disease, but the survival rates for those high-risk tumors carrying no MYCN amplification remain equally dismal. The paucity of mutations and molecular heterogeneity has hindered the development of targeted therapies for most advanced neuroblastomas. We use an alternative method to identify potential drugs that target nononcogene dependencies in high-risk neuroblastoma. Experimental Design: By using a gene expression–based integrative approach, we identified prognostic signatures and potentially effective single agents and drug combinations for high-risk neuroblastoma. Results: Among these predictions, we validated in vitro efficacies of some investigational and marketed drugs, of which niclosamide, an anthelmintic drug approved by the FDA, was further investigated in vivo. We also quantified the proteomic changes during niclosamide treatment to pinpoint nucleoside diphosphate kinase 3 (NME3) downregulation as a potential mechanism for its antitumor activity. Conclusions: Our results establish a gene expression–based strategy to interrogate cancer biology and inform drug discovery and repositioning for high-risk neuroblastoma.

Summary Cancer is a complex disease that relies on both oncogenic mutations and non-mutated genes for survival, and therefore coined as oncogene and non-oncogene addictions. The need for more effective combination therapies to overcome drug resistance in oncology has been increasingly recognized, but the identification of potentially synergistic drugs at scale remains challenging. Here we propose a gene-expression-based approach, which uses the recurrent perturbation-transcript regulatory relationships inferred from a large compendium of chemical and genetic perturbation experiments across multiple cell lines, to engender a testable hypothesis for combination therapies. These transcript-level recurrences were distinct from known compound-protein target counterparts, were reproducible in external datasets, and correlated with small-molecule sensitivity. We applied these recurrent relationships to predict synergistic drug pairs for cancer and experimentally confirmed two unexpected drug combinations in vitro. Our results corroborate a gene-expression-based strategy for combinatorial drug screening as a way to target non-mutated genes in complex diseases., Graphical Abstract, Highlights • Compound signatures targeting non-oncogene addiction for combinatorial drug discovery • These signatures are reproducible and linked to cancer hallmarks and drug sensitivity • Two synergistic drug combinations are experimentally confirmed in vitro, Bioinformatics; Cancer Systems Biology; Pharmacoinformatics

Neuroblastoma is a neural crest-derived embryonal tumor and accounts for about 15% of all cancer deaths in children. MYCN amplification is associated with aggressive and advanced stage of high-risk neuroblastoma, which remains difficult to treat and exhibits poor survival under current multimodality treatment. Here, we analyzed the transcriptomic profiles of neuroblastoma patients and showed that aurora kinases lead to poor survival and had positive correlation with MYCN amplification and high-risk disease. Further, pan-aurora kinase inhibitor (tozasertib) treatment not only induces cell-cycle arrest and suppresses cell proliferation, migration, and invasion ability in MYCN-amplified (MNA) neuroblastoma cell lines, but also inhibits tumor growth and prolongs animal survival in Th-MYCN transgenic mice. Moreover, we performed quantitative proteomics and identified 150 differentially expressed proteins after tozasertib treatment in the Th-MYCN mouse model. The functional and network-based enrichment revealed that tozasertib alters metabolic processes and identified a mitochondrial flavoenzyme in fatty acid β-oxidation, ACADM, which is correlated with aurora kinases and neuroblastoma patient survival. Our findings indicate that the aurora kinase inhibitor could cause metabolic imbalance, possibly by disturbing carbohydrate and fatty acid metabolic pathways, and ACADM may be a potential target in MNA neuroblastoma.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a causative agent for the outbreak of coronavirus disease 2019 (COVID-19). This global pandemic is now calling for efforts to develop more effective COVID-19 therapies. Here we use a host-directed approach, which focuses on cellular responses to diverse small-molecule treatments, to identify potentially effective drugs for COVID-19. This framework looks at the ability of compounds to elicit a similar transcriptional response to IFN-β, a type I interferon that fails to be induced at notable levels in response to SARS-CoV-2 infection. By correlating the perturbation profiles of ~3,000 small molecules with a high-quality signature of IFN-β-responsive genes in primary normal human bronchial epithelial cells, our analysis revealed four candidate COVID-19 compounds, namely homoharringtonine, narciclasine, anisomycin, and emetine. We experimentally confirmed that the predicted compounds significantly inhibited SARS-CoV-2 replication in Vero E6 cells at nanomolar, relatively non-toxic concentrations, with half-maximal inhibitory concentrations of 165.7 nM, 16.5 nM, and 31.4 nM for homoharringtonine, narciclasine, and anisomycin, respectively. Together, our results corroborate a host-centric strategy to inform protective antiviral therapies for COVID-19., Systems biology; COVID-19; Host-directed therapy; Type I interferon; Drug repurposing; Bioinformatics; Infectious disease; Transcriptomics; Pharmaceutical science

The interaction of long noncoding RNAs (lncRNAs) with one or more RNA-binding proteins (RBPs) is important to a plethora of cellular and physiological processes. The lncRNA SNHG1 was reported to be aberrantly expressed and associated with poor patient prognosis in several cancers including neuroblastoma. However, the interacting RBPs and biological functions associated with SNHG1 in neuroblastoma remain unknown. In this study, we identified 283, 31, and 164 SNHG1-interacting proteins in SK-N-BE(2)C, SK-N-DZ, and SK-N-AS neuroblastoma cells, respectively, using a RNA-protein pull-down assay coupled with liquid chromatography-tandem mass spectrometry (LC-MS/MS). Twenty-four SNHG1-interacting RBPs were identified in common from these three neuroblastoma cell lines. RBPs MATR3, YBX1, and HNRNPL have the binding sites for SNHG1 predicted by DeepBind motif analysis. Furthermore, the direct binding of MATR3 with SNHG1 was validated by Western blot and confirmed by RNA immunoprecipitation assay (RIP). Coexpression analysis revealed that the expression of SNHG1 is positively correlated with MATR3 ( P = 3.402 × 10-13). The high expression of MATR3 is associated with poor event-free survival ( P = 0.00711) and overall survival ( P = 0.00064). Biological functions such as ribonucleoprotein complex biogenesis, RNA processing, and RNA splicing are significantly enriched and in common between SNHG1 and MATR3. In conclusion, we identified MATR3 as binding to SNHG1 and the interaction might be involved in splicing events that enhance neuroblastoma progression.

Cancer is a complex disease that relies on both oncogenic mutations and non-mutated genes for survival, thereafter coined as oncogene and non-oncogene addictions. The need for more effective combination therapies to overcome drug resistance in oncology has been increasingly recognized, but the identification of potentially synergistic drugs at scale remains challenging. Here we propose a geneexpression- based approach, which uses the recurrent perturbation–transcript regulatory relationships inferred from a large compendium of chemical and genetic perturbation experiments across multiple cell lines, to engender a testable hypothesis for combination therapies. These transcript-level recurrences were distinct from known compound–protein target counterparts, reproducible in external datasets, and correlated with small-molecule sensitivity. Using this new approach, we predicted synergistic drug pairs for cancer, including most frequently the combinations of a topoisomerase inhibitor and an mTOR or PI3K inhibitor. We also experimentally confirmed the synergistic effects of one combination of CD-437 (a retinoid) and sirolimus (an mTOR inhibitor) and the other combination of narciclasine (a protein synthesis inhibitor) and purvalanol A (a CDK inhibitor) in two breast and two lung cancer cell lines. Our results corroborate a gene-expression-based strategy for combinatorial drug screening as a way to target non-mutated genes in complex diseases.

Terminators, which signal the end of transcription processes, are typically placed behind the last coding sequence of an operon to prevent interference between transcript units in most biologically synthetic systems. Here, we seek to extend the usability of terminators in genetic system design by using terminators as regulatory genetic parts. Terminators with different impacts on their upstream and downstream genes are characterized in detail via dynamic modeling to predict the behavior of the overall genetic system. Some nonlinear effects of terminators observed in our terminator measurements potentially facilitate regulation of gene expression. Through dynamic modeling in silico, we find that such genetic systems may behave like genetic filters. In agreement with the simulations, we successfully implement genetic high-pass and bandpass filters in vivo, demonstrating the potential of using terminators as regulatory parts. The genetic bandpass filter in this work is implemented through the interdependence between genetic parts, in which the termination efficiency of a terminator varies with the strength of the upstream promoter. This design strategy for a bandpass filter requires fewer base pairs than the conventional strategy of concatenating high-pass and low-pass filters. Our results show that this novel utilization of terminators as regulatory parts may provide a new perspective for efficient design of genetic circuits. We believe that further exploration of the complicated dynamics of terminators is important in the development of synthetic biology.

Metastasis of melanoma cells during the recurrence or the late stage of melanoma has been characterized as the dissemination of tumor cells under anchorage independency. The secreted interleukin-8 (IL-8) and its conical receptors from melanoma cells have been associated with melanoma malignancy. However, their correlations with melanoma cells under anchorage independency were unclear. Suspension of adherent melanoma cells generated the suspended melanoma cell model of anoikis resistance. The in-vivo xenograft experiment, in-vitro cell proliferation/migration assay, microarray, and bioinformatics analysis were used to compare the malignancy and gene expression profiling in adherent and suspended melanoma cells. PCR, enzyme-linked immunosorbent assay, immunohistochemistry, and kinase inhibition assay were adapted to validate the expression and regulation of IL-8 and CXCR1/2. Suspended melanoma cells were anoikis resistant and showed elevated malignancy in vivo and in vitro. Gene expression profiling of adherent and suspended melanoma cells showed extensive alteration associated with cell survival/death, cell signaling, and regulation of gene expression. Microarray and bioinformatics analysis on gene set enrichment analysis further showed elevated IL-8 expression in suspended melanoma cells. The upregulation of IL-8 and the effect on chemotaxis were mediated by MEK/ERK activation upon cell suspension. Change in JNK phosphorylation induced CXCR1 downregulation under cell suspension, but upregulation by cell reattachment. We suggest the possible roles of elevated IL-8 secretion and change in CXCR expression contributing toward elevated melanoma malignancy upon reattachment from cell suspension. We show that the suspension of melanoma cells is critical in promoting melanoma malignancy in vivo and in vitro.

Anchorage-independent survival is one of the key features for malignant tumor cells. Whether specific gene alterations contributed by anchorage independency would further affect metastatic phenotypes of melanoma cells was unclear. We adapted suspension culture of melanoma cells to establish anchorage independency. The suspended melanoma cells lost their invasive abilities in vitro. Specific loss of laminin-binding ability in suspended melanoma cells was observed, which was correlated with downregulation of syndecan-1 as revealed by microarray and validated by qPCR and Western blot. Modulation of syndecan-1 expression level affected laminin binding, transwell migration and matrix metalloproteinase-2 secretion in melanoma cells. SDC1 expression and transwell migration were correlated with activity or level of protein kinase Cδ as evidence by specific inhibitors and shRNA transfection. In this study, we compared metastatic phenotypes and gene expressions of adherent and suspended melanoma cells. The anchorage independency led to protein kinase Cδ-mediated syndecan-1 downregulation, which contributed to loss of laminin-binding ability, reduced metalloproteinase-2 secretion and loss of invasiveness.