Genomic Classification and Individualized Prognosis in Multiple Myeloma.

PURPOSE: Outcomes for patients with newly diagnosed multiple myeloma (NDMM) are heterogenous, with overall survival (OS) ranging from months to over 10 years. METHODS: To decipher and predict the molecular and clinical heterogeneity of NDMM, we assembled a series of 1,933 patients with available clinical, genomic, and therapeutic data. RESULTS: Leveraging a comprehensive catalog of genomic drivers, we identified 12 groups, expanding on previous gene expression-based molecular classifications. To build a model predicting individualized risk in NDMM (IRMMa), we integrated clinical, genomic, and treatment variables. To correct for time-dependent variables, including high-dose melphalan followed by autologous stem-cell transplantation (HDM-ASCT), and maintenance therapy, a multi-state model was designed. The IRMMa model accuracy was significantly higher than all comparator prognostic models, with a c-index for OS of 0.726, compared with International Staging System (ISS; 0.61), revised-ISS (0.572), and R2-ISS (0.625). Integral to model accuracy was 20 genomic features, including 1q21 gain/amp, del 1p, TP53 loss, NSD2 translocations, APOBEC mutational signatures, and copy-number signatures (reflecting the complex structural variant chromothripsis). IRMMa accuracy and superiority compared with other prognostic models were validated on 256 patients enrolled in the GMMG-HD6 (ClinicalTrials.gov identifier: NCT02495922) clinical trial. Individualized patient risks were significantly affected across the 12 genomic groups by different treatment strategies (ie, treatment variance), which was used to identify patients for whom HDM-ASCT is particularly effective versus patients for whom the impact is limited. CONCLUSION: Integrating clinical, demographic, genomic, and therapeutic data, to our knowledge, we have developed the first individualized risk-prediction model enabling personally tailored therapeutic decisions for patients with NDMM.

CK1δ and CK1ε Signaling Sustains Mitochondrial Metabolism and Cell Survival in Multiple Myeloma.

Multiple myeloma remains an incurable malignancy due to acquisition of intrinsic programs that drive therapy resistance. Here we report that casein kinase-1δ (CK1δ) and CK1ε are therapeutic targets in multiple myeloma that are necessary to sustain mitochondrial metabolism. Specifically, the dual CK1δ/CK1ε inhibitor SR-3029 had potent in vivo and ex vivo anti-multiple myeloma activity, including against primary multiple myeloma patient specimens. RNA sequencing (RNA-seq) and metabolic analyses revealed inhibiting CK1δ/CK1ε disables multiple myeloma metabolism by suppressing genes involved in oxidative phosphorylation (OxPhos), reducing citric acid cycle intermediates, and suppressing complexes I and IV of the electron transport chain. Finally, sensitivity of multiple myeloma patient specimens to SR-3029 correlated with elevated expression of mitochondrial genes, and RNA-seq from 687 multiple myeloma patient samples revealed that increased CSNK1D, CSNK1E, and OxPhos genes correlate with disease progression and inferior outcomes. Thus, increases in mitochondrial metabolism are a hallmark of multiple myeloma progression that can be disabled by targeting CK1δ/CK1ε. SIGNIFICANCE: CK1δ and CK1ε are attractive therapeutic targets in multiple myeloma whose expression increases with disease progression and connote poor outcomes, and that are necessary to sustain expression of genes directing OxPhos.

Whole-Exome sequencing analysis identified TMSB10/TRABD2A locus to be associated with carfilzomib-related cardiotoxicity among patients with multiple myeloma.

Background: Proteasome inhibitor Carfilzomib (CFZ) is effective in treating patients with refractory or relapsed multiple myeloma (MM) but has been associated with cardiovascular adverse events (CVAE) such as hypertension, cardiomyopathy, and heart failure. This study aimed to investigate the contribution of germline genetic variants in protein-coding genes in CFZ-CVAE among MM patients using whole-exome sequencing (WES) analysis. Methods: Exome-wide single-variant association analysis, gene-based analysis, and rare variant analyses were performed on 603,920 variants in 247 patients with MM who have been treated with CFZ and enrolled in the Oncology Research Information Exchange Network (ORIEN) at the Moffitt Cancer Center. Separate analyses were performed in European Americans and African Americans followed by a trans-ethnic meta-analysis. Results: The most significant variant in the exome-wide single variant analysis was a missense variant rs7148 in the thymosin beta-10/TraB Domain Containing 2A (TMSB10/TRABD2A) locus. The effect allele of rs7148 was associated with a higher risk of CVAE [odds ratio (OR) = 9.3 with a 95% confidence interval of 3.9-22.3, p = 5.42*10-7]. MM patients with rs7148 AG or AA genotype had a higher risk of CVAE (50%) than those with GG genotype (10%). rs7148 is an expression quantitative trait locus (eQTL) for TRABD2A and TMSB10. The gene-based analysis also showed TRABD2A as the most significant gene associated with CFZ-CVAE (p = 1.06*10-6). Conclusions: We identified a missense SNP rs7148 in the TMSB10/TRABD2A as associated with CFZ-CVAE in MM patients. More investigation is needed to understand the underlying mechanisms of these associations.

Replicative Instability Drives Cancer Progression.

In the past decade, defective DNA repair has been increasingly linked with cancer progression. Human tumors with markers of defective DNA repair and increased replication stress exhibit genomic instability and poor survival rates across tumor types. Seminal studies have demonstrated that genomic instability develops following inactivation of BRCA1, BRCA2, or BRCA-related genes. However, it is recognized that many tumors exhibit genomic instability but lack BRCA inactivation. We sought to identify a pan-cancer mechanism that underpins genomic instability and cancer progression in BRCA-wildtype tumors. Methods: Using multi-omics data from two independent consortia, we analyzed data from dozens of tumor types to identify patient cohorts characterized by poor outcomes, genomic instability, and wildtype BRCA genes. We developed several novel metrics to identify the genetic underpinnings of genomic instability in tumors with wildtype BRCA. Associated clinical data was mined to analyze patient responses to standard of care therapies and potential differences in metastatic dissemination. Results: Systematic analysis of the DNA repair landscape revealed that defective single-strand break repair, translesion synthesis, and non-homologous end-joining effectors drive genomic instability in tumors with wildtype BRCA and BRCA-related genes. Importantly, we find that loss of these effectors promotes replication stress, therapy resistance, and increased primary carcinoma to brain metastasis. Conclusions: Our results have defined a new pan-cancer class of tumors characterized by replicative instability (RIN). RIN is defined by the accumulation of intra-chromosomal, gene-level gain and loss events at replication stress sensitive (RSS) genome sites. We find that RIN accelerates cancer progression by driving copy number alterations and transcriptional program rewiring that promote tumor evolution. Clinically, we find that RIN drives therapy resistance and distant metastases across multiple tumor types.

Racial and ethnic differences in clonal hematopoiesis, tumor markers, and outcomes of patients with multiple myeloma.

Multiple myeloma (MM) incidence, mortality, and survival vary by race and ethnicity, but the causes of differences remain unclear. We investigated demographic, clinical, and molecular features of diverse MM patients to elucidate mechanisms driving clinical disparities. This study included 495 MM patients (self-reported Hispanic, n = 45; non-Hispanic Black, n = 52; non-Hispanic White, n = 398). Hispanic and non-Hispanic Black individuals had an earlier age of onset than non-Hispanic White individuals (53 and 57 vs 63 years, respectively, P < .001). There were no differences in treatment by race and ethnicity groups, but non-Hispanic Black patients had a longer time to hematopoietic cell transplant than non-Hispanic White patients (376 days vs 248 days; P = .01). Overall survival (OS) was improved for non-Hispanic Black compared with non-Hispanic White patients (HR, 0.50; 95% CI, 0.31-0.81; P = .005), although this association was attenuated after adjusting for clinical features (HR, 0.62; 95% CI, 0.37-1.03; P = .06). Tumor mutations in IRF4 were most common in Hispanic patients, and mutations in SP140, AUTS2, and SETD2 were most common in non-Hispanic Black patients. Differences in tumor expression of BCL7A, SPEF2, and ANKRD26 by race and ethnicity were observed. Clonal hematopoiesis was detected in 12% of patients and associated with inferior OS in non-Hispanic Black patients compared with patients without clonal hematopoiesis (HR, 4.36; 95% CI, 1.36-14.00). This study provides insight into differences in molecular features that may drive clinical disparities in MM patients receiving comparable treatment, with the novel inclusion of Hispanic individuals.

Plasma cell dependence on histone/protein deacetylase 11 reveals a therapeutic target in multiple myeloma.

The clinical utility of histone/protein deacetylase (HDAC) inhibitors in combinatorial regimens with proteasome inhibitors for patients with relapsed and refractory multiple myeloma (MM) is often limited by excessive toxicity due to HDAC inhibitor promiscuity with multiple HDACs. Therefore, more selective inhibition minimizing off-target toxicity may increase the clinical effectiveness of HDAC inhibitors. We demonstrated that plasma cell development and survival are dependent upon HDAC11, suggesting this enzyme is a promising therapeutic target in MM. Mice lacking HDAC11 exhibited markedly decreased plasma cell numbers. Accordingly, in vitro plasma cell differentiation was arrested in B cells lacking functional HDAC11. Mechanistically, we showed that HDAC11 is involved in the deacetylation of IRF4 at lysine103. Further, targeting HDAC11 led to IRF4 hyperacetylation, resulting in impaired IRF4 nuclear localization and target promoter binding. Importantly, transient HDAC11 knockdown or treatment with elevenostat, an HDAC11-selective inhibitor, induced cell death in MM cell lines. Elevenostat produced similar anti-MM activity in vivo, improving survival among mice inoculated with 5TGM1 MM cells. Elevenostat demonstrated nanomolar ex vivo activity in 34 MM patient specimens and synergistic activity when combined with bortezomib. Collectively, our data indicated that HDAC11 regulates an essential pathway in plasma cell biology establishing its potential as an emerging theraputic vulnerability in MM.

Metabolic Changes Are Associated with Melphalan Resistance in Multiple Myeloma.

Multiple myeloma is an incurable hematological malignancy that impacts tens of thousands of people every year in the United States. Treatment for eligible patients involves induction, consolidation with stem cell rescue, and maintenance. High-dose therapy with a DNA alkylating agent, melphalan, remains the primary drug for consolidation therapy in conjunction with autologous stem-cell transplantation; as such, melphalan resistance remains a relevant clinical challenge. Here, we describe a proteometabolomic approach to examine mechanisms of acquired melphalan resistance in two cell line models. Drug metabolism, steady-state metabolomics, activity-based protein profiling (ABPP, data available at PRIDE: PXD019725), acute-treatment metabolomics, and western blot analyses have allowed us to further elucidate metabolic processes associated with melphalan resistance. Proteometabolomic data indicate that drug-resistant cells have higher levels of pentose phosphate pathway metabolites. Purine, pyrimidine, and glutathione metabolisms were commonly altered, and cell-line-specific changes in metabolite levels were observed, which could be linked to the differences in steady-state metabolism of naïve cells. Inhibition of selected enzymes in purine synthesis and pentose phosphate pathways was evaluated to determine their potential to improve melphalan's efficacy. The clinical relevance of these proteometabolomic leads was confirmed by comparison of tumor cell transcriptomes from newly diagnosed MM patients and patients with relapsed disease after treatment with high-dose melphalan and autologous stem-cell transplantation. The observation of common and cell-line-specific changes in metabolite levels suggests that omic approaches will be needed to fully examine melphalan resistance in patient specimens and define personalized strategies to optimize the use of high-dose melphalan.

DACH1 mutation frequency in endometrial cancer is associated with high tumor mutation burden.

OBJECTIVE: DACH1 is a transcriptional repressor and tumor suppressor gene frequently mutated in melanoma, bladder, and prostate cancer. Loss of DACH1 expression is associated with poor prognostic features and reduced overall survival in uterine cancer. In this study, we utilized the Oncology Research Information Exchange Network (ORIEN) Avatar database to determine the frequency of DACH1 mutations in patients with endometrial cancer in our Kentucky population. METHODS: We obtained clinical and genomic data for 65 patients with endometrial cancer from the Markey Cancer Center (MCC). We examined the clinical attributes of the cancers by DACH1 status by comparing whole-exome sequencing (WES), RNA Sequencing (RNASeq), microsatellite instability (MSI), and tumor mutational burden (TMB). RESULTS: Kentucky women with endometrial cancer had an increased frequency of DACH1 mutations (12/65 patients, 18.5%) compared to The Cancer Genome Atlas (TCGA) endometrial cancer population (25/586 patients, 3.8%) with p-value = 1.04E-05. DACH1 mutations were associated with increased tumor mutation count in both TCGA (median 65 vs. 8972, p-value = 7.35E-09) and our Kentucky population (490 vs. 2160, p-value = 6.0E-04). DACH1 mutated patients have a higher tumor mutation burden compared to DACH1 wild-type (24 vs. 6.02, p-value = 4.29E-05). DACH1 mutations showed significant gene co-occurrence patterns with POLE, MLH1, and PMS2. DACH1 mutations were not associated with an increase in microsatellite instability at MCC (MSI-H) (p-value = 0.1342). CONCLUSIONS: DACH1 mutations are prevalent in Kentucky patients with endometrial cancer. These mutations are associated with high tumor mutational burden and co-occur with genome destabilizing gene mutations. These findings suggest DACH1 may be a candidate biomarker for future trials with immunotherapy, particularly in endometrial cancers.

NF-κB and STAT3 co-operation enhances high glucose induced aggressiveness of cholangiocarcinoma cells.

AIMS: The present report aimed to investigate the underlying genes and pathways of high glucose driving cholangiocarcinoma (CCA) aggressiveness. MAIN METHODS: We screened and compared the gene expression profiles obtained by RNA sequencing, of CCA cells cultured in high and normal glucose. Results from the transcriptomic analysis were confirmed in additional cell lines using in vitro migration-invasion assay, Western blotting and immunocytofluorescence. KEY FINDINGS: Data indicated that high glucose increased the expression of interleukin-1ß (IL-1ß), an upstream regulator of nuclear factor-κB (NF-κB) pathway, through the nuclear localization of NF-κB. High glucose-induced NF-κB increased the migration and invasion of CCA cells and the expression of downstream NF-κB targeted genes associated with aggressiveness, including interleukin-6, a potent triggering signal of the signal transducer and activator of transcription 3 (STAT3) pathway. Such effects were reversed by inhibiting NF-κB nuclear translocation which additionally reduced the phosphorylation of STAT3 at Y705. SIGNIFICANCE: These results indicate that NF-κB is activated by high glucose and they suggest that NF-κB interaction with STAT3 enhances CCA aggressiveness. Therefore, targeting multiple pathways such as STAT3 and NF-κB might improve CCA treatment outcome especially in condition such as hyperglycemia.

Functional Annotation of ESR1 Gene Fusions in Estrogen Receptor-Positive Breast Cancer.

RNA sequencing (RNA-seq) detects estrogen receptor alpha gene (ESR1) fusion transcripts in estrogen receptor-positive (ER+) breast cancer, but their role in disease pathogenesis remains unclear. We examined multiple ESR1 fusions and found that two, both identified in advanced endocrine treatment-resistant disease, encoded stable and functional fusion proteins. In both examples, ESR1-e6>YAP1 and ESR1-e6>PCDH11X, ESR1 exons 1-6 were fused in frame to C-terminal sequences from the partner gene. Functional properties include estrogen-independent growth, constitutive expression of ER target genes, and anti-estrogen resistance. Both fusions activate a metastasis-associated transcriptional program, induce cellular motility, and promote the development of lung metastasis. ESR1-e6>YAP1- and ESR1-e6>PCDH11X-induced growth remained sensitive to a CDK4/6 inhibitor, and a patient-derived xenograft (PDX) naturally expressing the ESR1-e6>YAP1 fusion was also responsive. Transcriptionally active ESR1 fusions therefore trigger both endocrine therapy resistance and metastatic progression, explaining the association with fatal disease progression, although CDK4/6 inhibitor treatment is predicted to be effective.

SVachra: a tool to identify genomic structural variation in mate pair sequencing data containing inward and outward facing reads.

BACKGROUND: Characterization of genomic structural variation (SV) is essential to expanding the research and clinical applications of genome sequencing. Reliance upon short DNA fragment paired end sequencing has yielded a wealth of single nucleotide variants and internal sequencing read insertions-deletions, at the cost of limited SV detection. Multi-kilobase DNA fragment mate pair sequencing has supplemented the void in SV detection, but introduced new analytic challenges requiring SV detection tools specifically designed for mate pair sequencing data. Here, we introduce SVachra - Structural Variation Assessment of CHRomosomal Aberrations, a breakpoint calling program that identifies large insertions-deletions, inversions, inter- and intra-chromosomal translocations utilizing both inward and outward facing read types generated by mate pair sequencing. RESULTS: We demonstrate SVachra's utility by executing the program on large-insert (Illumina Nextera) mate pair sequencing data from the personal genome of a single subject (HS1011). An additional data set of long-read (Pacific BioSciences RSII) was also generated to validate SV calls from SVachra and other comparison SV calling programs. SVachra exhibited the highest validation rate and reported the widest distribution of SV types and size ranges when compared to other SV callers. CONCLUSIONS: SVachra is a highly specific breakpoint calling program that exhibits a more unbiased SV detection methodology than other callers.

Therapy-related Acute Leukemia With Mixed Phenotype and Novel t(1: 6)(q25;p23) After Treatment for High-risk Neuroblastoma.

Neuroblastoma is the most common extracranial malignancy of childhood. Patients with high-risk disease receive multimodal treatment including chemotherapy combinations containing alkylating agents and topoisomerase inhibitors with potential for inducing therapy-related malignancy later in life. Most commonly, cytogenetic changes of pediatric therapy-related myelodysplastic syndrome/acute myeloid leukemia involve chromosome 5 or 7. Here we report a novel case of therapy-related myelodysplastic syndrome/acute myeloid leukemia 30 months after treatment for high-risk neuroblastoma with biphenotypic cell surface markers and a not yet described translocation t(1;6)(q25;p23).

A Children's Oncology Group and TARGET initiative exploring the genetic landscape of Wilms tumor.

We performed genome-wide sequencing and analyzed mRNA and miRNA expression, DNA copy number, and DNA methylation in 117 Wilms tumors, followed by targeted sequencing of 651 Wilms tumors. In addition to genes previously implicated in Wilms tumors (WT1, CTNNB1, AMER1, DROSHA, DGCR8, XPO5, DICER1, SIX1, SIX2, MLLT1, MYCN, and TP53), we identified mutations in genes not previously recognized as recurrently involved in Wilms tumors, the most frequent being BCOR, BCORL1, NONO, MAX, COL6A3, ASXL1, MAP3K4, and ARID1A. DNA copy number changes resulted in recurrent 1q gain, MYCN amplification, LIN28B gain, and MIRLET7A loss. Unexpected germline variants involved PALB2 and CHEK2. Integrated analyses support two major classes of genetic changes that preserve the progenitor state and/or interrupt normal development.

Activating MAPK1 (ERK2) mutation in an aggressive case of disseminated juvenile xanthogranuloma.

Juvenile xanthogranuloma (JXG) is a rare histiocytic disorder that is usually benign and self-limiting. We present a case of atypical, aggressive JXG harboring a novel mitogen-activated protein kinase (MAPK) pathway mutation in the MAPK1 gene, which encodes mitogen-activated protein kinase 1 or extracellular signal-regulated 2 (ERK2). Our analysis revealed that the mutation results in constitutive ERK activation that is resistant to BRAF or MEK inhibitors but susceptible to an ERK inhibitor. These data highlight the importance of identifying specific MAPK pathway alterations as part of the diagnostic workup for patients with histiocytic disorders rather than initiating empiric treatment with MEK inhibitors.

Non-malignant respiratory epithelial cells preferentially proliferate from resected non-small cell lung cancer specimens cultured under conditionally reprogrammed conditions.

The "conditionally reprogrammed cells" (CRC) method, using a Rho kinase inhibitor and irradiated mouse fibroblast cells has been described for the efficient growth of cells from malignant and non-malignant samples from primary tumor and non-malignant sites. Using the CRC method, four institutions independently cultured tumor tissues from 48 non-small cell lung cancers (NSCLC, mostly from primary resected tumors) and 22 non-malignant lungs. We found that epithelial cells could be cultured from tumor and non-malignant lung. However, epithelial cells cultured from tumors had features of non-malignant respiratory epithelial cells which include: 1) among 22 mutations found in the original tumors only two mutations were found in the CRC cultures with reduced frequency (31% to 13% and 92% to 15% from original tumor and CRC culture respectively); 2) copy number variation was analyzed in 9 tumor and their CRC cultures and only diploid patterns were found in CRC cultures; 3) mRNA expression profiles were similar to those of normal respiratory epithelial cells; and 4) co-culture of tumor and non-malignant lung epithelial cells resulted in mostly non-malignant cells. We conclude that CRC method is a highly selective and useful method for the growth of non-malignant respiratory epithelial cells from tumor specimens and only occasionally do such CRC cultures contain a small subpopulation of cancer cells marked by oncogenic mutations. While our findings are restricted to resected primary NSCLC, they indicated the necessity to fully characterize all CRC cultures and the need to develop culture technology that facilitates the growth of primary lung cancers.

Genomic analysis of hepatoblastoma identifies distinct molecular and prognostic subgroups.

Despite being the most common liver cancer in children, hepatoblastoma (HB) is a rare neoplasm. Consequently, few pretreatment tumors have been molecularly profiled, and there are no validated prognostic or therapeutic biomarkers for HB patients. We report on the first large-scale effort to profile pretreatment HBs at diagnosis. Our analysis of 88 clinically annotated HBs revealed three risk-stratifying molecular subtypes that are characterized by differential activation of hepatic progenitor cell markers and metabolic pathways: high-risk tumors were characterized by up-regulated nuclear factor, erythroid 2-like 2 activity; high lin-28 homolog B, high mobility group AT-hook 2, spalt-like transcription factor 4, and alpha-fetoprotein expression; and high coordinated expression of oncofetal proteins and stem-cell markers, while low-risk tumors had low lin-28 homolog B and lethal-7 expression and high hepatic nuclear factor 1 alpha activity. CONCLUSION: Analysis of immunohistochemical assays using antibodies targeting these genes in a prospective study of 35 HBs suggested that these candidate biomarkers have the potential to improve risk stratification and guide treatment decisions for HB patients at diagnosis; our results pave the way for clinical collaborative studies to validate candidate biomarkers and test their potential to improve outcome for HB patients. (Hepatology 2017;65:104-121).

Alternative genetic mechanisms of BRAF activation in Langerhans cell histiocytosis.

Langerhans cell histiocytosis (LCH) is characterized by inflammatory lesions containing pathologic CD207+ dendritic cells with constitutively activated ERK. Mutually exclusive somatic mutations in MAPK pathway genes have been identified in ∼75% of LCH cases, including recurrent BRAF-V600E and MAP2K1 mutations. To elucidate mechanisms of ERK activation in the remaining 25% of patients, we performed whole-exome sequencing (WES, n = 6), targeted BRAF sequencing (n = 19), and/or whole-transcriptome sequencing (RNA-seq, n = 6) on 24 LCH patient samples lacking BRAF-V600E or MAP2K1 mutations. WES and BRAF sequencing identified in-frame BRAF deletions in the ß3-αC loop in 6 lesions. RNA-seq revealed one case with an in-frame FAM73A-BRAF fusion lacking the BRAF autoinhibitory regulatory domain but retaining an intact kinase domain. High levels of phospho-ERK were detected in vitro in cells overexpressing either BRAF fusion or deletion constructs and ex vivo in CD207+ cells from lesions. ERK activation was resistant to BRAF-V600E inhibition, but responsive to both a second-generation BRAF inhibitor and a MEK inhibitor. These results support an emerging model of universal ERK-activating genetic alterations driving pathogenesis in LCH. A personalized approach in which patient-specific alterations are identified may be necessary to maximize benefit from targeted therapies for patients with LCH.

Significance of TP53 Mutation in Wilms Tumors with Diffuse Anaplasia: A Report from the Children's Oncology Group.

PURPOSE: To investigate the role and significance of TP53 mutation in diffusely anaplastic Wilms tumors (DAWTs). EXPERIMENTAL DESIGN: All DAWTs registered on National Wilms Tumor Study-5 (n = 118) with available samples were analyzed for TP53 mutations and copy loss. Integrative genomic analysis was performed on 39 selected DAWTs. RESULTS: Following analysis of a single random sample, 57 DAWTs (48%) demonstrated TP53 mutations, 13 (11%) copy loss without mutation, and 48 (41%) lacked both [defined as TP53-wild-type (wt)]. Patients with stage III/IV TP53-wt DAWTs (but not those with stage I/II disease) had significantly lower relapse and death rates than those with TP53 abnormalities. In-depth analysis of a subset of 39 DAWTs showed seven (18%) to be TP53-wt: These demonstrated gene expression evidence of an active p53 pathway. Retrospective pathology review of TP53-wt DAWT revealed no or very low volume of anaplasia in six of seven tumors. When samples from TP53-wt tumors known to contain anaplasia histologically were available, abnormal p53 protein accumulation was observed by immunohistochemistry. CONCLUSIONS: These data support the key role of TP53 loss in the development of anaplasia in WT, and support its significant clinical impact in patients with residual anaplastic tumor following surgery. These data also suggest that most DAWTs will show evidence of TP53 mutation when samples selected for the presence of anaplasia are analyzed. This suggests that modifications of the current criteria to also consider volume of anaplasia and documentation of TP53 aberrations may better reflect the risk of relapse and death and enable optimization of therapeutic stratification. Clin Cancer Res; 22(22); 5582-91. ©2016 AACR.

ITD assembler: an algorithm for internal tandem duplication discovery from short-read sequencing data.

BACKGROUND: Detection of tandem duplication within coding exons, referred to as internal tandem duplication (ITD), remains challenging due to inefficiencies in alignment of ITD-containing reads to the reference genome. There is a critical need to develop efficient methods to recover these important mutational events. RESULTS: In this paper we introduce ITD Assembler, a novel approach that rapidly evaluates all unmapped and partially mapped reads from whole exome NGS data using a De Bruijn graphs approach to select reads that harbor cycles of appropriate length, followed by assembly using overlap-layout-consensus. We tested ITD Assembler on The Cancer Genome Atlas AML dataset as a truth set. ITD Assembler identified the highest percentage of reported FLT3-ITDs when compared to other ITD detection algorithms, and discovered additional ITDs in FLT3, KIT, CEBPA, WT1 and other genes. Evidence of polymorphic ITDs in 54 genes were also found. Novel ITDs were validated by analyzing the corresponding RNA sequencing data. CONCLUSIONS: ITD Assembler is a very sensitive tool which can detect partial, large and complex tandem duplications. This study highlights the need to more effectively look for ITD's in other cancers and Mendelian diseases.

Novel patient-derived xenograft and cell line models for therapeutic testing of pediatric liver cancer.

BACKGROUND & AIMS: Pediatric liver cancer is a rare but serious disease whose incidence is rising, and for which the therapeutic options are limited. Development of more targeted, less toxic therapies is hindered by the lack of an experimental animal model that captures the heterogeneity and metastatic capability of these tumors. METHODS: Here we established an orthotopic engraftment technique to model a series of patient-derived tumor xenograft (PDTX) from pediatric liver cancers of all major histologic subtypes: hepatoblastoma, hepatocellular cancer and hepatocellular malignant neoplasm. We utilized standard (immuno) staining methods for histological characterization, RNA sequencing for gene expression profiling and genome sequencing for identification of druggable targets. We also adapted stem cell culturing techniques to derive two new pediatric cancer cell lines from the xenografted mice. RESULTS: The patient-derived tumor xenografts recapitulated the histologic, genetic, and biological characteristics-including the metastatic behavior-of the corresponding primary tumors. Furthermore, the gene expression profiles of the two new liver cancer cell lines closely resemble those of the primary tumors. Targeted therapy of PDTX from an aggressive hepatocellular malignant neoplasm with the MEK1 inhibitor trametinib and pan-class I PI3 kinase inhibitor NVP-BKM120 resulted in significant growth inhibition, thus confirming this PDTX model as a valuable tool to study tumor biology and patient-specific therapeutic responses. CONCLUSIONS: The novel metastatic xenograft model and the isogenic xenograft-derived cell lines described in this study provide reliable tools for developing mutation- and patient-specific therapies for pediatric liver cancer. LAY SUMMARY: Pediatric liver cancer is a rare but serious disease and no experimental animal model currently captures the complexity and metastatic capability of these tumors. We have established a novel animal model using human tumor tissue that recapitulates the genetic and biological characteristics of this cancer. We demonstrate that our patient-derived animal model, as well as two new cell lines, are useful tools for experimental therapies.