A Transcriptome-Based Precision Oncology Platform for Patient-Therapy Alignment in a Diverse Set of Treatment-Resistant Malignancies

Network-based approach elucidates critical genes in BRCA subtypes and chemotherapy response in triple negative breast cancer

Breast cancers (BRCA) exhibit substantial transcriptional heterogeneity, posing a significant clinical challenge. The global transcriptional changes in a disease context, however, are likely mediated by few key genes which reflect disease etiology better than the differentially expressed genes (DEGs). We apply our network-based tool PathExt to 1,059 BRCA tumors across 4 subtypes to identify key mediator genes in each subtype. Compared to conventional differential expression analysis, PathExt-identified genes exhibit greater concordance across tumors, revealing shared and subtype-specific biological processes; better recapitulate BRCA-associated genes in multiple benchmarks, and are more essential in BRCA subtype-specific cell lines. Single-cell transcriptomic analysis reveals a subtype-specific distribution of PathExt-identified genes in multiple cell types from the tumor microenvironment. Application of PathExt to a TNBC chemotherapy response dataset identified subtype-specific key genes and biological processes associated with resistance. We described putative drugs that target key genes potentially mediating drug resistance.

Network-based elucidation of colon cancer drug resistance mechanisms by phosphoproteomic time-series analysis

Aberrant signaling pathway activity is a hallmark of tumorigenesis and progression, which has guided targeted inhibitor design for over 30 years. Yet, adaptive resistance mechanisms, induced by rapid, context-specific signaling network rewiring, continue to challenge therapeutic efficacy. Leveraging progress in proteomic technologies and network-based methodologies, we introduce Virtual Enrichment-based Signaling Protein-activity Analysis (VESPA)-an algorithm designed to elucidate mechanisms of cell response and adaptation to drug perturbations-and use it to analyze 7-point phosphoproteomic time series from colorectal cancer cells treated with clinically-relevant inhibitors and control media. Interrogating tumor-specific enzyme/substrate interactions accurately infers kinase and phosphatase activity, based on their substrate phosphorylation state, effectively accounting for signal crosstalk and sparse phosphoproteome coverage. The analysis elucidates time-dependent signaling pathway response to each drug perturbation and, more importantly, cell adaptive response and rewiring, experimentally confirmed by CRISPR knock-out assays, suggesting broad applicability to cancer and other diseases.

Epigenetic targeting of PGBD5-dependent DNA damage in SMARCB1-deficient sarcomas

Despite the potential of targeted epigenetic therapies, most cancers do not respond to current epigenetic drugs. The Polycomb repressive complex EZH2 inhibitor tazemetostat was recently approved for the treatment of SMARCB1 -deficient epithelioid sarcomas, based on the functional antagonism between PRC2 and loss of SMARCB1. Through the analysis of tazemetostat-treated patient tumors, we recently defined key principles of their response and resistance to EZH2 epigenetic therapy. Here, using transcriptomic inference from SMARCB1 -deficient tumor cells, we nominate the DNA damage repair kinase ATR as a target for rational combination EZH2 epigenetic therapy. We show that EZH2 inhibition promotes DNA damage in epithelioid and rhabdoid tumor cells, at least in part via its induction of the transposase-derived PGBD5. We leverage this collateral synthetic lethal dependency to target PGBD5-dependent DNA damage by inhibition of ATR but not CHK1 using elimusertib. Consequently, combined EZH2 and ATR inhibition improves therapeutic responses in diverse patient-derived epithelioid and rhabdoid tumors in vivo . This advances a combination epigenetic therapy based on EZH2-PGBD5 synthetic lethal dependency suitable for immediate translation to clinical trials for patients.

Subtype-selective prenylated isoflavonoids disrupt regulatory drivers of MYCN-amplified cancers

Transcription factors have proven difficult to target with small molecules because they lack pockets necessary for potent binding. Disruption of protein expression can suppress targets and enable therapeutic intervention. To this end, we developed a drug discovery workflow that incorporates cell-line-selective screening and high-throughput expression profiling followed by regulatory network analysis to identify compounds that suppress regulatory drivers of disease. Applying this approach to neuroblastoma (NBL), we screened bioactive molecules in cell lines representing its MYC-dependent (MYCNA) and mesenchymal (MES) subtypes to identify selective compounds, followed by PLATESeq profiling of treated cells. This revealed compounds that disrupt a sub-network of MYCNA-specific regulatory proteins, resulting in MYCN degradation in vivo. The top hit was isopomiferin, a prenylated isoflavonoid that inhibited casein kinase 2 (CK2) in cells. Isopomiferin and its structural analogs inhibited MYC and MYCN in NBL and lung cancer cells, highlighting the general MYC-inhibiting potential of this unique scaffold.

Elucidation and Pharmacologic Targeting of Master Regulator Dependencies in Coexisting Diffuse Midline Glioma Subpopulations

Diffuse Midline Gliomas (DMGs) are universally fatal, primarily pediatric malignancies affecting the midline structures of the central nervous system. Despite decades of clinical trials, treatment remains limited to palliative radiation therapy. A major challenge is the coexistence of molecularly distinct malignant cell states with potentially orthogonal drug sensitivities. To address this challenge, we leveraged established network-based methodologies to elucidate Master Regulator (MR) proteins representing mechanistic, non-oncogene dependencies of seven coexisting subpopulations identified by single-cell analysis-whose enrichment in essential genes was validated by pooled CRISPR/Cas9 screens. Perturbational profiles of 372 clinically relevant drugs helped identify those able to invert the activity of subpopulation-specific MRs for follow-up in vivo validation. While individual drugs predicted to target individual subpopulations-including avapritinib, larotrectinib, and ruxolitinib-produced only modest tumor growth reduction in orthotopic models, systemic co-administration induced significant survival extension, making this approach a valuable contribution to the rational design of combination therapy.

Phenotypic maps for precision medicine: a promising systems biology tool for assessing therapy response and resistance at a personalized level

In this perspective we discuss how tumor heterogeneity and therapy resistance necessitate a focus on more personalized approaches, prompting a shift toward precision medicine. At the heart of the shift towards personalized medicine, omics-driven systems biology becomes a driving force as it leverages high-throughput technologies and novel bioinformatics tools. These enable the creation of systems-based maps, providing a comprehensive view of individual tumor's functional plasticity. We highlight the innovative PHENOSTAMP program, which leverages high-dimensional data to construct a visually intuitive and user-friendly map. This map was created to encapsulate complex transitional states in cancer cells, such as Epithelial-Mesenchymal Transition (EMT) and Mesenchymal-Epithelial Transition (MET), offering a visually intuitive way to understand disease progression and therapeutic responses at single-cell resolution in relation to EMT-related single-cell phenotypes. Most importantly, PHENOSTAMP functions as a reference map, which allows researchers and clinicians to assess one clinical specimen at a time in relation to their phenotypic heterogeneity, setting the foundation on constructing phenotypic maps for personalized medicine. This perspective argues that such dynamic predictive maps could also catalyze the development of personalized cancer treatment. They hold the potential to transform our understanding of cancer biology, providing a foundation for a future where therapy is tailored to each patient's unique molecular and cellular tumor profile. As our knowledge of cancer expands, these maps can be continually refined, ensuring they remain a valuable tool in precision oncology.

Elucidating Compound Mechanism of Action and Polypharmacology with a Large-scale Perturbational Profile Compendium

The Mechanism of Action (MoA) of a drug is generally represented as a small, non-tissue-specific repertoire of high-affinity binding targets. Yet, drug activity and polypharmacology are increasingly associated with a broad range of off-target and tissue-specific effector proteins. To address this challenge, we have implemented an efficient integrative experimental and computational framework leveraging the systematic generation and analysis of drug perturbational profiles representing >700 FDA-approved and experimental oncology drugs, in cell lines selected as high-fidelity models of 23 aggressive tumor subtypes. Protein activity-based analyses revealed highly reproducible, drug-mediated modulation of tissue-specific targets, leading to generation of a proteome-wide polypharmacology map, characterization of MoA-related drug clusters and off-target effects, and identification and experimental validation of novel, tissue-specific inhibitors of undruggable oncoproteins. The proposed framework, which is easily extended to elucidating the MoA of novel small-molecule libraries, could help support more systematic and quantitative approaches to precision oncology.

Computational drug prediction in hepatoblastoma by integrating pan-cancer transcriptomics with pharmacological response

BACKGROUND AND AIMS

Hepatoblastoma (HB) is the predominant form of pediatric liver cancer, though it remains exceptionally rare. While treatment outcomes for children with HB have improved, patients with advanced tumors face limited therapeutic choices. Additionally, survivors often suffer from long-term adverse effects due to treatment, including ototoxicity, cardiotoxicity, delayed growth, and secondary tumors. Consequently, there is a pressing need to identify new and effective therapeutic strategies for patients with HB. Computational methods to predict drug sensitivity from a tumor's transcriptome have been successfully applied for some common adult malignancies, but specific efforts in pediatric cancers are lacking because of the paucity of data.

APPROACH AND RESULTS

In this study, we used DrugSense to assess drug efficacy in patients with HB, particularly those with the aggressive C2 subtype associated with poor clinical outcomes. Our method relied on publicly available collections of pan-cancer transcriptional profiles and drug responses across 36 tumor types and 495 compounds. The drugs predicted to be most effective were experimentally validated using patient-derived xenograft models of HB grown in vitro and in vivo. We thus identified 2 cyclin-dependent kinase 9 inhibitors, alvocidib and dinaciclib as potent HB growth inhibitors for the high-risk C2 molecular subtype. We also found that in a cohort of 46 patients with HB, high cyclin-dependent kinase 9 tumor expression was significantly associated with poor prognosis.

CONCLUSIONS

Our work proves the usefulness of computational methods trained on pan-cancer data sets to reposition drugs in rare pediatric cancers such as HB, and to help clinicians in choosing the best treatment options for their patients.

Network-based approach elucidates critical genes in BRCA subtypes and chemotherapy response in Triple Negative Breast Cancer

Breast cancers exhibit substantial transcriptional heterogeneity, posing a significant challenge to the prediction of treatment response and prognostication of outcomes. Especially, translation of TNBC subtypes to the clinic remains a work in progress, in part because of a lack of clear transcriptional signatures distinguishing the subtypes. Our recent network-based approach, PathExt, demonstrates that global transcriptional changes in a disease context are likely mediated by a small number of key genes, and these mediators may better reflect functional or translationally relevant heterogeneity. We apply PathExt to 1059 BRCA tumors and 112 healthy control samples across 4 subtypes to identify frequent, key-mediator genes in each BRCA subtype. Compared to conventional differential expression analysis, PathExt-identified genes (1) exhibit greater concordance across tumors, revealing shared as well as BRCA subtype-specific biological processes, (2) better recapitulate BRCA-associated genes in multiple benchmarks, and (3) exhibit greater dependency scores in BRCA subtype-specific cancer cell lines. Single cell transcriptomes of BRCA subtype tumors reveal a subtype-specific distribution of PathExt-identified genes in multiple cell types from the tumor microenvironment. Application of PathExt to a TNBC chemotherapy response dataset identified TNBC subtype-specific key genes and biological processes associated with resistance. We described putative drugs that target top novel genes potentially mediating drug resistance. Overall, PathExt applied to breast cancer refines previous views of gene expression heterogeneity and identifies potential mediators of TNBC subtypes, including potential therapeutic targets.

NetBID2 provides comprehensive hidden driver analysis

Many signaling and other genes known as "hidden" drivers may not be genetically or epigenetically altered or differentially expressed at the mRNA or protein levels, but, rather, drive a phenotype such as tumorigenesis via post-translational modification or other mechanisms. However, conventional approaches based on genomics or differential expression are limited in exposing such hidden drivers. Here, we present a comprehensive algorithm and toolkit NetBID2 (data-driven network-based Bayesian inference of drivers, version 2), which reverse-engineers context-specific interactomes and integrates network activity inferred from large-scale multi-omics data, empowering the identification of hidden drivers that could not be detected by traditional analyses. NetBID2 has substantially re-engineered the previous prototype version by providing versatile data visualization and sophisticated statistical analyses, which strongly facilitate researchers for result interpretation through end-to-end multi-omics data analysis. We demonstrate the power of NetBID2 using three hidden driver examples. We deploy NetBID2 Viewer, Runner, and Cloud apps with 145 context-specific gene regulatory and signaling networks across normal tissues and paediatric and adult cancers to facilitate end-to-end analysis, real-time interactive visualization and cloud-based data sharing. NetBID2 is freely available at https://jyyulab.github.io/NetBID .