Functional genomics reveals an off-target dependency of drug synergy in gastric cancer therapy.

BACKGROUND: Integrating molecular-targeted agents into combination chemotherapy is transformative for enhancing treatment outcomes in cancer. However, realizing the full potential of this approach requires a clear comprehension of the genetic dependencies underlying drug synergy. While the interactions between conventional chemotherapeutics are well-explored, the interplay of molecular-targeted agents with conventional chemotherapeutics remains a frontier in cancer treatment. Hence, we leveraged a powerful functional genomics approach to decode genomic dependencies that drive synergy in molecular-targeted agent/chemotherapeutic combinations in gastric adenocarcinoma, addressing a critical need in gastric cancer therapy. METHODS: We screened pharmacological interactions between fifteen molecular-targeted agent/conventional chemotherapeutic pairs in gastric adenocarcinoma cells, and examined the genome-scale genetic dependencies of synergy integrating genome-wide CRISPR screening with the shRNA-based signature assay. We validated the synergy in cell death using fluorescence-based and lysis-dependent inference of cell death kinetics assay, and validated the genetic dependencies by single-gene knockout experiments. RESULTS: Our combination screen identified SN-38/erlotinib as the drug pair with the strongest synergism. Functional genomics assays unveiled a genetic dependency signature of SN-38/erlotinib identical to SN-38. Remarkably, the enhanced cell death with improved kinetics induced by SN-38/erlotinib was attributed to erlotinib's off-target effect, inhibiting ABCG2, rather than its on-target effect on EGFR. CONCLUSION: In the era of precision medicine, where emphasis on primary drug targets prevails, our research challenges this paradigm by showcasing a robust synergy underpinned by an off-target dependency. Further dissection of the intricate genetic dependencies that underlie synergy can pave the way to developing more effective combination strategies in gastric cancer therapy.

Simultaneous screening of overexpressed genes in breast cancer for oncogenic drivers and tumor dependencies.

There are hundreds of genes typically overexpressed in breast cancer cells and it's often assumed that their overexpression contributes to cancer progression. However, the precise proportion of these overexpressed genes contributing to tumorigenicity remains unclear. To address this gap, we undertook a comprehensive screening of a diverse set of seventy-two genes overexpressed in breast cancer. This systematic screening evaluated their potential for inducing malignant transformation and, concurrently, assessed their impact on breast cancer cell proliferation and viability. Select genes including ALDH3B1, CEACAM5, IL8, PYGO2, and WWTR1, exhibited pronounced activity in promoting tumor formation and establishing gene dependencies critical for tumorigenicity. Subsequent investigations revealed that CEACAM5 overexpression triggered the activation of signaling pathways involving ß-catenin, Cdk4, and mTOR. Additionally, it conferred a growth advantage independent of exogenous insulin in defined medium and facilitated spheroid expansion by inducing multiple layers of epithelial cells while preserving a hollow lumen. Furthermore, the silencing of CEACAM5 expression synergized with tamoxifen-induced growth inhibition in breast cancer cells. These findings underscore the potential of screening overexpressed genes for both oncogenic drivers and tumor dependencies to expand the repertoire of therapeutic targets for breast cancer treatment.

Leukemia circulation kinetics revealed through blood exchange method.

Leukemias and their bone marrow microenvironments undergo dynamic changes over the course of disease. However, little is known about the circulation kinetics of leukemia cells, nor the impact of specific factors on the clearance of circulating leukemia cells (CLCs) from the blood. To gain a basic understanding of CLC dynamics over the course of disease progression and therapeutic response, we apply a blood exchange method to mouse models of acute leukemia. We find that CLCs circulate in the blood for 1-2 orders of magnitude longer than solid tumor circulating tumor cells. We further observe that: (i) leukemia presence in the marrow can limit the clearance of CLCs in a model of acute lymphocytic leukemia (ALL), and (ii) CLCs in a model of relapsed acute myeloid leukemia (AML) can clear faster than their untreated counterparts. Our approach can also directly quantify the impact of microenvironmental factors on CLC clearance properties. For example, data from two leukemia models suggest that E-selectin, a vascular adhesion molecule, alters CLC clearance. Our research highlights that clearance rates of CLCs can vary in response to tumor and treatment status and provides a strategy for identifying basic processes and factors that govern the kinetics of circulating cells.

Inherent genome instability underlies trisomy 21-associated myeloid malignancies.

Constitutional trisomy 21 (T21) is a state of aneuploidy associated with high incidence of childhood acute myeloid leukemia (AML). T21-associated AML is preceded by transient abnormal myelopoiesis (TAM), which is triggered by truncating mutations in GATA1 generating a short GATA1 isoform (GATA1s). T21-associated AML emerges due to secondary mutations in hematopoietic clones bearing GATA1s. Since aneuploidy generally impairs cellular fitness, the paradoxically elevated risk of myeloid malignancy in T21 is not fully understood. We hypothesized that individuals with T21 bear inherent genome instability in hematopoietic lineages that promotes leukemogenic mutations driving the genesis of TAM and AML. We found that individuals with T21 show increased chromosomal copy number variations (CNVs) compared to euploid individuals, suggesting that genome instability could be underlying predisposition to TAM and AML. Acquisition of GATA1s enforces myeloid skewing and maintenance of the hematopoietic progenitor state independently of T21; however, GATA1s in T21 hematopoietic progenitor cells (HPCs) further augments genome instability. Increased dosage of the chromosome 21 (chr21) gene DYRK1A impairs homology-directed DNA repair as a mechanism of elevated mutagenesis. These results posit a model wherein inherent genome instability in T21 drives myeloid malignancy in concert with GATA1s mutations.

Leukemia-intrinsic determinants of CAR-T response revealed by iterative in vivo genome-wide CRISPR screening.

CAR-T therapy is a promising, novel treatment modality for B-cell malignancies and yet many patients relapse through a variety of means, including loss of CAR-T cells and antigen escape. To investigate leukemia-intrinsic CAR-T resistance mechanisms, we performed genome-wide CRISPR-Cas9 loss-of-function screens in an immunocompetent murine model of B-cell acute lymphoblastic leukemia (B-ALL) utilizing a modular guide RNA library. We identified IFNγR/JAK/STAT signaling and components of antigen processing and presentation pathway as key mediators of resistance to CAR-T therapy in vivo; intriguingly, loss of this pathway yielded the opposite effect in vitro (sensitized leukemia to CAR-T cells). Transcriptional characterization of this model demonstrated upregulation of these pathways in tumors relapsed after CAR-T treatment, and functional studies showed a surprising role for natural killer (NK) cells in engaging this resistance program. Finally, examination of data from B-ALL patients treated with CAR-T revealed an association between poor outcomes and increased expression of JAK/STAT and MHC-I in leukemia cells. Overall, our data identify an unexpected mechanism of resistance to CAR-T therapy in which tumor cell interaction with the in vivo tumor microenvironment, including NK cells, induces expression of an adaptive, therapy-induced, T-cell resistance program in tumor cells.

REV7 Monomer Is Unable to Participate in Double Strand Break Repair and Translesion Synthesis but Suppresses Mitotic Errors.

Rev7 is a regulatory protein with roles in translesion synthesis (TLS), double strand break (DSB) repair, replication fork protection, and cell cycle regulation. Rev7 forms a homodimer in vitro using its HORMA (Hop, Rev7, Mad2) domain; however, the functional importance of Rev7 dimerization has been incompletely understood. We analyzed the functional properties of cells expressing either wild-type mouse Rev7 or Rev7K44A/R124A/A135D, a mutant that cannot dimerize. The expression of wild-type Rev7, but not the mutant, rescued the sensitivity of Rev7-/- cells to X-rays and several alkylating agents and reversed the olaparib resistance phenotype of Rev7-/- cells. Using a novel fluorescent host-cell reactivation assay, we found that Rev7K44A/R124A/A135D is unable to promote gap-filling TLS opposite an abasic site analog. The Rev7 dimerization interface is also required for shieldin function, as both Rev7-/- cells and Rev7-/- cells expressing Rev7K44A/R124A/A135D exhibit decreased proficiency in rejoining some types of double strand breaks, as well as increased homologous recombination. Interestingly, Rev7K44A/R124A/A135D retains some function in cell cycle regulation, as it maintains an interaction with Ras-related nuclear protein (Ran) and partially rescues the formation of micronuclei. The mutant Rev7 also rescues the G2/M accumulation observed in Rev7-/- cells but does not affect progression through mitosis following nocodazole release. We conclude that while Rev7 dimerization is required for its roles in TLS, DSB repair, and regulation of the anaphase promoting complex, dimerization is at least partially dispensable for promoting mitotic spindle assembly through its interaction with Ran.

Functional genomics reveals an off-target dependency of drug synergy in gastric cancer therapy.

The rational combination of anticancer agents is critical to improving patient outcomes in cancer. Nonetheless, most combination regimens in the clinic result from empirical methodologies disregarding insight into the mechanism of action and missing the opportunity to improve therapy outcomes incrementally. Deciphering the genetic dependencies and vulnerabilities responsible for synergistic interactions is crucial for rationally developing effective anticancer drug combinations. Hence, we screened pairwise pharmacological interactions between molecular-targeted agents and conventional chemotherapeutics and examined the genome-scale genetic dependencies in gastric adenocarcinoma cell models. Since this type of cancer is mainly chemoresistant and incurable, clinical situations demand effective combination strategies. Our pairwise combination screen revealed SN38/erlotinib as the drug pair with the most robust synergism. Genome-wide CRISPR screening and a shRNA-based signature assay indicated that the genetic dependency/vulnerability signature of SN38/erlotinib is the same as SN38 alone. Additional investigation revealed that the enhanced cell death with improved death kinetics caused by the SN38/erlotinib combination is surprisingly due to erlotinib's off-target effect that inhibits ABCG2 but not its on-target effect on EGFR. Our results confirm that a genetic dependency signature different from the single-drug application may not be necessary for the synergistic interaction of molecular-targeted agents with conventional chemotherapeutics in gastric adenocarcinoma. The findings also demonstrated the efficacy of functional genomics approaches in unveiling biologically validated mechanisms of pharmacological interactions.

Leukemia circulation kinetics revealed through blood exchange method.

Leukemias and their bone marrow microenvironment are known to undergo dynamic changes over the course of disease. However, relatively little is known about the circulation kinetics of leukemia cells, nor the impact of specific factors on the clearance of circulating leukemia cells (CLCs) from the blood. To gain a basic understanding of leukemia cell dynamics over the course of disease progression and therapeutic response, we apply a blood exchange method to mouse models of acute leukemia. We find that CLCs circulate in the blood for 1-2 orders of magnitude longer than solid tumor circulating tumor cells. We further observe that: i) leukemia presence in the marrow can limit the clearance of CLCs in a model of acute lymphocytic leukemia (ALL), and ii) CLCs in a model of relapsed acute myeloid leukemia (AML) can clear faster than their untreated counterparts. Our approach can also directly quantify the impact of microenvironmental factors on CLC clearance properties. For example, data from two leukemia models suggest that E-selectin, a vascular adhesion molecule, alters CLC clearance. Our research highlights that clearance rates of CLCs can vary in response to tumor and treatment status and provides a strategy for identifying basic processes and factors that govern the kinetics of circulating cells.

Integrated multi-omics analyses reveal homology-directed repair pathway as a unique dependency in near-haploid leukemia.

Whole chromosome losses resulting in near-haploid karyotypes are found in a rare subgroup of treatment-refractory acute lymphoblastic leukemia. To systematically dissect the unique physiology and uncover susceptibilities that can be exploited in near-haploid leukemia, we leveraged single-cell RNA-Seq and computational inference of cell cycle stages to pinpoint key differences between near-haploid and diploid leukemia cells. Combining cell cycle stage-specific differential expression with gene essentiality scores from a genome-wide CRISPR-Cas9-mediated knockout screen, we identified the homologous recombination pathway component RAD51B as an essential gene in near-haploid leukemia. DNA damage analyses revealed significantly increased sensitivity of RAD51-mediated repair to RAD51B loss in the G2/M stage of near-haploid cells, suggesting a unique role of RAD51B in the homologous recombination pathway. Elevated G2/M and G1/S checkpoint signaling was part of a RAD51B signature expression program in response to chemotherapy in a xenograft model of human near-haploid B-ALL, and RAD51B and its associated programs were overexpressed in a large panel of near-haploid B-ALL patients. These data highlight a unique genetic dependency on DNA repair machinery in near-haploid leukemia and demarcate RAD51B as a promising candidate for targeted therapy in this treatment-resistant disease.

Screening for CD19-specific chimaeric antigen receptors with enhanced signalling via a barcoded library of intracellular domains.

The immunostimulatory intracellular domains (ICDs) of chimaeric antigen receptors (CARs) are essential for converting antigen recognition into antitumoural function. Although there are many possible combinations of ICDs, almost all current CARs rely on combinations of CD3𝛇, CD28 and 4-1BB. Here we show that a barcoded library of 700,000 unique CD19-specific CARs with diverse ICDs cloned into lentiviral vectors and transduced into Jurkat T cells can be screened at high throughput via cell sorting and next-generation sequencing to optimize CAR signalling for antitumoural functions. By using this screening approach, we identified CARs with new ICD combinations that, compared with clinically available CARs, endowed human primary T cells with comparable tumour control in mice and with improved proliferation, persistence, exhaustion and cytotoxicity after tumour rechallenge in vitro. The screening strategy can be adapted to other disease models, cell types and selection conditions, and could be used to improve adoptive cell therapies and to expand their utility to new disease indications.

Collateral responses to classical cytotoxic chemotherapies are heterogeneous and sensitivities are sparse.

Chemotherapy resistance is a major obstacle to curing cancer patients. Combination drug regimens have shown promise as a method to overcome resistance; however, to date only some cancers have been cured with this method. Collateral sensitivity-the phenomenon whereby resistance to one drug is co-occurrent with sensitivity to a second drug-has been gaining traction as a promising new concept to guide rational design of combination regimens. Here we evolved over 100 subclones of the Eµ-Myc; p19ARF-/- cell line to be resistant to one of four classical chemotherapy agents: doxorubicin, vincristine, paclitaxel, and cisplatin. We then surveyed collateral responses to acquisition of resistance to these agents. Although numerous collateral sensitivities have been documented for antibiotics and targeted cancer therapies, we observed only one collateral sensitivity: half of cell lines that acquired resistance to paclitaxel also acquired a collateral sensitivity to verapamil. However, we found that the mechanism of this collateral sensitivity was unrelated to the mechanism of paclitaxel resistance. Interestingly, we observed heterogeneity in the phenotypic response to acquisition of resistance to most of the drugs we tested, most notably for paclitaxel, suggesting the existence of multiple different states of resistance. Surprisingly, this phenotypic heterogeneity in paclitaxel resistant cell lines was unrelated to transcriptomic heterogeneity among those cell lines. These features of phenotypic and transcriptomic heterogeneity must be taken into account in future studies of treated tumor subclones and in design of chemotherapy combinations.

Microenvironmental IL-6 inhibits anti-cancer immune responses generated by cytotoxic chemotherapy.

Cytotoxic chemotherapeutics primarily function through DNA damage-induced tumor cell apoptosis, although the inflammation provoked by these agents can stimulate anti-cancer immune responses. The mechanisms that control these distinct effects and limit immunogenic responses to DNA-damage mediated cell death in vivo are currently unclear. Using a mouse model of BCR-ABL+ B-cell acute lymphoblastic leukemia, we show that chemotherapy-induced anti-cancer immunity is suppressed by the tumor microenvironment through production of the cytokine IL-6. The chemotherapeutic doxorubicin is curative in IL-6-deficient mice through the induction of CD8+ T-cell-mediated anti-cancer responses, while moderately extending lifespan in wild type tumor-bearing mice. We also show that IL-6 suppresses the effectiveness of immune-checkpoint inhibition with anti-PD-L1 blockade. Our results suggest that IL-6 is a key regulator of anti-cancer immune responses induced by genotoxic stress and that its inhibition can switch cancer cell clearance from primarily apoptotic to immunogenic, promoting and maintaining durable anti-tumor immune responses.

Acquired resistance to PRMT5 inhibition induces concomitant collateral sensitivity to paclitaxel.

Epigenetic regulators play key roles in cancer and are increasingly being targeted for treatment. However, for many, little is known about mechanisms of resistance to the inhibition of these regulators. We have generated a model of resistance to inhibitors of protein arginine methyltransferase 5 (PRMT5). This study was conducted in KrasG12D;Tp53-null lung adenocarcinoma (LUAD) cell lines. Resistance to PRMT5 inhibitors (PRMT5i) arose rapidly, and barcoding experiments showed that this resulted from a drug-induced transcriptional state switch, not selection of a preexisting population. This resistant state is both stable and conserved across variants arising from distinct LUAD lines. Moreover, it brought with it vulnerabilities to other chemotherapeutics, especially the taxane paclitaxel. This paclitaxel sensitivity depended on the presence of stathmin 2 (STMN2), a microtubule regulator that is specifically expressed in the resistant state. Remarkably, STMN2 was also essential for resistance to PRMT5 inhibition. Thus, a single gene is required for both acquisition of resistance to PRMT5i and collateral sensitivity to paclitaxel in our LUAD cells. Accordingly, the combination of PRMT5i and paclitaxel yielded potent and synergistic killing of the murine LUAD cells. Importantly, the synergy between PRMT5i and paclitaxel also extended to human cancer cell lines. Finally, analysis of The Cancer Genome Atlas patient data showed that high STMN2 levels correlate with complete regression of tumors in response to taxane treatment. Collectively, this study reveals a recurring mechanism of PRMT5i resistance in LUAD and identifies collateral sensitivities that have potential clinical relevance.

Mitochondrial apoptotic priming is a key determinant of cell fate upon p53 restoration.

Reactivation of p53 in established tumors typically results in one of two cell fates, cell cycle arrest or apoptosis, but it remains unclear how this cell fate is determined. We hypothesized that high mitochondrial priming prior to p53 reactivation would lead to apoptosis, while low priming would lead to survival and cell cycle arrest. Using a panel of Kras-driven, p53 restorable cell lines derived from genetically engineered mouse models of lung adenocarcinoma and sarcoma (both of which undergo cell cycle arrest upon p53 restoration), as well as lymphoma (which instead undergo apoptosis), we show that the level of mitochondrial apoptotic priming is a critical determinant of p53 reactivation outcome. Cells with high initial priming (e.g., lymphomas) lacked sufficient reserve antiapoptotic capacity and underwent apoptosis after p53 restoration. Forced BCL-2 or BCL-XL expression reduced priming and resulted in survival and cell cycle arrest. Cells with low initial priming (e.g., lung adenocarcinoma and sarcoma) survived and proceeded to arrest in the cell cycle. When primed by inhibition of their antiapoptotic proteins using genetic (BCL-2 or BCL-XL deletion or BAD overexpression) or pharmacologic (navitoclax) means, apoptosis resulted upon p53 restoration in vitro and in vivo. These data demonstrate that mitochondrial apoptotic priming is a key determining factor of cell fate upon p53 activation. Moreover, it is possible to enforce apoptotic cell fate following p53 activation in less primed cells using p53-independent drugs that increase apoptotic priming, including BH3 mimetic drugs.

Targeting acute myeloid leukemia dependency on VCP-mediated DNA repair through a selective second-generation small-molecule inhibitor.

The development and survival of cancer cells require adaptive mechanisms to stress. Such adaptations can confer intrinsic vulnerabilities, enabling the selective targeting of cancer cells. Through a pooled in vivo short hairpin RNA (shRNA) screen, we identified the adenosine triphosphatase associated with diverse cellular activities (AAA-ATPase) valosin-containing protein (VCP) as a top stress-related vulnerability in acute myeloid leukemia (AML). We established that AML was the most responsive disease to chemical inhibition of VCP across a panel of 16 cancer types. The sensitivity to VCP inhibition of human AML cell lines, primary patient samples, and syngeneic and xenograft mouse models of AML was validated using VCP-directed shRNAs, overexpression of a dominant-negative VCP mutant, and chemical inhibition. By combining mass spectrometry-based analysis of the VCP interactome and phospho-signaling studies, we determined that VCP is important for ataxia telangiectasia mutated (ATM) kinase activation and subsequent DNA repair through homologous recombination in AML. A second-generation VCP inhibitor, CB-5339, was then developed and characterized. Efficacy and safety of CB-5339 were validated in multiple AML models, including syngeneic and patient-derived xenograft murine models. We further demonstrated that combining DNA-damaging agents, such as anthracyclines, with CB-5339 treatment synergizes to impair leukemic growth in an MLL-AF9-driven AML murine model. These studies support the clinical testing of CB-5339 as a single agent or in combination with standard-of-care DNA-damaging chemotherapy for the treatment of AML.

Targeting the De Novo Purine Synthesis Pathway Through Adenylosuccinate Lyase Depletion Impairs Liver Cancer Growth by Perturbing Mitochondrial Function.

BACKGROUND AND AIMS: Hepatocellular carcinoma (HCC) is among the most common cancer types worldwide, yet patients with HCC have limited treatment options. There is an urgent need to identify drug targets that specifically inhibit the growth of HCC cells. APPROACH AND RESULTS: We used a CRISPR library targeting ~2,000 druggable genes to perform a high-throughput screen and identified adenylosuccinate lyase (ADSL), a key enzyme involved in the de novo purine synthesis pathway, as a potential drug target for HCC. ADSL has been implicated as a potential oncogenic driver in some cancers, but its role in liver cancer progression remains unknown. CRISPR-mediated knockout of ADSL impaired colony formation of liver cancer cells by affecting AMP production. In the absence of ADSL, the growth of liver tumors is retarded in vivo. Mechanistically, we found that ADSL knockout caused S-phase cell cycle arrest not by inducing DNA damage but by impairing mitochondrial function. Using data from patients with HCC, we also revealed that high ADSL expression occurs during tumorigenesis and is linked to poor survival rate. CONCLUSIONS: Our findings uncover the role of ADSL-mediated de novo purine synthesis in fueling mitochondrial ATP production to promote liver cancer cell growth. Targeting ADSL may be a therapeutic approach for patients with HCC.

REV1 inhibitor JH-RE-06 enhances tumor cell response to chemotherapy by triggering senescence hallmarks.

REV1/POLζ-dependent mutagenic translesion synthesis (TLS) promotes cell survival after DNA damage but is responsible for most of the resulting mutations. A novel inhibitor of this pathway, JH-RE-06, promotes cisplatin efficacy in cancer cells and mouse xenograft models, but the mechanism underlying this combinatorial effect is not known. We report that, unexpectedly, in two different mouse xenograft models and four human and mouse cell lines we examined in vitro cisplatin/JH-RE-06 treatment does not increase apoptosis. Rather, it increases hallmarks of senescence such as senescence-associated ß-galactosidase, increased p21 expression, micronuclei formation, reduced Lamin B1, and increased expression of the immune regulators IL6 and IL8 followed by cell death. Moreover, although p-γ-H2AX foci formation was elevated and ATR expression was low in single agent cisplatin-treated cells, the opposite was true in cells treated with cisplatin/JH-RE-06. These observations suggest that targeting REV1 with JH-RE-06 profoundly affects the nature of the persistent genomic damage after cisplatin treatment and also the resulting physiological responses. These data highlight the potential of REV1/POLζ inhibitors to alter the biological response to DNA-damaging chemotherapy and enhance the efficacy of chemotherapy.

Rev7 loss alters cisplatin response and increases drug efficacy in chemotherapy-resistant lung cancer.

Cisplatin is a standard of care for lung cancer, yet platinum therapy rarely results in substantial tumor regression or a dramatic extension in patient survival. Here, we examined whether targeting Rev7 (also referred to as Mad2B, Mad2L2, and FANCV), a component of the translesion synthesis (TLS) machinery, could potentiate the action of cisplatin in non-small cell lung cancer (NSCLC) treatment. Rev7 loss led to an enhanced tumor cell sensitivity to cisplatin and dramatically improved chemotherapeutic response in a highly drug-resistant mouse model of NSCLC. While cisplatin monotherapy resulted in tumor cell apoptosis, Rev7 deletion promoted a cisplatin-induced senescence phenotype. Moreover, Rev7 deficiency promoted greater cisplatin sensitivity than that previously shown following targeting of other Pol ζ-proteins, suggesting that Pol ζ-dependent and -independent roles of Rev7 are relevant to cisplatin response. Thus, targeting Rev7 may represent a unique strategy for altering and enhancing chemotherapeutic response.

Integrated regulatory models for inference of subtype-specific susceptibilities in glioblastoma.

Glioblastoma multiforme (GBM) is a highly malignant form of cancer that lacks effective treatment options or well-defined strategies for personalized cancer therapy. The disease has been stratified into distinct molecular subtypes; however, the underlying regulatory circuitry that gives rise to such heterogeneity and its implications for therapy remain unclear. We developed a modular computational pipeline, Integrative Modeling of Transcription Regulatory Interactions for Systematic Inference of Susceptibility in Cancer (inTRINSiC), to dissect subtype-specific regulatory programs and predict genetic dependencies in individual patient tumors. Using a multilayer network consisting of 518 transcription factors (TFs), 10,733 target genes, and a signaling layer of 3,132 proteins, we were able to accurately identify differential regulatory activity of TFs that shape subtype-specific expression landscapes. Our models also allowed inference of mechanisms for altered TF behavior in different GBM subtypes. Most importantly, we were able to use the multilayer models to perform an in silico perturbation analysis to infer differential genetic vulnerabilities across GBM subtypes and pinpoint the MYB family member MYBL2 as a drug target specific for the Proneural subtype.