RNA interference screening identifies lenalidomide sensitizers in multiple myeloma, including RSK2.

To identify molecular targets that modify sensitivity to lenalidomide, we measured proliferation in multiple myeloma (MM) cells transfected with 27 968 small interfering RNAs in the presence of increasing concentrations of drug and identified 63 genes that enhance activity of lenalidomide upon silencing. Ribosomal protein S6 kinase (RPS6KA3 or RSK2) was the most potent sensitizer. Other notable gene targets included 5 RAB family members, 3 potassium channel proteins, and 2 peroxisome family members. Single genes of interest included I-κ-B kinase-α (CHUK), and a phosphorylation dependent transcription factor (CREB1), which associate with RSK2 to regulate several signaling pathways. RSK2 knockdown induced cytotoxicity across a panel of MM cell lines and consistently increased sensitivity to lenalidomide. Accordingly, 3 small molecular inhibitors of RSK2 demonstrated synergy with lenalidomide cytotoxicity in MM cells even in the presence of stromal contact. Both RSK2 knockdown and small molecule inhibition downregulate interferon regulatory factor 4 and MYC, and provides an explanation for the synergy between lenalidomide and RSK2 inhibition. Interestingly, RSK2 inhibition also sensitized MM cells to bortezomib, melphalan, and dexamethasone, but did not downregulate Ikaros or influence lenalidomide-mediated downregulation of tumor necrosis factor-α or increase lenalidomide-induced IL-2 upregulation. In summary, inhibition of RSK2 may prove a broadly useful adjunct to MM therapy.

Use of mechanistic models to integrate and analyze multiple proteomic datasets.

Proteins in cell signaling networks tend to interact promiscuously through low-affinity interactions. Consequently, evaluating the physiological importance of mapped interactions can be difficult. Attempts to do so have tended to focus on single, measurable physicochemical factors, such as affinity or abundance. For example, interaction importance has been assessed on the basis of the relative affinities of binding partners for a protein of interest, such as a receptor. However, multiple factors can be expected to simultaneously influence the recruitment of proteins to a receptor (and the potential of these proteins to contribute to receptor signaling), including affinity, abundance, and competition, which is a network property. Here, we demonstrate that measurements of protein copy numbers and binding affinities can be integrated within the framework of a mechanistic, computational model that accounts for mass action and competition. We use cell line-specific models to rank the relative importance of protein-protein interactions in the epidermal growth factor receptor (EGFR) signaling network for 11 different cell lines. Each model accounts for experimentally characterized interactions of six autophosphorylation sites in EGFR with proteins containing a Src homology 2 and/or phosphotyrosine-binding domain. We measure importance as the predicted maximal extent of recruitment of a protein to EGFR following ligand-stimulated activation of EGFR signaling. We find that interactions ranked highly by this metric include experimentally detected interactions. Proteins with high importance rank in multiple cell lines include proteins with recognized, well-characterized roles in EGFR signaling, such as GRB2 and SHC1, as well as a protein with a less well-defined role, YES1. Our results reveal potential cell line-specific differences in recruitment.

In Vitro Transcription/Translation System: A Versatile Tool in the Search for Missing Proteins.

Approximately 18% of all human genes purported to encode proteins have not been directly evidenced at the protein level, according to the validation criteria established by neXtProt, and are considered to be "missing" proteins. One of the goals of the Chromosome-Centric Human Proteome Project (C-HPP) is to identify as many of these missing proteins as possible in human samples using mass spectrometry-based methods. To further this goal, a consortium of C-HPP teams (chromosomes 5, 10, 16, and 19) has joined forces to devise new strategies to identify missing proteins by use of a cell-free in vitro transcription/translation system (IVTT). The proposed strategy employs LC-MS/MS data-dependent acquisition (DDA) and targeted selective reaction monitoring (SRM) methods to scrutinize low-complexity samples derived from IVTT. The optimized assays are then applied to identify missing proteins in human cells and tissues. We describe the approach and show proof-of-concept results for development of LC-SRM assays for identification of 18 missing proteins. We believe that the IVTT system, when coupled with downstream mass spectrometric identification, can be applied to identify proteins that have eluded more traditional methods of detection.

Whole-genome sequencing of an aggressive BRAF wild-type papillary thyroid cancer identified EML4-ALK translocation as a therapeutic target.

BACKGROUND: Recent advances in the treatment of cancer have focused on targeting genomic aberrations with selective therapeutic agents. In radioiodine resistant aggressive papillary thyroid cancers, there remain few effective therapeutic options. A 62-year-old man who underwent multiple operations for papillary thyroid cancer and whose metastases progressed despite standard treatments provided tumor tissue. METHODS: We analyzed tumor and whole blood DNA by whole genome sequencing, achieving 80× or greater coverage over 94 % of the exome and 90 % of the genome. We determined somatic mutations and structural alterations. RESULTS: We found a total of 57 somatic mutations in 55 genes of the cancer genome. There was notably a lack of mutations in NRAS and BRAF, and no RET/PTC rearrangement. There was a mutation in the TRAPP oncogene and a loss of heterozygosity of the p16, p18, and RB1 tumor suppressor genes. The oncogenic driver for this tumor is a translocation involving the genes for anaplastic lymphoma receptor tyrosine kinase (ALK) and echinoderm microtubule associated protein like 4 (EML4). The EML4-ALK translocation has been reported in approximately 5 % of lung cancers, as well as in pediatric neuroblastoma, and is a therapeutic target for crizotinib. CONCLUSIONS: This is the first report of the whole genomic sequencing of a papillary thyroid cancer in which we identified an EML4-ALK translocation of a TRAPP oncogene mutation. These findings suggest that this tumor has a more distinct oncogenesis than BRAF mutant papillary thyroid cancer. Whole genome sequencing can elucidate an oncogenic context and expose potential therapeutic vulnerabilities in rare cancers.

NONO and RALY proteins are required for YB-1 oxaliplatin induced resistance in colon adenocarcinoma cell lines.

BACKGROUND: YB-1 is a multifunctional protein that affects transcription, splicing, and translation. Overexpression of YB-1 in breast cancers causes cisplatin resistance. Recent data have shown that YB-1 is also overexpress in colorectal cancer. In this study, we tested the hypothesis that YB-1 also confers oxaliplatin resistance in colorectal adenocarcinomas. RESULTS: We show for the first time that transfection of YB-1 cDNA confers oxaliplatin resistance in two colorectal cancer cell lines (SW480 and HT29 cell lines). Furthermore, we identified by mass spectrometry analyses important YB-1 interactors required for such oxaliplatin resistance in these colorectal cancer cell lines. A tagged YB-1 construct was used to identify proteins interacting directly to YB-1 in such cells. We then focused on proteins that are potentially involved in colorectal cancer progression based on the Oncomine microarray database. Genes encoding for these YB-1 interactors were also examined in the public NCBI comparative genomic hybridization database to determine whether these genes are localized to regions of chromosomes rearranged in colorectal cancer tissues. From these analyses, we obtained a list of proteins interacting with YB-1 and potentially involved in oxaliplatin resistance. Oxaliplatin dose response curves of SW480 and HT29 colorectal cancer cell lines transfected with several siRNAs corresponding to each of these YB-1 interactors were obtained to identify proteins significantly affecting oxaliplatin sensitivity upon gene silencing. Only the depletion of either NONO or RALY sensitized both colorectal cancer cell lines to oxaliplatin. Furthermore, depletion of NONO or RALY sensitized otherwise oxaliplatin resistant overexpressing YB-1 SW480 or HT29 cells. CONCLUSION: These results suggest knocking down NONO or RALY significant counteracts oxaliplatin resistance in colorectal cancers overexpressing the YB-1 protein.

Cancer gene mutation frequencies for the U.S. population.

Mutations play a fundamental role in the development of cancer, and many create targetable vulnerabilities. There are both public health and basic science benefits from the determination of the proportion of all cancer cases within a population that include a mutant form of a gene. Here, we provide the first such estimates by combining genomic and epidemiological data. We estimate KRAS is mutated in only 11% of all cancers, which is less than PIK3CA (13%) and marginally higher than BRAF (8%). TP53 is the most commonly mutated gene (35%), and KMT2C, KMT2D, and ARID1A are among the ten most commonly mutated driver genes, highlighting the role of epigenetic dysregulation in cancer. Analysis of major cancer subclassifications highlighted varying dependencies upon individual cancer drivers. Overall, we find that cancer genetics is less dominated by high-frequency, high-profile cancer driver genes than studies limited to a subset of cancer types have suggested.

A phase 2 trial of personalized cytotoxic therapy based on tumor immunohistochemistry in previously treated metastatic pancreatic cancer patients.

BACKGROUND: The choice of a regimen in metastatic pancreatic cancer patients following progression on 1st line therapy is empiric and outcomes are unsatisfactory. This phase II study was performed to evaluate the efficacy of therapy selected by immunohistochemistry (IHC) in these patients following progression after one or more therapies. METHODS: Eligible patients underwent a percutaneous biopsy of a metastatic lesion and treatment selection was determined by IHC. The study required 35 evaluable patients (power of 86%) for detecting a true 1-year survival rate of >20%. RESULTS: A tumor biopsy was performed in 48 of 49 accrued patients. Study therapy was not given (n=13) either due to insufficient tumor on biopsy (n=8) or due to worsening cancer related symptoms after biopsy (n=5). The demographics of evaluable patients (n=35) are male/female (59%/41%), with age range 34-78 years (median 63 years). Patients had 1-6 prior regimens (median of 2). The most common IHC targets were topoisomerase 1 or 2, thymidylate synthase, excision repair cross-complementation group 1 protein (ERCC1), and osteonectin secreted protein acidic and rich in cysteine (SPARC). Commercially available treatment regimens prescribed included FOLFIRI, FOLFOX, irinotecan, and doxorubicin. The response (RECIST) was 9%, the median survival was 5.6 months (94% CI, 3.8-8.2), and the 1-year survival was 20% (95% CI, 7-33%). CONCLUSIONS: In all patients, IHC assays resulted in identification of at least two targets for therapy and a non-cross resistant regimen could be prescribed for therapy with evidence of some benefit. An IHC based treatment strategy is feasible and needs validation in larger studies.

Long-Term ERK Inhibition in KRAS-Mutant Pancreatic Cancer Is Associated with MYC Degradation and Senescence-like Growth Suppression.

Induction of compensatory mechanisms and ERK reactivation has limited the effectiveness of Raf and MEK inhibitors in RAS-mutant cancers. We determined that direct pharmacologic inhibition of ERK suppressed the growth of a subset of KRAS-mutant pancreatic cancer cell lines and that concurrent phosphatidylinositol 3-kinase (PI3K) inhibition caused synergistic cell death. Additional combinations that enhanced ERK inhibitor action were also identified. Unexpectedly, long-term treatment of sensitive cell lines caused senescence, mediated in part by MYC degradation and p16 reactivation. Enhanced basal PI3K-AKT-mTOR signaling was associated with de novo resistance to ERK inhibitor, as were other protein kinases identified by kinome-wide siRNA screening and a genetic gain-of-function screen. Our findings reveal distinct consequences of inhibiting this kinase cascade at the level of ERK.

Integrated genomic approaches identify upregulation of SCRN1 as a novel mechanism associated with acquired resistance to erlotinib in PC9 cells harboring oncogenic EGFR mutation.

Therapies targeting the tyrosine kinase activity of Epidermal Growth Factor Receptor (EGFR) have been proven to be effective in treating a subset of non-small cell lung cancer (NSCLC) patients harboring activating EGFR mutations. Inevitably these patients develop resistance to the EGFR-targeted tyrosine kinase inhibitors (TKIs). Here, we performed integrated genomic analyses using an in vitro system to uncover alternative genomic mechanisms responsible for acquired resistance to EGFR-TKIs. Specifically, we identified 80 genes whose expression is significantly increased in the erlotinib-resistant clones. RNAi-based systematic synthetic lethal screening of these candidate genes revealed that suppression of one upregulated transcript, SCRN1, a secernin family member, restores sensitivity to erlotinib by enhancing inhibition of PI3K/AKT signaling pathway. Furthermore, immunohistochemical analysis revealed increased levels of SCRN1 in 5 of 11 lung tumor specimens from EGFR-TKIs resistant patients. Taken together, we propose that upregulation of SCRN1 is an additional mechanism associated with acquired resistance to EGFR-TKIs and that its suppression serves as a novel therapeutic strategy to overcome drug resistance in these patients.

Disseminated Medulloblastoma in a Child with Germline BRCA2 6174delT Mutation and without Fanconi Anemia.

Medulloblastoma, the most common malignant brain tumor in children, occurs with increased frequency in individuals with Fanconi anemia who have biallelic germline mutations in BRCA2. We describe an 8-year-old child who had disseminated anaplastic medulloblastoma and a deleterious heterozygous BRCA2 6174delT germline mutation. Molecular profiling was consistent with Group 4 medulloblastoma. The posterior fossa mass was resected and the patient received intensive chemotherapy and craniospinal irradiation. Despite this, the patient succumbed to a second recurrence of his medulloblastoma, which presented 8 months after diagnosis as malignant pleural and peritoneal effusions. Continuous medulloblastoma cell lines were isolated from the original tumor (CHLA-01-MED) and the malignant pleural effusion (CHLA-01R-MED). Here, we provide their analyses, including in vitro and in vivo growth, drug sensitivity, comparative genomic hybridization, and next generation sequencing analysis. In addition to the BRCA2 6174delT, the medulloblastoma cells had amplification of MYC, deletion at Xp11.2, and isochromosome 17, but no structural variations or overexpression of GFI1 or GFI1B. To our knowledge, this is the first pair of diagnosis/recurrence medulloblastoma cell lines, the only medulloblastoma cell lines with BRCA2 6174delT described to date, and the first reported case of a child with medulloblastoma associated with a germline BRCA2 6174delT who did not also have Fanconi anemia.

STAG2 is a clinically relevant tumor suppressor in pancreatic ductal adenocarcinoma.

BACKGROUND: Pancreatic ductal adenocarcinoma (PDA) is a highly lethal cancer characterized by complex aberrant genomes. A fundamental goal of current studies is to identify those somatic events arising in the variable landscape of PDA genomes that can be exploited for improved clinical outcomes. METHODS: We used DNA content flow sorting to identify and purify tumor nuclei of PDA samples from 50 patients. The genome of each sorted sample was profiled by oligonucleotide comparative genomic hybridization and targeted resequencing of STAG2. Transposon insertions within STAG2 in a KRAS (G12D)-driven genetically engineered mouse model of PDA were screened by RT-PCR. We then used a tissue microarray to survey STAG2 protein expression levels in 344 human PDA tumor samples and adjacent tissues. Univariate Kaplan Meier analysis and multivariate Cox Regression analysis were used to assess the association of STAG2 expression relative to overall survival and response to adjuvant therapy. Finally, RNAi-based assays with PDA cell lines were used to assess the potential therapeutic consequence of STAG2 expression in response to 18 therapeutic agents. RESULTS: STAG2 is targeted by somatic aberrations in a subset (4%) of human PDAs. Transposon-mediated disruption of STAG2 in a KRAS (G12D) genetically engineered mouse model promotes the development of PDA and its progression to metastatic disease. There was a statistically significant loss of STAG2 protein expression in human tumor tissue (Wilcoxon-Rank test) with complete absence of STAG2 staining observed in 15 (4.3%) patients. In univariate Kaplan Meier analysis nearly complete STAG2 positive staining (>95% of nuclei positive) was associated with a median survival benefit of 6.41 months (P = 0.031). The survival benefit of adjuvant chemotherapy was only seen in patients with a STAG2 staining of less than 95% (median survival benefit 7.65 months; P = 0.028). Multivariate Cox Regression analysis showed that STAG2 is an independent prognostic factor for survival in pancreatic cancer patients. Finally, we show that RNAi-mediated knockdown of STAG2 selectively sensitizes human PDA cell lines to platinum-based therapy. CONCLUSIONS: Based on these iterative findings we propose that STAG2 is a clinically significant tumor suppressor in PDA.

Hypoxia triggers hedgehog-mediated tumor-stromal interactions in pancreatic cancer.

Pancreatic cancer is characterized by a desmoplastic reaction that creates a dense fibroinflammatory microenvironment, promoting hypoxia and limiting cancer drug delivery due to decreased blood perfusion. Here, we describe a novel tumor-stroma interaction that may help explain the prevalence of desmoplasia in this cancer. Specifically, we found that activation of hypoxia-inducible factor-1α (HIF-1α) by tumor hypoxia strongly activates secretion of the sonic hedgehog (SHH) ligand by cancer cells, which in turn causes stromal fibroblasts to increase fibrous tissue deposition. In support of this finding, elevated levels of HIF-1α and SHH in pancreatic tumors were determined to be markers of decreased patient survival. Repeated cycles of hypoxia and desmoplasia amplified each other in a feed forward loop that made tumors more aggressive and resistant to therapy. This loop could be blocked by HIF-1α inhibition, which was sufficient to block SHH production and hedgehog signaling. Taken together, our findings suggest that increased HIF-1α produced by hypoxic tumors triggers the desmoplasic reaction in pancreatic cancer, which is then amplified by a feed forward loop involving cycles of decreased blood flow and increased hypoxia. Our findings strengthen the rationale for testing HIF inhibitors and may therefore represent a novel therapeutic option for pancreatic cancer.

SMG1 identified as a regulator of Parkinson's disease-associated alpha-synuclein through siRNA screening.

Synucleinopathies are a broad class of neurodegenerative disorders characterized by the presence of intracellular protein aggregates containing α-synuclein protein. The aggregated α-synuclein protein is hyperphosphorylated on serine 129 (S129) compared to the unaggregated form of the protein. While the precise functional consequences of S129 hyperphosphorylation are still being clarified, numerous in vitro and in vivo studies suggest that S129 phosphorylation is an early event in α-synuclein dysfunction and aggregation. Identifying the kinases and phosphatases that regulate this critical phosphorylation event may ultimately prove beneficial by allowing pharmacological mitigation of synuclein dysfunction and toxicity in Parkinson's disease and other synucleinopathies. We report here the development of a high-content, fluorescence-based assay to quantitate levels of total and S129 phosphorylated α-synuclein protein. We have applied this assay to conduct high-throughput loss-of-function screens with siRNA libraries targeting 711 known and predicted human kinases and 206 phosphatases. Specifically, knockdown of the phosphatidylinositol 3-kinase related kinase SMG1 resulted in significant increases in the expression of pS129 phosphorylated α-synuclein (p-syn). Moreover, SMG1 protein levels were significantly reduced in brain regions with high p-syn levels in both dementia with Lewy bodies (DLB) and Parkinson's disease with dementia (PDD). These findings suggest that SMG1 may play an important role in increased α-synuclein pathology during the course of PDD, DLB, and possibly other synucleinopathies.

Specification, annotation, visualization and simulation of a large rule-based model for ERBB receptor signaling.

BACKGROUND: Mathematical/computational models are needed to understand cell signaling networks, which are complex. Signaling proteins contain multiple functional components and multiple sites of post-translational modification. The multiplicity of components and sites of modification ensures that interactions among signaling proteins have the potential to generate myriad protein complexes and post-translational modification states. As a result, the number of chemical species that can be populated in a cell signaling network, and hence the number of equations in an ordinary differential equation model required to capture the dynamics of these species, is prohibitively large. To overcome this problem, the rule-based modeling approach has been developed for representing interactions within signaling networks efficiently and compactly through coarse-graining of the chemical kinetics of molecular interactions. RESULTS: Here, we provide a demonstration that the rule-based modeling approach can be used to specify and simulate a large model for ERBB receptor signaling that accounts for site-specific details of protein-protein interactions. The model is considered large because it corresponds to a reaction network containing more reactions than can be practically enumerated. The model encompasses activation of ERK and Akt, and it can be simulated using a network-free simulator, such as NFsim, to generate time courses of phosphorylation for 55 individual serine, threonine, and tyrosine residues. The model is annotated and visualized in the form of an extended contact map. CONCLUSIONS: With the development of software that implements novel computational methods for calculating the dynamics of large-scale rule-based representations of cellular signaling networks, it is now possible to build and analyze models that include a significant fraction of the protein interactions that comprise a signaling network, with incorporation of the site-specific details of the interactions. Modeling at this level of detail is important for understanding cellular signaling.

Functional genomics reveals diverse cellular processes that modulate tumor cell response to oxaliplatin.

Oxaliplatin is widely used to treat colorectal cancer, as both adjuvant therapy for resected disease and palliative treatment of metastatic disease. However, a significant number of patients experience serious side effects, including prolonged neurotoxicity, from oxaliplatin treatment creating an urgent need for biomarkers of oxaliplatin response or resistance to direct therapy to those most likely to benefit. As a first step to improve selection of patients for oxaliplatin-based chemotherapy, we have conducted an in vitro cell-based small interfering RNA (siRNA) screen of 500 genes aimed at identifying genes whose loss of expression alters tumor cell response to oxaliplatin. The siRNA screen identified twenty-seven genes, which when silenced, significantly altered colon tumor cell line sensitivity to oxaliplatin. Silencing of a group of putative resistance genes increased the extent of oxaliplatin-mediated DNA damage and inhibited cell-cycle progression in oxaliplatin-treated cells. The activity of several signaling nodes, including AKT1 and MEK1, was also altered. We used cDNA transfection to overexpress two genes (LTBR and TMEM30A) that were identified in the siRNA screen as mediators of oxaliplatin sensitivity. In both instances, overexpression conferred resistance to oxaliplatin. In summary, this study identified numerous putative predictive biomarkers of response to oxaliplatin that should be studied further in patient specimens for potential clinical application. Diverse gene networks seem to influence tumor survival in response to DNA damage by oxaliplatin. Finally, those genes whose loss of expression (or function) is related to oxaliplatin sensitivity may be promising therapeutic targets to increase patient response to oxaliplatin.