PitViper: a software for comparative meta-analysis and annotation of functional screening data.

Recent advancements in shRNA and Cas protein technologies have enabled functional screening methods targeting genes or non-coding regions using single or pooled shRNA and sgRNA. CRISPR-based systems have also been developed for modulating DNA accessibility, resulting in CRISPR-mediated interference (CRISPRi) or activation (CRISPRa) of targeted genes or genomic DNA elements. However, there is still a lack of software tools for integrating diverse array of functional genomics screening outputs that could offer a cohesive framework for comprehensive data integration. Here, we developed PitViper, a flexible and interactive open-source software designed to fill this gap, providing reliable results for the type of elements being screened. It is an end-to-end automated and reproducible bioinformatics pipeline integrating gold-standard methods for functional screening analysis. Our sensitivity analyses demonstrate that PitViper is a useful tool for identifying potential super-enhancer liabilities in a leukemia cell line through genome-wide CRISPRi-based screening. It offers a robust, flexible, and interactive solution for integrating data analysis and reanalysis from functional screening methods, making it a valuable resource for researchers in the field.

Integrative PTEN Enhancer Discovery Reveals a New Model of Enhancer Organization.

Enhancers possess both structural elements mediating promoter looping and functional elements mediating gene expression. Traditional models of enhancer-mediated gene regulation imply genomic overlap or immediate adjacency of these elements. We test this model by combining densely-tiled CRISPRa screening with nucleosome-resolution Region Capture Micro-C topology analysis. Using this integrated approach, we comprehensively define the cis-regulatory landscape for the tumor suppressor PTEN, identifying and validating 10 distinct enhancers and defining their 3D spatial organization. Unexpectedly, we identify several long-range functional enhancers whose promoter proximity is facilitated by chromatin loop anchors several kilobases away, and demonstrate that accounting for this spatial separation improves the computational prediction of validated enhancers. Thus, we propose a new model of enhancer organization incorporating spatial separation of essential functional and structural components.

MCB-613 exploits a collateral sensitivity in drug resistant EGFR-mutant non-small cell lung cancer through covalent inhibition of KEAP1.

Targeted therapies have revolutionized cancer chemotherapy. Unfortunately, most patients develop multifocal resistance to these drugs within a matter of months. Here, we used a high-throughput phenotypic small molecule screen to identify MCB-613 as a compound that selectively targets EGFR-mutant, EGFR inhibitor-resistant non-small cell lung cancer (NSCLC) cells harboring diverse resistance mechanisms. Subsequent proteomic and functional genomic screens involving MCB-613 identified its target in this context to be KEAP1, revealing that this gene is selectively essential in the setting of EGFR inhibitor resistance. In-depth molecular characterization demonstrated that (1) MCB-613 binds KEAP1 covalently; (2) a single molecule of MCB-613 is capable of bridging two KEAP1 monomers together; and, (3) this modification interferes with the degradation of canonical KEAP1 substrates such as NRF2. Surprisingly, NRF2 knockout sensitizes cells to MCB-613, suggesting that the drug functions through modulation of an alternative KEAP1 substrate. Together, these findings advance MCB-613 as a new tool for exploiting the selective essentiality of KEAP1 in drug-resistant, EGFR-mutant NSCLC cells.

Small-molecule targeted therapies induce dependence on DNA double-strand break repair in residual tumor cells.

Residual cancer cells that survive drug treatments with targeted therapies act as a reservoir from which eventual resistant disease emerges. Although there is great interest in therapeutically targeting residual cells, efforts are hampered by our limited knowledge of the vulnerabilities existing in this cell state. Here, we report that diverse oncogene-targeted therapies, including inhibitors of epidermal growth factor receptor (EGFR), anaplastic lymphoma kinase (ALK), KRAS, and BRAF, induce DNA double-strand breaks and, consequently, ataxia-telangiectasia mutated (ATM)-dependent DNA repair in oncogene-matched residual tumor cells. This DNA damage response, observed in cell lines, mouse xenograft models, and human patients, is driven by a pathway involving the activation of caspases 3 and 7 and the downstream caspase-activated deoxyribonuclease (CAD). CAD is, in turn, activated through caspase-mediated degradation of its endogenous inhibitor, ICAD. In models of EGFR mutant non-small cell lung cancer (NSCLC), tumor cells that survive treatment with small-molecule EGFR-targeted therapies are thus synthetically dependent on ATM, and combined treatment with an ATM kinase inhibitor eradicates these cells in vivo. This led to more penetrant and durable responses in EGFR mutant NSCLC mouse xenograft models, including those derived from both established cell lines and patient tumors. Last, we found that rare patients with EGFR mutant NSCLC harboring co-occurring, loss-of-function mutations in ATM exhibit extended progression-free survival on first generation EGFR inhibitor therapy relative to patients with EGFR mutant NSCLC lacking deleterious ATM mutations. Together, these findings establish a rationale for the mechanism-based integration of ATM inhibitors alongside existing targeted therapies.

P2RY2-AKT activation is a therapeutically actionable consequence of XPO1 inhibition in acute myeloid leukemia.

Selinexor is a first-in-class inhibitor of the nuclear exportin XPO1 that was recently approved by the US Food and Drug Administration for the treatment of multiple myeloma and diffuse large B-cell lymphoma. In relapsed/refractory acute myeloid leukemia (AML), selinexor has shown promising activity, suggesting that selinexor-based combination therapies may have clinical potential. Here, motivated by the hypothesis that selinexor's nuclear sequestration of diverse substrates imposes pleiotropic fitness effects on AML cells, we systematically catalog the pro- and anti-fitness consequences of selinexor treatment. We discover that selinexor activates PI3Kγ-dependent AKT signaling in AML by upregulating the purinergic receptor P2RY2. Inhibiting this axis potentiates the anti-leukemic effects of selinexor in AML cell lines, patient-derived primary cultures and multiple mouse models of AML. In a syngeneic, MLL-AF9-driven mouse model of AML, treatment with selinexor and ipatasertib outperforms both standard-of-care chemotherapy and chemotherapy with selinexor. Together, these findings establish drug-induced P2RY2-AKT signaling as an actionable consequence of XPO1 inhibition in AML.

The Folate Cycle Enzyme MTHFR Is a Critical Regulator of Cell Response to MYC-Targeting Therapies.

Deciphering the impact of metabolic intervention on response to anticancer therapy may elucidate a path toward improved clinical responses. Here, we identify amino acid-related pathways connected to the folate cycle whose activation predicts sensitivity to MYC-targeting therapies in acute myeloid leukemia (AML). We establish that folate restriction and deficiency of the rate-limiting folate cycle enzyme MTHFR, which exhibits reduced-function polymorphisms in about 10% of Caucasians, induce resistance to MYC targeting by BET and CDK7 inhibitors in cell lines, primary patient samples, and syngeneic mouse models of AML. Furthermore, this effect is abrogated by supplementation with the MTHFR enzymatic product CH3-THF. Mechanistically, folate cycle disturbance reduces H3K27/K9 histone methylation and activates a SPI1 transcriptional program counteracting the effect of BET inhibition. Our data provide a rationale for screening MTHFR polymorphisms and folate cycle status to nominate patients most likely to benefit from MYC-targeting therapies. SIGNIFICANCE: Although MYC-targeting therapies represent a promising strategy for cancer treatment, evidence of predictors of sensitivity to these agents is limited. We pinpoint that folate cycle disturbance and frequent polymorphisms associated with reduced MTHFR activity promote resistance to BET inhibitors. CH3-THF supplementation thus represents a low-risk intervention to enhance their effects.See related commentary by Marando and Huntly, p. 1791.This article is highlighted in the In This Issue feature, p. 1775.

The creatine kinase pathway is a metabolic vulnerability in EVI1-positive acute myeloid leukemia.

Expression of the MECOM (also known as EVI1) proto-oncogene is deregulated by chromosomal translocations in some cases of acute myeloid leukemia (AML) and is associated with poor clinical outcome. Here, through transcriptomic and metabolomic profiling of hematopoietic cells, we reveal that EVI1 overexpression alters cellular metabolism. A screen using pooled short hairpin RNAs (shRNAs) identified the ATP-buffering, mitochondrial creatine kinase CKMT1 as necessary for survival of EVI1-expressing cells in subjects with EVI1-positive AML. EVI1 promotes CKMT1 expression by repressing the myeloid differentiation regulator RUNX1. Suppression of arginine-creatine metabolism by CKMT1-directed shRNAs or by the small molecule cyclocreatine selectively decreased the viability, promoted the cell cycle arrest and apoptosis of human EVI1-positive cell lines, and prolonged survival in both orthotopic xenograft models and mouse models of primary AML. CKMT1 inhibition altered mitochondrial respiration and ATP production, an effect that was abrogated by phosphocreatine-mediated reactivation of the arginine-creatine pathway. Targeting CKMT1 is thus a promising therapeutic strategy for this EVI1-driven AML subtype that is highly resistant to current treatment regimens.

Targeting MTHFD2 in acute myeloid leukemia.

Drugs targeting metabolism have formed the backbone of therapy for some cancers. We sought to identify new such targets in acute myeloid leukemia (AML). The one-carbon folate pathway, specifically methylenetetrahydrofolate dehydrogenase-cyclohydrolase 2 (MTHFD2), emerged as a top candidate in our analyses. MTHFD2 is the most differentially expressed metabolic enzyme in cancer versus normal cells. Knockdown of MTHFD2 in AML cells decreased growth, induced differentiation, and impaired colony formation in primary AML blasts. In human xenograft and MLL-AF9 mouse leukemia models, MTHFD2 suppression decreased leukemia burden and prolonged survival. Based upon primary patient AML data and functional genomic screening, we determined that FLT3-ITD is a biomarker of response to MTHFD2 suppression. Mechanistically, MYC regulates the expression of MTHFD2, and MTHFD2 knockdown suppresses the TCA cycle. This study supports the therapeutic targeting of MTHFD2 in AML.

Increased SYK activity is associated with unfavorable outcome among patients with acute myeloid leukemia.

Recent discoveries have led to the testing of novel targeted therapies for the treatment of acute myeloid leukemia (AML). To better inform the results of clinical trials, there is a need to identify and systematically assess biomarkers of response and pharmacodynamic markers of successful target engagement. Spleen tyrosine kinase (SYK) is a candidate therapeutic target in AML. Small-molecule inhibitors of SYK induce AML differentiation and impair leukemia progression in preclinical studies. However, tools to predict response to SYK inhibition and to routinely evaluate SYK activation in primary patient samples have been lacking. In this study we quantified phosphorylated SYK (P-SYK) in AML cell lines and establish that increasing levels of baseline P-SYK are correlated with an increasing sensitivity to small-molecule inhibitors targeting SYK. In addition, we found that pharmacological inhibition of SYK activity extinguishes P-SYK expression as detected by an immunohistochemical (IHC) test. Quantitative analysis of P-SYK expression by the IHC test in a series of 70 primary bone marrow biopsy specimens revealed a spectrum of P-SYK expression across AML cases and that high P-SYK expression is associated with unfavourable outcome independent of age, cytogenetics, and white blood cell count. This study thus establishes P-SYK as a critical biomarker in AML that identifies tumors sensitive to SYK inhibition, identifies an at-risk patient population, and allows for the monitoring of target inhibition during treatment.

SYK is a critical regulator of FLT3 in acute myeloid leukemia.

Cooperative dependencies between mutant oncoproteins and wild-type proteins are critical in cancer pathogenesis and therapy resistance. Although spleen tyrosine kinase (SYK) has been implicated in hematologic malignancies, it is rarely mutated. We used kinase activity profiling to identify collaborators of SYK in acute myeloid leukemia (AML) and determined that FMS-like tyrosine kinase 3 (FLT3) is transactivated by SYK via direct binding. Highly activated SYK is predominantly found in FLT3-ITD positive AML and cooperates with FLT3-ITD to activate MYC transcriptional programs. FLT3-ITD AML cells are more vulnerable to SYK suppression than FLT3 wild-type counterparts. In a FLT3-ITD in vivo model, SYK is indispensable for myeloproliferative disease (MPD) development, and SYK overexpression promotes overt transformation to AML and resistance to FLT3-ITD-targeted therapy.

Main principles and outcomes of DNA methylation analysis.

Epigenetic modifications, including DNA methylation, are critically important mediators of normal cell function over the course of our lives. These modifications therefore also can play prominent roles in the development of disorders and diseases, including ovarian cancer. Genome-wide studies are now beginning to comprehensively decipher the methylome in normal and diseased tissues and cells, providing new insights into the distribution, specificity, and magnitude of modifications that occur and raising questions about these changes at specific loci. Further study of these alterations in specific tissues usually involves targeted approaches, of which there are a number available, all with distinct advantages and disadvantages. Here we provide a brief overview of DNA methylation and some of the methylation alterations that have been identified in ovarian cancer, as well as some of the technical approaches used to study these modifications.

Methylation-specific PCR.

Defining DNA methylation patterns in the genome has become essential for understanding diverse biological processes including the regulation of gene expression, imprinted genes, and X chromosome inactivation and how these patterns are deregulated in human diseases. Methylation-specific (MS)-PCR is a useful tool for qualitative DNA methylation analysis with multiple advantages, including ease of design and execution, sensitivity in the ability to detect small quantities of methylated DNA, and the ability to rapidly screen a large number of samples without the need for purchase of expensive laboratory equipment. This assay requires modification of the genomic DNA by sodium bisulfite and two independent primer sets for PCR amplification, one pair designed to recognize the methylated and the other pair the unmethylated versions of the bisulfite-modified sequence. The amplicons are visualized using ethidium bromide staining following agarose gel electrophoresis. Amplicons of the expected size produced from either primer pair are indicative of the presence of DNA in the original sample with the respective methylation status.

Bisulfite sequencing of cloned alleles.

Bisulfite sequencing of cloned alleles is a widely used method for capturing the methylation profiles of single alleles. This method combines PCR amplification of the bisulfite-modified DNA with the subcloning of the amplicons into plasmids followed by transformation into bacteria and plating on selective media. The resulting colony forming units are each comprised of bacterial clones containing the same plasmid reflecting a single allele in the original PCR reaction. Following whole cell PCR and sequencing, the results provide highly detailed information about the status of each CG site within an allele. Sequencing of a large number of individual clones can provide quantitative information, assuming unbiased PCR, subcloning and clone selection. The proportion of methylated cytosine at a particular position within the sequenced alleles can be determined by counting the number of alleles showing methylation at the position of interest and dividing this by the total number of clones sequenced.

Bisulfite pyrosequencing.

Bisulfite pyrosequencing is a sequencing-by-synthesis method used to quantitatively determine the methylation of individual CG cytosines from PCR amplicons of a region up to 115 bases in length. The procedure relies on prior bisulfite conversion of all potentially methylated CG cytosines to either cytosine (methylated) or thymine (unmethylated) and involves the stepwise incorporation of deoxynucleotide triphosphates into the growing strand of nascent DNA. The incorporation of these dNTPs results in the proportional release of pyrophosphate, which is converted into ATP to aid in a subsequent conversion of luciferin to oxyluciferin. The amount of light released in the process is proportional to the number of nucleotides incorporated, and the procedure provides a quantitative portrait of the methylation profile for the amplicon in question.

Targeting MYCN in neuroblastoma by BET bromodomain inhibition.

Bromodomain inhibition comprises a promising therapeutic strategy in cancer, particularly for hematologic malignancies. To date, however, genomic biomarkers to direct clinical translation have been lacking. We conducted a cell-based screen of genetically defined cancer cell lines using a prototypical inhibitor of BET bromodomains. Integration of genetic features with chemosensitivity data revealed a robust correlation between MYCN amplification and sensitivity to bromodomain inhibition. We characterized the mechanistic and translational significance of this finding in neuroblastoma, a childhood cancer with frequent amplification of MYCN. Genome-wide expression analysis showed downregulation of the MYCN transcriptional program accompanied by suppression of MYCN transcription. Functionally, bromodomain-mediated inhibition of MYCN impaired growth and induced apoptosis in neuroblastoma. BRD4 knockdown phenocopied these effects, establishing BET bromodomains as transcriptional regulators of MYCN. BET inhibition conferred a significant survival advantage in 3 in vivo neuroblastoma models, providing a compelling rationale for developing BET bromodomain inhibitors in patients with neuroblastoma.

Selective HDAC1/HDAC2 inhibitors induce neuroblastoma differentiation.

While cytotoxic chemotherapy remains the hallmark of cancer treatment, intensive regimens fall short in many malignancies, including high-risk neuroblastoma. One alternative strategy is to therapeutically promote tumor differentiation. We created a gene expression signature to measure neuroblast maturation, adapted it to a high-throughput platform, and screened a diversity oriented synthesis-generated small-molecule library for differentiation inducers. We identified BRD8430, containing a nine-membered lactam, an ortho-amino anilide functionality, and three chiral centers, as a selective class I histone deacetylase (HDAC) inhibitor (HDAC1 > 2 > 3). Further investigation demonstrated that selective HDAC1/HDAC2 inhibition using compounds or RNA interference induced differentiation and decreased viability in neuroblastoma cell lines. Combined treatment with 13-cis retinoic acid augmented these effects and enhanced activation of retinoic acid signaling. Therefore, by applying a chemical genomic screening approach, we identified selective HDAC1/HDAC2 inhibition as a strategy to induce neuroblastoma differentiation.