Genomic analysis across 53 canine cancer types reveals novel mutations and high clinical actionability potential.

A genomic understanding of the oncogenic processes and individual variability of human cancer has steadily fueled improvement in patient outcomes over the past 20 years. Mutations within tumour tissues are routinely assessed through clinical genomic diagnostic assays by academic and commercial laboratories to facilitate diagnosis, prognosis and effective treatment stratification. The application of genomics has unveiled a wealth of mutation-based biomarkers in canine cancers, suggesting that the transformative principles that have revolutionized human cancer medicine can be brought to bear in veterinary oncology. To advance clinical genomics and genomics-guided medicine in canine oncology, we have developed and validated a canine cancer next-generation sequencing gene panel for the identification of multiple mutation types in clinical specimens. With this panel, we examined the genomic landscapes of 828 tumours from 813 dogs, spanning 53 cancer types. We identified 7856 alterations, encompassing copy number variants, single nucleotide variants, indels and internal tandem duplications. Additionally, we evaluated the clinical utility of these alterations by incorporating a biomarker framework from comprehensive curation of primary canine literature and inferences from human cancer genomic biomarker literature and clinical diagnostics. Remarkably, nearly 90% of the cases exhibited mutations with diagnostic, prognostic or therapeutic implications. Our work represents a thorough assessment of genomic landscapes in a large cohort of canine cancers, the first of its kind for its comprehensive inclusion of multiple mutation types and structured annotation of biomarkers, demonstrating the clinical potential of leveraging mutation-based biomarkers in veterinary oncology.

Development and Validation of an Ultradeep Next-Generation Sequencing Assay for Testing of Plasma Cell-Free DNA from Patients with Advanced Cancer.

Purpose: Tumor-derived cell-free DNA (cfDNA) in plasma can be used for molecular testing and provide an attractive alternative to tumor tissue. Commonly used PCR-based technologies can test for limited number of alterations at the time. Therefore, novel ultrasensitive technologies capable of testing for a broad spectrum of molecular alterations are needed to further personalized cancer therapy.Experimental Design: We developed a highly sensitive ultradeep next-generation sequencing (NGS) assay using reagents from TruSeqNano library preparation and NexteraRapid Capture target enrichment kits to generate plasma cfDNA sequencing libraries for mutational analysis in 61 cancer-related genes using common bioinformatics tools. The results were retrospectively compared with molecular testing of archival primary or metastatic tumor tissue obtained at different points of clinical care.Results: In a study of 55 patients with advanced cancer, the ultradeep NGS assay detected 82% (complete detection) to 87% (complete and partial detection) of the aberrations identified in discordantly collected corresponding archival tumor tissue. Patients with a low variant allele frequency (VAF) of mutant cfDNA survived longer than those with a high VAF did (P = 0.018). In patients undergoing systemic therapy, radiological response was positively associated with changes in cfDNA VAF (P = 0.02), and compared with unchanged/increased mutant cfDNA VAF, decreased cfDNA VAF was associated with longer time to treatment failure (TTF; P = 0.03).Conclusions: Ultradeep NGS assay has good sensitivity compared with conventional clinical mutation testing of archival specimens. A high VAF in mutant cfDNA corresponded with shorter survival. Changes in VAF of mutated cfDNA were associated with TTF. Clin Cancer Res; 23(18); 5648-56. ©2017 AACR.

Open-access synthetic spike-in mRNA-seq data for cancer gene fusions.

BACKGROUND: Oncogenic fusion genes underlie the mechanism of several common cancers. Next-generation sequencing based RNA-seq analyses have revealed an increasing number of recurrent fusions in a variety of cancers. However, absence of a publicly available gene-fusion focused RNA-seq data impedes comparative assessment and collaborative development of novel gene fusions detection algorithms. We have generated nine synthetic poly-adenylated RNA transcripts that correspond to previously reported oncogenic gene fusions. These synthetic RNAs were spiked at known molarity over a wide range into total RNA prior to construction of next-generation sequencing mRNA libraries to generate RNA-seq data. RESULTS: Leveraging a priori knowledge about replicates and molarity of each synthetic fusion transcript, we demonstrate utility of this dataset to compare multiple gene fusion algorithms' detection ability. In general, more fusions are detected at higher molarity, indicating that our constructs performed as expected. However, systematic detection differences are observed based on molarity or algorithm-specific characteristics. Fusion-sequence specific detection differences indicate that for applications where specific sequences are being investigated, additional constructs may be added to provide quantitative data that is specific for the sequence of interest. CONCLUSIONS: To our knowledge, this is the first publicly available synthetic RNA-seq data that specifically leverages known cancer gene-fusions. The proposed method of designing multiple gene-fusion constructs over a wide range of molarity allows granular performance analyses of multiple fusion-detection algorithms. The community can leverage and augment this publicly available data to further collaborative development of analytical tools and performance assessment frameworks for gene fusions from next-generation sequencing data.

Proteome-wide discovery of mislocated proteins in cancer.

Several studies have sought systematically to identify protein subcellular locations, but an even larger task is to map which of these proteins conditionally relocates in disease (the mislocalizome). Here, we report an integrative computational framework for mapping conditional location and mislocation of proteins on a proteome-wide scale, called a conditional location predictor (CoLP). Using CoLP, we mapped the locations of over 10,000 proteins in normal human brain and in glioma. The prediction showed 0.9 accuracy using 100 location tests of 20 randomly selected proteins. Of the 10,000 proteins, over 150 have a strong likelihood of mislocation under glioma, which is striking considering that few mislocation events have been identified in this disease previously. Using immunofluorescence and Western blotting in both primary cells and tissues, we successfully experimentally confirmed 15 mislocations. The most common type of mislocation occurs between the endoplasmic reticulum and the nucleus; for example, for RNF138, TLX3, and NFRKB. In particular, we found that the gene for the mislocating protein GFRA4 had a nonsynonymous point mutation in exon 2. Moreover, redirection of GFRA4 to its normal location, the plasma membrane, led to marked reductions in phospho-STAT3 and proliferation of glioma cells. This framework has the potential to track changes in protein location in many human diseases.

Subnetwork-based analysis of chronic lymphocytic leukemia identifies pathways that associate with disease progression.

The clinical course of patients with chronic lymphocytic leukemia (CLL) is heterogeneous. Several prognostic factors have been identified that can stratify patients into groups that differ in their relative tendency for disease progression and/or survival. Here, we pursued a subnetwork-based analysis of gene expression profiles to discriminate between groups of patients with disparate risks for CLL progression. From an initial cohort of 130 patients, we identified 38 prognostic subnetworks that could predict the relative risk for disease progression requiring therapy from the time of sample collection, more accurately than established markers. The prognostic power of these subnetworks then was validated on 2 other cohorts of patients. We noted reduced divergence in gene expression between leukemia cells of CLL patients classified at diagnosis with aggressive versus indolent disease over time. The predictive subnetworks vary in levels of expression over time but exhibit increased similarity at later time points before therapy, suggesting that degenerate pathways apparently converge into common pathways that are associated with disease progression. As such, these results have implications for understanding cancer evolution and for the development of novel treatment strategies for patients with CLL.

ROR1 is expressed in human breast cancer and associated with enhanced tumor-cell growth.

Receptor-tyrosine-kinase-like orphan receptor 1 (ROR1) is expressed during embryogenesis and by certain leukemias, but not by normal adult tissues. Here we show that the neoplastic cells of many human breast cancers express the ROR1 protein and high-level expression of ROR1 in breast adenocarcinoma was associated with aggressive disease. Silencing expression of ROR1 in human breast cancer cell lines found to express this protein impaired their growth in vitro and also in immune-deficient mice. We found that ROR1 could interact with casein kinase 1 epsilon (CK1ε) to activate phosphoinositide 3-kinase-mediated AKT phosphorylation and cAMP-response-element-binding protein (CREB), which was associated with enhanced tumor-cell growth. Wnt5a, a ligand of ROR1, could induce ROR1-dependent signaling and enhance cell growth. This study demonstrates that ROR1 is expressed in human breast cancers and has biological and clinical significance, indicating that it may be a potential target for breast cancer therapy.

B-cell activating factor and v-Myc myelocytomatosis viral oncogene homolog (c-Myc) influence progression of chronic lymphocytic leukemia.

Mice bearing a v-Myc myelocytomatosis viral oncogene homolog (c-Myc) transgene controlled by an Ig-alpha heavy-chain enhancer (iMyc(Cα) mice) rarely develop lymphomas but instead have increased rates of memory B-cell turnover and impaired antibody responses to antigen. We found that male progeny of iMyc(Cα) mice mated with mice transgenic (Tg) for CD257 (B-cell activating factor, BAFF) developed CD5(+) B-cell leukemia resembling human chronic lymphocytic leukemia (CLL), which also displays a male gender bias. Surprisingly, leukemic cells of Myc/Baff Tg mice expressed higher levels of c-Myc than did B cells of iMyc(Cα) mice. We found that CLL cells of many patients with progressive disease also expressed high amounts of c-MYC, particularly CLL cells whose survival depends on nurse-like cells (NLC), which express high-levels of BAFF. We find that BAFF could enhance CLL-cell expression of c-MYC via activation the canonical IκB kinase (IKK)/NF-κB pathway. Inhibition of the IKK/NF-κB pathway in mouse or human leukemia cells blocked the capacity of BAFF to induce c-MYC or promote leukemia-cell survival and significantly impaired disease progression in Myc/Baff Tg mice. This study reveals an important relationship between BAFF and c-MYC in CLL which may affect disease development and progression, and suggests that inhibitors of the canonical NF-κB pathway may be effective in treatment of patients with this disease.

Chronic lymphocytic leukemia of Emu-TCL1 transgenic mice undergoes rapid cell turnover that can be offset by extrinsic CD257 to accelerate disease progression.

Results of heavy-water labeling studies have challenged the notion that chronic lymphocytic leukemia (CLL) represents an accumulation of noncycling B cells. We examined leukemia cell turnover in Emu-TCL1 transgenic (TCL1-Tg) mice, which develop a CLL-like disease at 8 to 12 months of age. We found that leukemia cells in these mice not only had higher proportions of proliferating cells but also apoptotic cells than did nonleukemic lymphocytes. We crossed TCL1-Tg with BAFF-Tg mice, which express high levels of CD257. TCL1 x BAFF-Tg mice developed CLL-like disease at a significantly younger age and had more rapid disease progression and shorter survival than TCL1-Tg mice. Leukemia cells of TCL1 x BAFF-Tg mice had similar proportions of proliferating cells, but fewer proportions of dying cells, than did the CLL cells of TCL1-Tg mice. Moreover, leukemia cells from either TCL1 x BAFF-Tg or TCL1-Tg mice produced more aggressive disease when transferred into BAFF-Tg mice than into wild-type (WT) mice. Neutralization of CD257 resulted in rapid reduction in circulating leukemia cells. These results indicate that the leukemia cells of TCL1-Tg mice undergo high levels of spontaneous apoptosis that is offset by relatively high rates of leukemia cell proliferation, which might allow for acquisition of mutations that contribute to disease evolution.

Inferring pathway activity toward precise disease classification.

The advent of microarray technology has made it possible to classify disease states based on gene expression profiles of patients. Typically, marker genes are selected by measuring the power of their expression profiles to discriminate among patients of different disease states. However, expression-based classification can be challenging in complex diseases due to factors such as cellular heterogeneity within a tissue sample and genetic heterogeneity across patients. A promising technique for coping with these challenges is to incorporate pathway information into the disease classification procedure in order to classify disease based on the activity of entire signaling pathways or protein complexes rather than on the expression levels of individual genes or proteins. We propose a new classification method based on pathway activities inferred for each patient. For each pathway, an activity level is summarized from the gene expression levels of its condition-responsive genes (CORGs), defined as the subset of genes in the pathway whose combined expression delivers optimal discriminative power for the disease phenotype. We show that classifiers using pathway activity achieve better performance than classifiers based on individual gene expression, for both simple and complex case-control studies including differentiation of perturbed from non-perturbed cells and subtyping of several different kinds of cancer. Moreover, the new method outperforms several previous approaches that use a static (i.e., non-conditional) definition of pathways. Within a pathway, the identified CORGs may facilitate the development of better diagnostic markers and the discovery of core alterations in human disease.

Network-based classification of breast cancer metastasis.

Mapping the pathways that give rise to metastasis is one of the key challenges of breast cancer research. Recently, several large-scale studies have shed light on this problem through analysis of gene expression profiles to identify markers correlated with metastasis. Here, we apply a protein-network-based approach that identifies markers not as individual genes but as subnetworks extracted from protein interaction databases. The resulting subnetworks provide novel hypotheses for pathways involved in tumor progression. Although genes with known breast cancer mutations are typically not detected through analysis of differential expression, they play a central role in the protein network by interconnecting many differentially expressed genes. We find that the subnetwork markers are more reproducible than individual marker genes selected without network information, and that they achieve higher accuracy in the classification of metastatic versus non-metastatic tumors.