Germline copy number variants (CNVs) and somatic copy number alterations (SCNAs) are of significant importance in syndromic conditions and cancer. Massively parallel sequencing is increasingly used to infer copy number information from variations in the read depth in sequencing data. However, this approach has limitations in the case of targeted re-sequencing, which leaves gaps in coverage between the regions chosen for enrichment and introduces biases related to the efficiency of target capture and library preparation. We present a method for copy number detection, implemented in the software package CNVkit, that uses both the targeted reads and the nonspecifically captured off-target reads to infer copy number evenly across the genome. This combination achieves both exon-level resolution in targeted regions and sufficient resolution in the larger intronic and intergenic regions to identify copy number changes. In particular, we successfully inferred copy number at equivalent to 100-kilobase resolution genome-wide from a platform targeting as few as 293 genes. After normalizing read counts to a pooled reference, we evaluated and corrected for three sources of bias that explain most of the extraneous variability in the sequencing read depth: GC content, target footprint size and spacing, and repetitive sequences. We compared the performance of CNVkit to copy number changes identified by array comparative genomic hybridization. We packaged the components of CNVkit so that it is straightforward to use and provides visualizations, detailed reporting of significant features, and export options for integration into existing analysis pipelines. CNVkit is freely available from https://github.com/etal/cnvkit.

Biliary tract cancers such as cholangiocarcinoma represent a heterogeneous group of cancers that can be difficult to diagnose. Recent comprehensive genomic analyses in large cholangiocarcinoma cohorts have defined important molecular subgroups within cholangiocarcinoma that may relate to anatomic location and etiology [1], [2], [3], [4] and may predict responsiveness to targeted therapies in development [5], [6], [7]. These emerging data highlight the potential for tumor genomics to inform diagnosis and treatment options in this challenging tumor type. We report the case of a patient with a germline BRCA1 mutation who presented with a cholangiocarcinoma driven by the novel YWHAZ-BRAF fusion. Hybrid capture-based DNA sequencing and copy number analysis performed as part of clinical care demonstrated that two later-occurring tumors were clonally derived from the primary cholangiocarcinoma rather than distinct new primaries, revealing an unusual pattern of late metachronous metastasis. We discuss the clinical significance of these genetic alterations and their relevance to therapeutic strategies.

Key points

Hybrid capture-based next-generation DNA sequencing assays can provide diagnostic clarity in patients with unusual patterns of metastasis and recurrence in which the pathologic diagnosis is ambiguous.To our knowledge, this is the first reported case of a YWHAZ-BRAF fusion in pancreaticobiliary cancer, and a very rare case of cholangiocarcinoma in the setting of a germline BRCA1 mutation.The patient's BRCA1 mutation and YWHAZ-BRAF fusion constitute potential targets for future therapy.

Oncogenic fusions in TRK family receptor tyrosine kinases have been identified in several cancers and can serve as therapeutic targets. We identified ETV6-NTRK3, MYO5A-NTRK3 and MYH9-NTRK3 fusions in Spitz tumours, and demonstrated that NTRK3 fusions constitutively activate the mitogen-activated protein kinase, phosphoinositide 3-kinase and phospholipase Cγ1 pathways in melanocytes. This signalling was inhibited by DS-6051a, a small-molecule inhibitor of NTRK1/2/3 and ROS1. NTRK3 fusions expand the range of oncogenic kinase fusions in melanocytic neoplasms and offer targets for a small subset of melanomas for which no targeted options currently exist. Copyright © 2016 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

BRAF is the most prevalent oncogene and an important therapeutic target in melanoma. In some cancers, BRAF is activated by rearrangements that fuse its kinase domain to 5' partner genes. We examined 848 comparative genomic hybridization profiles of melanocytic tumors and found copy number transitions within BRAF in 10 tumors, of which six could be further characterized by sequencing. In all, the BRAF kinase domain was fused in-frame to six N-terminal partners. No other mutations were identified in melanoma oncogenes. One of the seven melanoma cell lines without known oncogenic mutations harbored a similar BRAF fusion, which constitutively activated the MAP kinase pathway. Sorafenib, but not vemurafenib, could block MAP kinase pathway activation and proliferation of the cell line at clinically relevant concentrations, whereas BRAF(V) (600E) mutant melanoma cell lines were significantly more sensitive to vemurafenib. The patient from whom the cell line was derived showed a durable clinical response to sorafenib.

BRAF fusions are detected in numerous neoplasms, but their clinical management remains unresolved. We identified six melanoma lines harboring BRAF fusions representative of the clinical cases reported in the literature. Their unexpected heterogeneous responses to RAF and MEK inhibitors could be categorized upon specific features of the fusion kinases. Higher expression level correlated with resistance, and fusion partners containing a dimerization domain promoted paradoxical activation of the mitogen-activated protein kinase (MAPK) pathway and hyperproliferation in response to first- and second-generation RAF inhibitors. By contrast, next-generation αC-IN/DFG-OUT RAF inhibitors blunted paradoxical activation across all lines and had their therapeutic efficacy further increased in vitro and in vivo by combination with MEK inhibitors, opening perspectives in the clinical management of tumors harboring BRAF fusions.

Loss of the CDKN2A tumor suppressor is associated with melanoma metastasis, but the mechanisms connecting the phenomena are unknown. Using CRISPR-Cas9 to engineer a cellular model of melanoma initiation from primary human melanocytes, we discovered that a lineage-restricted transcription factor, BRN2, is downstream of CDKN2A and directly regulated by E2F1. In a cohort of melanocytic tumors that capture distinct progression stages, we observed that CDKN2A loss coincides with both the onset of invasive behavior and increased BRN2 expression. Loss of the CDKN2A protein product p16^INK4A permitted metastatic dissemination of human melanoma lines in mice, a phenotype rescued by inhibition of BRN2. These results demonstrate a mechanism by which CDKN2A suppresses the initiation of melanoma invasion through inhibition of BRN2.

Spitzoid neoplasms are a group of melanocytic tumours with distinctive histopathological features. They include benign tumours (Spitz naevi), malignant tumours (spitzoid melanomas) and tumours with borderline histopathological features and uncertain clinical outcome (atypical Spitz tumours). Their genetic underpinnings are poorly understood, and alterations in common melanoma-associated oncogenes are typically absent. Here we show that spitzoid neoplasms harbour kinase fusions of ROS1 (17%), NTRK1 (16%), ALK (10%), BRAF (5%) and RET (3%) in a mutually exclusive pattern. The chimeric proteins are constitutively active, stimulate oncogenic signalling pathways, are tumourigenic and are found in the entire biologic spectrum of spitzoid neoplasms, including 55% of Spitz naevi, 56% of atypical Spitz tumours and 39% of spitzoid melanomas. Kinase inhibitors suppress the oncogenic signalling of the fusion proteins in vitro. In summary, kinase fusions account for the majority of oncogenic aberrations in spitzoid neoplasms and may serve as therapeutic targets for metastatic spitzoid melanomas.