Search

Supplementary Data from Exquisite Sensitivity to Dual BRG1/BRM ATPase Inhibitors Reveals Broad SWI/SNF Dependencies in Acute Myeloid Leukemia

Various subunits of mammalian SWI/SNF chromatin remodeling complexes display loss-of-function mutations characteristic of tumor suppressors in different cancers, but an additional role for SWI/SNF supporting cell survival in distinct cancer contexts is emerging. In particular, genetic dependence on the catalytic subunit BRG1/SMARCA4 has been observed in acute myelogenous leukemia (AML), yet the feasibility of direct therapeutic targeting of SWI/SNF catalytic activity in leukemia remains unknown. Here, we evaluated the activity of dual BRG1/BRM ATPase inhibitors across a genetically diverse panel of cancer cell lines and observed that hematopoietic cancer cell lines were among the most sensitive compared with other lineages. This result was striking in comparison with data from pooled short hairpin RNA screens, which showed that only a subset of leukemia cell lines display sensitivity to BRG1 knockdown. We demonstrate that combined genetic knockdown of BRG1 and BRM is required to recapitulate the effects of dual inhibitors, suggesting that SWI/SNF dependency in human leukemia extends beyond a predominantly BRG1-driven mechanism. Through gene expression and chromatin accessibility studies, we show that the dual inhibitors act at genomic loci associated with oncogenic transcription factors, and observe a downregulation of leukemic pathway genes, including MYC, a well-established target of BRG1 activity in AML. Overall, small-molecule inhibition of BRG1/BRM induced common transcriptional responses across leukemia models resulting in a spectrum of cellular phenotypes. Implications: Our studies reveal the breadth of SWI/SNF dependency in leukemia and support targeting SWI/SNF catalytic function as a potential therapeutic strategy in AML.

Purpose: Targeting RAF for antitumor therapy in RAS-mutant tumors holds promise. Herein, we describe in detail novel properties of the type II RAF inhibitor, LXH254. Experimental Design: LXH254 was profiled in biochemical, in vitro, and in vivo assays, including examining the activities of the drug in a large panel of cancer-derived cell lines and a comprehensive set of in vivo models. In addition, activity of LXH254 was assessed in cells where different sets of RAF paralogs were ablated, or that expressed kinase-impaired and dimer-deficient variants of ARAF. Results: We describe an unexpected paralog selectivity of LXH254, which is able to potently inhibit BRAF and CRAF, but has less activity against ARAF. LXH254 was active in models harboring BRAF alterations, including atypical BRAF alterations coexpressed with mutant K/NRAS, and NRAS mutants, but had only modest activity in KRAS mutants. In RAS-mutant lines, loss of ARAF, but not BRAF or CRAF, sensitized cells to LXH254. ARAF-mediated resistance to LXH254 required both kinase function and dimerization. Higher concentrations of LXH254 were required to inhibit signaling in RAS-mutant cells expressing only ARAF relative to BRAF or CRAF. Moreover, specifically in cells expressing only ARAF, LXH254 caused paradoxical activation of MAPK signaling in a manner similar to dabrafenib. Finally, in vivo, LXH254 drove complete regressions of isogenic variants of RAS-mutant cells lacking ARAF expression, while parental lines were only modestly sensitive. Conclusions: LXH254 is a novel RAF inhibitor, which is able to inhibit dimerized BRAF and CRAF, as well as monomeric BRAF, while largely sparing ARAF.

Chemogenetic libraries, collections of well-defined chemical probes, provide tremendous value to biomedical research but require substantial effort to ensure diversity as well as quality of the contents. We have assembled a chemogenetic library by data mining and crowdsourcing institutional expertise. We are sharing our approach, lessons learned, and disclosing our current collection of 4,185 compounds with their primary annotated gene targets (https://github.com/Novartis/MoaBox). This physical collection is regularly updated and used broadly both within Novartis and in collaboration with external partners.

Three limonoid natural products with selective anti-proliferative activity against BRAF(V600E) and NRAS(Q61K)-mutation-dependent melanoma cell lines were identified. Differential transcriptome analysis revealed dependency of compound activity on expression of the mitochondrial cytochrome P450 oxidase CYP27A1, a transcriptional target of melanogenesis-associated transcription factor (MITF). We determined that CYP27A1 activity is necessary for the generation of a reactive metabolite that proceeds to inhibit cellular proliferation. A genome-wide small interfering RNA screen in combination with chemical proteomics experiments revealed gene-drug functional epistasis, suggesting that these compounds target mitochondrial biogenesis and inhibit tumor bioenergetics through a covalent mechanism. Our work suggests a strategy for melanoma-specific targeting by exploiting the expression of MITF target gene CYP27A1 and inhibiting mitochondrial oxidative phosphorylation in BRAF mutant melanomas.

Cellular senescence is a stress-induced state of stable growth arrest characterized by high expression of cell cycle inhibitors; a dramatic change in cell morphology, including an increase in lysosomal content; and secretion of large numbers of proteins involved in immune signaling and extracellular matrix remodeling. The physiological importance of cellular senescence has been attributed to prevention of carcinogenesis, aging, development, and tissue repair, and tumor cells can undergo senescence in response to therapeutic agents. In this work, we sought to validate the senescence-inducing activity of two known inducers (doxorubicin and TGFβ1) and a CDK4/6 inhibitor (ribociclib) and to identify proteins that can kill senescent tumor cells (senolytic targets) if knocked out and to identify downstream components of the tumor cell senescence pathway. Huh7 hepatocellular carcinoma cells and SK-MEL-28 melanoma cells were induced to senescence by treatment with three different agents: ribociclib, low doses of doxorubicin, or TGFβ1. The induction of senescence was confirmed by observing growth arrest, an increase in SA-β-gal staining, dramatic cell morphology changes, loss of c-Myc protein, increased expression of p15, and increased expression of senescence-associated secretory phenotype (SASP) proteins. Induction of SASP components was measured by RNAseq and SOMAscan. All three agents induced a senescent state, with blockage at different stages of the cell cycle observed. Induction of known immune factors, including IL-8 and IL-11, were identified in senescent cells (Huh7). A whole-genome CRISPR screen identified proteins required to enter senescence and those that were incompatible with the senescent state if knocked out. Expected hits were observed (eg, TGFBR1/TGFBR2 for TGFβ1, RB for ribociclib, and TOP2A for doxorubicin) for guide DNAs (gDNAs) that blocked entry into the senescent state. No gDNA candidates for common downstream senescence pathway components were observed, suggesting that these components are essential genes or that they do not exist. gDNA-induced knockouts that were incompatible with the senescent state dropped out of the screen and represent potential senolytic targets. The screen identified BCL2L1 as the only common senolytic hit across multiple senescence-inducing reagents, confirming published reports suggesting it is a senolytic target. These data show that ribociclib, doxorubicin, and TGFβ1 induced senescence in cancer cell lines. Whole-genome CRISPR screens identified senescence pathway components for each of these agents, as well as a common senolytic target. Citation Format: Pasupuleti Rao, Jennifer Tullai, Peter Aspesi, Felipa Mapa, Nadire Cochran, Frederic Sigoillot, Guglielmo Roma, Scott Gleim, Jaison Jacob, Jason Marchese, Jonathan Solomon. Characterization of cancer cell lines made senescent by exposure to ribociclib, doxorubicin, or TGFβ1, and identification of genes required for entry into senescence and senescent cell survival [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 2025.

We develop CellSIUS (Cell Subtype Identification from Upregulated gene Sets) to fill a methodology gap for rare cell population identification for scRNA-seq data. CellSIUS outperforms existing algorithms for specificity and selectivity for rare cell types and their transcriptomic signature identification in synthetic and complex biological data. Characterization of a human pluripotent cell differentiation protocol recapitulating deep-layer corticogenesis using CellSIUS reveals unrecognized complexity in human stem cell-derived cellular populations. CellSIUS enables identification of novel rare cell populations and their signature genes providing the means to study those populations in vitro in light of their role in health and disease., Genome Biology, 20 (1), ISSN:1474-760X

Comprehensive benchmarking of computational methods for single-cell RNA sequencing (scRNA-seq) analysis is scarce. Using a modular workflow and a large dataset with known cell composition, we benchmarked feature selection and clustering methodologies for scRNA-seq data. Results highlighted a methodology gap for rare cell population identification for which we developed CellSIUS (Cell Subtype Identification from Upregulated gene Sets). CellSIUS outperformed existing approaches, enabled the identification of rare cell populations and, in contrast to other methods, simultaneously revealed transcriptomic signatures indicative of the rare cells’ function. We exemplified the use of our workflow and CellSIUS for the characterization of a human pluripotent cell 3D spheroid differentiation protocol recapitulating deep-layer corticogenesis in vitro. Results revealed lineage bifurcation between Cajal-Retzius cells and layer V/VI neurons as well as rare cell populations that differ by migratory, metabolic, or cell cycle status, including a choroid plexus neuroepithelial subgroup, revealing previously unrecognized complexity in human stem cell-derived cellular populations.

The systematic translation of cancer genomic data into knowledge of tumour biology and therapeutic possibilities remains challenging. Such efforts should be greatly aided by robust preclinical model systems that reflect the genomic diversity of human cancers and for which detailed genetic and pharmacological annotation is available. Here we describe the Cancer Cell Line Encyclopedia (CCLE): a compilation of gene expression, chromosomal copy number and massively parallel sequencing data from 947 human cancer cell lines. When coupled with pharmacological profiles for 24 anticancer drugs across 479 of the cell lines, this collection allowed identification of genetic, lineage, and gene-expression-based predictors of drug sensitivity. In addition to known predictors, we found that plasma cell lineage correlated with sensitivity to IGF1 receptor inhibitors; AHR expression was associated with MEK inhibitor efficacy in NRAS-mutant lines; and SLFN11 expression predicted sensitivity to topoisomerase inhibitors. Together, our results indicate that large, annotated cell-line collections may help to enable preclinical stratification schemata for anticancer agents. The generation of genetic predictions of drug response in the preclinical setting and their incorporation into cancer clinical trial design could speed the emergence of 'personalized' therapeutic regimens.

Cell-based, phenotypic screening of small molecules often identifies compounds with provocative biological properties. However, determining the cellular target(s) and/or mechanism of action (MoA) of lead compounds remains an extremely challenging and time-consuming exercise. To provide insights into a compound's cellular action and greatly reduce the time required for MoA determination, we have developed a screening platform consisting of an extensive series of reporter gene assays (RGAs). A collection of > 11,000 compounds of known MoA (e.g., World Drug Index entries) were screened against the entire panel. The output provided evidence that an RGA signature could be ascribed to numerous, biologically diverse MoAs. The reference database generated suggested novel biological activity for particular compounds. For example, the profiling data led to the prediction that the cellular target of the natural product terprenin was dihydroorotate dehydrogenase (DHODH), which was confirmed experimentally. The screening methodology developed for this endeavor renders it amenable to the future examination of compounds with unknown MoA, in an automated, inexpensive, and time-efficient manner.

Cultivation of myxobacteria of the Nannocystis genus led to the isolation and structure elucidation of a class of novel cyclic lactone inhibitors of elongation factor 1. Whole genome sequence analysis and annotation enabled identification of the putative biosynthetic cluster and synthesis process. In biological assays the compounds displayed anti-fungal and cytotoxic activity. Combined genetic and proteomic approaches identified the eukaryotic translation elongation factor 1α (EF-1α) as the primary target for this compound class. Nannocystin A (1) displayed differential activity across various cancer cell lines and EEF1A1 expression levels appear to be the main differentiating factor. Biochemical and genetic evidence support an overlapping binding site of 1 with the anti-cancer compound didemnin B on EF-1α. This myxobacterial chemotype thus offers an interesting starting point for further investigations of the potential of therapeutics targeting elongation factor 1.

We recently reported the positional cloning of a candidate gene for hereditary hemochromatosis (HH), called HLA-H, which is a novel member of the major histocompatibility complex class I family. A mutation in this gene, cysteine 282 --tyrosine (C282Y), was found to be present in 83% of HH patient DNAs, while a second variant, histidine 63 --aspartate (H63D), was enriched in patients heterozygous for C282Y. The functional relevance of either mutation has not been described. Co-immunoprecipitation studies of cell lysates from human embryonic kidney cells transfected with wild-type or mutant HLA-H cDNA demonstrate that wild-type HLA-H binds beta2-microglobulin and that the C282Y mutation, but not the H63D mutation, completely abrogates this interaction. Immunofluorescence labeling and subcellular fractionations demonstrate that while the wild-type and H63D HLA-H proteins are expressed on the cell surface, the C282Y mutant protein is localized exclusively intracellularly. This report describes the first functional significance of the C282Y mutation by suggesting that an abnormality in protein trafficking and/or cell-surface expression of HLA-H leads to HH disease.

In the process of positionally cloning a candidate gene responsible for hereditary hemochromatosis (HH), we constructed a 1.1-Mb transcript map of the region of human chromosome 6p that lies 4.5 Mb telomeric to HLA-A. A combination of three gene-finding techniques, direct cDNA selection, exon trapping, and sample sequencing, were used initially for a saturation screening of the 1.1-Mb region for expressed sequence fragments. As genetic analysis further narrowed the HH candidate locus, we sequenced completely 0.25 Mb of genomic DNA as a final measure to identify all genes. Besides the novel MHC class 1-like HH candidate gene HLA-H, we identified a family of five butyrophilin-related sequences, two genes with structural similarity to a type 1 sodium phosphate transporter, 12 novel histone genes, and a gene we named RoRet based on its strong similarity to the 52-kD Ro/SSA lupus and Sjogren’s syndrome auto-antigen and the RET finger protein. Several members of the butyrophilin family and the RoRet gene share an exon of common evolutionary origin called B30-2. The B30-2 exon was originally isolated from the HLA class 1 region, yet has apparently “shuffled” into several genes along the chromosome telomeric to the MHC. The conservation of the B30-2 exon in several novel genes and the previously described amino acid homology of HLA-H to MHC class 1 molecules provide further support that this gene-rich region of 6p21.3 is related to the MHC. Finally, we performed an analysis of the four approaches for gene finding and conclude that direct selection provides the most effective probes for cDNA screening, and that as much as 30% of ESTs in this 1.1-Mb region may be derived from noncoding genomic DNA.[The sequence data described in this paper have been submitted to GenBank under accession nos. U90543–U90548,U90550–U90552, and U91328.]

Hereditary haemochromatosis (HH), which affects some 1 in 400 and has an estimated carrier frequency of 1 in 10 individuals of Northern European descent, results in multi-organ dysfunction caused by increased iron deposition, and is treatable if detected early. Using linkage-disequilibrium and full haplotype analysis, we have identified a 250-kilobase region more than 3 megabases telomeric of the major histocompatibility complex (MHC) that is identical-by-descent in 85% of patient chromosomes. Within this region, we have identified a gene related to the MHC class I family, termed HLA-H, containing two missense alterations. One of these is predicted to inactivate this class of proteins and was found homozygous in 83% of 178 patients. A role of this gene in haemochromatosis is supported by the frequency and nature of the major mutation and prior studies implicating MHC class I-like proteins in iron metabolism.

Death Receptor 5 (DR5) agonists demonstrate anti-tumor activity in preclinical models but have yet to demonstrate robust clinical responses. A key limitation may be the lack of patient selection strategies to identify those most likely to respond to treatment. To overcome this limitation, we screened a DR5 agonist Nanobody across >600 cell lines representing 21 tumor lineages and assessed molecular features associated with response. High expression of DR5 and Casp8 were significantly associated with sensitivity, but their expression thresholds were difficult to translate due to low dynamic ranges. To address the translational challenge of establishing thresholds of gene expression, we developed a classifier based on ratios of genes that predicted response across lineages. The ratio classifier outperformed the DR5+Casp8 classifier, as well as standard approaches for feature selection and classification using genes, instead of ratios. This classifier was independently validated using 11 primary patient-derived pancreatic xenograft models showing perfect predictions as well as a striking linearity between prediction probability and anti-tumor response. A network analysis of the genes in the ratio classifier captured important biological relationships mediating drug response, specifically identifying key positive and negative regulators of DR5 mediated apoptosis, including DR5, CASP8, BID, cFLIP, XIAP and PEA15. Importantly, the ratio classifier shows translatability across gene expression platforms (from Affymetrix microarrays to RNA-seq) and across model systems (in vitro to in vivo). Our approach of using gene expression ratios presents a robust and novel method for constructing translatable biomarkers of compound response, which can also probe the underlying biology of treatment response.

Presenilin is the enzymatic component of γ-secretase, a multisubunit intramembrane protease that processes several transmembrane receptors, such as the amyloid precursor protein (APP). Mutations in human Presenilins lead to altered APP cleavage and early-onset Alzheimer's disease. Presenilins also play an essential role in Notch receptor cleavage and signaling. The Notch pathway is a highly conserved signaling pathway that functions during the development of multicellular organisms, including vertebrates, Drosophila, and C. elegans. Recent studies have shown that Notch signaling is sensitive to perturbations in subcellular trafficking, although the specific mechanisms are largely unknown. To identify genes that regulate Notch pathway function, we have performed two genetic screens in Drosophila for modifiers of Presenilin-dependent Notch phenotypes. We describe here the cloning and identification of 19 modifiers, including nicastrin and several genes with previously undescribed involvement in Notch biology. The predicted functions of these newly identified genes are consistent with extracellular matrix and vesicular trafficking mechanisms in Presenilin and Notch pathway regulation and suggest a novel role for γ-tubulin in the pathway.

For nearly a century, genetic analysis in Drosophila melanogaster has been a powerful tool for analyzing gene function, yet Drosophila lacks the molecular genetic mapping tools that recently have revolutionized human, mouse, and plant genetics. Here, we describe the systematic characterization of a dense set of molecular markers in Drosophila by using a sequence tagged site-based physical map of the genome. We identify 474 biallelic markers in standard laboratory strains of Drosophila that span the genome. Most of these markers are single nucleotide polymorphisms and sequences for these variants are provided in an accessible format. The average density of the new markers is one per 225 kb on the autosomes and one per megabase on the X chromosome. We include in this survey a set of P-element strains that provide additional use for high-resolution mapping. We show one application of the new markers in a simple set of crosses to map a mutation in the hedgehog gene to an interval of Drosophila research community.