Synthetic lethality-mediated precision oncology via the tumor transcriptome

Harnessing genetic interactions for prediction of immune checkpoint inhibitors response signature in cancer cells

Genetic interactions (GIs) refer to two altered genes having a combined effect that is not seen individually. They play a crucial role in influencing drug efficacy. We utilized CGIdb 2.0 (http://www.medsysbio.org/CGIdb2/), an updated database of comprehensively published GIs information, encompassing synthetic lethality (SL), synthetic viability (SV), and chemical-genetic interactions. CGIdb 2.0 elucidates GIs relationships between or within protein complex models by integrating protein-protein physical interactions. Additionally, we introduced GENIUS (GENetic Interactions mediated drUg Signature) to leverage GIs for identifying the response signature of immune checkpoint inhibitors (ICIs). GENIUS identified high MAP4K4 expression as a resistant signature and high HERC4 expression as a sensitive signature for ICIs treatment. Melanoma patients with high expression of MAP4K4 were associated with decreased efficacy and poorer survival following ICIs treatment. Conversely, overexpression of HERC4 in melanoma patients correlated with a positive response to ICIs. Notably, HERC4 enhances sensitivity to immunotherapy by facilitating antigen presentation. Analyses of immune cell infiltration and single-cell data revealed that B cells expressing MAP4K4 may contribute to resistance to ICIs in melanoma. Overall, CGIdb 2.0, provides integrated GIs data, thus serving as a crucial tool for exploring drug effects.

Identifying actionable synthetically lethal cancer gene pairs using mutual exclusivity

Mutually exclusive loss-of-function alterations in gene pairs are those that occur together less frequently than may be expected and may denote a synthetically lethal relationship (SSL) between the genes. SSLs can be exploited therapeutically to selectively kill cancer cells. Here, we analysed mutation, copy number variation, and methylation levels in samples from The Cancer Genome Atlas, using the hypergeometric and the Poisson binomial tests to identify mutually exclusive inactivated genes. We focused on gene pairs where one is an inactivated tumour suppressor and the other a gene whose protein product can be inhibited by known drugs. This provided an abundance of potential targeted therapeutics and repositioning opportunities for several cancers. These data are available on the MexDrugs website, https://bioinformaticslab.sussex.ac.uk/mexdrugs.

TimiGP-Response: the pan-cancer immune landscape associated with response to immunotherapy

Accumulating evidence suggests that the tumor immune microenvironment (TIME) significantly influences the response to immunotherapy, yet this complex relationship remains elusive. To address this issue, we developed TimiGP-Response (TIME Illustration based on Gene Pairing designed for immunotherapy Response), a computational framework leveraging single-cell and bulk transcriptomic data, along with response information, to construct cell-cell interaction networks associated with responders and estimate the role of immune cells in treatment response. This framework was showcased in triple-negative breast cancer treated with immune checkpoint inhibitors targeting the PD-1:PD-L1 interaction, and orthogonally validated with imaging mass cytometry. As a result, we identified CD8+ GZMB+ T cells associated with responders and its interaction with regulatory T cells emerged as a potential feature for selecting patients who may benefit from these therapies. Subsequently, we analyzed 3,410 patients with seven cancer types (melanoma, non-small cell lung cancer, renal cell carcinoma, metastatic urothelial carcinoma, hepatocellular carcinoma, breast cancer, and esophageal cancer) treated with various immunotherapies and combination therapies, as well as several chemo- and targeted therapies as controls. Using TimiGP-Response, we depicted the pan-cancer immune landscape associated with immunotherapy response at different resolutions. At the TIME level, CD8 T cells and CD4 memory T cells were associated with responders, while anti-inflammatory (M2) macrophages and mast cells were linked to non-responders across most cancer types and datasets. Given that T cells are the primary targets of these immunotherapies and our TIME analysis highlights their importance in response to treatment, we portrayed the pan-caner landscape on 40 T cell subtypes. Notably, CD8+ and CD4+ GZMK+ effector memory T cells emerged as crucial across all cancer types and treatments, while IL-17-producing CD8+ T cells were top candidates associated with immunotherapy non-responders. In summary, this study provides a computational method to study the association between TIME and response across the pan-cancer immune landscape, offering resources and insights into immune cell interactions and their impact on treatment efficacy.

Advances in Integrated Multi-omics Analysis for Drug-Target Identification

As an essential component of modern drug discovery, the role of drug-target identification is growing increasingly prominent. Additionally, single-omics technologies have been widely utilized in the process of discovering drug targets. However, it is difficult for any single-omics level to clearly expound the causal connection between drugs and how they give rise to the emergence of complex phenotypes. With the progress of large-scale sequencing and the development of high-throughput technologies, the tendency in drug-target identification has shifted towards integrated multi-omics techniques, gradually replacing traditional single-omics techniques. Herein, this review centers on the recent advancements in the domain of integrated multi-omics techniques for target identification, highlights the common multi-omics analysis strategies, briefly summarizes the selection of multi-omics analysis tools, and explores the challenges of existing multi-omics analyses, as well as the applications of multi-omics technology in drug-target identification.

A muti-modal feature fusion method based on deep learning for predicting immunotherapy response

Immune checkpoint therapy (ICT) has greatly improved the survival of cancer patients in the past few years, but only a small number of patients respond to ICT. To predict ICT response, we developed a multi-modal feature fusion model based on deep learning (MFMDL). This model utilizes graph neural networks to map gene-gene relationships in gene networks to low dimensional vector spaces, and then fuses biological pathway features and immune cell infiltration features to make robust predictions of ICT. We used five datasets to validate the predictive performance of the MFMDL. These five datasets span multiple types of cancer, including melanoma, lung cancer, and gastric cancer. We found that the prediction performance of multi-modal feature fusion model based on deep learning is superior to other traditional ICT biomarkers, such as ICT targets or tumor microenvironment-associated markers. In addition, we also conducted ablation experiments to demonstrate the necessity of fusing different modal features, which can improve the prediction accuracy of the model.

Predicting transcriptional outcomes of novel multigene perturbations with GEARS

Understanding cellular responses to genetic perturbation is central to numerous biomedical applications, from identifying genetic interactions involved in cancer to developing methods for regenerative medicine. However, the combinatorial explosion in the number of possible multigene perturbations severely limits experimental interrogation. Here, we present graph-enhanced gene activation and repression simulator (GEARS), a method that integrates deep learning with a knowledge graph of gene-gene relationships to predict transcriptional responses to both single and multigene perturbations using single-cell RNA-sequencing data from perturbational screens. GEARS is able to predict outcomes of perturbing combinations consisting of genes that were never experimentally perturbed. GEARS exhibited 40% higher precision than existing approaches in predicting four distinct genetic interaction subtypes in a combinatorial perturbation screen and identified the strongest interactions twice as well as prior approaches. Overall, GEARS can predict phenotypically distinct effects of multigene perturbations and thus guide the design of perturbational experiments.

PERCEPTION predicts patient response and resistance to treatment using single-cell transcriptomics of their tumors

Tailoring optimal treatment for individual cancer patients remains a significant challenge. To address this issue, we developed PERCEPTION (PERsonalized Single-Cell Expression-Based Planning for Treatments In ONcology), a precision oncology computational pipeline. Our approach uses publicly available matched bulk and single-cell (sc) expression profiles from large-scale cell-line drug screens. These profiles help build treatment response models based on patients' sc-tumor transcriptomics. PERCEPTION demonstrates success in predicting responses to targeted therapies in cultured and patient-tumor-derived primary cells, as well as in two clinical trials for multiple myeloma and breast cancer. It also captures the resistance development in patients with lung cancer treated with tyrosine kinase inhibitors. PERCEPTION outperforms published state-of-the-art sc-based and bulk-based predictors in all clinical cohorts. PERCEPTION is accessible at https://github.com/ruppinlab/PERCEPTION . Our work, showcasing patient stratification using sc-expression profiles of their tumors, will encourage the adoption of sc-omics profiling in clinical settings, enhancing precision oncology tools based on sc-omics.

A Phase II, Open-Label, Randomized Trial of Durvalumab With Olaparib or Cediranib in Patients With Mismatch Repair-Proficient Colorectal or Pancreatic Cancer

BACKGROUND

The use of immunotherapy in mismatch repair proficient colorectal cancer (pMMR-CRC) or pancreatic adenocarcinoma (PDAC) is associated with limited efficacy. DAPPER (NCT03851614) is a phase 2, basket study randomizing patients with pMMR CRC or PDAC to durvalumab with olaparib (durvalumab + olaparib) or durvalumab with cediranib (durvalumab + cediranib).

METHODS

PDAC or pMMR-CRC patients were randomized to either durvalumab+olaparib (arm A), or durvalumab + cediranib (arm B). Co-primary endpoints included pharmacodynamic immune changes in the tumor microenvironment (TME) and safety. Objective response rate, progression-free survival (PFS) and overall survival (OS) were determined. Paired tumor samples were analyzed by multiplexed immunohistochemistry and RNA-sequencing.

RESULTS

A total of 31 metastatic pMMR-CRC patients were randomized to arm A (n = 16) or B (n = 15). In 28 evaluable patients, 3 patients had stable disease (SD) (2 patients treated with durvalumab + olaparib and 1 patient treated with durvalumab + cediranib) while 25 had progressive disease (PD). Among patients with PDAC (n = 19), 9 patients were randomized to arm A and 10 patients were randomized to arm B. In 18 evaluable patients, 1 patient had a partial response (unconfirmed) with durvalumab + cediranib, 1 patient had SD with durvalumab + olaparib while 16 had PD. Safety profile was manageable and no grade 4-5 treatment-related adverse events were observed in either arm A or B. No significant changes were observed for CD3+/CD8+ immune infiltration in on-treatment biopsies as compared to baseline for pMMR-CRC and PDAC independent of treatment arms. Increased tumor-infiltrating lymphocytes at baseline, low baseline CD68+ cells and different immune gene expression signatures at baseline were associated with outcomes.

CONCLUSIONS

In patients with pMMR-CRC or PDAC, durvalumab + olaparib and durvalumab + cediranib showed limited antitumor activity. Different immune components of the TME were associated with treatment outcomes.

Multicellular ecotypes shape progression of lung adenocarcinoma from ground-glass opacity toward advanced stages

Lung adenocarcinoma is a type of cancer that exhibits a wide range of clinical radiological manifestations, from ground-glass opacity (GGO) to pure solid nodules, which vary greatly in terms of their biological characteristics. Our current understanding of this heterogeneity is limited. To address this gap, we analyze 58 lung adenocarcinoma patients via machine learning, single-cell RNA sequencing (scRNA-seq), and whole-exome sequencing, and we identify six lung multicellular ecotypes (LMEs) correlating with distinct radiological patterns and cancer cell states. Notably, GGO-associated neoantigens in early-stage cancers are recognized by CD8+ T cells, indicating an immune-active environment, while solid nodules feature an immune-suppressive LME with exhausted CD8+ T cells, driven by specific stromal cells such as CTHCR1+ fibroblasts. This study also highlights EGFR(L858R) neoantigens in GGO samples, suggesting potential CD8+ T cell activation. Our findings offer valuable insights into lung adenocarcinoma heterogeneity, suggesting avenues for targeted therapies in early-stage disease.

A comprehensive review of machine learning techniques for multi-omics data integration: challenges and applications in precision oncology

Multi-omics data play a crucial role in precision medicine, mainly to understand the diverse biological interaction between different omics. Machine learning approaches have been extensively employed in this context over the years. This review aims to comprehensively summarize and categorize these advancements, focusing on the integration of multi-omics data, which includes genomics, transcriptomics, proteomics and metabolomics, alongside clinical data. We discuss various machine learning techniques and computational methodologies used for integrating distinct omics datasets and provide valuable insights into their application. The review emphasizes both the challenges and opportunities present in multi-omics data integration, precision medicine and patient stratification, offering practical recommendations for method selection in various scenarios. Recent advances in deep learning and network-based approaches are also explored, highlighting their potential to harmonize diverse biological information layers. Additionally, we present a roadmap for the integration of multi-omics data in precision oncology, outlining the advantages, challenges and implementation difficulties. Hence this review offers a thorough overview of current literature, providing researchers with insights into machine learning techniques for patient stratification, particularly in precision oncology. Contact:  anirban@klyuniv.ac.in.

Cell-cell communication network-based interpretable machine learning predicts cancer patient response to immune checkpoint inhibitors

Immune checkpoint inhibitors (ICIs) have revolutionized cancer treatment. However, only some patients respond to ICIs, and current biomarkers for ICI efficacy have limited performance. Here, we devised an interpretable machine learning (ML) model trained using patient-specific cell-cell communication networks (CCNs) decoded from the patient's bulk tumor transcriptome. The model could (i) predict ICI efficacy for patients across four cancer types (median AUROC: 0.79) and (ii) identify key communication pathways with crucial players responsible for patient response or resistance to ICIs by analyzing more than 700 ICI-treated patient samples from 11 cohorts. The model prioritized chemotaxis communication of immune-related cells and growth factor communication of structural cells as the key biological processes underlying response and resistance to ICIs, respectively. We confirmed the key communication pathways and players at the single-cell level in patients with melanoma. Our network-based ML approach can be used to expand ICIs' clinical benefits in cancer patients.

Chromosome 7 to the rescue: overcoming chromosome 10 loss in gliomas

The co-occurrence of chromosome 10 loss and chromosome 7 gain in gliomas is the most frequent loss-gain co-aneuploidy pair in human cancers, a phenomenon that has been investigated without resolution since the late 1980s. Expanding beyond previous gene-centric studies, we investigate the co-occurrence in a genome-wide manner taking an evolutionary perspective. First, by mining large tumor aneuploidy data, we predict that the more likely order is 10 loss followed by 7 gain. Second, by analyzing extensive genomic and transcriptomic data from both patients and cell lines, we find that this co-occurrence can be explained by functional rescue interactions that are highly enriched on 7, which can possibly compensate for any detrimental consequences arising from the loss of 10. Finally, by analyzing transcriptomic data from normal, non-cancerous, human brain tissues, we provide a plausible reason why this co-occurrence happens preferentially in cancers originating in certain regions of the brain.

Predicting patient outcomes after treatment with immune checkpoint blockade: A review of biomarkers derived from diverse data modalities

Immune checkpoint blockade (ICB) therapy targeting cytotoxic T-lymphocyte-associated protein 4, programmed death 1, and programmed death ligand 1 has shown durable remission and clinical success across different cancer types. However, patient outcomes vary among disease indications. Studies have identified prognostic biomarkers associated with immunotherapy response and patient outcomes derived from diverse data types, including next-generation bulk and single-cell DNA, RNA, T cell and B cell receptor sequencing data, liquid biopsies, and clinical imaging. Owing to inter- and intra-tumor heterogeneity and the immune system's complexity, these biomarkers have diverse efficacy in clinical trials of ICB. Here, we review the genetic and genomic signatures and image features of ICB studies for pan-cancer applications and specific indications. We discuss the advantages and disadvantages of computational approaches for predicting immunotherapy effectiveness and patient outcomes. We also elucidate the challenges of immunotherapy prognostication and the discovery of novel immunotherapy targets.

Predicting host-based, synthetic lethal antiviral targets from omics data

Traditional antiviral therapies often have limited effectiveness due to toxicity and the emergence of drug resistance. Host-based antivirals are an alternative, but can cause nonspecific effects. Recent evidence shows that virus-infected cells can be selectively eliminated by targeting synthetic lethal (SL) partners of proteins disrupted by viral infection. Thus, we hypothesized that genes depleted in CRISPR knockout (KO) screens of virus-infected cells may be enriched in SL partners of proteins altered by infection. To investigate this, we established a computational pipeline predicting antiviral SL drug targets. First, we identified SARS-CoV-2-induced changes in gene products via a large compendium of omics data. Second, we identified SL partners for each altered gene product. Last, we screened CRISPR KO data for SL partners required for cell viability in infected cells. Despite differences in virus-induced alterations detected by various omics data, they share many predicted SL targets, with significant enrichment in CRISPR KO-depleted datasets. Our comparison of SARS-CoV-2 and influenza infection data revealed potential broad-spectrum, host-based antiviral SL targets. This suggests that CRISPR KO data are replete with common antiviral targets due to their SL relationship with virus-altered states and that such targets can be revealed from analysis of omics datasets and SL predictions.

An immunosuppressive subtype of senescent tumor cells predicted worse immunotherapy response in lung adenocarcinoma

Senescent tumor cells (STCs) can induce immunosuppression, promoting tumor progression and therapy resistance. However, the specific characteristics of immunosuppressive STC have not been thoroughly investigated. This study aimed to characterize and elucidate the immunosuppressive phenotype of STC in lung adenocarcinoma by employing single-cell and bulk transcriptomics, as well as serum proteomics profiling. We identified senescence-related genes specific to tumors and identified Cluster10 of STC as the immunomodulatory subtype. Cluster10 exhibited a distinct secretome dominated by cytokines such as CXCL1, CXCL2, and CXCL8 and showed activation of transcription factors associated with cytokine secretion, including NFKB1, RELA, and STAT3. Notably, Cluster10 demonstrated the highest degree of intercellular communication among all cell types, with interactions as LGALS9-TIM3 and MIF-CD74. Furthermore, Cluster10 showed significant associations with poor prognosis and diminished response to immunotherapy. Analysis of serum proteomics data from our in-house cohort identified CXCL8 as a potential marker for predicting immunotherapeutic outcomes.

Synthetic viability induces resistance to immune checkpoint inhibitors in cancer cells

BACKGROUND

Immune checkpoint inhibitors (ICI) have revolutionized the treatment for multiple cancers. However, most of patients encounter resistance. Synthetic viability (SV) between genes could induce resistance. In this study, we established SV signature to predict the efficacy of ICI treatment for melanoma.

METHODS

We collected features and predicted SV gene pairs by random forest classifier. This work prioritized SV gene pairs based on CRISPR/Cas9 screens. SV gene pairs signature were constructed to predict the response to ICI for melanoma patients.

RESULTS

This study predicted robust SV gene pairs based on 14 features. Filtered by CRISPR/Cas9 screens, we identified 1,861 SV gene pairs, which were also related with prognosis across multiple cancer types. Next, we constructed the six SV pairs signature to predict resistance to ICI for melanoma patients. This study applied the six SV pairs signature to divide melanoma patients into high-risk and low-risk. High-risk melanoma patients were associated with worse response after ICI treatment. Immune landscape analysis revealed that high-risk melanoma patients had lower natural killer cells and CD8+ T cells infiltration.

CONCLUSIONS

In summary, the 14 features classifier accurately predicted robust SV gene pairs for cancer. The six SV pairs signature could predict resistance to ICI.

Increased interleukin-6/C-reactive protein levels are associated with the upregulation of the adenosine pathway and serve as potential markers of therapeutic resistance to immune checkpoint inhibitor-based therapies in non-small cell lung cancer

BACKGROUND

Systemic immune activation, hallmarked by C-reactive protein (CRP) and interleukin-6 (IL-6), can modulate antitumor immune responses. In this study, we evaluated the role of IL-6 and CRP in the stratification of patients with non-small cell lung cancer (NSCLC) treated with immune checkpoint inhibitors (ICIs). We also interrogated the underlying immunosuppressive mechanisms driven by the IL-6/CRP axis.

METHODS

In cohort A (n=308), we estimated the association of baseline CRP with objective response rate (ORR), progression-free survival (PFS), and overall survival (OS) in patients with NSCLC treated with ICIs alone or with chemo-immunotherapy (Chemo-ICI). Baseline tumor bulk RNA sequencing (RNA-seq) of lung adenocarcinomas (LUADs) treated with pembrolizumab (cohort B, n=59) was used to evaluate differential expression of purine metabolism, as well as correlate IL-6 expression with PFS. CODEFACS approach was applied to deconvolve cohort B to characterize the tumor microenvironment by reconstructing the cell-type-specific transcriptome from bulk expression. Using the LUAD cohort from The Cancer Genome Atlas (TCGA) we explored the correlation between IL-6 expression and adenosine gene signatures. In a third cohort (cohort C, n=18), plasma concentrations of CRP, adenosine 2a receptor (A2aR), and IL-6 were measured using ELISA.

RESULTS

In cohort A, 67.2% of patients had a baseline CRP≥10 mg/L (CRP-H). Patients with CRP-H achieved shorter OS (8.6 vs 14.8 months; p=0.006), shorter PFS (3.3 vs 6.6 months; p=0.013), and lower ORR (24.7% vs 46.3%; p=0.015). After adjusting for relevant clinical variables, CRP-H was confirmed as an independent predictor of increased risk of death (HR 1.51, 95% CI: 1.09 to 2.11) and lower probability of achieving disease response (OR 0.34, 95% CI: 0.13 to 0.89). In cohort B, RNA-seq analysis demonstrated higher IL-6 expression on tumor cells of non-responders, along with a shorter PFS (p<0.05) and enrichment of the purinergic pathway. Within the TCGA LUAD cohort, tumor IL-6 expression strongly correlated with the adenosine signature (R=0.65; p<2.2e-16). Plasma analysis in cohort C demonstrated that CRP-H patients had a greater median baseline level of A2aR (6.0 ng/mL vs 1.3 ng/mL; p=0.01).

CONCLUSIONS

This study demonstrates CRP as a readily available blood-based prognostic biomarker in ICI-treated NSCLC. Additionally, we elucidate a potential link of the CRP/IL-6 axis with the immunosuppressive adenosine signature pathway that could drive inferior outcomes to ICIs in NSCLC and also offer novel therapeutic avenues.

Functional Profiling of Soft Tissue Sarcoma Using Mechanistic Models

Soft tissue sarcoma is an umbrella term for a group of rare cancers that are difficult to treat. In addition to surgery, neoadjuvant chemotherapy has shown the potential to downstage tumors and prevent micrometastases. However, finding effective therapeutic targets remains a research challenge. Here, a previously developed computational approach called mechanistic models of signaling pathways has been employed to unravel the impact of observed changes at the gene expression level on the ultimate functional behavior of cells. In the context of such a mechanistic model, RNA-Seq counts sourced from the Recount3 resource, from The Cancer Genome Atlas (TCGA) Sarcoma project, and non-diseased sarcomagenic tissues from the Genotype-Tissue Expression (GTEx) project were utilized to investigate signal transduction activity through signaling pathways. This approach provides a precise view of the relationship between sarcoma patient survival and the signaling landscape in tumors and their environment. Despite the distinct regulatory alterations observed in each sarcoma subtype, this study identified 13 signaling circuits, or elementary sub-pathways triggering specific cell functions, present across all subtypes, belonging to eight signaling pathways, which served as predictors for patient survival. Additionally, nine signaling circuits from five signaling pathways that highlighted the modifications tumor samples underwent in comparison to normal tissues were found. These results describe the protective role of the immune system, suggesting an anti-tumorigenic effect in the tumor microenvironment, in the process of tumor cell detachment and migration, or the dysregulation of ion homeostasis. Also, the analysis of signaling circuit intermediary proteins suggests multiple strategies for therapy.

Prediction of cancer treatment response from histopathology images through imputed transcriptomics

Advances in artificial intelligence have paved the way for leveraging hematoxylin and eosin (H&E)-stained tumor slides for precision oncology. We present ENLIGHT-DeepPT, an approach for predicting response to multiple targeted and immunotherapies from H&E-slides. In difference from existing approaches that aim to predict treatment response directly from the slides, ENLIGHT-DeepPT is an indirect two-step approach consisting of (1) DeepPT, a new deep-learning framework that predicts genome-wide tumor mRNA expression from slides, and (2) ENLIGHT, which predicts response based on the DeepPT inferred expression values. DeepPT successfully predicts transcriptomics in all 16 TCGA cohorts tested and generalizes well to two independent datasets. Our key contribution is showing that ENLIGHT-DeepPT successfully predicts true responders in five independent patients' cohorts involving four different treatments spanning six cancer types with an overall odds ratio of 2.44, increasing the baseline response rate by 43.47% among predicted responders, without the need for any treatment data for training. Furthermore, its prediction accuracy on these datasets is comparable to a supervised approach predicting the response directly from the images, which needs to be trained and tested on the same cohort. ENLIGHT-DeepPT future application could provide clinicians with rapid treatment recommendations to an array of different therapies and importantly, may contribute to advancing precision oncology in developing countries.

Tumor aerobic glycolysis confers immune evasion through modulating sensitivity to T cell-mediated bystander killing via TNF-α

Metabolic reprogramming toward glycolysis is a hallmark of cancer malignancy. The molecular mechanisms by which the tumor glycolysis pathway promotes immune evasion remain to be elucidated. Here, by performing genome-wide CRISPR screens in murine tumor cells co-cultured with cytotoxic T cells (CTLs), we identified that deficiency of two important glycolysis enzymes, Glut1 (glucose transporter 1) and Gpi1 (glucose-6-phosphate isomerase 1), resulted in enhanced killing of tumor cells by CTLs. Mechanistically, Glut1 inactivation causes metabolic rewiring toward oxidative phosphorylation, which generates an excessive amount of reactive oxygen species (ROS). Accumulated ROS potentiate tumor cell death mediated by tumor necrosis factor alpha (TNF-α) in a caspase-8- and Fadd-dependent manner. Genetic and pharmacological inactivation of Glut1 sensitizes tumors to anti-tumor immunity and synergizes with anti-PD-1 therapy through the TNF-α pathway. The mechanistic interplay between tumor-intrinsic glycolysis and TNF-α-induced killing provides new therapeutic strategies to enhance anti-tumor immunity.

MAGCN: A Multiple Attention Graph Convolution Networks for Predicting Synthetic Lethality

Synthetic lethality (SL) is a potential cancer therapeutic strategy and drug discovery. Computational approaches to identify synthetic lethality genes have become an effective complement to wet experiments which are time consuming and costly. Graph convolutional networks (GCN) has been utilized to such prediction task as be good at capturing the neighborhood dependency in a graph. However, it is still a lack of the mechanism of aggregating the complementary neighboring information from various heterogeneous graphs. Here, we propose the Multiple Attention Graph Convolution Networks for predicting synthetic lethality (MAGCN). First, we obtain the functional similarity features and topological structure features of genes from different data sources respectively, such as Gene Ontology data and Protein-Protein Interaction. Then, graph convolutional network is utilized to accumulate the knowledge from neighbor nodes according to synthetic lethal associations. Meanwhile, we propose a multiple graphs attention model and construct a multiple graphs attention network to learn the contribution factors of different graphs to generate embedded representation by aggregating these graphs. Finally, the generated feature matrix is decoded to predict potential synthetic lethal interaction. Experimental results show that MAGCN is superior to other baseline methods. Case study demonstrates the ability of MAGCN to predict human SL gene pairs.

Predicting host-based, synthetic lethal antiviral targets from omics data

Traditional antiviral therapies often have limited effectiveness due to toxicity and development of drug resistance. Host-based antivirals, while an alternative, may lead to non-specific effects. Recent evidence shows that virus-infected cells can be selectively eliminated by targeting synthetic lethal (SL) partners of proteins disrupted by viral infection. Thus, we hypothesized that genes depleted in CRISPR KO screens of virus-infected cells may be enriched in SL partners of proteins altered by infection. To investigate this, we established a computational pipeline predicting SL drug targets of viral infections. First, we identified SARS-CoV-2-induced changes in gene products via a large compendium of omics data. Second, we identified SL partners for each altered gene product. Last, we screened CRISPR KO data for SL partners required for cell viability in infected cells. Despite differences in virus-induced alterations detected by various omics data, they share many predicted SL targets, with significant enrichment in CRISPR KO-depleted datasets. Comparing data from SARS-CoV-2 and influenza infections, we found possible broad-spectrum, host-based antiviral SL targets. This suggests that CRISPR KO data are replete with common antiviral targets due to their SL relationship with virus-altered states and that such targets can be revealed from analysis of omics datasets and SL predictions.

Prognostic and clinicopathological role of RACK1 for cancer patients: a systematic review and meta-analysis

Background

The receptor for activated C kinase 1 (RACK1) expression is associated with clinicopathological characteristics and the prognosis of various cancers; however, the conclusions are controversial. As a result, this study aimed to explore the clinicopathological and prognostic values of RACK1 expression in patients with cancer.

Methodology

PubMed, Embase, Web of Science, Cochrane Library, and Scopus were comprehensively explored from their inception to April 20, 2023, for selecting studies on the clinicopathological and prognostic role of RACK1 in patients with cancer that met the criteria for inclusion in this review. Pooled hazard ratios (HRs) and 95% confidence intervals (CIs) were used to assess the prognosis-predictive value of RACK1 expression, while pooled odds ratios (ORs) and 95% CIs were used to evaluate the correlation between RACK1 expression and the clinicopathological characteristics of patients with cancer. The quality of the included studies was evaluated using the Newcastle-Ottawa Scale.

Results

Twenty-two studies (13 on prognosis and 20 on clinicopathological characteristics) were included in this systematic review and meta-analysis. The findings indicated that high RACK1 expression was significantly associated with poor overall survival (HR = 1.62; 95% CI, 1.13-2.33; P = 0.009; I2 = 89%) and reversely correlated with disease-free survival/recurrence-free survival (HR = 1.87; 95% CI, 1.22-2.88; P = 0.004; I2 = 0%). Furthermore, increased RACK1 expression was significantly associated with lymphatic invasion/N+ stage (OR = 1.74; 95% CI, 1.04-2.90; P = 0.04; I2 = 79%) of tumors.

Conclusions

RACK1 may be a global predictive marker of poor prognosis in patients with cancer and unfavorable clinicopathological characteristics. However, further clinical studies are required to validate these findings.

TimiGP: Inferring cell-cell interactions and prognostic associations in the tumor immune microenvironment through gene pairs

Determining the prognostic association of different immune cell types in the tumor microenvironment is critical for understanding cancer biology and developing new therapeutic strategies. However, this is challenging in certain cancer types, where the abundance of different immune subsets is highly correlated. In this study, we develop a computational method named TimiGP to overcome this challenge. Based on bulk gene expression and survival data, TimiGP infers cell-cell interactions that reveal the association between immune cell relative abundance and prognosis. As demonstrated in metastatic melanoma, TimiGP prioritizes immune cells critical in prognosis based on the identified cell-cell interactions. Highly consistent results are obtained by TimiGP when applied to seven independent melanoma datasets and when different cell-type marker sets are used as inputs. Additionally, TimiGP can leverage single-cell RNA sequencing data to delineate the tumor immune microenvironment at high resolutions across a wide range of cancer types.

In Silico Discovery of 5'-Modified 7-Deoxy-7-ethynyl-4'-thioadenosine as a HASPIN Inhibitor and Its Synergistic Anticancer Effect with the PLK1 Inhibitor

Despite genetic perturbations resulting in embryo lethality for most mitotic kinases, loss of the histone H3 mitotic kinase HASPIN reveals no adverse effect in mice models, establishing HASPIN as a promising target for anticancer therapy. However, developing a HASPIN inhibitor from conventional pharmacophores poses a technical challenge as this atypical kinase shares slight similarities with eukaryotic protein kinases. Chemically modifying a cytotoxic 4'-thioadenosine analogue through high genotoxicity yielded several novel nongenotoxic kinase inhibitors. In silico apporoaches utilizing transcriptomic and chemical similarities with known compounds and KINOMEscan profiles unveiled the HASPIN inhibitor LJ4827. LJ4827's specificity and potency as a HASPIN inhibitor were verified through in vitro kinase assay and X-ray crystallography. HASPIN inhibition by LJ4827 reduced histone H3 phosphorylation and impeded Aurora B recruitment in cancer cell centromeres but not in noncancer cells. Through transcriptome analysis of lung cancer patients, PLK1 was determined as a druggable synergistic partner to complement HASPIN inhibition. Chemical or genetic PLK1 perturbation with LJ4827 effectuated pronounced lung cancer cytotoxicity in vitro and in vivo. Therefore, LJ4827 is a novel anticancer therapeutic for selectively impeding cancer mitosis through potent HASPIN inhibition, and simultaneous HASPIN and PLK1 interference is a promising therapeutic strategy for lung cancer.

Immune Resistance Mechanisms and the Road to Personalized Immunotherapy

What does the future of cancer immunotherapy look like and how do we get there? Find out where we've been and where we're headed in A Report on Resistance: The Road to personalized immunotherapy.

Genomic hallmarks and therapeutic implications of G0 cell cycle arrest in cancer

BACKGROUND

Therapy resistance in cancer is often driven by a subpopulation of cells that are temporarily arrested in a non-proliferative G0 state, which is difficult to capture and whose mutational drivers remain largely unknown.

RESULTS

We develop methodology to robustly identify this state from transcriptomic signals and characterise its prevalence and genomic constraints in solid primary tumours. We show that G0 arrest preferentially emerges in the context of more stable, less mutated genomes which maintain TP53 integrity and lack the hallmarks of DNA damage repair deficiency, while presenting increased APOBEC mutagenesis. We employ machine learning to uncover novel genomic dependencies of this process and validate the role of the centrosomal gene CEP89 as a modulator of proliferation and G0 arrest capacity. Lastly, we demonstrate that G0 arrest underlies unfavourable responses to various therapies exploiting cell cycle, kinase signalling and epigenetic mechanisms in single-cell data.

CONCLUSIONS

We propose a G0 arrest transcriptional signature that is linked with therapeutic resistance and can be used to further study and clinically track this state.

piTracer - Automatic reconstruction of molecular cascades for the identification of synergistic drug targets

Cancer cells frequently undergo metabolic reprogramming as a mechanism of resistance against chemotherapeutic drugs. Metabolomic profiling provides a direct readout of metabolic changes and can thus be used to identify these tumor escape mechanisms. Here, we introduce piTracer, a computational tool that uses multi-scale molecular networks to identify potential combination therapies from pre- and post-treatment metabolomics data. We first demonstrate piTracer’s core ability to reconstruct cellular cascades by inspecting well-characterized molecular pathways and previously studied associations between genetic variants and metabolite levels. We then apply a new gene ranking algorithm on differential metabolomic profiles from human breast cancer cells after glutaminase inhibition. Four of the automatically identified gene targets were experimentally tested by simultaneous inhibition of the respective targets and glutaminase. Of these combination treatments, two were be confirmed to induce synthetic lethality in the cell line. In summary, piTracer integrates the molecular monitoring of escape mechanisms into comprehensive pathway networks to accelerate drug target identification. The tool is open source and can be accessed at https://github.com/krumsieklab/pitracer .

Sex-biased molecular differences in lung adenocarcinoma are ethnic and smoking specific

BACKGROUND

Sex-related differences in cancer epidemiology, tumor biology, immune system activity, and pharmacogenomics have been suggested to be important considerations for precision cancer control. Here we elucidated systematically sex biases in genetic variants, gene expression profiles, and immunological landscapes of lung adenocarcinoma patients (LUADs) with different ancestry and smoking status.

METHODS

Somatic mutation and mRNA expression data of Asian and Non-Asian LUADs were obtained from public databases. Sex-biased genetic mutations, gene expression, biological pathways, and immune infiltration were identified in the context of smoking status and race.

RESULTS

Among nonsmokers, male-biased mutations were prevalent in Asian LUADs, while few sex-biased mutations were detected in Non-Asian LUADs. EGFR was the only mutation whose frequency was significantly higher in females than males in both Asian and Non-Asian nonsmokers. More genes exhibited sex-biased expression in Non-Asian LUADs compared to Asian LUADs. Moreover, genes distinctly expressed in females were mainly related to immune-related pathways, whereas those in males were more involved in activation of DNA repair, E2F_targets, and MYC_targets pathways. We also detected sex-specific immune infiltration in the context of genetic variation. In EGFR-mutant LUADs, males had a significantly increased infiltration of CD8 + T cells, whereas resting CD4 + memory T cells were more abundant in females. Additionally, in KRAS-mutant LUADs, CD8 + and CD4 + T cells were more abundant in females than males. In addition, we detected all female patients with high SCGB3A2 expression were exclusively sensitive to immunotherapy, while this phenomenon was not observed in male patients.

CONCLUSIONS

Our findings provided evidence that sex-related molecular and cellular components are involved in shaping tumor distinct genetic and immune features, which might have important impact on personalized targeted and immune therapy.

Comprehensive transcriptomic analysis to identify biological and clinical differences in cholangiocarcinoma

BACKGROUND

Cholangiocarcinoma (CC) is a rare and aggressive disease with limited therapeutic options and a poor prognosis. All available public records of cohorts reporting transcriptomic data on intrahepatic cholangiocarcinoma (ICC) and extrahepatic cholangiocarcinoma (ECC) were collected with the aim to provide a comprehensive gene expression-based classification with clinical relevance.

METHODS

A total of 543 patients with primary tumor tissues profiled by RNAseq and microarray platforms from seven public datasets were used as a discovery set to identify distinct biological subgroups. Group predictors developed on the discovery sets were applied to a single cohort of 131 patients profiled with RNAseq for validation and assessment of clinical relevance leveraging machine learning techniques.

RESULTS

By unsupervised clustering analysis of gene expression data we identified both in the ICC and ECC discovery datasets four subgroups characterized by a distinct type of immune infiltrate and signaling pathways. We next developed class predictors using short gene list signatures and identified in an independent dataset subgroups of ICC tumors at different prognosis.

CONCLUSIONS

The developed class-predictor allows identification of CC subgroups with specific biological features and clinical behavior at single-sample level. Such results represent the starting point for a complete molecular characterization of CC, including integration of genomics data to develop in clinical practice.

Biological knowledge graph-guided investigation of immune therapy response in cancer with graph neural network

The determination of transcriptome profiles that mediate immune therapy in cancer remains a major clinical and biological challenge. Despite responses induced by immune-check points inhibitors (ICIs) in diverse tumor types and all the big breakthroughs in cancer immunotherapy, most patients with solid tumors do not respond to ICI therapies. It still remains a big challenge to predict the ICI treatment response. Here, we propose a framework with multiple prior knowledge networks guided for immune checkpoints inhibitors prediction-DeepOmix-ICI (or ICInet for short). ICInet can predict the immune therapy response by leveraging geometric deep learning and prior biological knowledge graphs of gene-gene interactions. Here, we demonstrate more than 600 ICI-treated patients with ICI response data and gene expression profile to apply on ICInet. ICInet was used for ICI therapy responses prediciton across different cancer types-melanoma, gastric cancer and bladder cancer, which includes 7 cohorts from different data sources. ICInet is able to robustly generalize into multiple cancer types. Moreover, the performance of ICInet in those cancer types can outperform other ICI biomarkers in the clinic. Our model [area under the curve (AUC = 0.85)] generally outperformed other measures, including tumor mutational burden (AUC = 0.62) and programmed cell death ligand-1 score (AUC = 0.74). Therefore, our study presents a prior-knowledge guided deep learning method to effectively select immunotherapy-response-associated biomarkers, thereby improving the prediction of immunotherapy response for precision oncology.

Leveraging transcriptomics for precision diagnosis: Lessons learned from cancer and sepsis

Diagnostics require precision and predictive ability to be clinically useful. Integration of multi-omic with clinical data is crucial to our understanding of disease pathogenesis and diagnosis. However, interpretation of overwhelming amounts of information at the individual level requires sophisticated computational tools for extraction of clinically meaningful outputs. Moreover, evolution of technical and analytical methods often outpaces standardisation strategies. RNA is the most dynamic component of all -omics technologies carrying an abundance of regulatory information that is least harnessed for use in clinical diagnostics. Gene expression-based tests capture genetic and non-genetic heterogeneity and have been implemented in certain diseases. For example patients with early breast cancer are spared toxic unnecessary treatments with scores based on the expression of a set of genes (e.g., Oncotype DX). The ability of transcriptomics to portray the transcriptional status at a moment in time has also been used in diagnosis of dynamic diseases such as sepsis. Gene expression profiles identify endotypes in sepsis patients with prognostic value and a potential to discriminate between viral and bacterial infection. The application of transcriptomics for patient stratification in clinical environments and clinical trials thus holds promise. In this review, we discuss the current clinical application in the fields of cancer and infection. We use these paradigms to highlight the impediments in identifying useful diagnostic and prognostic biomarkers and propose approaches to overcome them and aid efforts towards clinical implementation.

Genomic and transcriptomic analyses identify a prognostic gene signature and predict response to therapy in pleural and peritoneal mesothelioma

Malignant mesothelioma is an aggressive cancer with limited treatment options and poor prognosis. A better understanding of mesothelioma genomics and transcriptomics could advance therapies. Here, we present a mesothelioma cohort of 122 patients along with their germline and tumor whole-exome and tumor RNA sequencing data as well as phenotypic and drug response information. We identify a 48-gene prognostic signature that is highly predictive of mesothelioma patient survival, including CCNB1, the expression of which is highly predictive of patient survival on its own. In addition, we analyze the transcriptomics data to study the tumor immune microenvironment and identify synthetic-lethality-based signatures predictive of response to therapy. This germline and somatic whole-exome sequencing as well as transcriptomics data from the same patient are a valuable resource to address important biological questions, including prognostic biomarkers and determinants of treatment response in mesothelioma.

Spatial multi-omics revealed the impact of tumor ecosystem heterogeneity on immunotherapy efficacy in patients with advanced non-small cell lung cancer treated with bispecific antibody

BACKGROUND

Immunotherapy for malignant tumors has made great progress, but many patients do not benefit from it. The complex intratumoral heterogeneity (ITH) hindered the in-depth exploration of immunotherapy. Conventional bulk sequencing has masked intratumor complexity, preventing a more detailed discovery of the impact of ITH on treatment efficacy. Hence, we initiated this study to explore ITH at the multi-omics spatial level and to seek prognostic biomarkers of immunotherapy efficacy considering the presence of ITH.

METHODS

Using the segmentation strategy of digital spatial profiling (DSP), we obtained differential information on tumor and stromal regions at the proteomic and transcriptomic levels. Based on the consideration of ITH, signatures constructed by candidate proteins in different regions were used to predict the efficacy of immunotherapy.

RESULTS

Eighteen patients treated with a bispecific antibody (bsAb)-KN046 were enrolled in this study. The tumor and stromal areas of the same samples exhibited distinct features. Signatures consisting of 11 and 18 differentially expressed DSP markers from the tumor and stromal areas, respectively, were associated with treatment response. Furthermore, the spatially resolved signature identified from the stromal areas showed greater predictive power for bsAb immunotherapy response (area under the curve=0.838). Subsequently, our stromal signature was validated in an independent cohort of patients with non-small cell lung cancer undergoing immunotherapy.

CONCLUSION

We deciphered ITH at the spatial level and demonstrated for the first time that genetic information in the stromal region can better predict the efficacy of bsAb treatment.

TRIAL REGISTRATION NUMBER

NCT03838848.

Clinically oriented prediction of patient response to targeted and immunotherapies from the tumor transcriptome

BACKGROUND

Precision oncology is gradually advancing into mainstream clinical practice, demonstrating significant survival benefits. However, eligibility and response rates remain limited in many cases, calling for better predictive biomarkers.

METHODS

We present ENLIGHT, a transcriptomics-based computational approach that identifies clinically relevant genetic interactions and uses them to predict a patient's response to a variety of therapies in multiple cancer types without training on previous treatment response data. We study ENLIGHT in two translationally oriented scenarios: personalized oncology (PO), aimed at prioritizing treatments for a single patient, and clinical trial design (CTD), selecting the most likely responders in a patient cohort.

FINDINGS

Evaluating ENLIGHT's performance on 21 blinded clinical trial datasets in the PO setting, we show that it can effectively predict a patient's treatment response across multiple therapies and cancer types. Its prediction accuracy is better than previously published transcriptomics-based signatures and is comparable with that of supervised predictors developed for specific indications and drugs. In combination with the interferon-γ signature, ENLIGHT achieves an odds ratio larger than 4 in predicting response to immune checkpoint therapy. In the CTD scenario, ENLIGHT can potentially enhance clinical trial success for immunotherapies and other monoclonal antibodies by excluding non-responders while overall achieving more than 90% of the response rate attainable under an optimal exclusion strategy.

CONCLUSIONS

ENLIGHT demonstrably enhances the ability to predict therapeutic response across multiple cancer types from the bulk tumor transcriptome.

FUNDING

This research was supported in part by the Intramural Research Program, NIH and by the Israeli Innovation Authority.

Oncosuppressive miRNAs loaded in lipid nanoparticles potentiate targeted therapies in BRAF-mutant melanoma by inhibiting core escape pathways of resistance

BRAF-mutated melanoma relapsing after targeted therapies is an aggressive disease with unmet clinical need. Hence the need to identify novel combination therapies able to overcome drug resistance. miRNAs have emerged as orchestrators of non-genetic mechanisms adopted by melanoma cells to challenge therapies. In this context we previously identified a subset of oncosuppressor miRNAs downregulated in drug-resistant melanomas. Here we demonstrate that lipid nanoparticles co-encapsulating two of them, miR-199-5p and miR-204-5p, inhibit tumor growth both in vitro and in vivo in combination with target therapy and block the development of drug resistance. Mechanistically they act by directly reducing melanoma cell growth and also indirectly by hampering the recruitment and reprogramming of pro-tumoral macrophages. Molecularly, we demonstrate that the effects on macrophages are mediated by the dysregulation of a newly identified miR-204-5p-miR-199b-5p/CCL5 axis. Finally, we unveiled that M2 macrophages programs are molecular signatures of resistance and predict response to therapy in patients. Overall, these findings have strong translational implications to propose new combination therapies making use of RNA therapeutics for metastatic melanoma patients.

Fibrolamellar carcinoma transcriptomic-based treatment prediction: complete response after nivolumab and ipilimumab

Fibrolamellar carcinoma (FLC) is a rare cancer of the liver that most commonly affects children and young adults. There is no clear standard of care for the disease, whose response to treatment seems to be very different from that of hepatocellular carcinoma. We present a case of FLC in a patient in her mid 30s that recurred and persisted despite resection and multiple lines of treatment. Following transcriptomic analysis, a combination of ipilimumab (anti-CTLA4) and nivolumab (anti-PD-1) led to complete remission, although common biomarkers for immune checkpoint blockade were all negative in this case. The patient is still in remission. Here, combined checkpoint blockade guided by novel transcriptomic analysis led to complete remission after failure of several lines of treatment.

ACAP1 Deficiency Predicts Inferior Immunotherapy Response in Solid Tumors

BACKGROUND

ACAP1 plays a key role in endocytic recycling, which is essential for the normal function of lymphocytes. However, the expression and function of ACAP1 in lymphocytes have rarely been studied.

METHODS

Large-scale genomic data, including multiple bulk RNA-sequencing datasets, single-cell sequencing datasets, and immunotherapy cohorts, were exploited to comprehensively characterize ACAP1 expression, regulation, and function. Gene set enrichment analysis (GSEA) was used to uncover the pathways associated with ACAP1 expression. Eight algorithms, including TIMER, CIBERSORT, CIBERSORT-ABS, QUANTISEQ, xCELL, MCPCOUNTER, EPIC, and TIDE, were applied to estimate the infiltrating level of immune cells. Western blotting, qPCR, and ChIP-PCR were used to validate the findings from bioinformatic analyses. A T-cell co-culture killing assay was used to investigate the function of ACAP1 in lymphocytes.

RESULTS

ACAP1 was highly expressed in immune-related tissues and cells and minimally in other tissues. Moreover, single-cell sequencing analysis in tumor samples revealed that ACAP1 is expressed primarily in tumor-infiltrating lymphocytes (TILs), including T, B, and NK cells. ACAP1 expression is negatively regulated by promoter DNA methylation, with its promoter hypo-methylated in immune cells but hyper-methylated in other cells. Furthermore, SPI1 binds to the ACAP1 promoter and positively regulates its expression in immune cells. ACAP1 levels positively correlate with the infiltrating levels of TILs, especially CD8+ T cells, across a broad range of solid cancer types. ACAP1 deficiency is associated with poor prognosis and immunotherapeutic response in multiple cancer types treated with checkpoint blockade therapy (ICT). Functionally, the depletion of ACAP1 by RNA interference significantly impairs the T cell-mediated killing of tumor cells.

CONCLUSIONS

Our study demonstrates that ACAP1 is essential for the normal function of TILs, and its deficiency indicates an immunologically "cold" status of tumors that are resistant to ICT.

Super enhancers: Pathogenic roles and potential therapeutic targets for acute myeloid leukemia (AML)

Acute myeloid leukemia (AML) is a malignant hematological tumor with disordered oncogenes/tumor suppressor genes and limited treatments. The potent anti-cancer effects of bromodomain and extra-terminal domain (BET) inhibitors, targeting the key component of super enhancers, in early clinical trials on AML patients, implies the critical role of super enhancers in AML. Here, we review the concept and characteristic of super enhancer, and then summarize the current researches about super enhancers in AML pathogenesis, diagnosis and classification, followed by illustrate the potential super enhancer-related targets and drugs, and propose the future directions of super enhancers in AML. This information provides integrated insight into the roles of super enhancers in this disease.

Big data in basic and translational cancer research

Historically, the primary focus of cancer research has been molecular and clinical studies of a few essential pathways and genes. Recent years have seen the rapid accumulation of large-scale cancer omics data catalysed by breakthroughs in high-throughput technologies. This fast data growth has given rise to an evolving concept of 'big data' in cancer, whose analysis demands large computational resources and can potentially bring novel insights into essential questions. Indeed, the combination of big data, bioinformatics and artificial intelligence has led to notable advances in our basic understanding of cancer biology and to translational advancements. Further advances will require a concerted effort among data scientists, clinicians, biologists and policymakers. Here, we review the current state of the art and future challenges for harnessing big data to advance cancer research and treatment.

Molecular pathways enhance drug response prediction using transfer learning from cell lines to tumors and patient-derived xenografts

Computational models have been successful in predicting drug sensitivity in cancer cell line data, creating an opportunity to guide precision medicine. However, translating these models to tumors remains challenging. We propose a new transfer learning workflow that transfers drug sensitivity predicting models from large-scale cancer cell lines to both tumors and patient derived xenografts based on molecular pathways derived from genomic features. We further compute feature importance to identify pathways most important to drug response prediction. We obtained good performance on tumors (AUROC = 0.77) and patient derived xenografts from triple negative breast cancers (RMSE = 0.11). Using feature importance, we highlight the association between ER-Golgi trafficking pathway in everolimus sensitivity within breast cancer patients and the role of class II histone deacetylases and interlukine-12 in response to drugs for triple-negative breast cancer. Pathway information support transfer of drug response prediction models from cell lines to tumors and can provide biological interpretation underlying the predictions, serving as a steppingstone towards usage in clinical setting.

Transcriptomic Harmonization as the Way for Suppressing Cross-Platform Bias and Batch Effect

(1) Background: Emergence of methods interrogating gene expression at high throughput gave birth to quantitative transcriptomics, but also posed a question of inter-comparison of expression profiles obtained using different equipment and protocols and/or in different series of experiments. Addressing this issue is challenging, because all of the above variables can dramatically influence gene expression signals and, therefore, cause a plethora of peculiar features in the transcriptomic profiles. Millions of transcriptomic profiles were obtained and deposited in public databases of which the usefulness is however strongly limited due to the inter-comparison issues; (2) Methods: Dozens of methods and software packages that can be generally classified as either flexible or predefined format harmonizers have been proposed, but none has become to the date the gold standard for unification of this type of Big Data; (3) Results: However, recent developments evidence that platform/protocol/batch bias can be efficiently reduced not only for the comparisons of limited transcriptomic datasets. Instead, instruments were proposed for transforming gene expression profiles into the universal, uniformly shaped format that can support multiple inter-comparisons for reasonable calculation costs. This forms a basement for universal indexing of all or most of all types of RNA sequencing and microarray hybridization profiles; (4) Conclusions: In this paper, we attempted to overview the landscape of modern approaches and methods in transcriptomic harmonization and focused on the practical aspects of their application.

Artificial intelligence-based prediction of clinical outcome in immunotherapy and targeted therapy of lung cancer

Lung cancer accounts for the main proportion of malignancy-related deaths and most patients are diagnosed at an advanced stage. Immunotherapy and targeted therapy have great advances in application in clinics to treat lung cancer patients, yet the efficacy is unstable. The response rate of these therapies varies among patients. Some biomarkers have been proposed to predict the outcomes of immunotherapy and targeted therapy, including programmed cell death-ligand 1 (PD-L1) expression and oncogene mutations. Nevertheless, the detection tests are invasive, time-consuming, and have high demands on tumor tissue. The predictive performance of conventional biomarkers is also unsatisfactory. Therefore, novel biomarkers are needed to effectively predict the outcomes of immunotherapy and targeted therapy. The application of artificial intelligence (AI) can be a possible solution, as it has several advantages. AI can help identify features that are unable to be used by humans and perform repetitive tasks. By combining AI methods with radiomics, pathology, genomics, transcriptomics, proteomics, and clinical data, the integrated model has shown predictive value in immunotherapy and targeted therapy, which significantly improves the precision treatment of lung cancer patients. Herein, we reviewed the application of AI in predicting the outcomes of immunotherapy and targeted therapy in lung cancer patients, and discussed the challenges and future directions in this field.

Exploring synthetic lethal network for the precision treatment of clear cell renal cell carcinoma

The emerging targeted therapies have revolutionized the treatment of advanced clear cell renal cell carcinoma (ccRCC) over the past 15 years. Nevertheless, lack of personalized treatment limits the development of effective clinical guidelines and improvement of patient prognosis. In this study, large-scale genomic profiles from ccRCC cohorts were explored for integrative analysis. A credible method was developed to identify synthetic lethality (SL) pairs and a list of 72 candidate pairs was determined, which might be utilized to selectively eliminate tumors with genetic aberrations using SL partners of specific mutations. Further analysis identified BRD4 and PRKDC as novel medical targets for patients with BAP1 mutations. After mapping these target genes to the comprehensive drug datasets, two agents (BI-2536 and PI-103) were found to have considerable therapeutic potentials in the BAP1 mutant tumors. Overall, our findings provided insight into the overview of ccRCC mutation patterns and offered novel opportunities for improving individualized cancer treatment.

Single-cell transcriptomic analysis reveals circadian rhythm disruption associated with poor prognosis and drug-resistance in lung adenocarcinoma

Circadian rhythm disruption (CRD) represents a major contributor to tumor proliferation. Nonetheless, the role of CRD in the clinical prediction of cancer outcomes has not been well studied. In this study, we developed a computational algorithm, which was implemented in an open-source R package CRDscore, to define the intratumoral status of circadian disruption in three representative single-cell RNA-seq data sets of lung adenocarcinoma. We found that the malignant cells with high CRDscore were characterized by activation of glycolysis and epithelial-mesenchymal transition pathways. Furthermore, cell communication analysis indicated that CRD played a pivotal role in T cell exhaustion, which may be responsible for the poor prognosis of the malignancy. We then validated the findings with public bulk transcriptome datasets involving 22 cancer types. Cox regression analysis revealed that the CRDscore was a valuable prognostic biomarker. A model containing 23 circadian-related genes performed well in predicting immunotherapeutic outcomes in 14 independent cohorts. Importantly, decreased CRDscore was detect by RNA sequencing on H1299 cells with melatonin treatment. Meanwhile, the cells downregulated the expression level of SNAIL and TWIST, which contributed to an invasive phenotype. In conclusion, this study provides a novel computational framework for characterizing CRD status using single-cell transcriptomic data and further confirmed the molecular mechanisms underlying metabolic reprogramming and T cell exhaustion under CRD. The better understanding of the mechanisms may provide new possibilities for incorporating "anticancer approaches based on circadian clocks" into the treatment protocols of precision medicine.

Functional Drug Screening in the Era of Precision Medicine

The focus of precision medicine is providing the right treatment to each unique patient. This scientific movement has incited monumental advances in oncology including the approval of effective, targeted agnostic therapies. Yet, precision oncology has focused largely on genomics in the treatment decision making process, and several recent clinical trials demonstrate that genomics is not the only variable to be considered. Drug screening in three dimensional (3D) models, including patient derived organoids, organs on a chip, xenografts, and 3D-bioprinted models provide a functional medicine perspective and necessary complement to genomic testing. In this review, we discuss the practicality of various 3D drug screening models and each model's ability to capture the patient's tumor microenvironment. We highlight the potential for enhancing precision medicine that personalized functional drug testing holds in combination with genomic testing and emerging mathematical models.

Network-based machine learning approach to predict immunotherapy response in cancer patients

Immune checkpoint inhibitors (ICIs) have substantially improved the survival of cancer patients over the past several years. However, only a minority of patients respond to ICI treatment (~30% in solid tumors), and current ICI-response-associated biomarkers often fail to predict the ICI treatment response. Here, we present a machine learning (ML) framework that leverages network-based analyses to identify ICI treatment biomarkers (NetBio) that can make robust predictions. We curate more than 700 ICI-treated patient samples with clinical outcomes and transcriptomic data, and observe that NetBio-based predictions accurately predict ICI treatment responses in three different cancer types—melanoma, gastric cancer, and bladder cancer. Moreover, the NetBio-based prediction is superior to predictions based on other conventional ICI treatment biomarkers, such as ICI targets or tumor microenvironment-associated markers. This work presents a network-based method to effectively select immunotherapy-response-associated biomarkers that can make robust ML-based predictions for precision oncology.

Identifying biomarkers for response to immunotherapy in cancer remains challenging. Here, the authors develop an approach based on network biology and machine learning -NetBio- to identify molecular biomarkers of response to immunotherapy across different cancer types and cohorts.

Phase II clinical trial of Eribulin-gemcitabine combination therapy in previously treated patients with advanced liposarcoma or leiomyosarcoma

PURPOSE

Monotherapy with eribulin or gemcitabine has been found to be moderately effective in treating soft-tissue sarcomas (STS). In this study, we evaluated the efficacy and safety of eribulin-gemcitabine combination therapy for the two most common histologic types of STS, liposarcoma and leiomyosarcoma.

PATIENTS AND METHODS

In this non-randomized, multicenter, phase II study, we included patients with progressive disease who had received one or two courses of chemotherapy that included doxorubicin. Patients were administered 1.4 mg/m2 eribulin and 1,000 mg/m2 gemcitabine on days 1 and 8 every 3 weeks. The primary endpoint was progression-free survival rate at 12 weeks (PFSR12wks), with null and alternative hypotheses of PFSR12wks {less than or equal to}20.0% and {greater than or equal to}40.0%, respectively. Exploratory biomarker analyses with next-generation sequencing (NGS) were performed on pretreatment tumor samples.

RESULTS

Among the 37 patients included, the overall PFSR12wks was 73.0%, achieving the primary endpoint. The objective response rate, disease control rate, median progression-free survival, and median overall survival were 16.2%, 78.4%, 5.6 months, and 31.9 months, respectively, without differences according to histologic type. New safety signals and treatment-related deaths were not documented. NGS-based transcriptome analysis revealed that functional enrichment in the transforming growth factor-ß pathway was mostly associated with a poor outcome, whereas single genetic alterations largely failed to predict treatment outcome.

CONCLUSIONS

Eribulin-gemcitabine combination therapy showed promising activity and an acceptable safety profile in patients with liposarcoma or leiomyosarcoma. Gene expression profiling with pathway enrichment analysis would have possibilities to have predictive value for survival outcome, necessitating further investigation to confirm.

Precision combination therapies based on recurrent oncogenic co-alterations

Cancer cells depend on multiple driver alterations whose oncogenic effects can be suppressed by drug combinations. Here, we provide a comprehensive resource of precision combination therapies tailored to oncogenic co-alterations that are recurrent across patient cohorts. To generate the resource, we developed Recurrent Features Leveraged for Combination Therapy (REFLECT), which integrates machine learning and cancer informatics algorithms. Using multi-omic data, the method maps recurrent co-alteration signatures in patient cohorts to combination therapies. We validated the REFLECT pipeline using data from patient-derived xenografts, in vitro drug screens, and a combination therapy clinical trial. These validations demonstrate that REFLECT-selected combination therapies have significantly improved efficacy, synergy, and survival outcomes. In patient cohorts with immunotherapy response markers, DNA repair aberrations, and HER2 activation, we have identified therapeutically actionable and recurrent co-alteration signatures. REFLECT provides a resource and framework to design combination therapies tailored to tumor cohorts in data-driven clinical trials and pre-clinical studies.

Computational methods, databases and tools for synthetic lethality prediction

Synthetic lethality (SL) occurs between two genes when the inactivation of either gene alone has no effect on cell survival but the inactivation of both genes results in cell death. SL-based therapy has become one of the most promising targeted cancer therapies in the last decade as PARP inhibitors achieve great success in the clinic. The key point to exploiting SL-based cancer therapy is the identification of robust SL pairs. Although many wet-lab-based methods have been developed to screen SL pairs, known SL pairs are less than 0.1% of all potential pairs due to large number of human gene combinations. Computational prediction methods complement wet-lab-based methods to effectively reduce the search space of SL pairs. In this paper, we review the recent applications of computational methods and commonly used databases for SL prediction. First, we introduce the concept of SL and its screening methods. Second, various SL-related data resources are summarized. Then, computational methods including statistical-based methods, network-based methods, classical machine learning methods and deep learning methods for SL prediction are summarized. In particular, we elaborate on the negative sampling methods applied in these models. Next, representative tools for SL prediction are introduced. Finally, the challenges and future work for SL prediction are discussed.