oncoPredict: an R package for predicting in vivo or cancer patient drug response and biomarkers from cell line screening data.

Cell line drug screening datasets can be utilized for a range of different drug discovery applications from drug biomarker discovery to building translational models of drug response. Previously, we described three separate methodologies to (1) correct for general levels of drug sensitivity to enable drug-specific biomarker discovery, (2) predict clinical drug response in patients and (3) associate these predictions with clinical features to perform in vivo drug biomarker discovery. Here, we unite and update these methodologies into one R package (oncoPredict) to facilitate the development and adoption of these tools. This new OncoPredict R package can be applied to various in vitro and in vivo contexts for drug and biomarker discovery.

A review of computational methods for predicting cancer drug response at the single-cell level through integration with bulk RNAseq data.

Cancer treatment failure is often attributed to tumor heterogeneity, where diverse malignant cell clones exist within a patient. Despite a growing understanding of heterogeneous tumor cells depicted by single-cell RNA sequencing (scRNA-seq), there is still a gap in the translation of such knowledge into treatment strategies tackling the pervasive issue of therapy resistance. In this review, we survey methods leveraging large-scale drug screens to generate cellular sensitivities to various therapeutics. These methods enable efficient drug screens in scRNA-seq data and serve as the bedrock of drug discovery for specific cancer cell groups. We envision that they will become an indispensable tool for tailoring patient care in the era of heterogeneity-aware precision medicine.

Sex differences in the molecular profile of adult diffuse glioma are shaped by IDH status and tumor microenvironment.

BACKGROUND: Sex differences in adult diffuse glioma (ADG) are well-established clinically, yet the underlying molecular mechanisms remain inadequately understood. Here, we aim to reveal molecular features and cellular compositions unique to each sex in ADG to comprehend the role of sex in disease etiology. METHODS: We quantified sex differences in transcriptome of ADG using multiple independent glioma patient datasets. Next, we delved into the single-cell landscape to examine sex differences in gene expression and cellular composition. To explore how sex influences disease progression, we analyzed paired samples from primary and recurrent ADG cases, aiming to identify sex-specific differences in molecular and cellular features. RESULTS: Our analysis revealed that mutations in isocitrate dehydrogenase (IDH) genes and the tumor microenvironment emerged as primary influencers of sex-differential molecular enrichments. In IDHwt tumors, genes in neuronal signaling pathway are found to be enriched in male tumors, while genes in hypoxia and inflammatory response pathways are enriched in female tumors. This pattern was reversed in IDHmut gliomas. We hypothesized that these distinctions could be attributed to heterogeneous cellular composition between sexes. Using single-cell data, we observed distinctive patterns of sex differences in cell states, cell composition and cell-cell interaction in IDHwt and IDHmut tumors separately. Further, by comparing molecular changes in paired primary and recurrent ADG samples, we identified sex-specific differences in molecular characteristics and cellular compositions of recurrent tumors. CONCLUSION: Our results provide a comprehensive multi-level characterization of sex differences in ADG, such findings provide novel insights into glioma disease progression in each sex.

Integration of Pan-Cancer Cell Line and Single-Cell Transcriptomic Profiles Enables Inference of Therapeutic Vulnerabilities in Heterogeneous Tumors.

Single-cell RNA sequencing (scRNA-seq) greatly advanced the understanding of intratumoral heterogeneity by identifying distinct cancer cell subpopulations. However, translating biological differences into treatment strategies is challenging due to a lack of tools to facilitate efficient drug discovery that tackles heterogeneous tumors. Developing such approaches requires accurate prediction of drug response at the single-cell level to offer therapeutic options to specific cell subpopulations. Here, we developed a transparent computational framework (nicknamed scIDUC) to predict therapeutic efficacies on an individual cell basis by integrating single-cell transcriptomic profiles with large, data-rich pan-cancer cell line screening data sets. This method achieved high accuracy in separating cells into their correct cellular drug response statuses. In three distinct prospective tests covering different diseases (rhabdomyosarcoma, pancreatic ductal adenocarcinoma, and castration-resistant prostate cancer), the predicted results using scIDUC were accurate and mirrored biological expectations. In the first two tests, the framework identified drugs for cell subpopulations that were resistant to standard-of-care (SOC) therapies due to intrinsic resistance or tumor microenvironmental effects, and the results showed high consistency with experimental findings from the original studies. In the third test using newly generated SOC therapy-resistant cell lines, scIDUC identified efficacious drugs for the resistant line, and the predictions were validated with in vitro experiments. Together, this study demonstrates the potential of scIDUC to quickly translate scRNA-seq data into drug responses for individual cells, displaying the potential as a tool to improve the treatment of heterogenous tumors. SIGNIFICANCE: A versatile method that infers cell-level drug response in scRNA-seq data facilitates the development of therapeutic strategies to target heterogeneous subpopulations within a tumor and address issues such as treatment failure and resistance.

Computational Modeling to Identify Drugs Targeting Metastatic Castration-Resistant Prostate Cancer Characterized by Heightened Glycolysis.

Metastatic castration-resistant prostate cancer (mCRPC) remains a deadly disease due to a lack of efficacious treatments. The reprogramming of cancer metabolism toward elevated glycolysis is a hallmark of mCRPC. Our goal is to identify therapeutics specifically associated with high glycolysis. Here, we established a computational framework to identify new pharmacological agents for mCRPC with heightened glycolysis activity under a tumor microenvironment, followed by in vitro validation. First, using our established computational tool, OncoPredict, we imputed the likelihood of drug responses to approximately 1900 agents in each mCRPC tumor from two large clinical patient cohorts. We selected drugs with predicted sensitivity highly correlated with glycolysis scores. In total, 77 drugs predicted to be more sensitive in high glycolysis mCRPC tumors were identified. These drugs represent diverse mechanisms of action. Three of the candidates, ivermectin, CNF2024, and P276-00, were selected for subsequent vitro validation based on the highest measured drug responses associated with glycolysis/OXPHOS in pan-cancer cell lines. By decreasing the input glucose level in culture media to mimic the mCRPC tumor microenvironments, we induced a high-glycolysis condition in PC3 cells and validated the projected higher sensitivity of all three drugs under this condition (p < 0.0001 for all drugs). For biomarker discovery, ivermectin and P276-00 were predicted to be more sensitive to mCRPC tumors with low androgen receptor activities and high glycolysis activities (AR(low)Gly(high)). In addition, we integrated a protein-protein interaction network and topological methods to identify biomarkers for these drug candidates. EEF1B2 and CCNA2 were identified as key biomarkers for ivermectin and CNF2024, respectively, through multiple independent biomarker nomination pipelines. In conclusion, this study offers new efficacious therapeutics beyond traditional androgen-deprivation therapies by precisely targeting mCRPC with high glycolysis.

Integration of Computational Pipeline to Streamline Efficacious Drug Nomination and Biomarker Discovery in Glioblastoma.

Glioblastoma multiforme (GBM) is the deadliest, most heterogeneous, and most common brain cancer in adults. Not only is there an urgent need to identify efficacious therapeutics, but there is also a great need to pair these therapeutics with biomarkers that can help tailor treatment to the right patient populations. We built patient drug response models by integrating patient tumor transcriptome data with high-throughput cell line drug screening data as well as Bayesian networks to infer relationships between patient gene expression and drug response. Through these discovery pipelines, we identified agents of interest for GBM to be effective across five independent patient cohorts and in a mouse avatar model: among them are a number of MEK inhibitors (MEKis). We also predicted phosphoglycerate dehydrogenase enzyme (PHGDH) gene expression levels to be causally associated with MEKi efficacy, where knockdown of this gene increased tumor sensitivity to MEKi and overexpression led to MEKi resistance. Overall, our work demonstrated the power of integrating computational approaches. In doing so, we quickly nominated several drugs with varying known mechanisms of action that can efficaciously target GBM. By simultaneously identifying biomarkers with these drugs, we also provide tools to select the right patient populations for subsequent evaluation.

Revealing Pan-Histology Immunomodulatory Targets in Pediatric Central Nervous System Tumors.

BACKGROUND: The application of immunotherapy for pediatric CNS malignancies has been limited by the poorly understood immune landscape in this context. The aim of this study was to uncover the mechanisms of immune suppression common among pediatric brain tumors. METHODS: We apply an immunologic clustering algorithm validated by The Cancer Genome Atlas Project to an independent pediatric CNS transcriptomic dataset. Within the clusters, the mechanisms of immunosuppression are explored via tumor microenvironment deconvolution and survival analyses to identify relevant immunosuppressive genes with translational relevance. RESULTS: High-grade diseases fall predominantly within an immunosuppressive subtype (C4) that independently lowers overall survival time and where common immune checkpoints (e.g., PDL1, CTLA4) are less relevant. Instead, we identify several alternative immunomodulatory targets with relevance across histologic diseases. Specifically, we show how the mechanism of EZH2 inhibition to enhance tumor immunogenicity in vitro via the upregulation of MHC class 1 is applicable to a pediatric CNS oncologic context. Meanwhile, we identify that the C3 (inflammatory) immune subtype is more common in low-grade diseases and find that immune checkpoint inhibition may be an effective way to curb progression for this subset. CONCLUSIONS: Three predominant immunologic clusters are identified across pediatric brain tumors. Among high-risk diseases, the predominant immune cluster is associated with recurrent immunomodulatory genes that influence immune infiltrate, including a subset that impacts survival across histologies.

Inferring therapeutic vulnerability within tumors through integration of pan-cancer cell line and single-cell transcriptomic profiles.

Single-cell RNA sequencing greatly advanced our understanding of intratumoral heterogeneity through identifying tumor subpopulations with distinct biologies. However, translating biological differences into treatment strategies is challenging, as we still lack tools to facilitate efficient drug discovery that tackles heterogeneous tumors. One key component of such approaches tackles accurate prediction of drug response at the single-cell level to offer therapeutic options to specific cell subpopulations. Here, we present a transparent computational framework (nicknamed scIDUC) to predict therapeutic efficacies on an individual-cell basis by integrating single-cell transcriptomic profiles with large, data-rich pan-cancer cell line screening datasets. Our method achieves high accuracy, with predicted sensitivities easily able to separate cells into their true cellular drug resistance status as measured by effect size (Cohen's d > 1.0). More importantly, we examine our method's utility with three distinct prospective tests covering different diseases (rhabdomyosarcoma, pancreatic ductal adenocarcinoma, and castration-resistant prostate cancer), and in each our predicted results are accurate and mirrored biological expectations. In the first two, we identified drugs for cell subpopulations that are resistant to standard-of-care (SOC) therapies due to intrinsic resistance or effects of tumor microenvironments. Our results showed high consistency with experimental findings from the original studies. In the third test, we generated SOC therapy resistant cell lines, used scIDUC to identify efficacious drugs for the resistant line, and validated the predictions with in-vitro experiments. Together, scIDUC quickly translates scRNA-seq data into drug response for individual cells, displaying the potential as a first-line tool for nuanced and heterogeneity-aware drug discovery.