Abstract We have conducted systematic small interfering RNA (siRNA) lethality screening of the kinome and “druggable genome” in human myeloma cell lines (HMCL) to functionally generate a comprehensive map of critical genes and vulnerabilities in human myeloma tumors. KMS11 human myeloma cells were screened with both an 1800 oligo (639 gene) siRNA library targeting the kinome and with a 13,984 oligo (6,791 gene) library targeting the druggable genome in duplicate using optimized conditions that resulted in >95% transfection efficiency and <5% background cytotoxicity. siRNA were used at low concentration (13nM) to minimise off-target cellular effects. Each gene was initially screened with at least 2 distinct siRNA. Universally lethal and non-silencing siRNA were tested in parallel as controls. Viability was measured at 96 hours by ATP-dependent luminescence and was normalized to the mean and distribution of controls by the B-score method. Approximately 14% of kinome siRNA and 5.8% of druggable genome siRNA caused important reductions in HMCL viability (greater than three standard deviations from control samples treated with non-silencing siRNA), compared with an anticipated rate due to chance of only 0.135%. Short-listed kinome vulnerabilities were validated by repeat targeting in a second lethality screening experiment in KMS11 using a total of four independent siRNA per gene. Forty-two kinases were confirmed as critical for myeloma cell survival. Unsupervised pathway analysis indicates that approximately half are concentrated within established pathways involved in cytokine (VEGFR, FGFR3, IGFR/IL6R) signaling and regulation of cell cycle, apoptosis or metabolism. The role of other survival kinases in myeloma are novel. The selectivity of kinome vulnerabilities identified in KMS11 for myeloma cells versus other tissues was examined by parallel silencing of these targets in KMS11 and JJN3 myeloma cells and in non-myeloma 293 and A549 epithelial cells. The most lethal kinase targets in myeloma that simultaneously show experimental evidence of tumor selectivity include aurora kinase A (AURKA), G-protein receptor coupled kinase 6 (GRK6) and protein kinase N1 (PKN1). By gene expression analyses, PKN1 and GRK6 are both preferentially expressed in plasma cells, compared with an analysis of >50 somatic tissues in which they are absent or only weakly expressed. The extent to which the kinome vulnerabilities identified in KMS11 occur in other HMCL was assessed by tertiary siRNA screening and by lentiviral shRNA techniques. More than 75% of kinases vulnerable in KMS11 cells are recurrently vulnerable in other myeloma tumor lines. Secondary validation studies of all genomic vulnerabilities identified from the much larger siRNA lethality screen of the druggable genome in myeloma cells is ongoing, however, to date 134 genes have been validated as essential myeloma survival factors during primary screening by concordant lethality of 2/2 independent siRNA (compared with an expected rate due to chance of n<1 genes). Significantly, proteosome subunits and specific ubiquitin ligases number amongst the most potently lethal genomic vulnerabilities in KMS11. Overall, this attempt to systematically and functionally document all genetic vulnerabilities in myeloma provides a strong rational impetus for the development of novel therapies that optimally target susceptibilities inherent within myeloma tumors.
The Multiple Myeloma Research Consortium Genomics Initiative.
Abstract The Multiple Myeloma Research Consortium (MMRC) Genomics Initiative is a three-year program to analyze tumor tissue from hundreds of multiple myeloma (MM) patients via gene expression profiling (GEP), comparative genomic hybridization (aCGH), and exon re-sequencing. In addition, RNAi knockdown of selected genes in MM tumor cell lines is being evaluated to identify potential new targets. All genomic data generated is scheduled for placement in an open-access Multiple Myeloma Genomics Portal pre-publication and in near real-time (www.broad.mit.edu/mmgp). Additionally, samples are also destined for drug validation and correlative science on clinical protocols as this study moves forward. This comprehensive project is spearheaded by the MMRC and conducted via collaboration with the Eli and Edythe L. Broad Institute of MIT and Harvard, the Translational Genomics Research Institute (TGen), Mayo Clinic Arizona, and The Dana-Farber Cancer Center. The study is supported by the collection from member institutions of the MMRC of bone marrow aspirates and matched peripheral blood samples from over 1000 patients. Specific genomic technologies that are currently being employed across this sample set include GEP using Affymetrix Human Genome U133A 2.0 Plus Arrays, and, in parallel, efforts to identify regions of genomic gain and loss are using Agilent Human Genome CGH arrays. In contrast to other large-scale genomic projects based on exon-sequencing of targeted gene sets, this project will be the first to perform genome-scale single molecule sequencing (SMS) of DNA from patient specimens. Results will be targeted against candidate classes of genes (e.g. kinases, phosphatases, known oncogenes and tumor suppressors), and genes from GEP or within candidate regions of copy gain or loss identified by the aCGH experiments. Mutations will be further validated in an independent set of patient specimens. Finally we will attempt to identify points of vulnerability of MM through systematic loss-of-function screens in myeloma cell lines using high-throughput RNA interference (using both shRNA and siRNA platforms). Importantly, data generated from this genomics initiative will ultimately be made public pre-publication through the established MMRC Multiple Myeloma Genomics Portal. Data from all aspects of this project (sample collection and analyte isolation, GEP, aCGH, SMS, RNAi and bioinformatics) will be described in this presentation. The power of this study is the comprehensive collection of gene expression, CGH, and genome sequencing on a single reference set of clinically annotated samples. The addition of RNAi screens makes this a very important and unique data resource, which we hope will help expedite the discovery of novel targeted agents for MM scientific community.
RNAi Screen of the Druggable Genome To Identify Modulators of Bortezomib in Myeloma.
Abstract Bortezomib and related proteasome inhibitors have profound pre-clinical and clinical activity in Multiple Myeloma (MM). However, the critical molecular target(s) modulated by proteosome inhibition remains unknown. The overall response rate to bortezomib in relapsed and refractory MM patients is consistently 35%–50%, suggesting specific molecular targets or resistance mechanisms are inherent in patients. None of the known effects of bortezomib has been shown to correlate with clinical responses in patients or to explain the unique sensitivity of MM (as opposed to other malignancies) to these reagents. Therefore, in this study, we have conducted a small interfering RNA (siRNA) screen of the druggable genome in human myeloma cells to identify genes that modulate sensitivity to bortezomib. KMS11 human myeloma cells were reverse transfected with 13,984 synthetic short interfering RNAs (siRNA) against each of 6992 druggable targets (2 siRNA/gene with a single-siRNA per well/384 well format). At 24 hours post transfection, bortezomib was added at IC10, 30, 60 and 90. Cell viability was measured by ATP-dependent luminescence assay 72 hours post drug treatment. The resulting data from primary screening was evaluated for multiple quality control metrics and found to exceed all expected performance parameters with >99% global transfection efficiency, <0.25 CV values across the dataset with minimal plate-to-plate and set-to-set variation observed. From analysis of the primary screening data, a core panel of 320 (4.6% of druggable targets in library) bortezomib sensitizers were identified. 187 (58%) of these were further confirmed in a secondary screening experiment by multiple (2 to 4) independent siRNA. Among these 187 genes (2.6% of the genes investigated), some known targets, such as proteasome subunits are evident. Other genes belonging to the same functional gene families (e.g. TNF receptor family, Ras family, kinases, oncogenes and transporter proteins) or functioning within the same network of signal transduction (e.g. NF-kB, Wnt and estrogen receptor pathways) are evident. Our primary screen also identified genes that attenuate bortezomib-induced cytotoxicity. Members in this group include genes implicated in the negative regulation of cell proliferation and survival, such as NF-kB inhibitors, JNK pathway activators, tumor suppressors, apoptosis inducers and phosphatases. These studies enable dissection of the molecular mechanism of bortezomib in myeloma cells and may serve to point to biomarkers for response prediction and novel chemotherapeutic combination therapies.
A Comprehensive Study of Chromosome 13 Abnormalities (Δ 13) in Multiple Myeloma Using High-Throughput, Array-Based Methods.
Abstract Chromosome 13 abnormalities (Δ13), monosomy or deletion, are seen in greater than 50% in patients with Multiple Myeloma (MM) and are associated with poor survival and reduced response to therapy. Δ13 is suggestive of a tumor suppressor model. In this study, we use a comprehensive, high-throughput approach involving array-based technologies to identify genes on chromosome 13 that are important in the pathogenesis and/or progression of MM. Array-based Comparative Genomic Hybridization (aCGH) is being conducted to identify a minimum region of 13q loss in a series of MM cell lines and patient samples. DNA is isolated, labeled and hybridized with a differentially labeled normal DNA reference to determine gene/genomic copy number changes. Arrays are analyzed to search for the minimum region of loss based upon single copy loss for a series of nearby mapping transcripts in each cell line. Data will be cross-referenced to determine the true minimum region of 13q loss in MM. Also, to discover genes harboring potential inactivating mutations, we are using inhibition of nonsense-mediated RNA decay (NMD) coupled with cDNA microarray. This modification of the Gene Identification by NMD Inhibition (GINI) method, first introduced by Noensie and Dietz (2001), is used to identify genes on chromosome 13 that may be inactivated due to premature truncating codons (PTCs). This method is based upon a cell’s natural ability to survey for and degrade RNA species that contain PTCs by activation of the nonsense mediated RNA decay (NMD) pathway. Treating cells with the drug Emetine can inhibit NMD and increase the expression of genes likely harboring truncating mutations. Preliminary data have revealed several candidate genes, with increased ratios after Emetine treatment, which map to chromosome 13. Upon defining a minimum region of loss based upon aCGH data, we will prioritize genes from our NMD study for mutational analysis. It our hope to also perform an array-based CpG island screen to simultaneously analyze all CpG islands on chromosome 13 to detect genes possibly inactivated by hypermethylation. This comprehensive, organized approach incorporating high-throughput techniques for the detection of genomic loss, mutation, and methylation is the first of its kind to our knowledge and may maximize the potential for the identification of a 13q tumor suppressor in MM.
RNAi Screen of the Druggable Genome Identifies Modulators of Proteasome Inhibitor Sensitivity in Myeloma Including CDK5.
Abstract Abstract 602 The molecular target(s) which co-operate with proteasome inhibition in inducing drug sensitivity or resistance in Multiple Myeloma (MM) remain unknown. We therefore conducted a genome scale small interfering RNA (siRNA) lethality study in KMS11 MM cells in the presence or absence of bortezomib without regard to pre-conceived mechanistic notions. Primary screening was performed in a single-siRNA-per-well format with the human druggable genome siRNA set V4 comprising 13,984 siRNA targeting 6,992 genes and comprising two RNAi per gene. siRNA were transfected at low concentration (13nM) to minimize off-target effects using conditions that resulted in transfection of >95% cells and <5% background cytotoxicity. Bortezomib was added (at the IC10, 25, 70 and 90) 24 hours post transfection. After 96 hours (72 hours after bortezomib), viability was measured by ATP-dependent luminescence. In primary screening 320 candidate bortezomib sensitizing genes were identified and rescreened with four siRNA oligos each gene. 57 of the top sensitizer hits were selected for which at least two distinct siRNA decreased the EC50 by 2 standard deviations from cells treated with control siRNAs. By cross referencing gene expression profile (GEP) data of KMS11, we furthered pared the list to 37 plausibly expressed targets (0.5% of those genes originally screened) as bortezomib sensitizers. After silencing, 50% of these 37 genes also sensitized the lung cancer cell line A549 to bortezomib. After further parsing of genes which also modulated the sensitivity of MM to Melphalan (non specific chemosensitizers), 34 genes remained: The strongest sensitizers to bortezomib were the proteasome subunits PSMA5, PSMB2, PSMB3, PSMB7 but included less obvious targets such as BAZ1B, CDK5, CDC42SE2, MDM4, NME7, TFE3, TNFAIP3, TNK1, TOP1, VAMP2 and YY1 were also identified. Of these, the most potent synergetic effects were observed with siRNAs against the proteasome subunits and against cyclin dependent kinase 5 (CDK5), which caused the greatest shift in EC50. CDK5 is of particular interest as a therapeutic target as it is expressed at high levels in MM and neural tissues but has low expression in other organs. Using viral shRNA expression, silencing CDK5 consistently increased the sensitivity of genetically variable MM cell lines (n=5) to all of the proteasome inhibitors tested: bortezomib, carfilzomib and PR047, and the effect could be at least partially rescued by overexpression of an RNAi resistant CDK5. To explore therapeutic relevance the small molecule CDK5 inhibitor, Roscovitin, was shown to be synergistic or additive with bortezomib in both MM cell lines and primary patient samples. Gene expression profiling was then performed to seek an explanation for the CDK5 sensitization effect and regulation of a proteasome subunit PSMB5 by CDK5 was identified as a probable route to sensitization. In summary inhibition of the existing proteasome either directly by suppression of proteasome subunits, or indirectly by suppression of modulators such as CDK5 appears to confer the greatest sensitization effect suggesting that combinations of bortezomib with other unique proteasome inhibitor drugs or combinations with inhibitors of CDK5 is a logical avenue for clinical exploration. Disclosures: No relevant conflicts of interest to declare.
Kinome-Wide RNAi Studies in Human Multiple Myeloma Identify a Lymphoid Restricted Kinase GRK6 as a Selectively Vulnerable Target That Regulates STAT3/MCL1.
Abstract Abstract 601 A paucity of validated kinase targets in human multiple myeloma (MM) has delayed clinical deployment of kinase inhibitors in treatment strategies. We therefore conducted a kinome-wide small interfering RNA (siRNA) lethality study in MM tumor lines bearing common t(4;14), t(14;16) and t(11;14) translocations to identify critically vulnerable kinases in MM tumor cells without regard to preconceived mechanistic notions. Primary screening was performed in duplicate using an 1800-oligo siRNA library in a single-siRNA-per-well format. siRNA were transfected at low concentration (13nM) to minimize off-target effects using conditions that resulted in transfection of >95% cells and <5% background cytotoxicity. After 96 hours, viability was measured by ATP-dependent luminescence. Fifteen kinases were consistently vulnerable in MM cells, including AKT1, AK3L1, AURKA, AURKB, CDC2L1, CDK5R2, FES, FLT4, GAK, GRK6, HK1, PKN1, PLK1, SMG1, and TNK2. While several kinases (PLK1, HK1) were equally vulnerable in epithelial cells, others and particularly the G-protein coupled receptor kinase, GRK6, appeared selectively vulnerable in MM. GRK6 inhibition is selectively lethal to KMS11, OPM1, H929, KMS18 and OCI-MY5 myeloma cells and has minimal effect on 293, MCF7, SF767, A375 or A549 epithelial cells. Persistent expression of FLAG-GRK6 via cDNA rescued KMS11 cells from the lethal effect of a 3'UTR-targeted GRK6 siRNA, but not from control siRNA, validating identification of GRK6 as an essential myeloma survival kinase. Furthermore, concordant results were obtained using four different exon-based GRK6 siRNA, all of which induced GRK6 silencing and similar inhibition of KMS11 proliferation and viability. Significantly, GRK6 is ubiquitously expressed in lymphoid tissues and myeloma, but by comparison appears absent or only weakly expressed in most primary human somatic tissues. From co-immunoprecipitation experiments we demonstrate that GRK6 is highly expressed in myeloma cells via direct association with the HSP90 chaperone. Inhibition of HSP90 with geldanamycin blocks GRK6 protein expression. Importantly, direct GRK6 silencing causes rapid and selective suppression of STAT3 phosphorylation that is associated with sustained reductions in total MCL1 protein levels and MCL1 phosphorylation (within 24 hours), providing a potent mechanism for the cytotoxicity of GRK6 inhibition in MM tumor cells. GF109203X is an inhibitor of both protein kinase C and of GRK6 that causes near total inhibition of these kinases in vitro at distinct concentrations of 0.1μM and 1-10μM respectively. Notably, GF109203X was substantially cytotoxic to 10/14 myeloma tumor lines at concentrations most consistent with GRK6 inhibition (5-20μM), and was selectively more cytotoxic to myeloma tumor cells than to non-myeloma cell lines (P=0.01), highlighting the potential of GRK6 as a pharmaceutical target for selective therapeutic intervention in myeloma. As mice that lack GRK6 are healthy, inhibition of GRK6 represents a uniquely targeted novel therapeutic strategy in human multiple myeloma. Disclosures: Perkins: MMRC: Employment. Reeder:Celgene: Research Funding; Millennium: Research Funding. Fonseca:Otsuka: Consultancy; BMS: Consultancy; Amgen: Consultancy; Medtronic: Consultancy; Genzyme: Consultancy.
Achilles Heel Molecular Vulnerabilities in Multiple Myeloma Identified From Genome-Scale RNA Interference Screening.
Abstract Abstract 2801 Poster Board II-777 To provide a rationale basis for targeted drug development for multiple myeloma patients, we have conducted high-throughput genome-scale RNA interference (RNAi) synthetic lethality studies in human myeloma cell lines (HMCL) to generate a comprehensive map of critical genes and molecular vulnerabilities in human myeloma tumor cells. KMS11 human myeloma cells were screened with a 13,982-oligo library targeting the ‘druggable' genome (6,991 genes) using optimized conditions that resulted in >95% transfection efficiency. Each gene was screened with 2 or more distinct oligos, in duplicate, using a single-siRNA-per-well format, testing >34,000 wells. Replicate high throughput experiments yielded highly reproducible results (R2=0.82). Viability was measured at 96h by ATP-dependent luminescence. Universally lethal siRNA and non-targeted siRNA were employed as controls. The specificity (FDR) of lethal RNAi was evaluated by custom-developed statistical methods based on RNAi result concordancy. From screening, 5.8% of druggable genome siRNA caused statistically relevant reductions in HMCL viability (greater than three standard deviations from control samples treated with non-silencing siRNA, compared with an anticipated rate due to chance of only 0.135%). Of these, two hundred and nineteen genes, targeted by the most lethal siRNA, were forwarded to validation studies. Validation high-throughput RNAi studies, using 4 oligos per gene, were conducted in order to verify target gene vulnerability. Ultimately, seventy two genes were validated as highly critical for myeloma cell survival (each with multiple concordant siRNA hits plus high reproducibility of R2 =0.94). Among top-ranked lethal molecular vulnerabilities in KMS11 myeloma cells we recurrently identified the proteasome (8 PSM subunits were independently identified as highly vulnerable); and BCL2 family member MCL1; both known therapeutic targets, validating a functional genomics approach to target discovery. In addition, a number of novel, equally lethal, molecular vulnerabilities were identified in KMS11 including WEE1, kinetochore complex component KNTC2, the aurora kinases, polo-like kinase 1, a previously uncharacterized DNA methyl transferase, ribonucleotide reductase, ribosomal protein L38, leptin receptor overlapping transcript LEPROT and SLC25A23 (a carrier responsible for mitochondrial adenine nucleotide flux). Various ubiquitous cellular proteins involved in RNA or protein processing (e.g. SF3a, SNW1, SNRPA1, EIF3s8, amongst others) were also identified as non redundant and critical for cellular viability. Thirty nine top-ranked molecular vulnerabilities in KMS11 were evaluated for comparative vulnerability in a second myeloma cell line, JJN3, and in A549 and 293 epithelial cells. A majority were vulnerable in JJN3. While 25 genes proved equally lethal on silencing in A549 and 293 epithelial cells, several, notably including MCL1 and three proteasome subunits, appear differentially susceptible in myeloma cells versus epithelial cells. From this genome scale study, MCL1 ranks third as a target capable of inducing cell death in myeloma cells, highlighting the potent vulnerability of myeloma cell viability to alterations in this anti-apoptosis regulator. Moreover, from active comparative studies of molecular vulnerability in myeloma versus epithelial cell lines, MCL1 as an RNAi target demonstrates perhaps the greatest cytotoxic selectivity between myeloma and epithelial cell lines of all targets examined. Gene expression studies indicate that MCL1 is highly expressed in myeloma tumor cells but is absent or only weakly expressed in human primary somatic tissues. Overall, therefore, from an unsupervised genome-scale screening approach, our data support MCL1 as an optimal Achilles heel molecular target in multiple myeloma, potentially offering greater specificity and deeper cytotoxic effect than proteasome inhibition. Disclosures: No relevant conflicts of interest to declare.
Characterization of Farnesyl Diphosphate Farnesyl Transferase 1 (
FDFT1
) Expression In Cancer
Abstract The molecular target(s) cooperating with proteasome inhibition in multiple myeloma (MM) remain unknown. We therefore measured proliferation in MM cells transfected with 13 984 small interfering RNAs in the absence or presence of increasing concentrations of bortezomib. We identified 37 genes, which when silenced, are not directly cytotoxic but do synergistically potentiate the growth inhibitory effects of bortezomib. To focus on bortezomib sensitizers, genes that also sensitized MM to melphalan were excluded. When suppressed, the strongest bortezomib sensitizers were the proteasome subunits PSMA5, PSMB2, PSMB3, and PSMB7 providing internal validation, but others included BAZ1B, CDK5, CDC42SE2, MDM4, NME7, RAB8B, TFE3, TNFAIP3, TNK1, TOP1, VAMP2, and YY1. The strongest hit CDK5 also featured prominently in pathway analysis of primary screen data. Cyclin-dependent kinase 5 (CDK5) is expressed at high levels in MM and neural tissues with relatively low expression in other organs. Viral shRNA knockdown of CDK5 consistently sensitized 5 genetically variable MM cell lines to proteasome inhibitors (bortezomib and carfilzomib). Small-molecule CDK5 inhibitors were demonstrated to synergize with bortezomib to induce cytotoxicity of primary myeloma cells and myeloma cell lines. CDK5 regulation of proteasome subunit PSMB5 was identified as a probable route to sensitization.