A Network of Conserved Synthetic Lethal Interactions for Exploration of Precision 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.

Unleashing the Power of Synthetic Lethality: Augmenting Treatment Efficacy through Synergistic Integration with Chemotherapy Drugs

Cancer is the second leading cause of death in the world, and chemotherapy is one of the main methods of cancer treatment. However, the resistance of cancer cells to chemotherapeutic drugs has always been the main reason affecting the therapeutic effect. Synthetic lethality has emerged as a promising approach to augment the sensitivity of cancer cells to chemotherapy agents. Synthetic lethality (SL) refers to the specific cell death resulting from the simultaneous mutation of two non-lethal genes, which individually allow cell survival. This comprehensive review explores the classification of SL, screening methods, and research advancements in SL inhibitors, including Poly (ADP-ribose) polymerase (PARP) inhibitors, Ataxia telangiectasia and Rad3-related (ATR) inhibitors, WEE1 G2 checkpoint kinase (WEE1) inhibitors, and protein arginine methyltransferase 5 (PRMT5) inhibitors. Emphasizing their combined use with chemotherapy drugs, we aim to unveil more effective treatment strategies for cancer patients.

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.

A landscape of response to drug combinations in non-small cell lung cancer

Combination of anti-cancer drugs is broadly seen as way to overcome the often-limited efficacy of single agents. The design and testing of combinations are however very challenging. Here we present a uniquely large dataset screening over 5000 targeted agent combinations across 81 non-small cell lung cancer cell lines. Our analysis reveals a profound heterogeneity of response across the tumor models. Notably, combinations very rarely result in a strong gain in efficacy over the range of response observable with single agents. Importantly, gain of activity over single agents is more often seen when co-targeting functionally proximal genes, offering a strategy for designing more efficient combinations. Because combinatorial effect is strongly context specific, tumor specificity should be achievable. The resource provided, together with an additional validation screen sheds light on major challenges and opportunities in building efficacious combinations against cancer and provides an opportunity for training computational models for synergy prediction.

Targeting Mutant Dicer Tumorigenesis in Pleuropulmonary Blastoma via Inhibition of RNA Polymerase I

DICER1 mutations predispose to increased risk for various cancers, particularly pleuropulmonary blastoma (PPB), the commonest lung malignancy of childhood. There is a paucity of directly actionable molecular targets as these tumors are driven by loss-of-function mutations of DICER1. Therapeutic development for PPB is further limited by a lack of biologically and physiologically-representative disease models. Given recent evidence of Dicer's role as a haploinsufficient tumor suppressor regulating RNA polymerase I (Pol I), Pol I inhibition could abrogate mutant Dicer-mediated accumulation of stalled polymerases to trigger apoptosis. Hence, we developed a novel sub-pleural orthotopic PPB patient-derived xenograft (PDX) model that retained both RNase IIIa and IIIb hotspot mutations and recapitulated the cardiorespiratory physiology of intra-thoracic disease, and with it evaluated the tolerability and efficacy of first-in-class Pol I inhibitor CX-5461. In PDX tumors, CX-5461 significantly reduced H3K9 di-methylation and increased nuclear p53 expression, within 24 hours' exposure. Following treatment at the maximum tolerated dosing regimen (12 doses, 30mg/kg), tumors were smaller and less hemorrhagic than controls, with significantly decreased cellular proliferation, and increased apoptosis. As demonstrated in a novel intra-thoracic tumor model of PPB, Pol I inhibition with CX-5461 could be a tolerable and clinically-feasible therapeutic strategy for mutant Dicer tumors, inducing anti-tumor effects by decreasing H3K9 methylation and enhancing p53-mediated apoptosis.

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.

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.

Genome-wide base editor screen identifies regulators of protein abundance in yeast

Proteins are key molecular players in a cell, and their abundance is extensively regulated not just at the level of gene expression but also post-transcriptionally. Here, we describe a genetic screen in yeast that enables systematic characterization of how protein abundance regulation is encoded in the genome. The screen combines a CRISPR/Cas9 base editor to introduce point mutations with fluorescent tagging of endogenous proteins to facilitate a flow-cytometric readout. We first benchmarked base editor performance in yeast with individual gRNAs as well as in positive and negative selection screens. We then examined the effects of 16,452 genetic perturbations on the abundance of eleven proteins representing a variety of cellular functions. We uncovered hundreds of regulatory relationships, including a novel link between the GAPDH isoenzymes Tdh1/2/3 and the Ras/PKA pathway. Many of the identified regulators are specific to one of the eleven proteins, but we also found genes that, upon perturbation, affected the abundance of most of the tested proteins. While the more specific regulators usually act transcriptionally, broad regulators often have roles in protein translation. Overall, our novel screening approach provides unprecedented insights into the components, scale and connectedness of the protein regulatory network.

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/

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.

Predicting and explaining the impact of genetic disruptions and interactions on organismal viability

MOTIVATION

Existing computational models can predict single- and double-mutant fitness but they do have limitations. First, they are often tested via evaluation metrics that are inappropriate for imbalanced datasets. Second, all of them only predict a binary outcome (viable or not, and negatively interacting or not). Third, most are uninterpretable black box machine learning models.

RESULTS

Budding yeast datasets were used to develop high-performance Multinomial Regression (MN) models capable of predicting the impact of single, double and triple genetic disruptions on viability. These models are interpretable and give realistic non-binary predictions and can predict negative genetic interactions (GIs) in triple-gene knockouts. They are based on a limited set of gene features and their predictions are influenced by the probability of target gene participating in molecular complexes or pathways. Furthermore, the MN models have utility in other organisms such as fission yeast, fruit flies and humans, with the single gene fitness MN model being able to distinguish essential genes necessary for cell-autonomous viability from those required for multicellular survival. Finally, our models exceed the performance of previous models, without sacrificing interpretability.

AVAILABILITY AND IMPLEMENTATION

All code and processed datasets used to generate results and figures in this manuscript are available at our Github repository at https://github.com/KISRDevelopment/cell_viability_paper. The repository also contains a link to the GI prediction website that lets users search for GIs using the MN models.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Impairment of RAD17 Functions by miR-506-3p as a Novel Synthetic Lethal Approach Targeting DNA Repair Pathways in Ovarian Cancer

Epithelial ovarian cancer (EOC) remains the most lethal gynecological cancer and development of chemo-resistance is a major factor in disease relapse. Homologous recombination (HR) is a critical pathway for DNA double strand break repair and its deficiency is associated to a better response to DNA damage-inducing agents. Strategies to inhibit HR-mediated DNA repair is a clinical need to improve patients’ outcome. MicroRNA (miRNAs) affect most of cellular processes including response to cancer treatment. We previously showed that miR-506-3p targets RAD51, an essential HR component. In this study we demonstrated that: i) another HR component, RAD17, is also a direct target of miR-506-3p and that it is involved in mediating miR-506-3p phenotypic effects; ii) the impairment of miR-506-3p binding to RAD17 3’ UTR reverted the miR-506-3p induced platinum sensitization; iii) miR-506-3p/RAD17 axis reduces the ability of EOC cell to sense DNA damage, abrogates the G2/M cell cycle checkpoint thus delaying the G2/M cell cycle arrest likely allowing the entry into mitosis of heavily DNA-damaged cells with a consequent mitotic catastrophe; iv) RAD17 expression, regulated by miR-506-3p, is synthetically lethal with inhibitors of cell cycle checkpoint kinases Chk1 and Wee1 in platinum resistant cell line. Overall miR-506-3p expression may recapitulate a BRCAness phenotype sensitizing EOC cells to chemotherapy and helping in selecting patients susceptible to DNA damaging drugs in combination with new small molecules targeting DNA-damage repair pathway.

The PICLS high-throughput screening method for agents extending cellular longevity identifies 2,5-anhydro-D-mannitol as novel anti-aging compound

Although aging is the biggest risk factor for human chronic (cancer, diabetic, cardiovascular, and neurodegenerative) diseases, few interventions are known besides caloric restriction and a small number of drugs (with substantial side effects) that directly address aging. Thus, there is an urgent need for new options that can generally delay aging processes and prevent age-related diseases. Cellular aging is at the basis of aging processes. Chronological lifespan (CLS) of yeast Saccharomyces cerevisiae is the well-established model system for investigating the interventions of human post-mitotic cellular aging. CLS is defined as the number of days cells remain viable in a stationary phase. We developed a new, cheap, and fast quantitative method for measuring CLS in cell cultures incubated together with various chemical agents and controls on 96-well plates. Our PICLS protocol with (1) the use of propidium iodide for fluorescent-based cell survival reading in a microplate reader and (2) total cell count measurement via OD600nm absorption from the same plate provides real high-throughput capacity. Depending on logistics, large numbers of plates can be processed in parallel so that the screening of thousands of compounds becomes feasible in a short time. The method was validated by measuring the effect of rapamycin and calorie restriction on yeast CLS. We utilized this approach for chemical agent screening. We discovered the anti-aging/geroprotective potential of 2,5-anhydro-D-mannitol (2,5-AM) and suggest its usage individually or in combination with other anti-aging interventions.

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.

Identification of Kic1p and Cdc42p as Novel Targets to Engineer Yeast Acetic Acid Stress Tolerance

Robust yeast strains that are tolerant to multiple stress environments are desired for an efficient biorefinery. Our previous studies revealed that zinc sulfate serves as an important nutrient for stress tolerance of budding yeast Saccharomyces cerevisiae. Acetic acid is a common inhibitor in cellulosic hydrolysate, and the development of acetic acid-tolerant strains is beneficial for lignocellulosic biorefineries. In this study, comparative proteomic studies were performed using S. cerevisiae cultured under acetic acid stress with or without zinc sulfate addition, and novel zinc-responsive proteins were identified. Among the differentially expressed proteins, the protein kinase Kic1p and the small rho-like GTPase Cdc42p, which is required for cell integrity and regulation of cell polarity, respectively, were selected for further studies. Overexpression of KIC1 and CDC42 endowed S. cerevisiae with faster growth and ethanol fermentation under the stresses of acetic acid and mixed inhibitors, as well as in corncob hydrolysate. Notably, the engineered yeast strains showed a 12 h shorter lag phase under the three tested conditions, leading to up to 52.99% higher ethanol productivity than that of the control strain. Further studies showed that the transcription of genes related to stress response was significantly upregulated in the engineered strains under the stress condition. Our results in this study provide novel insights in exploring zinc-responsive proteins for applications of synthetic biology in developing a robust industrial yeast.

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.

Harnessing Synthetic Lethal Interactions for Personalized Medicine

Two genes are said to have synthetic lethal (SL) interactions if the simultaneous mutations in a cell lead to lethality, but each individual mutation does not. Targeting SL partners of mutated cancer genes can kill cancer cells but leave normal cells intact. The applicability of translating this concept into clinics has been demonstrated by three drugs that have been approved by the FDA to target PARP for tumors bearing mutations in BRCA1/2. This article reviews applications of the SL concept to translational cancer medicine over the past five years. Topics are (1) exploiting the SL concept for drug combinations to circumvent tumor resistance, (2) using synthetic lethality to identify prognostic and predictive biomarkers, (3) applying SL interactions to stratify patients for targeted and immunotherapy, and (4) discussions on challenges and future directions.

Identifying digenic disease genes via machine learning in the Undiagnosed Diseases Network

Rare diseases affect millions of people worldwide, and discovering their genetic causes is challenging. More than half of the individuals analyzed by the Undiagnosed Diseases Network (UDN) remain undiagnosed. The central hypothesis of this work is that many of these rare genetic disorders are caused by multiple variants in more than one gene. However, given the large number of variants in each individual genome, experimentally evaluating combinations of variants for potential to cause disease is currently infeasible. To address this challenge, we developed the digenic predictor (DiGePred), a random forest classifier for identifying candidate digenic disease gene pairs by features derived from biological networks, genomics, evolutionary history, and functional annotations. We trained the DiGePred classifier by using DIDA, the largest available database of known digenic-disease-causing gene pairs, and several sets of non-digenic gene pairs, including variant pairs derived from unaffected relatives of UDN individuals. DiGePred achieved high precision and recall in cross-validation and on a held-out test set (PR area under the curve > 77%), and we further demonstrate its utility by using digenic pairs from the recent literature. In contrast to other approaches, DiGePred also appropriately controls the number of false positives when applied in realistic clinical settings. Finally, to enable the rapid screening of variant gene pairs for digenic disease potential, we freely provide the predictions of DiGePred on all human gene pairs. Our work enables the discovery of genetic causes for rare non-monogenic diseases by providing a means to rapidly evaluate variant gene pairs for the potential to cause digenic disease.

A systematic analysis of genetic interactions and their underlying biology in childhood cancer

Childhood cancer is a major cause of child death in developed countries. Genetic interactions between mutated genes play an important role in cancer development. They can be detected by searching for pairs of mutated genes that co-occur more (or less) often than expected. Co-occurrence suggests a cooperative role in cancer development, while mutual exclusivity points to synthetic lethality, a phenomenon of interest in cancer treatment research. Little is known about genetic interactions in childhood cancer. We apply a statistical pipeline to detect genetic interactions in a combined dataset comprising over 2,500 tumors from 23 cancer types. The resulting genetic interaction map of childhood cancers comprises 15 co-occurring and 27 mutually exclusive candidates. The biological explanation of most candidates points to either tumor subtype, pathway epistasis or cooperation while synthetic lethality plays a much smaller role. Thus, other explanations beyond synthetic lethality should be considered when interpreting genetic interaction test results.

Folliculin: A Regulator of Transcription Through AMPK and mTOR Signaling Pathways

Folliculin (FLCN) is a tumor suppressor gene responsible for the inherited Birt-Hogg-Dubé (BHD) syndrome, which affects kidneys, skin and lungs. FLCN is a highly conserved protein that forms a complex with folliculin interacting proteins 1 and 2 (FNIP1/2). Although its sequence does not show homology to known functional domains, structural studies have determined a role of FLCN as a GTPase activating protein (GAP) for small GTPases such as Rag GTPases. FLCN GAP activity on the Rags is required for the recruitment of mTORC1 and the transcriptional factors TFEB and TFE3 on the lysosome, where mTORC1 phosphorylates and inactivates these factors. TFEB/TFE3 are master regulators of lysosomal biogenesis and function, and autophagy. By this mechanism, FLCN/FNIP complex participates in the control of metabolic processes. AMPK, a key regulator of catabolism, interacts with FLCN/FNIP complex. FLCN loss results in constitutive activation of AMPK, which suggests an additional mechanism by which FLCN/FNIP may control metabolism. AMPK regulates the expression and activity of the transcriptional cofactors PGC1α/β, implicated in the control of mitochondrial biogenesis and oxidative metabolism. In this review, we summarize our current knowledge of the interplay between mTORC1, FLCN/FNIP, and AMPK and their implications in the control of cellular homeostasis through the transcriptional activity of TFEB/TFE3 and PGC1α/β. Other pathways and cellular processes regulated by FLCN will be briefly discussed.

A Perspective on Synthetic Biology in Drug Discovery and Development-Current Impact and Future Opportunities

The global impact of synthetic biology has been accelerating, because of the plummeting cost of DNA synthesis, advances in genetic engineering, growing understanding of genome organization, and explosion in data science. However, much of the discipline's application in the pharmaceutical industry remains enigmatic. In this review, we highlight recent examples of the impact of synthetic biology on target validation, assay development, hit finding, lead optimization, and chemical synthesis, through to the development of cellular therapeutics. We also highlight the availability of tools and technologies driving the discipline. Synthetic biology is certainly impacting all stages of drug discovery and development, and the recognition of the discipline's contribution can further enhance the opportunities for the drug discovery and development value chain.

New Insights into CDK Regulators: Novel Opportunities for Cancer Therapy

Cyclins and their catalytic partners, the cyclin-dependent kinases (CDKs), control the transition between different phases of the cell cycle. CDK/cyclin activity is regulated by CDK inhibitors (CKIs), currently comprising the CDK-interacting protein/kinase inhibitory protein (CIP/KIP) family and the inhibitor of kinase (INK) family. Recent studies have identified a third group of CKIs, called ribosomal protein-inhibiting CDKs (RPICs). RPICs were discovered in the context of cellular senescence, a stable cell cycle arrest with tumor-suppressing abilities. RPICs accumulate in the nonribosomal fraction of senescent cells due to a decrease in rRNA biogenesis. Accordingly, RPICs are often downregulated in human cancers together with other ribosomal proteins, the tumor-suppressor functions of which are still under study. In this review, we discuss unique therapies that have been developed to target CDK activity in the context of cancer treatment or senescence-associated pathologies, providing novel tools for precision medicine.

Molecularly Targeted Clinical Trials

Glioblastoma remains incurable despite advances in surgery, radiation, and chemotherapy, underscoring the need for new therapies. The genetic heterogenicity, presence of redundant molecular pathways, and the blood-brain barrier have limited the applicability of molecularly targeted agents. The therapeutic benefit seen with a small subset of patients suggests, however, that patient selection is critical. Recent investigations show that molecularly targeted synthetic lethality is a promising complementary approach. The article provides an overview of the challenges of molecularly targeted therapy in adults with glioblastoma, including current trials and future therapeutic directions.

τ-SGA: synthetic genetic array analysis for systematically screening and quantifying trigenic interactions in yeast

Systematic complex genetic interaction studies have provided insight into high-order functional redundancies and genetic network wiring of the cell. Here, we describe a method for screening and quantifying trigenic interactions from ordered arrays of yeast strains grown on agar plates as individual colonies. The protocol instructs users on the trigenic synthetic genetic array analysis technique, τ-SGA, for high-throughput screens. The steps describe construction of the double-mutant query strains and the corresponding single-mutant control query strains, which are screened in parallel in two replicates. The screening experimental set-up consists of sequential replica-pinning steps that enable automated mating, meiotic recombination and successive haploid selection steps for the generation of triple mutants, which are scored for colony size as a proxy for fitness, which enables the calculation of trigenic interactions. The procedure described here was used to conduct 422 trigenic interaction screens, which generated ~460,000 yeast triple mutants for trigenic interaction analysis. Users should be familiar with robotic equipment required for high-throughput genetic interaction screens and be proficient at the command line to execute the scoring pipeline. Large-scale screen computational analysis is achieved by using MATLAB pipelines that score raw colony size data to produce τ-SGA interaction scores. Additional recommendations are included for optimizing experimental design and analysis of smaller-scale trigenic interaction screens by using a web-based analysis system, SGAtools. This protocol provides a resource for those who would like to gain a deeper, more practical understanding of trigenic interaction screening and quantification methodology.

Differential Co-Expression Analyses Allow the Identification of Critical Signalling Pathways Altered during Tumour Transformation and Progression

Biological systems respond to perturbations through the rewiring of molecular interactions, organised in gene regulatory networks (GRNs). Among these, the increasingly high availability of transcriptomic data makes gene co-expression networks the most exploited ones. Differential co-expression networks are useful tools to identify changes in response to an external perturbation, such as mutations predisposing to cancer development, and leading to changes in the activity of gene expression regulators or signalling. They can help explain the robustness of cancer cells to perturbations and identify promising candidates for targeted therapy, moreover providing higher specificity with respect to standard co-expression methods. Here, we comprehensively review the literature about the methods developed to assess differential co-expression and their applications to cancer biology. Via the comparison of normal and diseased conditions and of different tumour stages, studies based on these methods led to the definition of pathways involved in gene network reorganisation upon oncogenes’ mutations and tumour progression, often converging on immune system signalling. A relevant implementation still lagging behind is the integration of different data types, which would greatly improve network interpretability. Most importantly, performance and predictivity evaluation of the large variety of mathematical models proposed would urgently require experimental validations and systematic comparisons. We believe that future work on differential gene co-expression networks, complemented with additional omics data and experimentally tested, will considerably improve our insights into the biology of tumours.

Recruitment of Xrn1 to stress-induced genes allows efficient transcription by controlling RNA polymerase II backtracking

A new paradigm has emerged proposing that the crosstalk between nuclear transcription and cytoplasmic mRNA stability keeps robust mRNA levels in cells under steady-state conditions. A key piece in this crosstalk is the highly conserved 5'-3' RNA exonuclease Xrn1, which degrades most cytoplasmic mRNAs but also associates with nuclear chromatin to activate transcription by not well-understood mechanisms. Here, we investigated the role of Xrn1 in the transcriptional response of Saccharomyces cerevisiae cells to osmotic stress. We show that a lack of Xrn1 results in much lower transcriptional induction of the upregulated genes but in similar high levels of their transcripts because of parallel mRNA stabilization. Unexpectedly, lower transcription in xrn1 occurs with a higher accumulation of RNA polymerase II (RNAPII) at stress-inducible genes, suggesting that this polymerase remains inactive backtracked. Xrn1 seems to be directly implicated in the formation of a competent elongation complex because Xrn1 is recruited to the osmotic stress-upregulated genes in parallel with the RNAPII complex, and both are dependent on the mitogen-activated protein kinase Hog1. Our findings extend the role of Xrn1 in preventing the accumulation of inactive RNAPII at highly induced genes to other situations of rapid and strong transcriptional upregulation.

Quantitative and multiplexed chemical-genetic phenotyping in mammalian cells with QMAP-Seq

Chemical-genetic interaction profiling in model organisms has proven powerful in providing insights into compound mechanism of action and gene function. However, identifying chemical-genetic interactions in mammalian systems has been limited to low-throughput or computational methods. Here, we develop Quantitative and Multiplexed Analysis of Phenotype by Sequencing (QMAP-Seq), which leverages next-generation sequencing for pooled high-throughput chemical-genetic profiling. We apply QMAP-Seq to investigate how cellular stress response factors affect therapeutic response in cancer. Using minimal automation, we treat pools of 60 cell types—comprising 12 genetic perturbations in five cell lines—with 1440 compound-dose combinations, generating 86,400 chemical-genetic measurements. QMAP-Seq produces precise and accurate quantitative measures of acute drug response comparable to gold standard assays, but with increased throughput at lower cost. Moreover, QMAP-Seq reveals clinically actionable drug vulnerabilities and functional relationships involving these stress response factors, many of which are activated in cancer. Thus, QMAP-Seq provides a broadly accessible and scalable strategy for chemical-genetic profiling in mammalian cells.

Identifying chemical-genetic interactions in mammalian cells is limited to low-throughput or computational methods. Here, the authors present QMAP-Seq, a broadly accessible and scalable approach that uses NGS for pooled high-throughput chemical-genetic profiling in mammalian cells.

Integrative omics analysis reveals relationships of genes with synthetic lethal interactions through a pan-cancer analysis

Synthetic lethality is thought to play an important role in anticancer therapies. Herein, to understand the potential distributions and relationships between synthetic lethal interactions between genes, especially for pairs deriving from different sources, we performed an integrative analysis of genes at multiple molecular levels. Based on inter-species phylogenetic conservation of synthetic lethal interactions, gene pairs from yeast and humans were analyzed; a total of 37,588 candidate gene pairs containing 7,816 genes were collected. Of these, 49.74% of genes had 2–10 interactions, 22.93% were involved in hallmarks of cancer, and 21.61% were identified as core essential genes. Many genes were shown to have important biological roles via functional enrichment analysis, and 65 were identified as potentially crucial in the pathophysiology of cancer. Gene pairs with dysregulated expression patterns had higher prognostic values. Further screening based on mutation and expression levels showed that remaining gene pairs were mainly derived from human predicted or validated pairs, while most predicted pairs from yeast were filtered from analysis. Genes with synthetic lethality were further analyzed with their interactive microRNAs (miRNAs) at the isomiR level which have been widely studied as negatively regulatory molecules. The miRNA–mRNA interaction network revealed that many synthetic lethal genes contributed to the cell cycle (seven of 12 genes), cancer pathways (five of 12 genes), oocyte meiosis, the p53 signaling pathway, and hallmarks of cancer. Our study contributes to the understanding of synthetic lethal interactions and promotes the application of genetic interactions in further cancer precision medicine.

Predicting synthetic lethal interactions using heterogeneous data sources

MOTIVATION

A synthetic lethal (SL) interaction is a relationship between two functional entities where the loss of either one of the entities is viable but the loss of both entities is lethal to the cell. Such pairs can be used as drug targets in targeted anticancer therapies, and so, many methods have been developed to identify potential candidate SL pairs. However, these methods use only a subset of available data from multiple platforms, at genomic, epigenomic and transcriptomic levels; and hence are limited in their ability to learn from complex associations in heterogeneous data sources.

RESULTS

In this article, we develop techniques that can seamlessly integrate multiple heterogeneous data sources to predict SL interactions. Our approach obtains latent representations by collective matrix factorization-based techniques, which in turn are used for prediction through matrix completion. Our experiments, on a variety of biological datasets, illustrate the efficacy and versatility of our approach, that outperforms state-of-the-art methods for predicting SL interactions and can be used with heterogeneous data sources with minimal feature engineering.

AVAILABILITY AND IMPLEMENTATION

Software available at https://github.com/lianyh.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A cancer drug atlas enables synergistic targeting of independent drug vulnerabilities

Personalized cancer treatments using combinations of drugs with a synergistic effect is attractive but proves to be highly challenging. Here we present an approach to uncover the efficacy of drug combinations based on the analysis of mono-drug effects. For this we used dose-response data from pharmacogenomic encyclopedias and represent these as a drug atlas. The drug atlas represents the relations between drug effects and allows to identify independent processes for which the tumor might be particularly vulnerable when attacked by two drugs. Our approach enables the prediction of combination-therapy which can be linked to tumor-driving mutations. By using this strategy, we can uncover potential effective drug combinations on a pan-cancer scale. Predicted synergies are provided and have been validated in glioblastoma, breast cancer, melanoma and leukemia mouse-models, resulting in therapeutic synergy in 75% of the tested models. This indicates that we can accurately predict effective drug combinations with translational value.

The yeast exoribonuclease Xrn1 and associated factors modulate RNA polymerase II processivity in 5' and 3' gene regions

mRNA levels are determined by the balance between mRNA synthesis and decay. Protein factors that mediate both processes, including the 5'-3' exonuclease Xrn1, are responsible for a cross-talk between the two processes that buffers steady-state mRNA levels. However, the roles of these proteins in transcription remain elusive and controversial. Applying native elongating transcript sequencing (NET-seq) to yeast cells, we show that Xrn1 functions mainly as a transcriptional activator and that its disruption manifests as a reduction of RNA polymerase II (Pol II) occupancy downstream of transcription start sites. By combining our sequencing data and mathematical modeling of transcription, we found that Xrn1 modulates transcription initiation and elongation of its target genes. Furthermore, Pol II occupancy markedly increased near cleavage and polyadenylation sites in xrn1Δ cells, whereas its activity decreased, a characteristic feature of backtracked Pol II. We also provide indirect evidence that Xrn1 is involved in transcription termination downstream of polyadenylation sites. We noted that two additional decay factors, Dhh1 and Lsm1, seem to function similarly to Xrn1 in transcription, perhaps as a complex, and that the decay factors Ccr4 and Rpb4 also perturb transcription in other ways. Interestingly, the decay factors could differentiate between SAGA- and TFIID-dominated promoters. These two classes of genes responded differently to XRN1 deletion in mRNA synthesis and were differentially regulated by mRNA decay pathways, raising the possibility that one distinction between these two gene classes lies in the mechanisms that balance mRNA synthesis with mRNA decay.

Integrative analysis of large-scale loss-of-function screens identifies robust cancer-associated genetic interactions

Genetic interactions, including synthetic lethal effects, can now be systematically identified in cancer cell lines using high-throughput genetic perturbation screens. Despite this advance, few genetic interactions have been reproduced across multiple studies and many appear highly context-specific. Here, by developing a new computational approach, we identified 220 robust driver-gene associated genetic interactions that can be reproduced across independent experiments and across non-overlapping cell line panels. Analysis of these interactions demonstrated that: (i) oncogene addiction effects are more robust than oncogene-related synthetic lethal effects; and (ii) robust genetic interactions are enriched among gene pairs whose protein products physically interact. Exploiting the latter observation, we used a protein-protein interaction network to identify robust synthetic lethal effects associated with passenger gene alterations and validated two new synthetic lethal effects. Our results suggest that protein-protein interaction networks can be used to prioritise therapeutic targets that will be more robust to tumour heterogeneity.

G2G: A web-server for the prediction of human synthetic lethal interactions

Genetic interactions (GIs) are fundamental to our understanding of biological processes in the cell. While GIs have been systematically mapped in yeast, there is scarce information about them in humans. Recently, we have suggested a state-of-the-art hierarchical method that leverages gene ontology information for predicting GIs in yeast. Here, we adapt this method and apply it for the first time to predict GIs in human. We introduce a web service called G2G for this task that is available at http://bnet.cs.tau.ac.il/g2g/.

A Multimodal Genotoxic Anticancer Drug Characterized by Pharmacogenetic Analysis in Caenorhabditis elegans

New anticancer therapeutics require extensive in vivo characterization to identify endogenous and exogenous factors affecting efficacy, to measure toxicity and mutagenicity, and to determine genotypes that result in therapeutic sensitivity or resistance. We used Caenorhabditis elegans as a platform with which to characterize properties of the anticancer therapeutic CX-5461. To understand the processes that respond to CX-5461-induced damage, we generated pharmacogenetic profiles for a panel of C. elegans DNA replication and repair mutants with common DNA-damaging agents for comparison with the profile of CX-5461. We found that multiple repair pathways, including homology-directed repair, microhomology-mediated end joining, nucleotide excision repair, and translesion synthesis, were needed for CX-5461 tolerance. To determine the frequency and spectrum of CX-5461-induced mutations, we used a genetic balancer to capture CX-5461-induced mutations. We found that CX-5461 is mutagenic, resulting in both large copy number variations and a high frequency of single-nucleotide variations (SNVs), which are consistent with the pharmacogenetic profile for CX-5461. Whole-genome sequencing of CX-5461-exposed animals found that CX-5461-induced SNVs exhibited a distinct mutational signature. We also phenocopied the CX-5461 photoreactivity observed in clinical trials and demonstrated that CX-5461 generates reactive oxygen species when exposed to UVA radiation. Together, the data from C. elegans demonstrate that CX-5461 is a multimodal DNA-damaging anticancer agent.

Link synthetic lethality to drug sensitivity of cancer cells

Synthetic lethal (SL) interactions occur when alterations in two genes lead to cell death but alteration in only one of them is not lethal. SL interactions provide a new strategy for molecular-targeted cancer therapy. Currently, there are few drugs targeting SL interactions that entered into clinical trials. Therefore, it is necessary to investigate the link between SL interactions and drug sensitivity of cancer cells systematically for drug development purpose. We identified SL interactions by integrating the high-throughput data from The Cancer Genome Atlas, small hairpin RNA data and genetic interactions of yeast. By integrating SL interactions from other studies, we tested whether the SL pairs that consist of drug target genes and the genes with genomic alterations are related with drug sensitivity of cancer cells. We found that only 6.26%∼34.61% of SL interactions showed the expected significant drug sensitivity using the pooled cancer cell line data from different tissues, but the proportion increased significantly to approximately 90% using the cancer cell line data for each specific tissue. From an independent pharmacogenomics data of 41 breast cancer cell lines, we found three SL interactions (ABL1-IFI16, ABL1-SLC50A1 and ABL1-SYT11) showed significantly better prognosis for the patients with both genes being altered than the patients with only one gene being altered, which partially supports the SL effect between the gene pairs. Our study not only provides a new way for unraveling the complex mechanisms of drug sensitivity but also suggests numerous potentially important drug targets for cancer therapy.

A deep learning model based on sparse auto-encoder for prioritizing cancer-related genes and drug target combinations

Prioritization of cancer-related genes from gene expression profiles and proteomic data is vital to improve the targeted therapies research. Although computational approaches have been complementing high-throughput biological experiments on the understanding of human diseases, it still remains a big challenge to accurately discover cancer-related proteins/genes via automatic learning from large-scale protein/gene expression data and protein-protein interaction data. Most of the existing methods are based on network construction combined with gene expression profiles, which ignore the diversity between normal samples and disease cell lines. In this study, we introduced a deep learning model based on a sparse auto-encoder to learn the specific characteristics of protein interactions in cancer cell lines integrated with protein expression data. The model showed learning ability to identify cancer-related proteins/genes from the input of different protein expression profiles by extracting the characteristics of protein interaction information, which could also predict cancer-related protein combinations. Comparing with other reported methods including differential expression and network-based methods, our model got the highest area under the curve value (>0.8) in predicting cancer-related genes. Our study prioritized ~500 high-confidence cancer-related genes; among these genes, 211 already known cancer drug targets were found, which supported the accuracy of our method. The above results indicated that the proposed auto-encoder model could computationally prioritize candidate proteins/genes involved in cancer and improve the targeted therapies research.

EXP2SL: A Machine Learning Framework for Cell-Line-Specific Synthetic Lethality Prediction

Synthetic lethality (SL), an important type of genetic interaction, can provide useful insight into the target identification process for the development of anticancer therapeutics. Although several well-established SL gene pairs have been verified to be conserved in humans, most SL interactions remain cell-line specific. Here, we demonstrated that the cell-line-specific gene expression profiles derived from the shRNA perturbation experiments performed in the LINCS L1000 project can provide useful features for predicting SL interactions in human. In this paper, we developed a semi-supervised neural network-based method called EXP2SL to accurately identify SL interactions from the L1000 gene expression profiles. Through a systematic evaluation on the SL datasets of three different cell lines, we demonstrated that our model achieved better performance than the baseline methods and verified the effectiveness of using the L1000 gene expression features and the semi-supervise training technique in SL prediction.

Individualized genetic network analysis reveals new therapeutic vulnerabilities in 6,700 cancer genomes

Tumor-specific genomic alterations allow systematic identification of genetic interactions that promote tumorigenesis and tumor vulnerabilities, offering novel strategies for development of targeted therapies for individual patients. We develop an Individualized Network-based Co-Mutation (INCM) methodology by inspecting over 2.5 million nonsynonymous somatic mutations derived from 6,789 tumor exomes across 14 cancer types from The Cancer Genome Atlas. Our INCM analysis reveals a higher genetic interaction burden on the significantly mutated genes, experimentally validated cancer genes, chromosome regulatory factors, and DNA damage repair genes, as compared to human pan-cancer essential genes identified by CRISPR-Cas9 screenings on 324 cancer cell lines. We find that genes involved in the cancer type-specific genetic subnetworks identified by INCM are significantly enriched in established cancer pathways, and the INCM-inferred putative genetic interactions are correlated with patient survival. By analyzing drug pharmacogenomics profiles from the Genomics of Drug Sensitivity in Cancer database, we show that the network-predicted putative genetic interactions (e.g., BRCA2-TP53) are significantly correlated with sensitivity/resistance of multiple therapeutic agents. We experimentally validated that afatinib has the strongest cytotoxic activity on BT474 (IC₅₀ = 55.5 nM, BRCA2 and TP53 co-mutant) compared to MCF7 (IC₅₀ = 7.7 μM, both BRCA2 and TP53 wild type) and MDA-MB-231 (IC₅₀ = 7.9 μM, BRCA2 wild type but TP53 mutant). Finally, drug-target network analysis reveals several potential druggable genetic interactions by targeting tumor vulnerabilities. This study offers a powerful network-based methodology for identification of candidate therapeutic pathways that target tumor vulnerabilities and prioritization of potential pharmacogenomics biomarkers for development of personalized cancer medicine.

Recent efforts to map genetic interactions in tumor cells have suggested that tumor vulnerabilities can be exploited for development of novel targeted therapies. Tumor-specific genomic alterations derived from multi-center cancer genome projects allow identification of genetic interactions that promote tumor vulnerabilities, offering novel strategies for development of targeted cancer therapies. This study develops a novel Individualized Network-based Co-Mutation (termed INCM) methodology for quantifying the putative genetic interactions in cancer. Trained on over 2.5 million nonsynonymous somatic mutations derived from 6,789 tumor exomes across 14 cancer type, we found that genes identified in the cancer type-specific genetic subnetworks were significantly enriched in established cancer pathways. The network-predicted putative genetic interactions are correlated with patient survival. By analyzing drug pharmacogenomics profiles, we showed that the network-predicted putative genetic interactions (e.g., BRCA2-TP53) were significantly correlated with sensitivity/resistance of anticancer drugs (e.g., afatinib) and we experimentally validated it in breast cancer cell lines. Finally, drug-target network analysis reveals several potential druggable genetic interactions (e.g., PIK3CA-PTEN) by targeting tumor vulnerabilities. This study offers a generalizable network-based approach for comprehensive identification of candidate therapeutic pathways that target tumor vulnerabilities and prioritization of potential prognostic and pharmacogenomics biomarkers for development of personalized cancer medicine.

Synthetic lethality: A promising therapeutic strategy for hepatocellular carcinoma

Hepatocellular carcinoma (HCC), the main cause of liver cancer-related death, is one of the main cancers in terms of incidence and mortality. However, HCC is difficult to target and develops strong drug resistance. Therefore, a new treatment strategy is urgently needed. The clinical application of the concept of synthetic lethality in recent years provides a new therapeutic direction for the accurate treatment of HCC. Here, we introduce the concept of synthetic lethality, the screening used to study synthetic lethality, and the identified and potential genetic interactions that induce synthetic lethality in HCC. In addition, we propose opportunities and challenges for translating synthetic lethal interactions to the clinical treatment of HCC.

The mRNA degradation factor Xrn1 regulates transcription elongation in parallel to Ccr4

Co-transcriptional imprinting of mRNA by Rpb4 and Rpb7 subunits of RNA polymerase II (RNAPII) and by the Ccr4-Not complex conditions its post-transcriptional fate. In turn, mRNA degradation factors like Xrn1 are able to influence RNAPII-dependent transcription, making a feedback loop that contributes to mRNA homeostasis. In this work, we have used repressible yeast GAL genes to perform accurate measurements of transcription and mRNA degradation in a set of mutants. This genetic analysis uncovered a link from mRNA decay to transcription elongation. We combined this experimental approach with computational multi-agent modelling and tested different possibilities of Xrn1 and Ccr4 action in gene transcription. This double strategy brought us to conclude that both Xrn1-decaysome and Ccr4-Not regulate RNAPII elongation, and that they do it in parallel. We validated this conclusion measuring TFIIS genome-wide recruitment to elongating RNAPII. We found that xrn1Δ and ccr4Δ exhibited very different patterns of TFIIS versus RNAPII occupancy, which confirmed their distinct role in controlling transcription elongation. We also found that the relative influence of Xrn1 and Ccr4 is different in the genes encoding ribosomal proteins as compared to the rest of the genome.

Synergism Through WEE1 and CHK1 Inhibition in Acute Lymphoblastic Leukemia

Introduction: Screening for synthetic lethality markers has demonstrated that the inhibition of the cell cycle checkpoint kinases WEE1 together with CHK1 drastically affects stability of the cell cycle and induces cell death in rapidly proliferating cells. Exploiting this finding for a possible therapeutic approach has showed efficacy in various solid and hematologic tumors, though not specifically tested in acute lymphoblastic leukemia. Methods: The efficacy of the combination between WEE1 and CHK1 inhibitors in B and T cell precursor acute lymphoblastic leukemia (B/T-ALL) was evaluated in vitro and ex vivo studies. The efficacy of the therapeutic strategy was tested in terms of cytotoxicity, induction of apoptosis, and changes in cell cycle profile and protein expression using B/T-ALL cell lines. In addition, the efficacy of the drug combination was studied in primary B-ALL blasts using clonogenic assays. Results: This study reports, for the first time, the efficacy of the concomitant inhibition of CHK1/CHK2 and WEE1 in ALL cell lines and primary leukemic B-ALL cells using two selective inhibitors: PF-0047736 (CHK1/CHK2 inhibitor) and AZD-1775 (WEE1 inhibitor). We showed strong synergism in the reduction of cell viability, proliferation and induction of apoptosis. The efficacy of the combination was related to the induction of early S-phase arrest and to the induction of DNA damage, ultimately triggering cell death. We reported evidence that the efficacy of the combination treatment is independent from the activation of the p53-p21 pathway. Moreover, gene expression analysis on B-ALL primary samples showed that Chek1 and Wee1 are significantly co-expressed in samples at diagnosis (Pearson r = 0.5770, p = 0.0001) and relapse (Pearson r= 0.8919; p = 0.0001). Finally, the efficacy of the combination was confirmed by the reduction in clonogenic survival of primary leukemic B-ALL cells. Conclusion: Our findings suggest that the combination of CHK1 and WEE1 inhibitors may be a promising therapeutic strategy to be tested in clinical trials for adult ALL.

Mapping genetic interactions in cancer: a road to rational combination therapies

The discovery of synthetic lethal interactions between poly (ADP-ribose) polymerase (PARP) inhibitors and BRCA genes, which are involved in homologous recombination, led to the approval of PARP inhibition as a monotherapy for patients with BRCA1/2-mutated breast or ovarian cancer. Studies following the initial observation of synthetic lethality demonstrated that the reach of PARP inhibitors is well beyond just BRCA1/2 mutants. Insights into the mechanisms of action of anticancer drugs are fundamental for the development of targeted monotherapies or rational combination treatments that will synergize to promote cancer cell death and overcome mechanisms of resistance. The development of targeted therapeutic agents is premised on mapping the physical and functional dependencies of mutated genes in cancer. An important part of this effort is the systematic screening of genetic interactions in a variety of cancer types. Until recently, genetic-interaction screens have relied either on the pairwise perturbations of two genes or on the perturbation of genes of interest combined with inhibition by commonly used anticancer drugs. Here, we summarize recent advances in mapping genetic interactions using targeted, genome-wide, and high-throughput genetic screens, and we discuss the therapeutic insights obtained through such screens. We further focus on factors that should be considered in order to develop a robust analysis pipeline. Finally, we discuss the integration of functional interaction data with orthogonal methods and suggest that such approaches will increase the reach of genetic-interaction screens for the development of rational combination therapies.

A Humanized Yeast Phenomic Model of Deoxycytidine Kinase to Predict Genetic Buffering of Nucleoside Analog Cytotoxicity

Knowledge about synthetic lethality can be applied to enhance the efficacy of anticancer therapies in individual patients harboring genetic alterations in their cancer that specifically render it vulnerable. We investigated the potential for high-resolution phenomic analysis in yeast to predict such genetic vulnerabilities by systematic, comprehensive, and quantitative assessment of drug–gene interaction for gemcitabine and cytarabine, substrates of deoxycytidine kinase that have similar molecular structures yet distinct antitumor efficacy. Human deoxycytidine kinase (dCK) was conditionally expressed in the Saccharomyces cerevisiae genomic library of knockout and knockdown (YKO/KD) strains, to globally and quantitatively characterize differential drug–gene interaction for gemcitabine and cytarabine. Pathway enrichment analysis revealed that autophagy, histone modification, chromatin remodeling, and apoptosis-related processes influence gemcitabine specifically, while drug–gene interaction specific to cytarabine was less enriched in gene ontology. Processes having influence over both drugs were DNA repair and integrity checkpoints and vesicle transport and fusion. Non-gene ontology (GO)-enriched genes were also informative. Yeast phenomic and cancer cell line pharmacogenomics data were integrated to identify yeast–human homologs with correlated differential gene expression and drug efficacy, thus providing a unique resource to predict whether differential gene expression observed in cancer genetic profiles are causal in tumor-specific responses to cytotoxic agents.

Functional RNAi Screens Define Distinct Protein Kinase Vulnerabilities in EGFR-Dependent HNSCC Cell Lines

The inhibitory epidermal growth factor receptor (EGFR) antibody, cetuximab, is an approved therapy for head and neck squamous cell carcinoma (HNSCC). Despite tumor response observed in some HNSCC patients, cetuximab alone or combined with radio- or chemotherapy fails to yield long-term control or cures. We hypothesize that a flexible receptor tyrosine kinase coactivation signaling network supports HNSCC survival in the setting of EGFR blockade, and that drugs disrupting this network will provide superior tumor control when combined with EGFR inhibitors. In this work, we submitted EGFR-dependent HNSCC cell lines to RNA interference-based functional genomics screens to identify, in an unbiased fashion, essential protein kinases for growth and survival as well as synthetic lethal targets for combined inhibition with EGFR antagonists. Mechanistic target of rapamycin kinase (MTOR) and erythroblastosis oncogene B (ERBB)3 were identified as high-ranking essential kinase hits in the HNSCC cell lines. MTOR dependency was confirmed by distinct short hairpin RNAs (shRNAs) and high sensitivity of the cell lines to AZD8055, whereas ERBB3 dependency was validated by shRNA-mediated silencing. Furthermore, a synthetic lethal kinome shRNA screen with a pan-ERBB inhibitor, AZD8931, identified multiple components of the extracellular signal-regulated kinase (ERK) mitogen-activated protein kinase pathway, consistent with ERK reactivation and/or incomplete ERK pathway inhibition in response to EGFR inhibitor monotherapy. As validation, distinct mitogen-activated protein kinase kinase (MEK) inhibitors yielded synergistic growth inhibition when combined with the EGFR inhibitors, gefitinib and AZD8931. The findings identify ERBB3 and MTOR as important pharmacological vulnerabilities in HNSCC and support combining MEK and EGFR inhibitors to enhance clinical efficacy in HNSCC. SIGNIFICANCE STATEMENT: Many cancers are driven by nonmutated receptor tyrosine kinase coactivation networks that defy full inhibition with single targeted drugs. This study identifies erythroblastosis oncogene B (ERBB)3 as an essential protein kinase in epidermal growth factor receptor-dependent head and neck squamous cell cancer (HNSCC) cell lines and a synthetic lethal interaction with the extracellular signal-regulated kinase mitogen-activated protein kinase pathway that provides a rationale for combining pan-ERBB and mitogen-activated protein kinase inhibitors as a therapeutic approach in subsets of HNSCC.

Humanized yeast genetic interaction mapping predicts synthetic lethal interactions of FBXW7 in breast cancer

Background

Synthetic lethal interactions (SLIs) that occur between gene pairs are exploited for cancer therapeutics. Studies in the model eukaryote yeast have identified ~ 550,000 negative genetic interactions that have been extensively studied, leading to characterization of novel pathways and gene functions. This resource can be used to predict SLIs that can be relevant to cancer therapeutics.

Methods

We used patient data to identify genes that are down-regulated in breast cancer. InParanoid orthology mapping was performed to identify yeast orthologs of the down-regulated genes and predict their corresponding SLIs in humans. The predicted network graphs were drawn with Cytoscape. CancerRXgene database was used to predict drug response.

Results

Harnessing the vast available knowledge of yeast genetics, we generated a Humanized Yeast Genetic Interaction Network (HYGIN) for 1009 human genes with 10,419 interactions. Through the addition of patient-data from The Cancer Genome Atlas (TCGA), we generated a breast cancer specific subnetwork. Specifically, by comparing 1009 genes in HYGIN to genes that were down-regulated in breast cancer, we identified 15 breast cancer genes with 130 potential SLIs. Interestingly, 32 of the 130 predicted SLIs occurred with FBXW7, a well-known tumor suppressor that functions as a substrate-recognition protein within a SKP/CUL1/F-Box ubiquitin ligase complex for proteasome degradation. Efforts to validate these SLIs using chemical genetic data predicted that patients with loss of FBXW7 may respond to treatment with drugs like Selumitinib or Cabozantinib.

Conclusions

This study provides a patient-data driven interpretation of yeast SLI data. HYGIN represents a novel strategy to uncover therapeutically relevant cancer drug targets and the yeast SLI data offers a major opportunity to mine these interactions.

Electronic supplementary material

The online version of this article (10.1186/s12920-019-0554-z) contains supplementary material, which is available to authorized users.