Mia D. Champion

Integrated genomic characterization reveals novel, therapeutically relevant drug targets in FGFR and EGFR pathways in sporadic intrahepatic cholangiocarcinoma.

Advanced cholangiocarcinoma continues to harbor a difficult prognosis and therapeutic options have been limited. During the course of a clinical trial of whole genomic sequencing seeking druggable targets, we examined six patients with advanced cholangiocarcinoma. Integrated genome-wide and whole transcriptome sequence analyses were performed on tumors from six patients with advanced, sporadic intrahepatic cholangiocarcinoma (SIC) to identify potential therapeutically actionable events. Among the somatic events captured in our analysis, we uncovered two novel therapeutically relevant genomic contexts that when acted upon, resulted in preliminary evidence of anti-tumor activity. Genome-wide structural analysis of sequence data revealed recurrent translocation events involving the FGFR2 locus in three of six assessed patients. These observations and supporting evidence triggered the use of FGFR inhibitors in these patients. In one example, preliminary anti-tumor activity of pazopanib (in vitro FGFR2 IC50≈350 nM) was noted in a patient with an FGFR2-TACC3 fusion. After progression on pazopanib, the same patient also had stable disease on ponatinib, a pan-FGFR inhibitor (in vitro, FGFR2 IC50≈8 nM). In an independent non-FGFR2 translocation patient, exome and transcriptome analysis revealed an allele specific somatic nonsense mutation (E384X) in ERRFI1, a direct negative regulator of EGFR activation. Rapid and robust disease regression was noted in this ERRFI1 inactivated tumor when treated with erlotinib, an EGFR kinase inhibitor. FGFR2 fusions and ERRFI mutations may represent novel targets in sporadic intrahepatic cholangiocarcinoma and trials should be characterized in larger cohorts of patients with these aberrations.

Identification of cereblon-binding proteins and relationship with response and survival after IMiDs in multiple myeloma

Cereblon (CRBN) mediates immunomodulatory drug (IMiD) action in multiple myeloma (MM). Using 2 different methodologies, we identified 244 CRBN binding proteins and established relevance to MM biology by changes in their abundance after exposure to lenalidomide. Proteins most reproducibly binding CRBN (>fourfold vs controls) included DDB1, CUL4A, IKZF1, KPNA2, LTF, PFKL, PRKAR2A, RANGAP1, and SHMT2. After lenalidomide treatment, the abundance of 46 CRBN binding proteins decreased. We focused attention on 2 of these-IKZF1 and IKZF3. IZKF expression is similar across all MM stages or subtypes; however, IKZF1 is substantially lower in 3 of 5 IMiD-resistant MM cell lines. The cell line (FR4) with the lowest IKZF1 levels also harbors a damaging mutation and a translocation that upregulates IRF4, an IKZF target. Clinical relevance of CRBN-binding proteins was demonstrated in 44 refractory MM patients treated with pomalidomide and dexamethasone therapy in whom low IKZF1 gene expression predicted lack of response (0/11 responses in the lowest expression quartile). CRBN, IKZF1, and KPNA2 levels also correlate with significant differences in overall survival. Our study identifies CRBN-binding proteins and demonstrates that in addition to CRBN, IKZF1, and KPNA2, expression can predict survival outcomes.

IL-33 facilitates oncogene-induced cholangiocarcinoma in mice by an interleukin-6-sensitive mechanism.

Cholangiocarcinoma (CCA) is a lethal hepatobiliary neoplasm originating from the biliary apparatus. In humans, CCA risk factors include hepatobiliary inflammation and fibrosis. The recently identified interleukin (IL)−1 family member, IL‐33, has been shown to be a biliary mitogen which also promotes liver inflammation and fibrosis. Our aim was to generate a mouse model of CCA mimicking the human disease. Ectopic oncogene expression in the biliary tract was accomplished by the Sleeping Beauty transposon transfection system with transduction of constitutively active AKT (myr‐AKT) and Yes‐associated protein. Intrabiliary instillation of the transposon–transposase complex was coupled with lobar bile duct ligation in C57BL/6 mice, followed by administration of IL‐33 for 3 consecutive days. Tumors developed in 72% of the male mice receiving both oncogenes plus IL‐33 by 10 weeks but in only 20% of the male mice transduced with the oncogenes alone. Tumors expressed SOX9 and pancytokeratin (features of CCA) but were negative for HepPar1 (a marker of hepatocellular carcinoma). Substantive overlap with human CCA specimens was revealed by RNA profiling. Not only did IL‐33 induce IL‐6 expression by human cholangiocytes but it likely facilitated tumor development in vivo by an IL‐6–sensitive process as tumor development was significantly attenuated in Il‐6–/– male animals. Furthermore, tumor formation occurred at a similar rate when IL‐6 was substituted for IL‐33 in this model. Conclusion: The transposase‐mediated transduction of constitutively active AKT and Yes‐associated protein in the biliary epithelium coupled with lobar obstruction and IL‐33 administration results in the development of CCA with morphological and biochemical features of the human disease; this model highlights the role of inflammatory cytokines in CCA oncogenesis. (Hepatology 2015;61:1627–1642)

The clinical significance of cereblon expression in multiple myeloma

Cereblon (CRBN) mediates immunomodulatory drug (IMiD) action in multiple myeloma (MM). We demonstrate here that no patient with very low CRBN expression responded to IMiD plus dexamethasone therapy. In 53 refractory MM patients treated with pomalidomide and dexamethasone, CRBN levels predict for decreased response rates and significant differences in PFS (3.0 vs. 8.9 months, p<0.001) and OS (9.1 vs. 27.2 months, p=0.01) (lowest quartile vs. highest three quartiles). While higher CRBN levels can serve as a surrogate for low risk disease, our study demonstrates that low CRBN expression can predict resistance to IMiD monotherapy and is a predictive biomarker for survival outcomes.

High-resolution profiling of histone h3 lysine 36 trimethylation in metastatic renal cell carcinoma

Mutations in SETD2, a histone H3 lysine trimethyltransferase, have been identified in clear cell renal cell carcinoma (ccRCC); however it is unclear if loss of SETD2 function alters the genomic distribution of histone 3 lysine 36 trimethylation (H3K36me3) in ccRCC. Furthermore, published epigenomic profiles are not specific to H3K36me3 or metastatic tumors. To determine if progressive SETD2 and H3K36me3 dysregulation occurs in metastatic tumors, H3K36me3, SETD2 copy number (CN) or SETD2 mRNA abundance was assessed in two independent cohorts: metastatic ccRCC (n=71) and the Cancer Genome Atlas Kidney Renal Clear Cell Carcinoma data set (n=413). Although SETD2 CN loss occurs with high frequency (>90%), H3K36me3 is not significantly impacted by monoallelic loss of SETD2. H3K36me3-positive nuclei were reduced an average of ~20% in primary ccRCC (90% positive nuclei in uninvolved vs 70% positive nuclei in ccRCC) and reduced by ~60% in metastases (90% positive in uninvolved kidney vs 30% positive in metastases) (P<0.001). To define a kidney-specific H3K36me3 profile, we generated genome-wide H3K36me3 profiles from four cytoreductive nephrectomies and SETD2 isogenic renal cell carcinoma (RCC) cell lines using chromatin immunoprecipitation coupled with high-throughput DNA sequencing and RNA sequencing. SETD2 loss of methyltransferase activity leads to regional alterations of H3K36me3 associated with aberrant RNA splicing in a SETD2 mutant RCC and SETD2 knockout cell line. These data suggest that during progression of ccRCC, a decline in H3K36me3 is observed in distant metastases, and regional H3K36me3 alterations influence alternative splicing in ccRCC.

Targeted sequencing using a 47 gene multiple myeloma mutation panel (M3P) in -17p high risk disease

We constructed a multiple myeloma (MM)‐specific gene panel for targeted sequencing and investigated 72 untreated high‐risk (del17p) MM patients. Mutations were identified in 78% of the patients. While the majority of studied genes were mutated at similar frequency to published literature, the prevalence of TP53 mutation was increased (28%) and no mutations were found in FAM46C. This study provides a comprehensive insight into the mutational landscape of del17p high‐risk MM. Additionally, our work demonstrates the practical use of a customized sequencing panel, as an easy, cheap and fast approach to characterize the mutational profile of MM.

Clinical Implementation of Integrated Genomic Profiling in Patients with Advanced Cancers

DNA focused panel sequencing has been rapidly adopted to assess therapeutic targets in advanced/refractory cancer. Integrated Genomic Profiling (IGP) utilising DNA/RNA with tumour/normal comparisons in a Clinical Laboratory Improvement Amendments (CLIA) compliant setting enables a single assay to provide: therapeutic target prioritisation, novel target discovery/application and comprehensive germline assessment. A prospective study in 35 advanced/refractory cancer patients was conducted using CLIA-compliant IGP. Feasibility was assessed by estimating time to results (TTR), prioritising/assigning putative therapeutic targets, assessing drug access, ascertaining germline alterations, and assessing patient preferences/perspectives on data use/reporting. Therapeutic targets were identified using biointelligence/pathway analyses and interpreted by a Genomic Tumour Board. Seventy-five percent of cases harboured 1–3 therapeutically targetable mutations/case (median 79 mutations of potential functional significance/case). Median time to CLIA-validated results was 116 days with CLIA-validation of targets achieved in 21/22 patients. IGP directed treatment was instituted in 13 patients utilising on/off label FDA approved drugs (n = 9), clinical trials (n = 3) and single patient IND (n = 1). Preliminary clinical efficacy was noted in five patients (two partial response, three stable disease). Although barriers to broader application exist, including the need for wider availability of therapies, IGP in a CLIA-framework is feasible and valuable in selection/prioritisation of anti-cancer therapeutic targets.

A Multidisciplinary Biospecimen Bank of Renal Cell Carcinomas Compatible with Discovery Platforms at Mayo Clinic, Scottsdale, Arizona

To address the need to study frozen clinical specimens using next-generation RNA, DNA, chromatin immunoprecipitation (ChIP) sequencing and protein analyses, we developed a biobank work flow to prospectively collect biospecimens from patients with renal cell carcinoma (RCC). We describe our standard operating procedures and work flow to annotate pathologic results and clinical outcomes. We report quality control outcomes and nucleic acid yields of our RCC submissions (N=16) to The Cancer Genome Atlas (TCGA) project, as well as newer discovery platforms, by describing mass spectrometry analysis of albumin oxidation in plasma and 6 ChIP sequencing libraries generated from nephrectomy specimens after histone H3 lysine 36 trimethylation (H3K36me3) immunoprecipitation. From June 1, 2010, through January 1, 2013, we enrolled 328 patients with RCC. Our mean (SD) TCGA RNA integrity numbers (RINs) were 8.1 (0.8) for papillary RCC, with a 12.5% overall rate of sample disqualification for RIN <7. Banked plasma had significantly less albumin oxidation (by mass spectrometry analysis) than plasma kept at 25°C (P<.001). For ChIP sequencing, the FastQC score for average read quality was at least 30 for 91% to 95% of paired-end reads. In parallel, we analyzed frozen tissue by RNA sequencing; after genome alignment, only 0.2% to 0.4% of total reads failed the default quality check steps of Bowtie2, which was comparable to the disqualification ratio (0.1%) of the 786-O RCC cell line that was prepared under optimal RNA isolation conditions. The overall correlation coefficients for gene expression between Mayo Clinic vs TCGA tissues ranged from 0.75 to 0.82. These data support the generation of high-quality nucleic acids for genomic analyses from banked RCC. Importantly, the protocol does not interfere with routine clinical care. Collections over defined time points during disease treatment further enhance collaborative efforts to integrate genomic information with outcomes.

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.

Whole Genome Analyses of a Well-Differentiated Liposarcoma Reveals Novel SYT1 and DDR2 Rearrangements

Liposarcoma is the most common soft tissue sarcoma, but little is known about the genomic basis of this disease. Given the low cell content of this tumor type, we utilized flow cytometry to isolate the diploid normal and aneuploid tumor populations from a well-differentiated liposarcoma prior to array comparative genomic hybridization and whole genome sequencing. This work revealed massive highly focal amplifications throughout the aneuploid tumor genome including MDM2, a gene that has previously been found to be amplified in well-differentiated liposarcoma. Structural analysis revealed massive rearrangement of chromosome 12 and 11 gene fusions, some of which may be part of double minute chromosomes commonly present in well-differentiated liposarcoma. We identified a hotspot of genomic instability localized to a region of chromosome 12 that includes a highly conserved, putative L1 retrotransposon element, LOC100507498 which resides within a gene cluster (NAV3, SYT1, PAWR) where 6 of the 11 fusion events occurred. Interestingly, a potential gene fusion was also identified in amplified DDR2, which is a potential therapeutic target of kinase inhibitors such as dastinib, that are not routinely used in the treatment of patients with liposarcoma. Furthermore, 7 somatic, damaging single nucleotide variants have also been identified, including D125N in the PTPRQ protein. In conclusion, this work is the first to report the entire genome of a well-differentiated liposarcoma with novel chromosomal rearrangements associated with amplification of therapeutically targetable genes such as MDM2 and DDR2.