Predicting cancer-specific vulnerability via data-driven detection of synthetic lethality

Synergistic epistasis among cancer drivers can rescue early tumors from the accumulation of deleterious passengers

Epistasis among driver mutations is pervasive and explains relevant features of cancer, such as differential therapy response and convergence towards well-characterized molecular subtypes. Furthermore, a growing body of evidence suggests that tumor development could be hampered by the accumulation of slightly deleterious passenger mutations. In this work, we combined empirical epistasis networks, computer simulations, and mathematical models to explore how synergistic interactions among driver mutations affect cancer progression under the burden of slightly deleterious passengers. We found that epistasis plays a crucial role in tumor development by promoting the transformation of precancerous clones into rapidly growing tumors through a process that is analogous to evolutionary rescue. The triggering of epistasis-driven rescue is strongly dependent on the intensity of epistasis and could be a key rate-limiting step in many tumors, contributing to their unpredictability. As a result, central genes in cancer epistasis networks appear as key intervention targets for cancer therapy.

CSSLdb: Discovery of cancer-specific synthetic lethal interactions based on machine learning and statistic inference

Synthetic lethality (SL) occurs when the inactivation of two genes results in cell death while the inactivation of either gene alone is non-lethal. SL-based therapy has become a promising anti-cancer treatment option with the increasing researches and applications in clinical practice, while the specific therapeutic opportunities for various cancers have not yet been comprehensively investigated. Herein, we described a computational approach based on machine learning and statistical inference to discover the cancer-specific synthetic lethal interactions. First, Random Forest and One-Class SVM were used to perform cancer unbiased prediction of synthetic lethality. Then, two strategies, including mutual exclusivity and differential expression, were used to screen cancer-specific synthetic lethal interactions, resulting in 14,582 SL gene pairs in 33 cancer types. Finally, we developed a freely available database of CSSLdb (Cancer Specific Synthetic Lethality Database, http://www.tmliang.cn/CSSL/) to present cancer-specific synthetic lethal genetic interactions, which would enrich the relevant research and contribute to underlying therapy strategies based on synthetic lethality.

SL-Miner: a web server for mining evidence and prioritization of cancer-specific synthetic lethality

SUMMARY

Synthetic lethality (SL) refers to a type of genetic interaction in which the simultaneous inactivation of two genes leads to cell death, while the inactivation of a single gene does not affect cell viability. It significantly expands the range of potential therapeutic targets for anti-cancer treatments. SL interactions are primarily identified through experimental screening and computational prediction. Although various computational methods have been proposed, they tend to ignore providing evidence to support their predictions of SL. Besides, they are rarely user-friendly for biologists who likely have limited programming skills. Moreover, the genetic context specificity of SL interactions is often not taken into consideration. Here, we introduce a web server called SL-Miner, which is designed to mine the evidence of SL relationships between a primary gene and a few candidate SL partner genes in a specific type of cancer, and to prioritize these candidate genes by integrating various types of evidence. For intuitive data visualization, SL-Miner provides a range of charts (e.g. volcano plot and box plot) to help users get insights from the data.

AVAILABILITY AND IMPLEMENTATION

SL-Miner is available at https://slminer.sist.shanghaitech.edu.cn.

Leveraging synthetic lethality to uncover potential therapeutic target in gastric cancer

Since trastuzumab was approved in 2012 for the first-line treatment of gastric cancer (GC), no significant advancement in GC targeted therapies has occurred. Synthetic lethality refers to the concept that simultaneous dysfunction of a pair of genes results in a lethal effect on cells, while the loss of an individual gene does not cause this effect. Through exploiting synthetic lethality, novel targeted therapies can be developed for the individualized treatment of GC. In this study, we proposed a computational strategy named Gastric cancer Specific Synthetic Lethality inference (GSSL) to identify synthetic lethal interactions in GC. GSSL analysis was used to infer probable synthetic lethality in GC using four accessible clinical datasets. In addition, prediction results were confirmed by experiments. GSSL analysis identified a total of 34 candidate synthetic lethal pairs, which included 33 unique targets. Among the synthetic lethal gene pairs, TP53-CHEK1 was selected for further experimental validation. Both computational and experimental results indicated that inhibiting CHEK1 could be a potential therapeutic strategy for GC patients with TP53 mutation. Meanwhile, in vitro experimental validation of two novel synthetic lethal pairs TP53-AURKB and ARID1A-EP300 further proved the universality and reliability of GSSL. Collectively, GSSL has been shown to be a reliable and feasible method for comprehensive analysis of inferring synthetic lethal interactions of GC, which may offer novel insight into the precision medicine and individualized treatment of GC.

SLKB: synthetic lethality knowledge base

Emerging CRISPR-Cas9 technology permits synthetic lethality (SL) screening of large number of gene pairs from gene combination double knockout (CDKO) experiments. However, the poor integration and annotation of CDKO SL data in current SL databases limit their utility, and diverse methods of calculating SL scores prohibit their comparison. To overcome these shortcomings, we have developed SL knowledge base (SLKB) that incorporates data of 11 CDKO experiments in 22 cell lines, 16,059 SL gene pairs and 264,424 non-SL gene pairs. Additionally, within SLKB, we have implemented five SL calculation methods: median score with and without background control normalization (Median-B/NB), sgRNA-derived score (sgRNA-B/NB), Horlbeck score, GEMINI score and MAGeCK score. The five scores have demonstrated a mere 1.21% overlap among their top 10% SL gene pairs, reflecting high diversity. Users can browse SL networks and assess the impact of scoring methods using Venn diagrams. The SL network generated from all data in SLKB shows a greater likelihood of SL gene pair connectivity with other SL gene pairs than non-SL pairs. Comparison of SL networks between two cell lines demonstrated greater likelihood to share SL hub genes than SL gene pairs. SLKB website and pipeline can be freely accessed at https://slkb.osubmi.org and https://slkb.docs.osubmi.org/, respectively.

ELISL: early-late integrated synthetic lethality prediction in cancer

MOTIVATION

Anti-cancer therapies based on synthetic lethality (SL) exploit tumour vulnerabilities for treatment with reduced side effects, by targeting a gene that is jointly essential with another whose function is lost. Computational prediction is key to expedite SL screening, yet existing methods are vulnerable to prevalent selection bias in SL data and reliant on cancer or tissue type-specific omics, which can be scarce. Notably, sequence similarity remains underexplored as a proxy for related gene function and joint essentiality.

RESULTS

We propose ELISL, Early-Late Integrated SL prediction with forest ensembles, using context-free protein sequence embeddings and context-specific omics from cell lines and tissue. Across eight cancer types, ELISL showed superior robustness to selection bias and recovery of known SL genes, as well as promising cross-cancer predictions. Co-occurring mutations in a BRCA gene and ELISL-predicted pairs from the HH, FGF, WNT, or NEIL gene families were associated with longer patient survival times, revealing therapeutic potential.

AVAILABILITY AND IMPLEMENTATION

Data: 10.6084/m9.figshare.23607558 & Code: github.com/joanagoncalveslab/ELISL.

Pathway-driven analysis of synthetic lethal interactions in cancer using perturbation screens

Synthetic lethality offers a promising approach for developing effective therapeutic interventions in cancer when direct targeting of driver genes is impractical. In this study, we comprehensively analyzed large-scale CRISPR, shRNA, and PRISM screens to identify potential synthetic lethal (SL) interactions in pan-cancer and 12 individual cancer types, using a new computational framework that leverages the biological function and signaling pathway information of key driver genes to mitigate the confounding effects of background genetic alterations in different cancer cell lines. This approach has successfully identified several putative SL interactions, including KRAS-MAP3K2 and APC-TCF7L2 in pan cancer, and CCND1-METTL1, TP53-FRS3, SMO-MDM2, and CCNE1-MTOR in liver, blood, skin, and gastric cancers, respectively. In addition, we proposed several FDA-approved cancer-targeted drugs for various cancer types through PRISM drug screens, such as cabazitaxel for VHL-mutated kidney cancer and alectinib for lung cancer with NRAS or KRAS mutations. Leveraging pathway information can enhance the concordance of shRNA and CRISPR screens and provide clinically relevant findings such as the potential efficacy of dasatinib, an inhibitor of SRC, for colorectal cancer patients with mutations in the WNT signaling pathway. These analyses revealed that taking signaling pathway information into account results in the identification of more promising SL interactions.

Synthetic Lethality Screening with Recursive Feature Machines

Synthetic lethality refers to a genetic interaction where the simultaneous perturbation of gene pairs leads to cell death. Synthetically lethal gene pairs (SL pairs) provide a potential avenue for selectively targeting cancer cells based on genetic vulnerabilities. The rise of large-scale gene perturbation screens such as the Cancer Dependency Map (DepMap) offers the opportunity to identify SL pairs automatically using machine learning. We build on a recently developed class of feature learning kernel machines known as Recursive Feature Machines (RFMs) to develop a pipeline for identifying SL pairs based on CRISPR viability data from DepMap. In particular, we first train RFMs to predict viability scores for a given CRISPR gene knockout from cell line embeddings consisting of gene expression and mutation features. After training, RFMs use a statistical operator known as average gradient outer product to provide weights for each feature indicating the importance of each feature in predicting cellular viability. We subsequently apply correlation-based filters to re-weight RFM feature importances and identify those features that are most indicative of low cellular viability. Our resulting pipeline is computationally efficient, taking under 3 minutes for analyzing all 17, 453 knockouts from DepMap for candidate SL pairs. We show that our pipeline more accurately recovers experimentally verified SL pairs than prior approaches. Moreover, our pipeline finds new candidate SL pairs, thereby opening novel avenues for identifying genetic vulnerabilities in cancer.

Using graph-based model to identify cell specific synthetic lethal effects

Synthetic lethal (SL) pairs are pairs of genes whose simultaneous loss-of-function results in cell death, while a damaging mutation of either gene alone does not affect the cell's survival. This makes SL pairs attractive targets for precision cancer therapies, as targeting the unimpaired gene of the SL pair can selectively kill cancer cells that already harbor the impaired gene. Limited by the difficulty of finding true SL pairs, especially on specific cell types, current computational approaches provide only limited insights because of overlooking the crucial aspects of cellular context dependency and mechanistic understanding of SL pairs. As a result, the identification of SL targets still relies on expensive, time-consuming experimental approaches. In this work, we applied cell-line specific multi-omics data to a specially designed deep learning model to predict cell-line specific SL pairs. Through incorporating multiple types of cell-specific omics data with a self-attention module, we represent gene relationships as graphs. Our approach achieves the prediction of SL pairs in a cell-specific manner and demonstrates the potential to facilitate the discovery of cell-specific SL targets for cancer therapeutics, providing a tool to unearth mechanisms underlying the origin of SL in cancer biology. The code and data of our approach can be found at https://github.com/promethiume/SLwise.

The first multi-tissue genome-scale metabolic model of a woody plant highlights suberin biosynthesis pathways in Quercus suber

Over the last decade, genome-scale metabolic models have been increasingly used to study plant metabolic behaviour at the tissue and multi-tissue level under different environmental conditions. Quercus suber, also known as the cork oak tree, is one of the most important forest communities of the Mediterranean/Iberian region. In this work, we present the genome-scale metabolic model of the Q. suber (iEC7871). The metabolic model comprises 7871 genes, 6231 reactions, and 6481 metabolites across eight compartments. Transcriptomics data was integrated into the model to obtain tissue-specific models for the leaf, inner bark, and phellogen, with specific biomass compositions. The tissue-specific models were merged into a diel multi-tissue metabolic model to predict interactions among the three tissues at the light and dark phases. The metabolic models were also used to analyse the pathways associated with the synthesis of suberin monomers, namely the acyl-lipids, phenylpropanoids, isoprenoids, and flavonoids production. The models developed in this work provide a systematic overview of the metabolism of Q. suber, including its secondary metabolism pathways and cork formation.

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.

Synthetic lethality in large-scale integrated metabolic and regulatory network models of human cells

Synthetic lethality (SL) is a promising concept in cancer research. A wide array of computational tools has been developed to predict and exploit synthetic lethality for the identification of tumour-specific vulnerabilities. Previously, we introduced the concept of genetic Minimal Cut Sets (gMCSs), a theoretical approach to SL developed for genome-scale metabolic networks. The major challenge in our gMCS framework is to go beyond metabolic networks and extend existing algorithms to more complex protein-protein interactions. In this article, we take a step further and incorporate linear regulatory pathways into our gMCS approach. Extensive algorithmic modifications to compute gMCSs in integrated metabolic and regulatory models are presented in detail. Our extended approach is applied to calculate gMCSs in integrated models of human cells. In particular, we integrate the most recent genome-scale metabolic network, Human1, with 3 different regulatory network databases: Omnipath, Dorothea and TRRUST. Based on the computed gMCSs and transcriptomic data, we discovered new essential genes and their associated synthetic lethal for different cancer cell lines. The performance of the different integrated models is assessed with available large-scale in-vitro gene silencing data. Finally, we discuss the most relevant gene essentiality predictions based on published literature in cancer research.

B-cell-specific checkpoint molecules that regulate anti-tumour immunity

The role of B cells in anti-tumour immunity is still debated and, accordingly, immunotherapies have focused on targeting T and natural killer cells to inhibit tumour growth1,2. Here, using high-throughput flow cytometry as well as bulk and single-cell RNA-sequencing and B-cell-receptor-sequencing analysis of B cells temporally during B16F10 melanoma growth, we identified a subset of B cells that expands specifically in the draining lymph node over time in tumour-bearing mice. The expanding B cell subset expresses the cell surface molecule T cell immunoglobulin and mucin domain 1 (TIM-1, encoded by Havcr1) and a unique transcriptional signature, including multiple co-inhibitory molecules such as PD-1, TIM-3, TIGIT and LAG-3. Although conditional deletion of these co-inhibitory molecules on B cells had little or no effect on tumour burden, selective deletion of Havcr1 in B cells both substantially inhibited tumour growth and enhanced effector T cell responses. Loss of TIM-1 enhanced the type 1 interferon response in B cells, which augmented B cell activation and increased antigen presentation and co-stimulation, resulting in increased expansion of tumour-specific effector T cells. Our results demonstrate that manipulation of TIM-1-expressing B cells enables engagement of the second arm of adaptive immunity to promote anti-tumour immunity and inhibit tumour growth.

KR4SL: knowledge graph reasoning for explainable prediction of synthetic lethality

MOTIVATION

Synthetic lethality (SL) is a promising strategy for anticancer therapy, as inhibiting SL partners of genes with cancer-specific mutations can selectively kill the cancer cells without harming the normal cells. Wet-lab techniques for SL screening have issues like high cost and off-target effects. Computational methods can help address these issues. Previous machine learning methods leverage known SL pairs, and the use of knowledge graphs (KGs) can significantly enhance the prediction performance. However, the subgraph structures of KG have not been fully explored. Besides, most machine learning methods lack interpretability, which is an obstacle for wide applications of machine learning to SL identification.

RESULTS

We present a model named KR4SL to predict SL partners for a given primary gene. It captures the structural semantics of a KG by efficiently constructing and learning from relational digraphs in the KG. To encode the semantic information of the relational digraphs, we fuse textual semantics of entities into propagated messages and enhance the sequential semantics of paths using a recurrent neural network. Moreover, we design an attentive aggregator to identify critical subgraph structures that contribute the most to the SL prediction as explanations. Extensive experiments under different settings show that KR4SL significantly outperforms all the baselines. The explanatory subgraphs for the predicted gene pairs can unveil prediction process and mechanisms underlying synthetic lethality. The improved predictive power and interpretability indicate that deep learning is practically useful for SL-based cancer drug target discovery.

AVAILABILITY AND IMPLEMENTATION

The source code is freely available at https://github.com/JieZheng-ShanghaiTech/KR4SL.

Integration of cancer-related genetic landscape of Eph receptors and ephrins with proteomics identifies a crosstalk between EPHB6 and EGFR

Eph receptors and their ephrin ligands are viewed as promising targets for cancer treatment; however, targeting them is hindered by their context-dependent functionalities. To circumvent this, we explore molecular landscapes underlying their pro- and anti-malignant activities. Using unbiased bioinformatics approaches, we construct a cancer-related network of genetic interactions (GIs) of all Ephs and ephrins to assist in their therapeutic manipulation. We also apply genetic screening and BioID proteomics and integrate them with machine learning approaches to select the most relevant GIs of one Eph receptor, EPHB6. This identifies a crosstalk between EPHB6 and EGFR, and further experiments confirm the ability of EPHB6 to modulate EGFR signaling, enhancing the proliferation of cancer cells and tumor development. Taken together, our observations show EPHB6 involvement in EGFR action, suggesting its targeting might be beneficial in EGFR-dependent tumors, and confirm that the Eph family genetic interactome presented here can be effectively exploited in developing cancer treatment approaches.

Synthetic lethality prediction in DNA damage repair, chromatin remodeling and the cell cycle using multi-omics data from cell lines and patients

Discovering synthetic lethal (SL) gene partners of cancer genes is an important step in developing cancer therapies. However, identification of SL interactions is challenging, due to a large number of possible gene pairs, inherent noise and confounding factors in the observed signal. To discover robust SL interactions, we devised SLIDE-VIP, a novel framework combining eight statistical tests, including a new patient data-based test iSurvLRT. SLIDE-VIP leverages multi-omics data from four different sources: gene inactivation cell line screens, cancer patient data, drug screens and gene pathways. We applied SLIDE-VIP to discover SL interactions between genes involved in DNA damage repair, chromatin remodeling and cell cycle, and their potentially druggable partners. The top 883 ranking SL candidates had strong evidence in cell line and patient data, 250-fold reducing the initial space of 200K pairs. Drug screen and pathway tests provided additional corroboration and insights into these interactions. We rediscovered well-known SL pairs such as RB1 and E2F3 or PRKDC and ATM, and in addition, proposed strong novel SL candidates such as PTEN and PIK3CB. In summary, SLIDE-VIP opens the door to the discovery of SL interactions with clinical potential. All analysis and visualizations are available via the online SLIDE-VIP WebApp.

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 .

Machine learning based combination of multi-omics data for subgroup identification in non-small cell lung cancer

Non-small Cell Lung Cancer (NSCLC) is a heterogeneous disease with a poor prognosis. Identifying novel subtypes in cancer can help classify patients with similar molecular and clinical phenotypes. This work proposes an end-to-end pipeline for subgroup identification in NSCLC. Here, we used a machine learning (ML) based approach to compress the multi-omics NSCLC data to a lower dimensional space. This data is subjected to consensus K-means clustering to identify the five novel clusters (C1-C5). Survival analysis of the resulting clusters revealed a significant difference in the overall survival of clusters (p-value: 0.019). Each cluster was then molecularly characterized to identify specific molecular characteristics. We found that cluster C3 showed minimal genetic aberration with a high prognosis. Next, classification models were developed using data from each omic level to predict the subgroup of unseen patients. Decision‑level fused classification models were then built using these classifiers, which were used to classify unseen patients into five novel clusters. We also showed that the multi-omics-based classification model outperformed single-omic-based models, and the combination of classifiers proved to be a more accurate prediction model than the individual classifiers. In summary, we have used ML models to develop a classification method and identified five novel NSCLC clusters with different genetic and clinical characteristics.

Multi-omics characterization of synthetic lethality-related molecular features: implications for SL-based therapeutic target screening

Synthetic lethality (SL) represents the co-occurrence of two or more non-lethal disordered genes that could lead to cell death. SL-based anticancer therapeutics could specifically kill the cancer cells carrying the targeted mutated gene while leaving normal cells alive. Recent large-scale computational and experimental screenings provide rich resources of SL information while lacking systematic research on molecular features of SL genes. Combined with comprehensive multi-omics data analysis and experimental validation of one SL gene pair, Guo et al. portrayed a systematic layout of cancer-specific SL interactions that could improve understanding of carcinogenesis and potentially assist the subsequent screening of anticancer therapeutic targets.

Systematic analysis of cancer-specific synthetic lethal interactions provides insight into personalized anticancer therapy

The concept of synthetic lethality has great potential for anticancer therapy as a new strategy to specifically kill cancer cells while sparing normal cells. To further understand the potential molecular interactions and gene characteristics involved in synthetic lethality, we performed a comprehensive analysis of predicted cancer-specific genetic interactions. Many genes were identified as cancer-associated genes that contributed to multiple biological processes and pathways, and the gene features were not random, indicating their potential roles in human carcinogenesis. Some relevant genes detected in multiple cancers were prone to be enriched in specific biological progresses and pathways, especially processes associated with DNA damage, chromosome-related functions and cancer pathways. These findings strongly implicated potential roles for these genes in cancer pathophysiology and functional relationships, as well as applications for future anticancer drug discovery. Further experimental validation indicated that the synthetic lethal interaction of APC and GFER may provide a potential anticancer strategy for patients with APC-mutant colon cancer. These results will contribute to further exploration of synthetic lethal interactions and broader application of the concept of synthetic lethality in anticancer therapeutics.

SLGNN: synthetic lethality prediction in human cancers based on factor-aware knowledge graph neural network

MOTIVATION

Synthetic lethality (SL) is a form of genetic interaction that can selectively kill cancer cells without damaging normal cells. Exploiting this mechanism is gaining popularity in the field of targeted cancer therapy and anticancer drug development. Due to the limitations of identifying SL interactions from laboratory experiments, an increasing number of research groups are devising computational prediction methods to guide the discovery of potential SL pairs. Although existing methods have attempted to capture the underlying mechanisms of SL interactions, methods that have a deeper understanding of and attempt to explain SL mechanisms still need to be developed.

RESULTS

In this work, we propose a novel SL prediction method, SLGNN. This method is based on the following assumption: SL interactions are caused by different molecular events or biological processes, which we define as SL-related factors that lead to SL interactions. SLGNN, apart from identifying SL interaction pairs, also models the preferences of genes for different SL-related factors, making the results more interpretable for biologists and clinicians. SLGNN consists of three steps: first, we model the combinations of relationships in the gene-related knowledge graph as the SL-related factors. Next, we derive initial embeddings of genes through an explicit message aggregation process of the knowledge graph. Finally, we derive the final gene embeddings through an SL graph, constructed using known SL gene pairs, utilizing factor-based message aggregation. At this stage, a supervised end-to-end training model is used for SL interaction prediction. Based on experimental results, the proposed SLGNN model outperforms all current state-of-the-art SL prediction methods and provides better interpretability.

AVAILABILITY AND IMPLEMENTATION

SLGNN is freely available at https://github.com/zy972014452/SLGNN.

Multi-view graph convolutional network for cancer cell-specific synthetic lethality prediction

Synthetic lethal (SL) genetic interactions have been regarded as a promising focus for investigating potential targeted therapeutics to tackle cancer. However, the costly investment of time and labor associated with wet-lab experimental screenings to discover potential SL relationships motivates the development of computational methods. Although graph neural network (GNN) models have performed well in the prediction of SL gene pairs, existing GNN-based models are not designed for predicting cancer cell-specific SL interactions that are more relevant to experimental validation in vitro. Besides, neither have existing methods fully utilized diverse graph representations of biological features to improve prediction performance. In this work, we propose MVGCN-iSL, a novel multi-view graph convolutional network (GCN) model to predict cancer cell-specific SL gene pairs, by incorporating five biological graph features and multi-omics data. Max pooling operation is applied to integrate five graph-specific representations obtained from GCN models. Afterwards, a deep neural network (DNN) model serves as the prediction module to predict the SL interactions in individual cancer cells (iSL). Extensive experiments have validated the model's successful integration of the multiple graph features and state-of-the-art performance in the prediction of potential SL gene pairs as well as generalization ability to novel genes.

Discovery of synthetic lethal interactions from large-scale pan-cancer perturbation screens

The development of cancer therapies is limited by the availability of suitable drug targets. Potential candidate drug targets can be identified based on the concept of synthetic lethality (SL), which refers to pairs of genes for which an aberration in either gene alone is non-lethal, but co-occurrence of the aberrations is lethal to the cell. Here, we present SLIdR (Synthetic Lethal Identification in R), a statistical framework for identifying SL pairs from large-scale perturbation screens. SLIdR successfully predicts SL pairs even with small sample sizes while minimizing the number of false positive targets. We apply SLIdR to Project DRIVE data and find both established and potential pan-cancer and cancer type-specific SL pairs consistent with findings from literature and drug response screening data. We experimentally validate two predicted SL interactions (ARID1A-TEAD1 and AXIN1-URI1) in hepatocellular carcinoma, thus corroborating the ability of SLIdR to identify potential drug targets.

ToMExO: A probabilistic tree-structured model for cancer progression

Identifying the interrelations among cancer driver genes and the patterns in which the driver genes get mutated is critical for understanding cancer. In this paper, we study cross-sectional data from cohorts of tumors to identify the cancer-type (or subtype) specific process in which the cancer driver genes accumulate critical mutations. We model this mutation accumulation process using a tree, where each node includes a driver gene or a set of driver genes. A mutation in each node enables its children to have a chance of mutating. This model simultaneously explains the mutual exclusivity patterns observed in mutations in specific cancer genes (by its nodes) and the temporal order of events (by its edges). We introduce a computationally efficient dynamic programming procedure for calculating the likelihood of our noisy datasets and use it to build our Markov Chain Monte Carlo (MCMC) inference algorithm, ToMExO. Together with a set of engineered MCMC moves, our fast likelihood calculations enable us to work with datasets with hundreds of genes and thousands of tumors, which cannot be dealt with using available cancer progression analysis methods. We demonstrate our method’s performance on several synthetic datasets covering various scenarios for cancer progression dynamics. Then, a comparison against two state-of-the-art methods on a moderate-size biological dataset shows the merits of our algorithm in identifying significant and valid patterns. Finally, we present our analyses of several large biological datasets, including colorectal cancer, glioblastoma, and pancreatic cancer. In all the analyses, we validate the results using a set of method-independent metrics testing the causality and significance of the relations identified by ToMExO or competing methods.

Synthetic lethal gene pairs: Experimental approaches and predictive models

Synthetic lethality (SL) refers to a genetic interaction in which the simultaneous perturbation of two genes leads to cell or organism death, whereas viability is maintained when only one of the pair is altered. The experimental exploration of these pairs and predictive modeling in computational biology contribute to our understanding of cancer biology and the development of cancer therapies. We extensively reviewed experimental technologies, public data sources, and predictive models in the study of synthetic lethal gene pairs and herein detail biological assumptions, experimental data, statistical models, and computational schemes of various predictive models, speculate regarding their influence on individual sample- and population-based synthetic lethal interactions, discuss the pros and cons of existing SL data and models, and highlight potential research directions in SL discovery.

Epistasis and evolutionary dependencies in human cancers

Cancer evolution is driven by the concerted action of multiple molecular alterations, which emerge and are selected during tumor progression. An alteration is selected when it provides an advantage to the tumor cell. However, the advantage provided by a specific alteration depends on the tumor lineage, cell epigenetic state, and presence of additional alterations. In this case, we say that an evolutionary dependency exists between an alteration and what influences its selection. Epistatic interactions between altered genes lead to evolutionary dependencies (EDs), by favoring or vetoing specific combinations of events. Large-scale cancer genomics studies have discovered examples of such dependencies, and showed that they influence tumor progression, disease phenotypes, and therapeutic response. In the past decade, several algorithmic approaches have been proposed to infer EDs from large-scale genomics datasets. These methods adopt diverse strategies to address common challenges and shed new light on cancer evolutionary trajectories. Here, we review these efforts starting from a simple conceptualization of the problem, presenting the tackled and still unmet needs in the field, and discussing the implications of EDs in cancer biology and precision oncology.

The role of drug-metabolizing enzymes in synthetic lethality of cancer

Drug-metabolizing enzymes (DMEs) have shown increasing importance in anticancer therapy. It is not only due to their effect on activation or deactivation of anticancer drugs, but also because of their extensive connections with pathological and biochemistry changes during tumorigenesis. Meanwhile, it has become more accessible to discovery anticancer drugs that selectively targeted cancer cells with the development of synthetic lethal screen technology. Synthetic lethal strategy makes use of unique genetic markers that different cancer cells from normal tissues to discovery anticancer agents. Dysregulation of DMEs has been found in various cancers, making them promising candidates for synthetic lethal strategy. In this review, we will systematically discuss about the role of DMEs in tumor progression, the application of synthetic lethality strategy in drug discovery, and a link between DMEs and synthetic lethal of cancer.

DNA Repair Mechanisms, Protein Interactions and Therapeutic Targeting of the MRN Complex

Repair of a DNA double-strand break relies upon a pathway of proteins to identify damage, regulate cell cycle checkpoints, and repair the damage. This process is initiated by a sensor protein complex, the MRN complex, comprised of three proteins-MRE11, RAD50, and NBS1. After a double-stranded break, the MRN complex recruits and activates ATM, in-turn activating other proteins such as BRCA1/2, ATR, CHEK1/2, PALB2 and RAD51. These proteins have been the focus of many studies for their individual roles in hereditary cancer syndromes and are included on several genetic testing panels. These panels have enabled us to acquire large amounts of genetic data, much of which remains a challenge to interpret due to the presence of variants of uncertain significance (VUS). While the primary aim of clinical testing is to accurately and confidently classify variants in order to inform medical management, the presence of VUSs has led to ambiguity in genetic counseling. Pathogenic variants within MRN complex genes have been implicated in breast, ovarian, prostate, colon cancers and gliomas; however, the hundreds of VUSs within MRE11, RAD50, and NBS1 precludes the application of these data in genetic guidance of carriers. In this review, we discuss the MRN complex's role in DNA double-strand break repair, its interactions with other cancer predisposing genes, the variants that can be found within the three MRN complex genes, and the MRN complex's potential as an anti-cancer therapeutic target.

NSF4SL: negative-sample-free contrastive learning for ranking synthetic lethal partner genes in human cancers

MOTIVATION

Detecting synthetic lethality (SL) is a promising strategy for identifying anti-cancer drug targets. Targeting SL partners of a primary gene mutated in cancer is selectively lethal to cancer cells. Due to high cost of wet-lab experiments and availability of gold standard SL data, supervised machine learning for SL prediction has been popular. However, most of the methods are based on binary classification and thus limited by the lack of reliable negative data. Contrastive learning can train models without any negative sample and is thus promising for finding novel SLs.

RESULTS

We propose NSF4SL, a negative-sample-free SL prediction model based on a contrastive learning framework. It captures the characteristics of positive SL samples by using two branches of neural networks that interact with each other to learn SL-related gene representations. Moreover, a feature-wise data augmentation strategy is used to mitigate the sparsity of SL data. NSF4SL significantly outperforms all baselines which require negative samples, even in challenging experimental settings. To the best of our knowledge, this is the first time that SL prediction is formulated as a gene ranking problem, which is more practical than the current formulation as binary classification. NSF4SL is the first contrastive learning method for SL prediction and its success points to a new direction of machine-learning methods for identifying novel SLs.

AVAILABILITY AND IMPLEMENTATION

Our source code is available at https://github.com/JieZheng-ShanghaiTech/NSF4SL.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

PiLSL: pairwise interaction learning-based graph neural network for synthetic lethality prediction in human cancers

MOTIVATION

Synthetic lethality (SL) is a type of genetic interaction in which the simultaneous inactivation of two genes leads to cell death, while the inactivation of a single gene does not affect the cell viability. It can effectively expand the range of anti-cancer therapeutic targets. SL interactions are identified mainly by experimental screening and computational prediction. Recent machine-learning methods mostly learn the representation of each gene individually, ignoring the representation of the pairwise interaction between two genes. In addition, the mechanisms of SL, the key to translating SL into cancer therapeutics, are often unclear.

RESULTS

To fill the gaps, we propose a pairwise interaction learning-based graph neural network (GNN) named PiLSL to learn the representation of pairwise interaction between two genes for SL prediction. First, we construct an enclosing graph for each pair of genes from a knowledge graph. Secondly, we design an attentive embedding propagation layer in a GNN to discriminate the importance among the edges in the enclosing graph and to learn the latent features of the pairwise interaction from the weighted enclosing graph. Finally, we further fuse the latent features with explicit features extracted from multi-omics data to obtain powerful gene representations for SL prediction. Extensive experimental results demonstrate that PiLSL outperforms the best baseline by a large margin and generalizes well under three realistic scenarios. Besides, PiLSL provides an explanation of SL mechanisms via the weighted paths in the enclosing graphs by attention mechanism.

AVAILABILITY AND IMPLEMENTATION

Our source code is available at https://github.com/JieZheng-ShanghaiTech/PiLSL.

Overcoming selection bias in synthetic lethality prediction

MOTIVATION

Synthetic lethality (SL) between two genes occurs when simultaneous loss of function leads to cell death. This holds great promise for developing anti-cancer therapeutics that target synthetic lethal pairs of endogenously disrupted genes. Identifying novel SL relationships through exhaustive experimental screens is challenging, due to the vast number of candidate pairs. Computational SL prediction is therefore sought to identify promising SL gene pairs for further experimentation. However, current SL prediction methods lack consideration for generalizability in the presence of selection bias in SL data.

RESULTS

We show that SL data exhibit considerable gene selection bias. Our experiments designed to assess the robustness of SL prediction reveal that models driven by the topology of known SL interactions (e.g. graph, matrix factorization) are especially sensitive to selection bias. We introduce selection bias-resilient synthetic lethality (SBSL) prediction using regularized logistic regression or random forests. Each gene pair is described by 27 molecular features derived from cancer cell line, cancer patient tissue and healthy donor tissue samples. SBSL models are built and tested using approximately 8000 experimentally derived SL pairs across breast, colon, lung and ovarian cancers. Compared to other SL prediction methods, SBSL showed higher predictive performance, better generalizability and robustness to selection bias. Gene dependency, quantifying the essentiality of a gene for cell survival, contributed most to SBSL predictions. Random forests were superior to linear models in the absence of dependency features, highlighting the relevance of mutual exclusivity of somatic mutations, co-expression in healthy tissue and differential expression in tumour samples.

AVAILABILITY AND IMPLEMENTATION

https://github.com/joanagoncalveslab/sbsl.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

SLOAD: a comprehensive database of cancer-specific synthetic lethal interactions for precision cancer therapy via multi-omics analysis

Synthetic lethality has been widely concerned because of its potential role in cancer treatment, which can be harnessed to selectively kill cancer cells via identifying inactive genes in a specific cancer type and further targeting the corresponding synthetic lethal partners. Herein, to obtain cancer-specific synthetic lethal interactions, we aimed to predict genetic interactions via a pan-cancer analysis from multiple molecular levels using random forest and then develop a user-friendly database. First, based on collected public gene pairs with synthetic lethal interactions, candidate gene pairs were analyzed via integrating multi-omics data, mainly including DNA mutation, copy number variation, methylation and mRNA expression data. Then, integrated features were used to predict cancer-specific synthetic lethal interactions using random forest. Finally, SLOAD (http://www.tmliang.cn/SLOAD) was constructed via integrating these findings, which was a user-friendly database for data searching, browsing, downloading and analyzing. These results can provide candidate cancer-specific synthetic lethal interactions, which will contribute to drug designing in cancer treatment that can promote therapy strategies based on the principle of synthetic lethality.

Database URL http://www.tmliang.cn/SLOAD/

Pervasive conditional selection of driver mutations and modular epistasis networks in cancer

Cancer driver mutations often display mutual exclusion or co-occurrence, underscoring the key role of epistasis in carcinogenesis. However, estimating the magnitude of epistasis and quantifying its effect on tumor evolution remains a challenge. We develop a method (Coselens) to quantify conditional selection on the excess of nonsynonymous substitutions in cancer genes. Coselens infers the number of drivers per gene in different partitions of a cancer genomics dataset using covariance-based mutation models and determines whether coding mutations in a gene affect selection for drivers in any other gene. Using Coselens, we identify 296 conditionally selected gene pairs across 16 cancer types in the TCGA dataset. Conditional selection affects 25%-50% of driver substitutions in tumors with >2 drivers. Conditionally co-selected genes form modular networks, whose structures challenge the traditional interpretation of within-pathway mutual exclusivity and across-pathway synergy, suggesting a more complex scenario where gene-specific across-pathway epistasis shapes differentiated cancer subtypes.

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.

A comprehensive analysis of ncRNA-mediated interactions reveals potential prognostic biomarkers in prostate adenocarcinoma

As one of common malignancies, prostate adenocarcinoma (PRAD) has been a growing health problem and a leading cause of cancer-related death. To obtain expression and functional relevant RNAs, we firstly screened candidate hub mRNAs and characterized their associations with cancer. Eight deregulated genes were identified and used to build a risk model (AUC was 0.972 at 10 years) that may be a specific biomarker for cancer prognosis. Then, relevant miRNAs and lncRNAs were screened, and the constructed primarily interaction networks showed the potential cross-talks among diverse RNAs. IsomiR landscapes were surveyed to understand the detailed isomiRs in relevant homologous miRNA loci, which largely enriched RNA interaction network due to diversities of sequence and expression. We finally characterized TK1, miR-222-3p and SNHG3 as crucial RNAs, and the abnormal expression patterns of them were correlated with poor survival outcomes. TK1 was found synthetic lethal interactions with other genes, implicating potential therapeutic target in precision medicine. LncRNA SNHG3 can sponge miR-222-3p to perturb RNA regulatory network and TK1 expression. These results demonstrate that TK1:miR-222-3p:SNHG3 axis may be a potential prognostic biomarker, which will contribute to further understanding cancer pathophysiology and providing potential therapeutic targets in precision medicine.

SynLethDB 2.0: a web-based knowledge graph database on synthetic lethality for novel anticancer drug discovery

Two genes are synthetic lethal if mutations in both genes result in impaired cell viability, while mutation of either gene does not affect the cell survival. The potential usage of synthetic lethality (SL) in anticancer therapeutics has attracted many researchers to identify synthetic lethal gene pairs. To include newly identified SLs and more related knowledge, we present a new version of the SynLethDB database to facilitate the discovery of clinically relevant SLs. We extended the first version of SynLethDB database significantly by including new SLs identified through Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) screening, a knowledge graph about human SLs, a new web interface, etc. Over 16 000 new SLs and 26 types of other relationships have been added, encompassing relationships among 14 100 genes, 53 cancers, 1898 drugs, etc. Moreover, a brand-new web interface has been developed to include modules such as SL query by disease or compound, SL partner gene set enrichment analysis and knowledge graph browsing through a dynamic graph viewer. The data can be downloaded directly from the website or through the RESTful Application Programming Interfaces (APIs). Database URL:  https://synlethdb.sist.shanghaitech.edu.cn/v2.

SL-Cloud: A Cloud-based resource to support synthetic lethal interaction discovery

Synthetic lethal interactions (SLIs), genetic interactions in which the simultaneous inactivation of two genes leads to a lethal phenotype, are promising targets for therapeutic intervention in cancer, as exemplified by the recent success of PARP inhibitors in treating BRCA1/2-deficient tumors. We present SL-Cloud, a new component of the Institute for Systems Biology Cancer Gateway in the Cloud (ISB-CGC), that provides an integrated framework of cloud-hosted data resources and curated workflows to enable facile prediction of SLIs. This resource addresses two main challenges related to SLI inference: the need to wrangle and preprocess large multi-omic datasets and the availability of multiple comparable prediction approaches. SL-Cloud enables customizable computational inference of SLIs and testing of prediction approaches across multiple datasets. We anticipate that cancer researchers will find utility in this tool for discovery of SLIs to support further investigation into potential drug targets for anticancer therapies.

SL-Cloud: A Cloud-based resource to support synthetic lethal interaction discovery

Synthetic lethal interactions (SLIs), genetic interactions in which the simultaneous inactivation of two genes leads to a lethal phenotype, are promising targets for therapeutic intervention in cancer, as exemplified by the recent success of PARP inhibitors in treating BRCA1/2-deficient tumors. We present SL-Cloud, a new component of the Institute for Systems Biology Cancer Gateway in the Cloud (ISB-CGC), that provides an integrated framework of cloud-hosted data resources and curated workflows to enable facile prediction of SLIs. This resource addresses two main challenges related to SLI inference: the need to wrangle and preprocess large multi-omic datasets and the availability of multiple comparable prediction approaches. SL-Cloud enables customizable computational inference of SLIs and testing of prediction approaches across multiple datasets. We anticipate that cancer researchers will find utility in this tool for discovery of SLIs to support further investigation into potential drug targets for anticancer therapies.

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.

Functional buffering via cell-specific gene expression promotes tissue homeostasis and cancer robustness

Functional buffering that ensures biological robustness is critical for maintaining tissue homeostasis, organismal survival, and evolution of novelty. However, the mechanism underlying functional buffering, particularly in multicellular organisms, remains largely elusive. Here, we proposed that functional buffering can be mediated via expression of buffering genes in specific cells and tissues, by which we named Cell-specific Expression-BUffering (CEBU). We developed an inference index (C-score) for CEBU by computing C-scores across 684 human cell lines using genome-wide CRISPR screens and transcriptomic RNA-seq. We report that C-score-identified putative buffering gene pairs are enriched for members of the same duplicated gene family, pathway, and protein complex. Furthermore, CEBU is especially prevalent in tissues of low regenerative capacity (e.g., bone and neuronal tissues) and is weakest in highly regenerative blood cells, linking functional buffering to tissue regeneration. Clinically, the buffering capacity enabled by CEBU can help predict patient survival for multiple cancers. Our results suggest CEBU as a potential buffering mechanism contributing to tissue homeostasis and cancer robustness in humans.

An Algorithm to Mine Therapeutic Motifs for Cancer from Networks of Genetic Interactions

Study of pairwise genetic interactions, such as mutually exclusive mutations, has led to understanding of underlying mechanisms in cancer. Investigation of various combinatorial motifs within networks of such interactions can lead to deeper insights into its mutational landscape and inform therapy development. One such motif called the Between-Pathway Model (BPM) represents redundant or compensatory pathways that can be therapeutically exploited. Finding such BPM motifs is challenging since most formulations require solving variants of the NP-complete maximum weight bipartite subgraph problem. In this paper we design an algorithm based on Integer Linear Programming (ILP) to solve this problem. In our experiments, our approach outperforms the best previous method to mine BPM motifs. Further, our ILP-based approach allows us to easily model additional application-specific constraints. We illustrate this advantage through a new application of BPM motifs that can potentially aid in finding combination therapies to combat cancer.

Literature-based translation from synthetic lethality screening into therapeutics targets: CD82 is a novel target for KRAS mutation in colon cancer

Synthetic lethality (SL) is an emerging therapeutic paradigm in cancer. We introduced a different approach to prioritize SL gene pairs through literature mining and RAS-mutant high-throughput screening (HTS) data. We matched essential genes from text-mining and mutant genes from the COSMIC and CCLE HTS datasets to build a prediction model of SL gene pairs. CCLE gene expression data were used to enrich the essential-mutant SL gene pairs using Spearman’s correlation coefficient and literature mining. In total, 223 essential trigger terms were extracted and ranked. The threshold of the essential gene score (Sg) was set to 10. We identified 586 genes essential for the SL prediction model of colon cancer. Seven essential RAS-mutant SL gene pairs were identified in our model, including CD82-KRAS/NRAS, PEBP1-NRAS, MT-CO2-HRAS, IFI27-NRAS/KRAS, and SUMO1-HRAS gene pairs. Using RAS-mutant HTS data validation, we identified two potential SL gene pairs, including the CD82 (essential gene)–KRAS (mutant gene) pair and CD82–NRAS pair in the DLD-1 colon cancer cell line (Spearman’s correlation p-values = 0.004786 and 0.00249, respectively). Based on further annotations by PubChem, we observed that digitonin targeted the complex comprising CD82, especially in KRAS-mutated HCT116 cancer cells. Moreover, we experimentally demonstrated that CD82 exhibited selective vulnerability in KRAS-mutant colorectal cancer. We used literature mining and HTS data to identify candidates for SL targets for RAS-mutant colon cancer.

Statistical tests for intra-tumour clonal co-occurrence and exclusivity

Tumour progression is an evolutionary process in which different clones evolve over time, leading to intra-tumour heterogeneity. Interactions between clones can affect tumour evolution and hence disease progression and treatment outcome. Intra-tumoural pairs of mutations that are overrepresented in a co-occurring or clonally exclusive fashion over a cohort of patient samples may be suggestive of a synergistic effect between the different clones carrying these mutations. We therefore developed a novel statistical testing framework, called GeneAccord, to identify such gene pairs that are altered in distinct subclones of the same tumour. We analysed our framework for calibration and power. By comparing its performance to baseline methods, we demonstrate that to control type I errors, it is essential to account for the evolutionary dependencies among clones. In applying GeneAccord to the single-cell sequencing of a cohort of 123 acute myeloid leukaemia patients, we find 1 clonally co-occurring and 8 clonally exclusive gene pairs. The clonally exclusive pairs mostly involve genes of the key signalling pathways.

Tumours typically display high levels of heterogeneity, not only between different tumours but also within a single one. Intra-tumour heterogeneity results from an evolutionary process, giving rise to different populations of cancer cells known as clones. How clones interact may affect tumour evolution, which in turn determines disease progression and treatment outcome. In practice, we may observe pairs of mutations that co-occur in clones or exclude each other more often than we would expect for a given cohort of patient samples. Exclusive pairs are suggestive that clones carrying one or the other mutation may cooperate in the evolutionary process. Targeting only one of them may then suffice to alter the tumour evolution. Therefore it is critical to have statistical methods which allow us to identify such pairs. GeneAccord is a novel statistical testing framework we developed especially to identify pairs of genes altered in distinct clones of the same tumour. Accounting for the evolutionary dependencies among clones emerged as critical to adequately control testing errors. In a cohort of 123 acute myeloid leukaemia patients, GeneAccord identified one clonally co-occurring and eight clonally exclusive gene pairs. The latter predominantly involved genes of key signalling pathways.

Precise Characterization of Genetic Interactions in Cancer via Molecular Network Refining Processes

Genetic interactions (GIs), such as the synthetic lethal interaction, are promising therapeutic targets in precision medicine. However, despite extensive efforts to characterize GIs by large-scale perturbation screening, considerable false positives have been reported in multiple studies. We propose a new computational approach for improved precision in GI identification by applying constraints that consider actual biological phenomena. In this study, GIs were characterized by assessing mutation, loss of function, and expression profiles in the DEPMAP database. The expression profiles were used to exclude loss-of-function data for nonexpressed genes in GI characterization. More importantly, the characterized GIs were refined based on Kyoto Encyclopedia of Genes and Genomes (KEGG) or protein–protein interaction (PPI) networks, under the assumption that genes genetically interacting with a certain mutated gene are adjacent in the networks. As a result, the initial GIs characterized with CRISPR and RNAi screenings were refined to 65 and 23 GIs based on KEGG networks and to 183 and 142 GIs based on PPI networks. The evaluation of refined GIs showed improved precision with respect to known synthetic lethal interactions. The refining process also yielded a synthetic partner network (SPN) for each mutated gene, which provides insight into therapeutic strategies for the mutated genes; specifically, exploring the SPN of mutated BRAF revealed ELAVL1 as a potential target for treating BRAF-mutated cancer, as validated by previous research. We expect that this work will advance cancer therapeutic research.

Revealing protein-protein interactions at the transcriptome scale by sequencing

We describe PROPER-seq (protein-protein interaction sequencing) to map protein-protein interactions (PPIs) en masse. PROPER-seq first converts transcriptomes of input cells into RNA-barcoded protein libraries, in which all interacting protein pairs are captured through nucleotide barcode ligation, recorded as chimeric DNA sequences, and decoded at once by sequencing and mapping. We applied PROPER-seq to human embryonic kidney cells, T lymphocytes, and endothelial cells and identified 210,518 human PPIs (collected in the PROPER v.1.0 database). Among these, 1,365 and 2,480 PPIs are supported by published co-immunoprecipitation (coIP) and affinity purification-mass spectrometry (AP-MS) data, 17,638 PPIs are predicted by the prePPI algorithm without previous experimental validation, and 100 PPIs overlap human synthetic lethal gene pairs. In addition, four previously uncharacterized interaction partners with poly(ADP-ribose) polymerase 1 (PARP1) (a critical protein in DNA repair) known as XPO1, MATR3, IPO5, and LEO1 are validated in vivo. PROPER-seq presents a time-effective technology to map PPIs at the transcriptome scale, and PROPER v.1.0 provides a rich resource for studying PPIs.

Synthetic lethality theory approaches to effective substance discovery and functional mechanisms elucidation of anti-cancer phytomedicine

BACKGROUND

Longstanding, successful use of combinations of phytopharmaceuticals in traditional Chinese medicine (TCM) has caught the attention of several pharmacologists to natural medicines. However, the development and popularisation of TCM is mainly limited because of the unavailability of reports clarifying the mechanisms of action and pharmacologically active ingredients in such formulations. Previous studies on natural medicines have mostly focused on their dominant components using forward pharmacology which often neglects trace components. It is necessary to assess the pharmacological and therapeutic superiority of many such trace components in comparison with single constituents.

PURPOSE

In this study, we aimed to propose a new pharmacological research strategy for TCM. In particular, we presented the possibility that the effective mechanism of action of trace components of TCM is based on synthetic lethality. We sincerely hope to explore this theory further.

METHOD

We obtained retrieve published research information related to synthetic lethality, phytochemicals and Chinese medicine from PubMed and Google scholar. Based on the inclusion criteria, 71 studies were selected and discussed in this review.

RESULTS

As an interaction among genes, synthetic lethality can amplify co-regulatory biological effects exponentially. Synthetic strategies have been successfully applied for research and development of anti-tumour agents, including poly ADP-ribose polymerase inhibitors and clinical combination of chemotherapeutic agents for efficacy enhancement and toxicity reduction. TCM drugs contain several secondary metabolites to combat environmental stresses, providing a multi-component basis for corresponding synergistic targets. Therefore, we aimed to study whether this method could be used to identify active components present in trace amounts in TCM drugs. Based on a reverse concept of target-component-effect and identified synergistic targets, we explored the mechanisms of action of weakly active components present in trace amounts in TCM drugs to assess combinations of potential synergistic components.

CONCLUSION

This pattern of synthetic lethality not only elucidated the mechanisms of action of TCM drugs from a new perspective but also inspired future studies on discovering naturally occurring active components.

Systematic identification of novel cancer genes through analysis of deep shRNA perturbation screens

Systematic perturbation screens provide comprehensive resources for the elucidation of cancer driver genes. The perturbation of many genes in relatively few cell lines in such functional screens necessitates the development of specialized computational tools with sufficient statistical power. Here we developed APSiC (Analysis of Perturbation Screens for identifying novel Cancer genes) to identify genetic drivers and effectors in perturbation screens even with few samples. Applying APSiC to the shRNA screen Project DRIVE, APSiC identified well-known and novel putative mutational and amplified cancer genes across all cancer types and in specific cancer types. Additionally, APSiC discovered tumor-promoting and tumor-suppressive effectors, respectively, for individual cancer types, including genes involved in cell cycle control, Wnt/β-catenin and hippo signalling pathways. We functionally demonstrated that LRRC4B, a putative novel tumor-suppressive effector, suppresses proliferation by delaying cell cycle and modulates apoptosis in breast cancer. We demonstrate APSiC is a robust statistical framework for discovery of novel cancer genes through analysis of large-scale perturbation screens. The analysis of DRIVE using APSiC is provided as a web portal and represents a valuable resource for the discovery of novel cancer genes.

Uncovering cancer vulnerabilities by machine learning prediction of synthetic lethality

BACKGROUND

Synthetic lethality describes a genetic interaction between two perturbations, leading to cell death, whereas neither event alone has a significant effect on cell viability. This concept can be exploited to specifically target tumor cells. CRISPR viability screens have been widely employed to identify cancer vulnerabilities. However, an approach to systematically infer genetic interactions from viability screens is missing.

METHODS

Here we describe PAn-canceR Inferred Synthetic lethalities (PARIS), a machine learning approach to identify cancer vulnerabilities. PARIS predicts synthetic lethal (SL) interactions by combining CRISPR viability screens with genomics and transcriptomics data across hundreds of cancer cell lines profiled within the Cancer Dependency Map.

RESULTS

Using PARIS, we predicted 15 high confidence SL interactions within 549 DNA damage repair (DDR) genes. We show experimental validation of an SL interaction between the tumor suppressor CDKN2A, thymidine phosphorylase (TYMP) and the thymidylate synthase (TYMS), which may allow stratifying patients for treatment with TYMS inhibitors. Using genome-wide mapping of SL interactions for DDR genes, we unraveled a dependency between the aldehyde dehydrogenase ALDH2 and the BRCA-interacting protein BRIP1. Our results suggest BRIP1 as a potential therapeutic target in ~ 30% of all tumors, which express low levels of ALDH2.

CONCLUSIONS

PARIS is an unbiased, scalable and easy to adapt platform to identify SL interactions that should aid in improving cancer therapy with increased availability of cancer genomics data.

Mapping the landscape of synthetic lethal interactions in liver cancer

Almost all the current therapies against liver cancer are based on the "one size fits all" principle and offer only limited survival benefit. Fortunately, synthetic lethality (SL) may provide an alternate route towards individualized therapy in liver cancer. The concept that simultaneous losses of two genes are lethal to a cell while a single loss is non-lethal can be utilized to selectively eliminate tumors with genetic aberrations. Methods: To infer liver cancer-specific SL interactions, we propose a computational pipeline termed SiLi (statistical inference-based synthetic lethality identification) that incorporates five inference procedures. Based on large-scale sequencing datasets, SiLi analysis was performed to identify SL interactions in liver cancer. Results: By SiLi analysis, a total of 272 SL pairs were discerned, which included 209 unique target candidates. Among these, polo-like kinase 1 (PLK1) was considered to have considerable therapeutic potential. Further computational and experimental validation of the SL pair TP53-PLK1 demonstrated that inhibition of PLK1 could be a novel therapeutic strategy specifically targeting those patients with TP53-mutant liver tumors. Conclusions: In this study, we report a comprehensive analysis of synthetic lethal interactions of liver cancer. Our findings may open new possibilities for patient-tailored therapeutic interventions in liver cancer.