Search

Supplementary Tables 1-2 from Genetic Polymorphisms in Vitamin D Receptor VDR/RXRA Influence the Likelihood of Colon Adenoma Recurrence

Genome-wide fragmentation patterns in cell-free DNA (cfDNA) in plasma are strongly influenced by cellular origin due to variation in chromatin accessibility across cell types. Such differences between healthy and cancer cells provide the opportunity for development of novel cancer diagnostics. Here, we investigated whether analysis of cfDNA fragment end positions and their surrounding DNA sequences reveals the presence of tumor-derived DNA in blood. We performed genome-wide analysis of cfDNA from 521 samples and analyzed sequencing data from an additional 2147 samples, including healthy individuals and patients with 11 different cancer types. We developed a metric based on genome-wide differences in fragment positioning, weighted by fragment length and GC content [information-weighted fraction of aberrant fragments (iwFAF)]. We observed that iwFAF strongly correlated with tumor fraction, was higher for DNA fragments carrying somatic mutations, and was higher within genomic regions affected by copy number amplifications. We also calculated sample-level means of nucleotide frequencies observed at genomic positions spanning fragment ends. Using a combination of iwFAF and nine nucleotide frequencies from three positions surrounding fragment ends, we developed a machine learning model to differentiate healthy individuals from patients with cancer. We observed an area under the receiver operative characteristic curve (AUC) of 0.91 for detection of cancer at any stage and an AUC of 0.87 for detection of stage I cancer. Our findings remained robust with as few as 1 million fragments analyzed per sample, demonstrating that analysis of fragment ends can become a cost-effective and accessible approach for cancer detection and monitoring.

Background Altered tumor suppressor genes (TSG-alt) in prostate cancer are associated with worse outcomes. The prognostic value of TSG-alt in metastatic, hormone-sensitive prostate cancer (M1-HSPC) is unknown. We evaluated the effects of TSG-alt on outcomes in M1-HSPC and their prognostic impact by first-line treatment. Methods We retrospectively identified patients with M1-HSPC at our institution treated with first-line androgen deprivation therapy plus docetaxel (ADT + D) or abiraterone acetate (ADT + A). TSG-alt was defined as any alteration in one or more TSG. The main outcomes were Kaplan–Meier-estimated progression-free survival (PFS) and overall survival, analyzed with the log-rank test. Clinical characteristics were compared with the χ2 test and Kruskal–Wallis rank sum test. Cox regression was used for univariate and multivariable analyses. Results We identified 97 patients with M1-HSPC: 48 (49%) with ADT + A and 49 (51%) with ADT + D. Of 96 patients with data available, 33 (34%) had 1 TSG-alt, 16 (17%) had 2 TSG-alt, and 2 (2%) had 3 TSG-alt. The most common alterations were in TP53 (36%) and PTEN (31%); 6% had RB1 alterations. Median PFS was 13.1 (95% CI, 10.3–26.0) months for patients with normal TSGs (TSG-normal) vs. 7.8 (95% CI, 5.8–10.5) months for TSG-alt (P = 0.005). Median PFS was lower for patients with TSG-alt vs TSG-normal for those with ADT + A (TSG-alt: 8.0 [95% CI, 5.8–13.8] months vs. TSG-normal: 23.2 [95% CI, 13.1–not estimated] months), but not with ADT + D (TSG-alt: 7.8 [95% CI, 5.7–12.9] months vs. TSG-normal: 9.5 [95% CI, 4.8–24.7] months). On multivariable analysis, only TSG-alt predicted worse PFS (hazard ratio, 2.37; 95% CI, 1.42–3.96; P Conclusions The presence of TSG-alt outperforms clinical criteria for predicting early progression during first-line treatment of M1-HSPC. ADT + A was less effective in patients with than without TSG-alt. Confirmation of these findings may establish the need for inclusion of molecular stratification in treatment algorithms.

Background The goal of this study is to evaluate germline genetic variants in African American men with metastatic prostate cancer as compared to those in Caucasian men with metastatic prostate cancer in an effort to understand the role of genetic factors in these populations. Methods African American and Caucasian men with metastatic prostate cancer who had germline testing using multigene panels were used to generate comparisons. Germline genetic results, clinical parameters, and family histories between the two populations were analyzed. Results A total of 867 patients were included in this retrospective study, including 188 African American and 669 Caucasian patients. There was no significant difference in the likelihood of a pathogenic or likely‐pathogenic variants (PV/LPVs) between African American and Caucasian patients (p = .09). African American patients were more likely to have a variant of unknown significance than Caucasians (odds ratio [OR] = 1.95; p

NUC1031 is a first-in-class ProTide, that is a gemcitabine pro-drug designed to overcome putative mechanisms of resistance, including decreased expression of hENT/hCNT transporters, absence of activating enzymes such as deoxycytidine kinase (dCK) and presence of degrading enzymes such as cytidine deaminase (CDA). We undertook comprehensive pre-clinical evaluation of NUC1031 in biliary tract cancer (BTC) models, given that gemcitabine/cisplatin is a standard first-line therapy in advanced BTC. Here, we compared the in vitro activity of NUC1031 in comparison to gemcitabine, validate putative mechanism(s) of action, assessed potential biomarkers of sensitivity or resistance, and performed combination studies with cisplatin. We also evaluated the in vivo efficacy of NUC1031 and gemcitabine using a CDA-high cholangiocarcinoma patient-derived xenograft (PDX) model. In a panel of BTC cell lines (N = 10), NUC1031 had less potency than gemcitabine in multiple cellular assays. NUC1031 did not demonstrate evidence of greater synergy over gemcitabine in combination with cisplatin. Surprisingly, efficacy of both gemcitabine and NUC1031 was not found to be correlated with hENT/hCTN, dCK or CDA transcript levels. Gemcitabine and NUC1031 showed equivalent efficacy in a CDA-high PDX model in vivo contradicting the primary rationale of NUC1031 design. NUC1031 did not exhibit evidence of superior activity over gemcitabine, as a single-agent, or in combination with cisplatin, in either our in vivo or in vitro BTC models. Given that the largest Phase 3 study (ClinicalTrials.gov: NCT0314666) to date in BTC is underway (N = 828) comparing NUC1031/cisplatin to gemcitabine/cisplatin, our results suggest that a more conservative clinical evaluation path would be more appropriate.

Biliary tract cancers (BTCs) are uncommon, but highly lethal, gastrointestinal malignancies. Gemcitabine/cisplatin is a standard-of-care systemic therapy, but has a modest impact on survival and harbors toxicities, including myelosuppression, nephropathy, neuropathy, and ototoxicity. Whereas BTCs are characterized by aberrations activating the cyclinD1/cyclin-dependent kinase (CDK)4/6/CDK inhibitor 2a/retinoblastoma pathway, clinical use of CDK4/6 inhibitors as monotherapy is limited by lack of validated biomarkers, diffident preclinical efficacy, and development of acquired drug resistance. Emerging studies have explored therapeutic strategies to enhance the antitumor efficacy of CDK4/6 inhibitors by the combination with chemotherapy regimens, but their mechanism of action remains elusive.Here, we report in vitro and in vivo synergy in BTC models, showing enhanced efficacy, reduced toxicity, and better survival with a combination comprising gemcitabine/cisplatin and CDK4/6 inhibitors. Furthermore, we demonstrated that abemaciclib monotherapy had only modest efficacy attributable to autophagy-induced resistance. Notably, triplet therapy was able to potentiate efficacy through elimination of the autophagic flux. Correspondingly, abemaciclib potentiated ribonucleotide reductase catalytic subunit M1 reduction, resulting in sensitization to gemcitabine.As such, these data provide robust preclinical mechanistic evidence of synergy between gemcitabine/cisplatin and CDK4/6 inhibitors and delineate a path forward for translation of these findings to preliminary clinical studies in advanced BTC patients.

Fragmentation patterns observed in plasma DNA reflect chromatin accessibility in contributing cells. Since DNA shed from cancer cells and blood cells may differ in fragmentation patterns, we investigated whether analysis of genomic positioning and nucleotide sequence at fragment ends can reveal the presence of tumor DNA in blood and aid cancer diagnostics. We analyzed whole genome sequencing data from >2700 plasma DNA samples including healthy individuals and patients with 11 different cancer types. We observed higher fractions of fragments with aberrantly positioned ends in patients with cancer, driven by contribution of tumor DNA into plasma. Genomewide analysis of fragment ends using machine learning showed overall area under the receiver operative characteristic curve of 0.96 for detection of cancer. Our findings remained robust with as few as 1 million fragments analyzed per sample, suggesting that analysis of fragment ends can become a cost-effective and accessible approach for cancer detection and monitoring.One-sentence summaryAnalyzing the positioning and nucleotide sequence at fragment ends in plasma DNA may enable cancer diagnostics.

Post-pubertal testicular germ-cell tumours (TGCTs) can present with a variety of distinct histologies which are nevertheless lineage related and often co-occurring. The exact lineage relationships and developmental pathways leading to the different histologies is debated. In order to investigate the relationship of histologic populations, mate-pair sequencing (MPseq) and exome sequencing (ExomeSeq) were conducted on different histological populations within the same tumour. Ten TGCTs with 1–3 histologic types/tumour were sequenced. Junctions of somatic chromosomal rearrangements were identified on a per genome basis, with germ cell neoplasia in situ possessing the least (median 1, range 0–4) and embryonal carcinoma the most (median 8.5, range 6–12). Copy number variation revealed gains and losses, including isoform 12p (i12p) (10/10 samples), and chromosomes 7, 8, and 21 gains (7/10 samples). Mapping of shared junctions within a tumour revealed lineage relationships, but only i12p was shared between patients. ExomeSeq from two cases demonstrated a high level of copy-neutral loss of heterozygosity. Parallel assessment of separate histologies within a single TGCT demonstrated cumulative and divergent changes, suggesting the importance of parallel sequencing for detection of relevant biomarkers.

About 15% of intrahepatic cholangiocarcinomas (ICCs) express constitutively active fibroblast growth factor receptor 2 (FGFR2) fusion proteins (FFs) generated by chromosomal translocations. FFs have been nominated as oncogenic drivers because administration of FGFR tyrosine kinase inhibitors (F-TKIs) can elicit meaningful objective clinical responses in patients carrying FF-positive ICC. Thus, optimization of FF targeting is a pressing clinical need. Herein, we report that three different FFs, previously isolated from ICC samples, are heat shock protein 90 (HSP90) clients and undergo rapid degradation upon HSP90 pharmacological blockade by the clinically advanced HSP90 inhibitor ganetespib. Combining catalytic suppression by the F-TKI BGJ398 with HSP90 blockade by ganetespib suppressed FGFR2-TACC3 (transforming acidic coiled-coil containing protein 3) signaling in cultured cells more effectively than either BGJ398 or ganetespib in isolation. The BGJ398 + ganetespib combo was also superior to single agents when tested in mice carrying subcutaneous tumors generated by transplantation of FGFR2-TACC3 NIH3T3 transformants. Of note, FF mutants known to enforce clinical resistance to BGJ398 in ICC patients retained full sensitivity to ganetespib in cultured cells. Conclusion: Our data provide a proof of principle that upfront treatment with the BGJ398 + ganetespib combo improves therapeutic targeting of FGFR2 fusions in an experimental setting, which may be relevant to precision medicine approaches to FF-driven ICC.

Gestational trophoblastic disease (GTD) is a heterogeneous group of lesions arising from placental tissue. Epithelioid trophoblastic tumor (ETT), derived from chorionic-type trophoblast, is the rarest form of GTD with only approximately 130 cases described in the literature. Due to its morphologic mimicry of epithelioid smooth muscle tumors and carcinoma, ETT can be misdiagnosed. To date, molecular characterization of ETTs is lacking. Furthermore, ETT is difficult to treat when disease spreads beyond the uterus. Here using RNA-Seq analysis in a cohort of ETTs and other gestational trophoblastic lesions we describe the discovery of LPCAT1-TERT fusion transcripts that occur in ETTs and coincide with underlying genomic deletions. Through cell-growth assays we demonstrate that LPCAT1-TERT fusion proteins can positively modulate cell proliferation and therefore may represent future treatment targets. Furthermore, we demonstrate that TERT upregulation appears to be a characteristic of ETTs, even in the absence of LPCAT1-TERT fusions, and that it appears linked to copy number gains of chromosome 5. No evidence of TERT upregulation was identified in other trophoblastic lesions tested, including placental site trophoblastic tumors and placental site nodules, which are thought to be the benign chorionic-type trophoblast counterpart to ETT. These findings indicate that LPCAT1-TERT fusions and copy-number driven TERT activation may represent novel markers for ETT, with the potential to improve the diagnosis, treatment, and outcome for women with this rare form of GTD.

Biliary tract cancers (BTCs) are uncommon but highly lethal gastrointestinal malignancies. Gemcitabine/cisplatin is a standard-of-care (SOC) systemic therapy, but has a modest impact on survival and harbor toxicities including myelosuppression, nephropathy, neuropathy and ototoxicity. While BTCs are characterized by aberrations activating the cyclinD1-CDK4/6-CDKN2A-RB pathway, clinical use of CDK4/6 inhibitors as monotherapy is limited by lack of validated biomarkers, diffident pre-clinical efficacy and development of acquired drug resistance. Emerging studies have explored therapeutic strategies to enhance the anti-tumor efficacy of CDK4/6 inhibitors by combination with chemotherapy-regimens but their mechanism of action remains elusive. Here, we report in vitro and in vivo synergy in BTC models, showing enhanced efficacy, reduced toxicity and better survival with a combination comprising gemcitabine/cisplatin and CDK4/6 inhibitors. Furthermore, we demonstrated that abemaciclib monotherapy had only modest efficacy due to autophagy induced resistance. Notably, triplettherapy was able to potentiate efficacy through elimination of the autophagic flux. Correspondingly, abemaciclib potentiated RRM1 reduction, resulting in sensitization to gemcitabine. Conclusions: As such, these data provide robust pre-clinical mechanistic evidence of synergy between gemcitabine/cisplatin and CDK4/6 inhibitors, and delineate a path forward for translation of these findings to preliminary clinical studies in advanced BTC patients.

NUC1031 is a first-in-class ProTide, that is a gemcitabine pro-drug designed to overcome putative mechanisms of resistance, including decreased expression of hENT/hCNT transporters, absence of activating enzymes such as deoxycytidine kinase (dCK) and presence of degrading enzymes such as cytidine deaminase (CDA). We undertook comprehensive pre-clinical evaluation of NUC1031 in biliary tract cancer (BTC) models, given that gemcitabine/cisplatin is a standard first-line therapy in advanced BTC.Here, we compared the in vitro activity of NUC1031 in comparison to gemcitabine, validate putative mechanism(s) of action, assessed potential biomarkers of sensitivity or resistance, and performed combination studies with cisplatin. We also evaluated the in vivo efficacy of NUC1031 and gemcitabine using a CDA-high cholangiocarcinoma patient-derived xenograft (PDX) model.In a panel of BTC cell lines (N = 10), NUC1031 had less potency than gemcitabine in multiple cellular assays. NUC1031 did not demonstrate evidence of greater synergy over gemcitabine in combination with cisplatin. Surprisingly, efficacy of both gemcitabine and NUC1031 was not found to be correlated with hENT/hCTN, dCK or CDA transcript levels. Gemcitabine and NUC1031 showed equivalent efficacy in a CDA-high PDX model in vivo contradicting the primary rationale of NUC1031 design.NUC1031 did not exhibit evidence of superior activity over gemcitabine, as a single-agent, or in combination with cisplatin, in either our in vivo or in vitro BTC models. Given that the largest Phase 3 study (ClinicalTrials.gov: NCT0314666) to date in BTC is underway (N = 828) comparing NUC1031/cisplatin to gemcitabine/cisplatin, our results suggest that a more conservative clinical evaluation path would be more appropriate.

Upper gastrointestinal tract malignancies are among the most challenging cancers with regard to response to treatment and prognosis. Cancers of the esophagus, stomach, pancreas, liver, and biliary tree have dismal 5-year survival, and very modest improvements in this rate have been made in recent times. Oncolytic viruses are being developed to address these malignancies, with a focus on high safety profiles and low off-target toxicities. Each viral platform has evolved to enhance oncolytic potency and the clinical response to either single-agent viral therapy or combined viral treatment with radiotherapy and chemotherapy. A panel of genomic alterations, chimeric proteins, and pseudotyped capsids are the breakthroughs for vector success. This article revisits developments for each viral platform to each tumor type, in an attempt to achieve maximum tumor selectivity. From the bench to clinical trials, the scope of this review is to highlight the beginnings of translational oncolytic virotherapy research in upper gastrointestinal tract malignancies and provide a bioengineering perspective of the most promising platforms.

With the advancement of a growing number of oncolytic viruses (OVs) to clinical development, drug delivery is becoming an important barrier to overcome for optimal therapeutic benefits. Host immunity, tumor microenvironment and abnormal vascularity contribute to inefficient vector delivery. A number of novel approaches for enhanced OV delivery are under evaluation, including use of nanoparticles, immunomodulatory agents and complex viral-particle ligands along with manipulations of the tumor microenvironment. This field of OV delivery has quickly evolved to bioengineering of complex nanoparticles that could be deposited within the tumor using minimal invasive image-guided delivery. Some of the strategies include ultrasound (US)-mediated cavitation-enhanced extravasation, magnetic viral complexes delivery, image-guided infusions with focused US and targeting photodynamic virotherapy. In addition, strategies that modulate tumor microenvironment to decrease extracellular matrix deposition and increase viral propagation are being used to improve tumor penetration by OVs. Some involve modification of the viral genome to enhance their tumoral penetration potential. Here, we highlight the barriers to oncolytic viral delivery, and discuss the challenges to improving it and the perspectives of establishing new modes of active delivery to achieve enhanced oncolytic effects.

e16185 Background: Comprehensive sequencing efforts have identified fibroblast growth factor receptor (FGFR) translocations in numerous cancers. FGFR2 fusions were recently identified as an actionable therapeutic target in ̃17% of cholangiocarcinoma patients, predominantly found in patients with intrahepatic cholangiocarcinoma. Genomic partners for FGFR fusions are variable. In addition, obtaining high yield tissue biopsies for CCA remains challenging and tumor tissue available for molecular analyses is scarce. When combined, these two factors make detection of FGFR fusions in CCA challenging. An accurate molecular assay to detect FGFR fusions in CCA could improve utilization of FGFR-targeted therapies. Methods: The objective of our study was to develop an approach for targeted RNA analysis to improve sensitivity for FGFR fusion detection from limited tumor material. We adapted the sensitive targeted digital sequencing (TARDIS) approach recently developed in our lab for analysis of RNA. This approach utilizes a combination of open-ended targeted amplification followed by ligation, enabling detection of tyrosine kinase fusions without prior knowledge of the precise sequence of the fusion breakpoint or identity of the fusion partner. For evaluation of analytical performance, we analyzed RNA from a CCA organoid model with a known FGFR2 fusion, and Seraseq Fusion RNA Mix v3 Reference Material with two known FGFR3 fusions. Results: Without prior knowledge the partner or coordinates of the fusion breakpoint, we detected an FGFR2-KIF5C fusion in RNA from a CCA organoid model. The fusion breakpoint coordinates predicted by TARDIS-RNA were validated by comparison with RNA-Seq. To assess sensitivity and quantitative performance of the assay, we analyzed a serial dilution of the reference RNA sample from Seraseq with concentration of RNA molecules representing candidate fusions verified by digital PCR. Using TARDIS-RNA to analyze replicates with 5 ng RNA input, we were able to detect candidate fusions represented by as few as 2.7 fusion molecules on average. Concentrations of fusion molecules measured using the two methods were highly correlated (R = 0.96). Conclusions: These results demonstrate sensitivity and quantitative performance of a targeted RNA sequencing assay to detect and quantify FGFR fusions in tumor specimen from intrahepatic cholangiocarcinoma. On-going work is focused on further evaluating assay performance and characterizing FGFR fusions using additional tissue specimen.

// Alan H. Bryce 1, 2, 3, * , Jan B. Egan 3, * , Mitesh J. Borad 1, 2, 3 , A. Keith Stewart 1, 2, 3 , Grzegorz S. Nowakowski 3, 4 , Asher Chanan-Khan 3, 5 , Mrinal M. Patnaik 3, 4 , Stephen M. Ansell 3, 4 , Michaela S. Banck 3, 6 , Steven I. Robinson 3, 6 , Aaron S. Mansfield 3, 6 , Eric W. Klee 3, 7 , Gavin R. Oliver 3, 7 , Jennifer B. McCormick 7 , Norine E. Huneke 3 , Colleen M. Tagtow 3 , Robert B. Jenkins 8 , Kandelaria M. Rumilla 8 , Sarah E. Kerr 9 , Jean-Pierre A. Kocher 3, 7 , Scott A. Beck 3 , Martin E. Fernandez-Zapico 10 , Gianrico Farrugia 3, 11 , Konstantinos N. Lazaridis 3, 11 , Robert R. McWilliams 3, 6, 12 1 Hematology/Oncology, Mayo Clinic, Phoenix, AZ, U.S.A 2 Mayo Clinic Cancer Center, Phoenix, AZ, U.S.A 3 Center for Individualized Medicine, Mayo Clinic, Rochester, MN, U.S.A 4 Hematology, Mayo Clinic, Rochester, MN, U.S.A 5 Hematology/Oncology, Mayo Clinic, Jacksonville, FL, U.S.A 6 Medical Oncology, Mayo Clinic, Rochester, MN, U.S.A 7 Health Sciences Research, Mayo Clinic, Rochester, MN, U.S.A 8 Laboratory Medicine and Pathology, Mayo Clinic, Rochester, MN, U.S.A 9 Anatomic Pathology, Mayo Clinic, Rochester, MN, U.S.A 10 Schulze Center for Novel Therapeutics, Division of Oncology Research, Medical Oncology, Mayo Clinic, Rochester, MN, U.S.A 11 Gastroenterology, Mayo Clinic, Rochester, MN, U.S.A 12 Mayo Clinic Cancer Center, Rochester, MN, U.S.A * These authors contributed equally to this work Correspondence to: Robert R. McWilliams, email: mcwilliams.robert@mayo.edu Keywords: precision medicine, cancer genomics, targeted therapeutics, genomic tumor board Received: November 09, 2016 Accepted: February 15, 2017 Published: March 09, 2017 ABSTRACT Background: The ability to analyze the genomics of malignancies has opened up new possibilities for off-label targeted therapy in cancers that are refractory to standard therapy. At Mayo Clinic these efforts are organized through the Center for Individualized Medicine (CIM). Results: Prior to GTB, datasets were analyzed and integrated by a team of bioinformaticians and cancer biologists. Therapeutically actionable mutations were identified in 65% (92/141) of the patients tested with 32% (29/92) receiving genomically targeted therapy with FDA approved drugs or in an independent clinical trial with 45% (13/29) responding. Standard of care (SOC) options were continued by 15% (14/92) of patients tested before exhausting SOC options, with 71% (10/14) responding to treatment. Over 35% (34/92) of patients with actionable targets were not treated with 65% (22/34) choosing comfort measures or passing away. Materials and Methods: Patients ( N = 165) were referred to the CIM Clinic between October 2012 and December 2015. All patients received clinical genomic panel testing with selected subsets receiving array comparative genomic hybridization and clinical whole exome sequencing to complement and validate panel findings. A genomic tumor board (GTB) reviewed results and, when possible, developed treatment recommendations. Conclusions: Treatment decisions driven by tumor genomic analysis can lead to significant clinical benefit in a minority of patients. The success of genomically driven therapy depends both on access to drugs and robustness of bioinformatics analysis. While novel clinical trial designs are increasing the utility of genomic testing, robust data sharing of outcomes is needed to optimize clinical benefit for all patients.

// Cheryl E. Myers 1 , Dominique B. Hoelzinger 1 , Tiffany N. Truong 1 , Lindsey A. Chew 1 , Arpita Myles 2 , Leena Chaudhuri 2 , Jan B. Egan 3 , Jun Liu 1 , Sandra J. Gendler 1, 2, 4 , Peter A. Cohen 2, 4 1 Department of Biochemistry and Molecular Biology, Mayo Clinic, Scottsdale, AZ, USA 2 Department of Immunology, Mayo Clinic, Scottsdale, AZ, USA 3 Center for Individualized Medicine, Mayo Clinic, Scottsdale, AZ, USA 4 Department of Hematology and Oncology, Mayo Clinic, Scottsdale, AZ, USA Correspondence to: Cheryl E. Myers, email: myers.cheryl@mayo.edu Peter A. Cohen, email: cohen.peter@mayo.edu Keywords: obesity, chemotherapy, proton leak, preadipocyte, lipogenesis Received: July 27, 2016 Accepted: December 05, 2016 Published: January 10, 2017 ABSTRACT Morbidly obese patients who accomplish substantial weight loss often display a long-term decline in their resting metabolism, causing even relatively restrained caloric intake to trigger a relapse to the obese state. Paradoxically, we observed that morbidly obese mice receiving chemotherapy for cancer experienced spontaneous weight reduction despite unabated ingestion of their high fat diet (HFD). This response to chemotherapy could also be achieved in morbidly obese mice without cancer. Optimally dosed methotrexate (MTX) or cyclophosphamide (CY) enabled the mice to completely and safely normalize their body weight despite continued consumption of obesogenic quantities of HFD. Weight reduction was not attributable to decreased HFD intake, enhanced energy expenditure or malabsorption. MTX or CY dosing significantly depleted both adipose tissue and preadipocyte progenitors. Remarkably, however, despite continued high fat feeding, a compensatory increase in hepatocyte lipid storage was not observed, but rather the opposite. Gene microarray liver analyses demonstrated that HFD mice receiving MTX or CY experienced significantly inhibited lipogenesis and lipid storage, whereas Enho (energy homeostasis) gene expression was significantly upregulated. Further metabolic studies employing a human hepatocellular line revealed that MTX treatment preserved robust oxidative phosphorylation, but also promoted mitochondrial uncoupling with a surge in proton leak. This is the first report that certain optimally dosed chemotherapeutic agents can induce weight loss in morbidly obese mice without reduced dietary intake, apparently by depleting stores of adipocytes and their progenitors, curtailment of lipogenesis, and inconspicuous disposal of incoming dietary lipid via a steady state partial uncoupling of mitochondrial oxidative phosphorylation.

Background Loss of tumour suppressor genes (TSGs), TP53, PTEN and Rb are common in advanced prostate cancer. They are associated with worse outcomes in hormone-resistant prostate cancer. However, the prognostic value of TSGs loss in metastatic hormone-sensitive prostate cancer (mHSPC) is unknown. Therefore we sought to evaluate the prognostic value of TSGs loss on progression-free survival (PFS) in mHSPC. Methods We performed a retrospective analysis of patients with mHSPC treated at Mayo Clinic Arizona between 2015 and 2018. Clinical characteristics, first-line treatment, and TSGs alterations were recorded. TSG-altered (TSG-alt) was defined as an alteration in ≥ 1 TSG (TP53, PTEN or Rb). Demographics and clinical characteristics were compared between groups by chi-square test for frequency data or Kruskal Wallis rank-sum test for continuous measures. PFS was estimated based on the Kaplan-Meier method and compared between groups by the log-rank test. Results 59 patients with mHSPC were included. 41 (69%) patients received first-line treatment with abiraterone and 18 (31%) patients with docetaxel. 32 patients (54%) had TSG-alt. 20 patients (39%) had damage to 1 TSG, 11 patients (19%) had damage to 2 TSGs, and 1 patient had damage to all 3. Clinical characteristics, prior locoregional treatment, and metastatic disease on diagnosis were similar between TSG-alt and TSG-normal groups. The rate of damaging alterations to the genes was: TP53 21 patients (35.6%), PTEN 21 patients (35.6%), and Rb 3 patients (5.1%). Median PFS was 33.2 mos (95% CI 23.1-NE) for TSG-normal versus 19.7 mos (95% CI 14.4-30.0) for TSG-alt, p = 0.005. There was no difference in PFS between having 1 TSG alteration (18.4 mos; 95% CI 14.7 - NE) vs ≥ 2 TSG alterations (21.8 mos; 95% CI 10.6 - NE) or between first-line treatments with abiraterone (20.2 mos; 95% CI 14.7-35.1) vs docetaxel (19.6 mos; 95% CI 8.0 - NE) in those with TSG-alt. Conclusions Loss of TSGs function is common in patients with mHSPC, and is prognostic for early progression during first-line treatment. The presence of ≥ 1 TSG-alt indicates an adverse prognostic impact in mHSPC. Confirmation of these findings may establish the need for inclusion of molecular stratification when defining high risk mHSPC patients. Legal entity responsible for the study Mayo Clinic Cancer Center. Funding Has not received any funding. Disclosure All authors have declared no conflicts of interest.

Motivation Functions of cancer driver genes vary substantially across tissues and organs. Distinguishing passenger genes, oncogenes (OGs) and tumor-suppressor genes (TSGs) for each cancer type is critical for understanding tumor biology and identifying clinically actionable targets. Although many computational tools are available to predict putative cancer driver genes, resources for context-aware classifications of OGs and TSGs are limited. Results We show that the direction and magnitude of somatic selection of protein-coding mutations are significantly different for passenger genes, OGs and TSGs. Based on these patterns, we develop a new method (genes under selection in tumors) to discover OGs and TSGs in a cancer-type specific manner. Genes under selection in tumors shows a high accuracy (92%) when evaluated via strict cross-validations. Its application to 10 172 tumor exomes found known and novel cancer drivers with high tissue-specificities. In 11 out of 13 OGs shared among multiple cancer types, we found functional domains selectively engaged in different cancers, suggesting differences in disease mechanisms. Availability and implementation An R implementation of the GUST algorithm is available at https://github.com/liliulab/gust. A database with pre-computed results is available at https://liliulab.shinyapps.io/gust. Supplementary information Supplementary data are available at Bioinformatics online.

Functions of cancer driver genes depend on cellular contexts that vary substantially across tissues and organs. Distinguishing oncogenes (OGs) and tumor suppressor genes (TSGs) for each cancer type is critical to identifying clinically actionable targets. However, current resources for context-aware classifications of cancer drivers are limited. In this study, we show that the direction and magnitude of somatic selection of missense and truncating mutations of a gene are suggestive of its contextual activities. By integrating these features with ratiometric and conservation measures, we developed a computational method to categorize OGs and TSGs using exome sequencing data. This new method, named genes under selection in tumors (GUST) shows an overall accuracy of 0.94 when tested on manually curated benchmarks. Application of GUST to 10,172 tumor exomes of 33 cancer types identified 98 OGs and 179 TSGs, >70% of which promote tumorigenesis in only one cancer type. In broad-spectrum drivers shared across multiple cancer types, we found heterogeneous mutational hotspots modifying distinct functional domains, implicating the synchrony of convergent and divergent disease mechanisms. We further discovered two novel OGs and 28 novel TSGs with high confidence. The GUST program is available at https://github.com/liliulab/gust. A database with pre-computed classifications is available at https://liliulab.shinyapps.io/gust

Purpose: Primary central nervous system lymphoma (PCNSL) is an aggressive non-Hodgkin lymphoma confined to the central nervous system. Whether there is a PCNSL-specific genomic signature and, if so, how it differs from systemic diffuse large B-cell lymphoma (DLBCL) is uncertain. Experimental Design: We performed a comprehensive genomic study of tumor samples from 19 immunocompetent PCNSL patients. Testing comprised array-comparative genomic hybridization and whole exome sequencing. Results: Biallelic inactivation of TOX and PRKCD was recurrently found in PCNSL but not in systemic DLBCL, suggesting a specific role in PCNSL pathogenesis. In addition, we found a high prevalence of MYD88 mutations (79%) and CDKN2A biallelic loss (60%). Several genes recurrently affected in PCNSL were common with systemic DLBCL, including loss of TNFAIP3, PRDM1, GNA13, TMEM30A, TBL1XR1, B2M, CD58, activating mutations of CD79B, CARD11, and translocations IgH-BCL6. Overall, B-cell receptor/Toll-like receptor/NF-κB pathways were altered in >90% of PNCSL, highlighting its value for targeted therapeutic approaches. Furthermore, integrated analysis showed enrichment of pathways associated with immune response, proliferation, apoptosis, and lymphocyte differentiation. Conclusions: In summary, genome-wide analysis uncovered novel recurrent alterations, including TOX and PRKCD, helping to differentiate PCNSL from systemic DLBCL and related lymphomas. Clin Cancer Res; 21(17); 3986–94. ©2015 AACR.

Recent advances in genomic sequencing technologies now allow results from deep next-generation sequencing to be obtained within clinically meaningful timeframes, making this an attractive approach to better guide personalized treatment strategies. No multiple myeloma-specific gene panel has been established so far; we therefore designed a 47-gene-targeting gene panel, containing 39 genes known to be mutated in ≥3 % of multiple myeloma cases and eight genes in pathways therapeutically targeted in multiple myeloma (MM). We performed targeted sequencing on tumor/germline DNA of 25 MM patients in which we also had a sequential sample post treatment. Mutation analysis revealed KRAS as the most commonly mutated gene (36 % in each time point), followed by NRAS (20 and 16 %), TP53 (16 and 16 %), DIS3 (16 and 16 %), FAM46C (12 and 16 %), and SP140 (12 and 12 %). We successfully tracked clonal evolution and identified mutation acquisition and/or loss in FAM46C, FAT1, KRAS, NRAS, SPEN, PRDM1, NEB, and TP53 as well as two mutations in XBP1, a gene associated with bortezomib resistance. Thus, we present the first longitudinal analysis of a MM-specific targeted sequencing gene panel that can be used for individual tumor characterization and for tracking clonal evolution over time.

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.

Purpose Genomic testing has increased the quantity of information available to oncologists. Unfortunately, many identified sequence alterations are variants of unknown significance (VUSs), which thus limit the clinician’s ability to use these findings to inform treatment. We applied a combination of in silico prediction and molecular modeling tools and laboratory techniques to rapidly define actionable VUSs. Materials and Methods Exome sequencing was conducted on 308 tumors from various origins. Most single nucleotide alterations within gene coding regions were VUSs. These VUSs were filtered to identify a subset of therapeutically targetable genes that were predicted with in silico tools to be altered in function by their variant sequence. A subset of receptor tyrosine kinase VUSs was characterized by laboratory comparison of each VUS versus its wild-type counterpart in terms of expression and signaling activity. Results The study identified 4,327 point mutations of which 3,833 were VUSs. Filtering for mutations in genes that were therapeutically targetable and predicted to affect protein function reduced these to 522 VUSs of interest, including a large number of kinases. Ten receptor tyrosine kinase VUSs were selected to explore in the laboratory. Of these, seven were found to be functionally altered. Three VUSs (FGFR2 F276C, FGFR4 R78H, and KDR G539R) showed increased basal or ligand-stimulated ERK phosphorylation compared with their wild-type counterparts, which suggests that they support transformation. Treatment of a patient who carried FGFR2 F276C with an FGFR inhibitor resulted in significant and sustained tumor response with clinical benefit. Conclusion The findings demonstrate the feasibility of rapid identification of the biologic relevance of somatic mutations, which thus advances clinicians’ ability to make informed treatment decisions.

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.

Purpose: Cholangiocarcinoma (CCA) is a highly lethal malignancy that often presents at advanced stages. Gemcitabine/cisplatin (GP) represents the standard-of-care but provides only a modest response rate and survival benefit. Loss of the CDK4/6 inhibitor CDKN2A is a prevalent genetic event in CCA, thus CDK4/6 inhibitors, FDA-approved in some breast cancer regimens, are hypothesized to have a role in treating CCA. Clinical challenges to target cell cycle dysregulation and exploit cytotoxic chemotherapy include predictive biomarkers for selecting patients most likely to respond. We hypothesized autophagy is involved in therapeutic resistance to CDK4/6 inhibitors, and plays a role in the activity of GP in CCA. Thus, we investigated CDK4/6 inhibitors as monotherapy, and combined with GP in CCA, measuring multiple readouts pertaining to drug synergy, putative pathway biomarkers and autophagy. Methods: The activity of 3 FDA approved CDK4/6 inhibitors, palbociclib, ribociclib and abemaciclib, as monotherapy (N=24), and combined with GP (N=10), was assessed in CCA cell lines by measuring relative cell number with Cell TiterGlo, and apoptosis with Annexin/PI and cleaved caspase3 quantification via flow cytometry. To characterize RB status for biomarker potential, protein levels of CCND1, CDK4, CDK6, RB, and CDKN2A were measured by western blot. Autophagy was measured using an autophagy flux assay, LysoTracker confocal staining, and by quantifying protein levels of LC3B, p62 and Beclin1A. Results: 22/24 (92%) of CCA cell lines were found to have at least one defect likely relevant to CDK4/6 inhibition: CDKN2A (92%), Rb1 (63%), CDK4 (23%) and cyclin D1 (32%). Protein levels of CDKN2A and RB1 showed moderate correlation with CDK4/6 inhibitor activity. Abemaciclib was the most potent CDK inhibitor in vitro (IC50≤0.5uM; palbociclib IC50>1uM; ribociclib IC50>5uM). Combined abemaciclib and GP showed higher cytotoxicity and synergy at lower doses compared to chemotherapy alone or single-agent abemaciclib. Compared to CDK4/6 inhibitors or chemotherapy alone, combination therapy enhanced growth inhibition, cell cycle arrest and caspase-dependent apoptotic cell death. Abemaciclib triggered autophagy mediated by significant increases in reactive oxygen species (ROS). Combination therapy reduced ROS production and thus autophagy, as seen by reduced autophagic flux, reduced LysoTracker mean fluorescence intensity, a decrease in the expression of LC3B-II, Beclin1A and an increase in p62. Conclusion: These results suggest CDK4/6 inhibitors interact synergistically with GP in CCA and mechanistically disrupt autophagy. This suggests potential to further explore CDKN2A and RB as biomarkers for preclinical sensitivity in CCA. Pre-clinical data supports ongoing clinical studies of this concept, as well as planned clinical studies of GP and alternative CDK4/6 inhibitors. Citation Format: Mansi Arora, James Bogenberger, Nagalo Bolni, Yumei Zhou, Latha Kannan, Jan Egan, Raquel Yokoda, Daniel H. Ahn, Mitesh J. Borad. Synergy of cyclin dependentkinase 4/6 inhibitors with cytotoxic chemotherapy in cholangiocarcinoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 4417.

Extramedullary disease (EMD) in Multiple Myeloma (MM) is characterized by the detection of monoclonal plasma cells outside the bone marrow niche, and is frequently associated with poor prognosis. Here we describe novel genomic events leading to drug refractory disease in a heavily pretreated 37-year-old IgG-kappa MM patient presenting with progressive, multi-drug refractory EMD. For the first time we report an acquired truncating mutation of Cereblon (CRBN) as well as point mutations in proteasome subunit G2 and the glucocorticoid receptor as an explanation for drug resistance. Initial myeloma treatment for the patient occurred over multiple years and included the immunomodulatory drugs (IMiDs) thalidomide and lenalidomide, the proteasome inhibitor bortezomib, cortiosteroids, radiation, one autologous and two allogeneic transplantations. She experienced extramedullary relapse, presenting as an extensive neck mass and smaller soft tissue nodules in the upper left triceps. The most recent therapy immediately prior to genomic sequencing was hyper-CVAD (hyperfractionated cyclophosphamide, vincristine, doxorubicin and dexamethasone) incorporating alkylating agent cyclophosphamide with transient minor response. The patient was then enrolled in a pilot study utilizing next generation sequencing (NGS) to identify novel markers and potential therapeutic targets. Samples were acquired with patient consent in compliance with Mayo Clinic Institutional Review Board. For this study we completed array comparative genomic hybridization, whole exome, whole genome (insert = 1 kb) and RNA sequencing (RNASeq) of a biopsy taken from the neck mass to thoroughly interrogate the tumour genome of this patient. The presence of mutations of interest was evaluated by capillary sequencing in an expanded cohort of 25 CD138+ MM samples with low CRBN expression. The neck mass pathology confirmed sheets of atypical plasma cells, kappa light chain restriction, CD138+, CD20− and CD45−. Array comparative genomic hybridization revealed multiple copy number abnormalities, most notably del1(p13.2–34.2), monosomy 13 and monosomy X. Therapy subsequent to biopsy for genome sequencing was with pomalidomide and dexamethasone without response. Unfortunately, the patient succumbed to her disease in less than the 12 weeks required at the time for sequencing and data analysis. Sequencing revealed a highly disturbed genome (Figure 1) consisting of 4 somatic insertions/deletions, 38 intra-chromosomal rearrangements, and 35 translocations, including the high risk marker and initiating tumour event t(14;16). Furthermore, 271 nonsynonymous, somatic point mutations were detected in genes including KRAS, PIK3CA, ATM, and NFKB2 (Table I). Importantly, a Q99* truncating mutation as well as a R283K point mutation were observed in CRBN, that we recently demonstrated as essential for the anti-MM action of IMiDs (Zhu, et al 2011). To our knowledge this is the first documented mutation of Cereblon in a primary myeloma sample. Additional sequencing of CRBN in the expanded cohort of 25 patients revealed a synonymous mutation in only one sample. Figure 1 Circos plot depicting genome wide somatic variants, rearrangements and copy number changes derived from next generation sequencing. Numbers with circles around them indicate the following: 1) somatic single nucleotide variation (SNV), 2) location of SNV ... Table I Summary of clinically relevant single nucleotide variations We also observed in the patient biopsy a potentially clinically relevant nonsynonymous point mutation in proteasome assembly chaperone 2, PSMG2 (E171K). PSMG2 is a proteasome assembly protein involved in mammalian 20S proteasome maturation (Hirano, et al 2005). Mutations in proteasome assembly components contribute to proteasome inhibitor resistance (Keats, et al 2007), possibly explaining this patient’s bortezomib-refractory disease. Capillary sequencing of PSMG2 in our expanded cohort revealed no mutations, although exonic deletion of PSMG2 has also been reported in myeloma (Walker, et al 2012). The last nonsynonymous point mutation associated with drug resistance was identified in NR3C1 (G369A), a glucocorticoid receptor. Mutation of NR3C1 has been associated with resistance to steroid therapy (Bray and Cotton 2003), which this patient received and proved refractory. No NR3C1 mutations were identified in our expanded cohort and none have been previously reported in other myeloma genomes (Chapman, et al 2011, Walker, et al 2012). Mutations in NR3C1 have however been described in the glucocorticoid resistant MM.1R cell line (Moalli, et al 1992). Patients with low NR3C1 expression levels who received thalidomide demonstrated better progression-free survival and overall survival than those with low NR3C1 who did not receive thalidomide (Heuck, et al 2012). While these mutations suggest causality of drug-refractory disease, they do not identify pathways that can be exploited with targeted therapies. Additional mutations were observed in pathways for which targeted therapies are available. This patient had mutations in KRAS (G12C) and in ATM (T1985I), both of which affect the signalling of MEK downstream, thus making MEK a therapeutic target of interest in this patient. While there are no approved MEK inhibitors available for MM treatment, more than 100 trials are currently investigating MEK inhibitors, of which three of these trials are being conducted in MM patients (www.clinicaltrials.gov). The patient tumour also contained canonical, activating mutations in PIK3CA (E542K). Interestingly, one study demonstrated that 64% of PIK3CA mutations occur in exon 9, where codon 542 is located. Moreover, 19% of patients with PIK3CA mutations also presented with KRAS mutations, of which 9% are G12C (Janku, et al 2012), found in our patient. The PI3K pathway is vitally important as it regulates downstream targets, such as AKT and MTOR, which are responsible for cell proliferation, growth, survival and metastasis (Bartholomeusz and Gonzalez-Angulo 2012). In addition, a number of clinical trials are currently investigating PIK3 inhibitors (www.clinicaltrials.gov). In summary, this is the first description of CRBN mutations in a primary myeloma sample and furthermore of a “triple negative” MM patient possessing mutations probably contributing to resistance to all three major drug classes utilized in MM therapy. These mutations were not replicated in our validation cohort of 25 patients with low level CRBN expression and functional data have not yet been obtained, thus further investigation is necessary to better understand the mutation frequency and the functional significance of mutation in these genes. In summary, our approach utilizing comprehensive next generation sequencing not only identified mutations suggestive of the patient’s refractory disease, but also revealed unforeseen therapeutic options highlighting the importance of this technology in advancing individualized medicine.

The longitudinal evolution of a myeloma genome from diagnosis to plasma cell leukemia has not previously been reported. We used whole-genome sequencing (WGS) on 4 purified tumor samples and patient germline DNA drawn over a 5-year period in a t(4;14) multiple myeloma patient. Tumor samples were acquired at diagnosis, first relapse, second relapse, and end-stage secondary plasma cell leukemia (sPCL). In addition to the t(4;14), all tumor time points also shared 10 common single-nucleotide variants (SNVs) on WGS comprising shared initiating events. Interestingly, we observed genomic sequence variants that waxed and waned with time in progressive tumors, suggesting the presence of multiple independent, yet related, clones at diagnosis that rose and fell in dominance. Five newly acquired SNVs, including truncating mutations of RB1 and ZKSCAN3, were observed only in the final sPCL sample suggesting leukemic transformation events. This longitudinal WGS characterization of the natural history of a high-risk myeloma patient demonstrated tumor heterogeneity at diagnosis with shifting dominance of tumor clones over time and has also identified potential mutations contributing to myelomagenesis as well as transformation from myeloma to overt extramedullary disease such as sPCL.

The vitamin D metabolite 1,25(OH) 2 D is the bioactive ligand of the vitamin D receptor (VDR). VDR forms a heterodimer with the retinoid X receptors (RXRs) that when bound to ligand influences the transcriptional control of genes that regulate circulating levels of vitamin D metabolites. Whether genetic variation in VDR or RXRA affects circulating levels of 1,25(OH) 2 D or 25(OH)D has not been established. We used a single nucleotide polymorphism (SNP) tagging approach to evaluate the association between SNPs in VDR and RXRA and serum levels of 1,25(OH) 2 D and 25(OH)D. A total of 42 tagSNPs in VDR and 32 in RXRA were analyzed in a sample of 415 participants. Principal components analyses revealed a gene-level association between RXRA and serum 1,25(OH) 2 D concentrations ( P = 0.01), but not 25(OH)D. No gene-level association was found for VDR with either serum biomarker. At the single-SNP level, a significant positive trend was observed for increasing 1,25(OH) 2 D levels with each additional copy of the A allele for RXRA SNP rs9409929 ( P -trend = 0.003). After a multiple comparisons adjustment, no individual SNP in VDR or RXRA was significantly associated with either outcome. These results demonstrate an association between genetic variation in RXRA and 1,25(OH) 2 D serum concentrations.

293 Background: Cholangiocarcinoma (CCA) is a deadly cancer with ~20% five year survival. Obtaining tissue to confirm diagnosis of CCA is often challenging. This also limits the use of tumor genotyping that could guide potential molecularly targeted treatment (such as FGFR inhibitors for a subset of cases). We investigated whether circulating tumor DNA (ctDNA) analysis in plasma is feasible for noninvasive tumor genotyping in patients with advanced CCA. Methods: We collected and processed plasma samples for ctDNA analysis from 17 patients with CCA at presentation and from 11 healthy volunteers. Following cell-free DNA extraction and quality assessment using digital PCR, we performed low-pass, shallow whole genome sequencing (sWGS). We analyzed sequencing data to identify copy number aberrations and quantify tumor fraction in plasma using recently published computational approaches. Results: Median total cfDNA concentration in patients with CCA was 7,003 genome equivalents (GEs)/mL plasma (range: 1,512- 107,707 GEs/mL) and significantly higher than healthy volunteers (two tailed t-test p < 0.04, median: 1,658 GEs/mL, range: 661-4,938 GEs/mL). We successfully sequenced all tumor libraries, generating a median of 173 million paired-end reads. ctDNA was detectable in 10/17 samples using sWGS. In these 10 samples, median tumor fraction in plasma DNA was 7.1% (range: 3.9%-44.6%). We found recurrent copy number aberrations expected in CCA including 1q gain in 80%, 7p gain in 50%, 1p loss in 70% and 3p loss in 50% of cases. Conclusions: ctDNA analysis for noninvasive tumor genotyping is feasible in patients with advanced cholangiocarcinoma. On-going work is focused on expanding these results using targeted and whole-exome sequencing.

Low circulating levels of vitamin D affect colorectal cancer risk. The biological actions of the hormonal form of vitamin D, 1,25(OH)2D3, are mediated by the vitamin D receptor (VDR), which heterodimerizes with retinoid X receptors (RXR). Using a single nucleotide polymorphism (SNP) tagging approach, we assessed the association between genetic variations in RXRA and VDR and odds of recurrent (metachronous) colorectal neoplasia in a pooled population of two studies. A total of 32 tag SNPs in RXRA and 42 in VDR were analyzed in 1,439 participants. A gene-level association was observed for RXRA and any (P = 0.04) or proximal (P = 0.03) metachronous neoplasia. No gene-level associations were observed for VDR, nor was any single SNP in VDR related to any metachronous adenoma after correction for multiple comparisons. In contrast, the association between RXRA SNP rs7861779 and proximal metachronous neoplasia was of borderline statistical significance [odds ratio (OR), 0.68; 95% confidence interval (95% CI), 0.53-0.86; unadjusted P = 0.001; adjusted P = 0.06], including when observed independently in each individual study. Haplotypes within linkage blocks of RXRA support an ∼30% reduction in odds of metachronous neoplasia arising in the proximal colon among carriers of specific haplotypes, which was strongest (ORproximal, 0.67; 95% CI, 0.52–0.86) for carriers of a CGGGCA haplotype (rs1805352, rs3132297, rs3132296, rs3118529, rs3118536, and rs7861779). Our results indicate that allelic variation in RXRA affects metachronous colorectal neoplasia, perhaps of particular importance in the development of proximal lesions. Cancer Res; 70(4); 1496–504

The activity of beta-catenin, commonly dysregulated in human colon cancers, is inhibited by the vitamin D receptor (VDR), and this mechanism is postulated to explain the putative anti-cancer activity of vitamin D metabolites in the colon. We investigated the effect of a common FokI restriction site polymorphism (F/f) in the human VDR gene as well as the effect of anti-tumorigenic 1,25-dihydroxyvitamin D(3) (1,25D) and pro-tumorigenic lithocholic acid (LCA) VDR ligands on beta-catenin transcriptional activity. Furthermore, the influence of a major regulatory protein of beta-catenin, the APC tumor suppressor gene, on VDR-dependent inhibition of beta-catenin activity was examined. We report herein that beta-catenin-mediated transcription is most effectively suppressed by the VDR FokI variant F/M4 when 1,25D is limiting. Using Caco-2 colorectal cancer (CRC) cells, it was observed that VDR ligands, 1,25D and LCA, both suppress beta-catenin transcriptional activity, though 1,25D exhibited significantly greater inhibition. Moreover, 1,25D, but not LCA, suppressed endogenous expression of the beta-catenin target gene DKK-4 independent of VDR DNA-binding activity. These results support beta-catenin sequestration away from endogenous gene targets by 1,25D-VDR. This activity is most efficiently mediated by the FokI gene variant F/M4, a VDR allele previously associated with protection against CRC. Interestingly, we found the inhibition of beta-catenin activity by 1,25D-VDR was significantly enhanced by wild-type APC. These results reveal a previously unrecognized role for 1,25D-VDR in APC/beta-catenin cross talk. Collectively, these findings strengthen evidence favoring a direct effect on the Wnt-signaling molecule beta-catenin as one anti-cancer target of 1,25D-VDR action in the colorectum.

Objective: Bringing together genomics, transcriptomics, proteomics, and other –omics technologies is an important step towards developing highly personalized medicine. However, instrumentation has advances far beyond expectations and now we are able to generate data faster than it can be interpreted. Materials and Methods: We have developed PANDA ( P athway AND A nnotation) Explorer, a visualization tool that integrates gene-level annotation in the context of biological pathways to help interpret complex data from disparate sources. PANDA is a web-based application that displays data in the context of well-studied pathways like KEGG, BioCarta, and PharmGKB. PANDA represents data/annotations as icons in the graph while maintaining the other data elements (i.e. other columns for the table of annotations). Custom pathways from underrepresented diseases can be imported when existing data sources are inadequate. PANDA also allows sharing annotations among collaborators. Results : In our first use case, we show how easy it is to view supplemental data from a manuscript in the context of a user’s own data. Another use-case is provided describing how PANDA was leveraged to design a treatment strategy from the somatic variants found in the tumor of a patient with metastatic sarcomatoid renal cell carcinoma. Conclusion : PANDA facilitates the interpretation of gene-centric annotations by visually integrating this information with context of biological pathways. The application can be downloaded or used directly from our website http://bioinformaticstools.mayo.edu/research/panda-viewer/.

11571 Background: Relapsed/refractory (R/R) non-Hodgkin lymphomas (NHL) have a poor prognosis with limited treatment options. Our expanding knowledge of molecular alterations seen in R/R NHL allows identification of patients that potentially may benefit from a precision medicine approach. However, experience in routine clinical implementation of precision medicine has been limited. Here, we summarize our clinical experience in molecular characterization of RR NHL targeted therapy (TT) using next-generation sequencing (NGS), and selection of targeted therapy (TT) based on molecular profile. Methods: We conducted a prospective study in RR NHL through the Center for Individualized Medicine at Mayo Clinic. Consenting patients underwent NGS using FoundationOne Heme panel from biopsies done at time of relapse. Results of NGS were discussed at the Genomic Tumor Board and recommendations for TT were given based on matching specific molecular alteration(s) with potential agent(s) predicted to be active based on NGS. The agents could include FDA-approved, off-label use and clinical trial therapies. Results: 28 cases were enrolled: 18 aggressive NHL, 10 follicular lymphoma (FL). Molecular alterations were present in all cases. In aggressive B-cell NHL, CDKN2A/B gene cluster alterations were seen in 73% (8/11), while seen in only 1/7 T-cell lymphomas (TCL), and 1/10 FL. TP53 deletions were second most common genomic alterations in DLBCL (57%) and seen in 40% FL. JAK-STAT and ERBB pathways were altered in TCL (2/7 each). IGH-BCl-2 gene rearrangement were common in FL (70%), followed by MLL gene alterations (50%). Targetable mutations were present in 86% (24/28) of cases. A TT was recommended in all 24 cases, but received by 2 patients only. Remaining patients did not due to benefit from current therapy (10/24), ineligibility or lack of clinical trial (7/24) or interim clinical deterioration (5/24). Conclusions: Targetable mutations were identified in most cases of RR NHL with TT recommended for all cases. However, access to TT limits potential clinical benefit of molecular-based matching strategy. More studies are needed to assess impact on clinical outcomes.

BACKGROUND: Primary central nervous system lymphoma (PCNSL) is an aggressive non-Hodgkin lymphoma (NHL) confined to the CNS. Whether there is a PCNSL-specific genomic signature and, if so, how it differs from systemic DLBCL is uncertain. To address this gap in critical knowledge we performed a comprehensive genomic study in a cohort of 19 immunocompetent PCNSL patients by combining aCGH, whole exome sequencing and mate-pair whole genome sequencing. METHODS: 1) Tumor samples from 19 immunocompetent (HIV-, EBV-) PCNSL. 2) DNA from frozen tissue or from formalin-fixed paraffin embedded (FFPE) tissues. 3) All samples analyzed using a combination of aCGH, whole exome sequencing, mate-pair whole genome sequencing, targeted sequencing and FISH. 4) Overall survival estimated using Kaplan Meier method. RESULTS: We found a complex karyotype with a median of 21 copy-number abnormalities, 6 structural abnormalities and 3.3 non-synonymous mutations per Mb. We found a remarkably high prevalence of MYD88 activating mutations (79%) and CDKN2A biallelic loss (60%). Biallelic losses of TOX (a regulator of T-cell development), PRKCD (protein kinase C delta) and chromosomal breakpoints in DLGAP1 (discs large-associated protein 1) were recurrently found in PCNSL but not in systemic DLBCL, suggesting a specific role in PCNSL pathogenesis. Several other genes recurrently affected in PCNSL are common with systemic DLBCL of the ABC-type, such as loss of TNFAIP3, PRDM1, GNA13, TMEM30A, TBL1XR1, B2M and CD58, activating mutations of CD79B, CARD11, BRAF3 and translocations of IgH-BCL6. Overall, components of the NF-□B signaling pathway were altered in >90% of PNCSL, highlighting the potential value of this pathway for targeted therapeutic approaches. Furthermore, integrated analysis also showed an enrichment of networks associated with immune response, proliferation, regulation of apoptosis, and lymphocyte differentiation. CONCLUSIONS: We report a PCNSL genomic landscape closely related to ABC-DLBCL but with critical elements that support the existence of a subset of PCNSL-specific abnormalities, helping to genetically differentiate PCNSL from systemic DLBCL and related lymphomas. SECONDARY CATEGORY: Immunobiology & Immunotherapy.

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.

Supplemental digital content is available in the text., Objective Integration of carcinogenic human papillomaviruses (HPVs) into the host genome is a significant tumorigenic factor in specific cancers including cervical carcinoma. Although major strides have been made with respect to HPV diagnosis and prevention, identification and development of efficacious treatments for cervical cancer patients remains a goal and thus requires additional detailed characterization of both somatic events and HPV integration. Given this need, the goal of this study was to use the next generation sequencing to simultaneously evaluate somatic alterations and expression changes in a patient’s cervical squamous carcinoma lesion metastatic to the lung and to detect and analyze HPV infection in the same sample. Materials and Methods We performed tumor and normal exome, tumor and normal shallow whole-genome sequencing, and RNA sequencing of the patient’s lung metastasis. Results We generated over 1.2 billion mapped reads and identified 130 somatic point mutations and indels, 21 genic translocations, 16 coding regions demonstrating copy number changes, and over 36 genes demonstrating altered expression in the tumor (corrected P < 0.05). Sequencing also revealed the HPV type 18 (HPV-18) integration in the metastasis. Using both DNA and RNA reads, we pinpointed 3 major events indicating HPV-18 integration into an intronic region of chromosome 6p25.1 in the patient’s tumor and validated these events with Sanger sequencing. This integration site has not been reported for HPV-18. Conclusions We demonstrate that DNA and RNA sequencing can be used to concurrently characterize somatic alterations and expression changes in a biopsy and delineate HPV integration at base resolution in cervical cancer. Further sequencing will allow us to better understand the molecular basis of cervical cancer pathogenesis.

RNA interference screening identified XPO1 (exportin 1) among the 55 most vulnerable targets in multiple myeloma (MM). XPO1 encodes CRM1, a nuclear export protein. XPO1 expression increases with MM disease progression. Patients with MM have a higher expression of XPO1 compared with normal plasma cells (P

Next-generation sequencing has led to a revolution in the study of hematological malignancies with a substantial number of publications and discoveries in the last few years. Significant discoveries associated with disease diagnosis, risk stratification, clonal evolution and therapeutic intervention have been generated by this powerful technology. As part of the post-genomic era, sequencing analysis will likely become part of routine clinical testing and the challenge will ultimately be successfully transitioning from gene discovery to preventive and therapeutic intervention as part of individualized medicine strategies. In this report, we review recent advances in the understanding of hematological malignancies derived through genome-wide sequence analysis.

Pancreatic adenocarcinoma (PAC) is among the most lethal malignancies. While research has implicated multiple genes in disease pathogenesis, identification of therapeutic leads has been difficult and the majority of currently available therapies provide only marginal benefit. To address this issue, our goal was to genomically characterize individual PAC patients to understand the range of aberrations that are occurring in each tumor. Because our understanding of PAC tumorigenesis is limited, evaluation of separate cases may reveal aberrations, that are less common but may provide relevant information on the disease, or that may represent viable therapeutic targets for the patient. We used next generation sequencing to assess global somatic events across 3 PAC patients to characterize each patient and to identify potential targets. This study is the first to report whole genome sequencing (WGS) findings in paired tumor/normal samples collected from 3 separate PAC patients. We generated on average 132 billion mappable bases across all patients using WGS, and identified 142 somatic coding events including point mutations, insertion/deletions, and chromosomal copy number variants. We did not identify any significant somatic translocation events. We also performed RNA sequencing on 2 of these patients' tumors for which tumor RNA was available to evaluate expression changes that may be associated with somatic events, and generated over 100 million mapped reads for each patient. We further performed pathway analysis of all sequencing data to identify processes that may be the most heavily impacted from somatic and expression alterations. As expected, the KRAS signaling pathway was the most heavily impacted pathway (P

Emerging evidence indicates that tumors can follow several evolutionary paths over a patient's disease course. With the use of serial genomic analysis of samples collected at different points during the disease course of 28 patients with multiple myeloma, we found that the genomes of standard-risk patients show few changes over time, whereas those of cytogenetically high-risk patients show significantly more changes over time. The results indicate the existence of 3 temporal tumor types, which can either be genetically stable, linearly evolving, or heterogeneous clonal mixtures with shifting predominant clones. A detailed analysis of one high-risk patient sampled at 7 time points over the entire disease course identified 2 competing subclones that alternate in a back and forth manner for dominance with therapy until one clone underwent a dramatic linear evolution. With the use of the Vk*MYC genetically engineered mouse model of myeloma we modeled this competition between subclones for predominance occurring spontaneously and with therapeutic selection.

Acquired or inherited mutations in the adenomatous polyposis coli (APC) tumor suppressor gene are causally linked to colorectal cancer. Given the significance of APC in colorectal cancer, we investigated the association between common single-nucleotide polymorphisms (SNP) in the APC gene and the odds of developing metachronous colorectal adenomas as a surrogate measure of colorectal cancer risk. Coding SNPs at codons 486, 1678, 1822, 1960, and 2502 were analyzed in a total of 1,399 subjects who participated in two randomized clinical trials for the prevention of colorectal adenomas. No association was found for any single SNP and the odds of metachronous adenoma. In contrast, a TA haplotype (codons 486 and 1822) was associated with a statistically significant 27% and 26% reduction in the odds of any and nonadvanced metachronous adenoma after adjustment for baseline adenoma characteristics [odds ratio (OR), 0.73; 95% confidence interval (95% CI), 0.59–0.91 and OR, 0.74; 95% CI, 0.57–0.94], respectively. No significant reduction in odds was observed for advanced metachronous lesions. Diplotype analysis revealed a strong gene dose effect with carriers of two alleles containing TT-AA (codons 486 and 1822, respectively) having an 89% lower odds for advanced metachronous adenomas (OR, 0.11; 95% CI, 0.01–0.80) when compared with the common CC-AA diplotype (codons 486 and 1822, respectively). Our findings support an important role for germ-line allele sequence in the APC gene and individual risk of metachronous adenomatous polyps. [Cancer Res 2008;68(14):6006–13]

Background: Genomic assessment of cancer has been revolutionized by Next-Generation sequencing and is increasingly being applied to guide clinical prognostic and therapeutic decision-making. Initial clinical applications have been limited to gene panels and whole exome based strategies due to challenges with time to reporting of results, specimen quantity, analyte quality, and ethical-legal-social issues (ELSI). Methods: Excisional or core tumor biopsies or bone marrow biopsies were obtained from consenting participants. Nucleic acid was extracted from fresh or frozen samples followed by WGS/WES/RNASeq on the Illumina HiSeq2000 or HiSeq2500 and bioinformatics analysis. Therapeutic targets were then prioritized by a multi-disciplinary Genomic Tumor Board (GTB) and CLIA validated using Sanger sequencing, RT-qPCR, FISH and/or IHC. Treatment was delivered using on/off-label FDA approved drugs, available clinical trials and single patient INDs. Results: We consented 64 and enrolled 35 patients with advanced, treatment-refractory cancers. Median age was 59 (range 27-91) with 62% male. The majority had ECOG scores of 1 (94%). A median of 79 (range 28-8891) potentially functional somatic point mutations were identified in each case. One to two mutations/case were further identified as therapeutically targetable in 60% of the cases. The median time from tissue acquisition to CLIA validated results was 116 days (range 42-282) with CLIA validation of targets achieved in 21 of 22 patients. Genomic, target directed treatment was ultimately instituted in 13 patients utilizing: on/off label FDA approved drugs (n = 9), clinical trial (n = 3) and single patient IND (n = 1). Preliminary clinical efficacy was noted in 5 patients (2 PR, 3 SD). Integration of WES, long-insert WGS and RNA-Seq identified a case with an ERRFI1 mutant allele present in only 11% of the DNA reads while 82% of the RNA-Seq reads had the mutant transcript. The patient was treated with erlotinib and achieved a partial response by RECIST criteria. In another case, a fusion between FGFR2-MGEA5 was observed in WES and RNA-Seq data leading to prioritization as a drug target. Reasons for patients not receiving targeted therapy included: inability to access treatment (n = 1), death prior to intended treatment (n = 3), results returned after death (n = 3), no targets identified (n = 2) and decision to pursue alternate therapy (n = 5). Conclusions: Integrating whole genome analysis in a CLIA setting is not only feasible, but also valuable, in the prioritization and selection of potential targeted therapies for patients with advanced tumors. Continued barriers to broad application include the need for shorter time to reporting as well as broad availability of therapies through basket studies. Citation Format: Mitesh J. Borad, Jan Egan, Mia Champion, Katherine Hunt, Robert McWilliams, Ann McCullough, Jessica Aldrich, Sara Nasser, Winnie Liang, Michael Barrett, David Craig, Ramesh Ramanathan, John Carpten, A. Keith Stewart, Alan Bryce. Implementation of CLIA enabled integrated whole genome (WGS)/exome (WES)/transcriptome (RNAseq) next-gen sequencing to identify therapeutically relevant targets in advanced cancer patients. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr CT112. doi:10.1158/1538-7445.AM2015-CT112

Background: Whole genome/exome/RNA sequencing has revolutionized the ability to assess the genomic landscape of cancer and is increasingly being utilized for clinical decision-making. Initial clinical applications have been constrained by specimen quantity, analyte quality and the time from sample acquisition to results report. Methods: Patients with advanced cancers underwent surgical resection, excisional/core biopsies, or bone marrow biopsy. Samples were analyzed by whole genome or exome sequencing in addition to RNA sequencing, bioinformatics analysis, and therapeutic target prioritization by a multi-disciplinary Clinical Genomics Board. All prioritized targets were CLIA validated using Sanger sequencing, RT-qPCR, FISH, or IHC as appropriate. Treatment was delivered using off-label FDA approved drugs, clinical trials, or single patient INDs. Results: We have enrolled 40 patients for whom sequencing data is available on 33. The initial 6 patients were evaluated in a non-CLIA pilot phase and 27 in a CLIA-enabled phase. Tumor types in the CLIA-enabled phase with the highest enrollment were pancreatic cancer (n=8) and cholangiocarcinoma (n=8). We sought to quantify the targets identified along with clinical benefit, defining these as the “Actionable Index” (AI) (proportion of patients with ≥ 1 putative drug target) and “Utility Index” (UI) (proportion of patients who derive clinical benefit). Putative therapeutic targets were identified in 7/8 (AI=0.88) cholangiocarcinoma (CC) patients and in 5/8 (AI=0.63) pancreatic cancer (PC) patients. All 3 CC patients who received target directed treatment achieved a partial response (UI=0.38). In contrast, none of the 4 PC patients who received target directed therapy had treatment response (UI=0.0). Interestingly no actionable targets were identified in 1 CC and in 2 PCs. One CC with an identified target was unable to access the drug and subsequently died. A CC patient and a PC patient, each with identified targets, expired prior to the initiation of therapy. Conclusions: While whole genome/exome/RNA sequencing is providing unparalleled detail of tumor genomes, the application to the clinic must be carefully considered. Actionability of targets will eventually need to be defined in close relation to eventual clinical utility and appropriate refinements to disease-gene-drug databases implemented. Preliminary observations in pancreatic cancer and cholangiocarcinoma demonstrate disparity in correlation between utility indices and actionable indices. Application of these tools in larger cohorts and types of tumors will need to be conducted to ascertain more precise estimates. Additional measures that are organ-site agnostic but pertain to specific targets (e.g. BRAF) will also need to be developed in order to facilitate more judicious application of sequencing in the clinical setting. Citation Format: Jan B. Egan, Alan H. Bryce, Mia D. Champion, Winnie S. Liang, Rafael Fonseca, Ann E. McCullough, Michael T. Barrett, Katherine Hunt, Rachel M. Condjella, Robert R. McWilliams, Stephen D. Mastrian, Janine LoBello, Daniel Von Hoff, David W. Craig, A. Keith Stewart, John D. Carpten, Mitesh J. Borad. Indices of actionability and clinical utility in a CLIA-enabled study of whole genome/exome/RNA sequencing in 33 cancer patients: Actionable vs. utility. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 4694. doi:10.1158/1538-7445.AM2014-4694