Search

Cutaneous squamous cell carcinoma is a form of skin cancer originating from keratinocytes in the skin. It is the second most common type of cancer and is responsible for an estimated 8000 deaths per year in the United States. Compared to other cancer subtypes with similar incidences and death tolls, our understanding of the somatic mutations driving cutaneous squamous cell carcinoma is limited. The main challenge is that these tumors have high mutation burdens, primarily a consequence of UV-radiation-induced DNA damage from sunlight, making it difficult to distinguish driver mutations from passenger mutations. We overcame this challenge by performing a meta-analysis of publicly available sequencing data covering 105 tumors from 10 different studies. Moreover, we eliminated tumors with issues, such as low neoplastic cell content, and from the tumors that passed quality control, we utilized multiple strategies to reveal genes under selection. In total, we nominated 30 cancer genes. Among the more novel genes, mutations frequently affected EP300, PBRM1, USP28, and CHUK. Collectively, mutations in the NOTCH and p53 pathways were ubiquitous, and to a lesser extent, mutations affected genes in the Hippo pathway, genes in the Ras/MAPK/PI3K pathway, genes critical for cell-cycle checkpoint control, and genes encoding chromatin remodeling factors. Taken together, our study provides a catalogue of driver genes in cutaneous squamous cell carcinoma, offering points of therapeutic intervention and insights into the biology of cutaneous squamous cell carcinoma.

Every cell in the human body has a unique set of somatic mutations, yet it remains difficult to comprehensively genotype an individual cell1. Here, we developed solutions to overcome this obstacle in the context of normal human skin, thus offering the first glimpse into the genomic landscapes of individual melanocytes from human skin. As expected, sun-shielded melanocytes had fewer mutations than sun-exposed melanocytes. However, within sun-exposed sites, melanocytes on chronically sun-exposed skin (e.g. the face) displayed a lower mutation burden than melanocytes on intermittently sun-exposed skin (e.g. the back). Melanocytes located adjacent to a skin cancer had higher mutation burdens than melanocytes from donors without skin cancer, implying that the mutation burden of normal skin can be harnessed to measure cumulative sun damage and skin cancer risk. Moreover, melanocytes from healthy skin commonly harbor pathogenic mutations, though these mutations tended to be weakly oncogenic, likely explaining why they did not give rise to discernible lesions. Phylogenetic analyses identified groups of related melanocytes, suggesting that melanocytes spread throughout skin as fields of clonally related cells, invisible to the naked eye. Overall, our study offers an unprecedented view into the genomic landscapes of individual melanocytes, revealing key insights into the causes and origins of melanoma.

Epigenetic landscapes can shape physiologic and disease phenotypes. We used integrative, high resolution multi-omics methods to delineate the methylome landscape and characterize the oncogenic drivers of esophageal squamous cell carcinoma (ESCC). We found 98% of CpGs are hypomethylated across the ESCC genome. Hypo-methylated regions are enriched in areas with heterochromatin binding markers (H3K9me3, H3K27me3), while hyper-methylated regions are enriched in polycomb repressive complex (EZH2/SUZ12) recognizing regions. Altered methylation in promoters, enhancers, and gene bodies, as well as in polycomb repressive complex occupancy and CTCF binding sites are associated with cancer-specific gene dysregulation. Epigenetic-mediated activation of non-canonical WNT/β-catenin/MMP signaling and a YY1/lncRNA ESCCAL-1/ribosomal protein network are uncovered and validated as potential novel ESCC driver alterations. This study advances our understanding of how epigenetic landscapes shape cancer pathogenesis and provides a resource for biomarker and target discovery., The epigenetic landscape of esophageal squamous cell carcinoma (ESCC) at genome-wide high resolution is incompletely studied. Here, the authors performed an integrated multi-omics analysis of ESCC and non-tumor tissues to define the genome-wide methylome landscape and epigenetic alterations to uncover oncogenic drivers of ESCC.

To date, over 1000 melanocytic neoplasms, spanning all stages of tumorigenesis, have been sequenced, offering detailed views into their -omic landscapes. This has coincided with advances in genetic engineering technologies that allow molecular biologists to edit the human genome with extreme precision and new mouse models to simulate disease progression. In this review, we describe how these technologies are being harnessed to provide insights into the evolution of melanoma at an unprecedented resolution, revealing that prior models of melanoma evolution, in which pathways are turned ‘on’ or ‘off’ in a binary fashion during the run-up to melanoma, are oversimplified.

Summary: 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. : Botton et al. shed light on the heterogeneity of BRAF fusions encountered in melanocytic tumors and characterize features influencing their signaling and drug response. These findings unveil the singularities of BRAF fusions and establish general principles to guide their clinical management in melanoma and other malignancies. Keywords: BRAF fusion, melanoma, paradoxical activation, RAF inhibitor, MEK inhibitor, sequencing, rearrangement, translocation, kinase, pre-clinical

Combined hepatocellular-cholangiocarcinomas (CHC) are mixed tumours with both hepatocellular carcinoma (HCC) and cholangiocarcinoma (CC) components. CHC prognosis is similar to intrahepatic CC (ICC) and worse than HCC; staging and treatment generally follow ICC algorithms. However, the molecular biology of CHC remains poorly characterised. We performed capture-based next-generation sequencing of 20 CHC and, for comparison, 10 ICC arising in cirrhosis. Intratumour heterogeneity was assessed by separately sequencing the HCC and CC components of nine CHC. CHC demonstrated molecular profiles similar to HCC, even in the CC component. CHC harboured recurrent alterations in TERT (80%), TP53 (80%), cell cycle genes (40%; CCND1, CCNE1, CDKN2A), receptor tyrosine kinase/Ras/PI3-kinase pathway genes (55%; MET, ERBB2, KRAS, PTEN), chromatin regulators (20%; ARID1A, ARID2) and Wnt pathway genes (20%; CTNNB1, AXIN, APC). No CHC had alterations in IDH1, IDH2, FGFR2 or BAP1, genes typically mutated in ICC. TERT promoter mutations were consistently identified in both HCC and CC components, supporting TERT alteration as an early event in CHC evolution. TP53 mutations were present in both components in slightly over half the TP53-altered cases. By contrast, focal amplifications of CCND1, MET and ERRB2, as well as Wnt pathway alterations, were most often exclusive to one component, suggesting that these are late events in CHC evolution. ICC in cirrhosis demonstrated alterations similar to ICC in non-cirrhotic liver, including in IDH1 or IDH2 (30%), CDKN2A (40%), FGFR2 (20%), PBRM1 (20%), ARID1A (10%) and BAP1 (10%). TERT promoter and TP53 mutation were present in only one ICC each. Our data demonstrate that CHC genetics are distinct from ICC (even in cirrhosis) and similar to HCC, which has diagnostic utility and implications for treatment. Copyright © 2019 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

A single, neonatal UVB exposure drives greater UV-signature production in Braf- than Nras-mutant murine melanomas., BRAF-mutant melanomas are more likely than NRAS-mutant melanomas to arise in anatomical locations protected from chronic sun damage. We hypothesized that this discrepancy in tumor location is a consequence of the differential sensitivity of BRAF and NRAS-mutant melanocytes to ultraviolet light (UV)-mediated carcinogenesis. We tested this hypothesis by comparing the mutagenic consequences of a single neonatal, ultraviolet-AI (UVA; 340–400 nm) or ultraviolet-B (UVB; 280–390 nm) exposure in mouse models heterozygous for mutant Braf or homozygous for mutant Nras. Tumor onset was accelerated by UVB, but not UVA, and the resulting melanomas contained recurrent mutations affecting the RING domain of MAP3K1 and Actin-binding domain of Filamin A. Melanomas from UVB-irradiated, Braf-mutant mice averaged twice as many single-nucleotide variants and five times as many dipyrimidine variants than tumors from similarly irradiated Nras-mutant mice. A mutational signature discovered in UVB-accelerated tumors mirrored COSMIC signatures associated with human skin cancer and was more prominent in Braf- than Nras-mutant murine melanomas. These data show that a single UVB exposure yields a greater burden of mutations in murine tumors driven by oncogenic Braf.

Background Acral melanoma is a rare type of melanoma that affects world populations irrespective of skin color and has worse survival than other cutaneous melanomas. It has relatively few single nucleotide mutations without the UV signature of cutaneous melanomas, but instead has a genetic landscape characterized by structural rearrangements and amplifications. BRAF mutations are less common than in other cutaneous melanomas, and knowledge about alternative therapeutic targets is incomplete. Methods To identify alternative therapeutic targets, we performed targeted deep-sequencing on 122 acral melanomas. We confirmed the loss of the tumor suppressors p16 and NF1 by immunohistochemistry in select cases. Results In addition to BRAF (21.3%), NRAS (27.9%), and KIT (11.5%) mutations, we identified a broad array of MAPK pathway activating alterations, including fusions of BRAF (2.5%), NTRK3 (2.5%), ALK (0.8%), and PRKCA (0.8%), which can be targeted by available inhibitors. Inactivation of NF1 occurred in 18 cases (14.8%). Inactivation of the NF1 cooperating factor SPRED1 occurred in eight cases (6.6%) as an alternative mechanism of disrupting the negative regulation of RAS. Amplifications recurrently affected narrow loci containing PAK1 and GAB2 (n = 27, 22.1%), CDK4 (n = 27, 22.1%), CCND1 (n = 24, 19.7%), EP300 (n = 20, 16.4%), YAP1 (n = 15, 12.3%), MDM2 (n = 13, 10.7%), and TERT (n = 13, 10.7%) providing additional and possibly complementary therapeutic targets. Acral melanomas with BRAFV600E mutations harbored fewer genomic amplifications and were more common in patients with European ancestry. Conclusion Our findings support a new, molecularly based subclassification of acral melanoma with potential therapeutic implications: BRAFV600E mutant acral melanomas with characteristics similar to nonacral melanomas that could benefit from BRAF inhibitor therapy, and non-BRAFV600E mutant acral melanomas. Acral melanomas without BRAFV600E mutations harbor a broad array of therapeutically relevant alterations. Expanded molecular profiling would increase the detection of potentially targetable alterations for this subtype of acral melanoma.

Every cell in the human body has a unique set of somatic mutations, yet it remains difficult to comprehensively genotype an individual cell. Here, we developed solutions to overcome this obstacle in the context of normal human skin, thus offering the first glimpse into the genomic landscapes of individual melanocytes from human skin. We comprehensively genotyped 133 melanocytes from 19 sites across 6 donors. As expected, sun-shielded melanocytes had fewer mutations than sun-exposed melanocytes. However, within sun-exposed sites, melanocytes on chronically sun-exposed skin (e.g. the face) displayed a lower mutation burden than melanocytes on intermittently sun-exposed skin (e.g. the back). Melanocytes located adjacent to a skin cancer had higher mutation burdens than melanocytes from donors without skin cancer, implying that the mutation burden of normal skin can be harnessed to measure cumulative sun damage and skin cancer risk. Moreover, melanocytes from healthy skin commonly harbor pathogenic mutations, likely explaining the origins of the melanomas that arise in the absence of a pre-existing nevus. Phylogenetic analyses identified groups of related melanocytes, suggesting that melanocytes spread throughout skin as fields of clonally related cells, invisible to the naked eye. Overall, our study offers an unprecedented view into the genomic landscapes of individual melanocytes, revealing key insights into the causes and origins of melanoma.

The deubiquitinating enzyme BAP1 is mutated in a hereditary cancer syndrome with a high risk of mesothelioma and melanocytic tumors. Here, we show that Bap1 deletion in melanocytes cooperates with the constitutively active, oncogenic form of BRAF (BRAFV600E ) and UV to cause melanoma in mice, albeit at very low frequency. In addition, Bap1-null melanoma cells derived from mouse tumors are more aggressive and colonize and grow at distant sites more than their wild-type counterparts. Molecularly, Bap1-null melanoma cell lines have increased DNA damage measured by γH2aX and hyperubiquitination of histone H2a. Therapeutically, these Bap1-null tumors are completely responsive to BRAF- and MEK-targeted therapies. Therefore, BAP1 functions as a tumor suppressor and limits tumor progression in melanoma.

Two-step role for mutant TERT promoters Telomeres preserve genomic stability by preventing chromosomal fusions. The recent discovery that human tumors harbor mutations in the promoter region of the telomerase gene ( TERT ) produced a flurry of research aimed at elucidating the role of these mutations in cancer development. Chiba et al. present data that reconcile many of the conflicting results reported to date (see the Perspective by Shay). In human melanoma samples and a fibroblast model, TERT promoter mutations acted in two steps. First, the mutations allowed cells to proliferate with short telomeres. This fueled genomic instability and up-regulation of telomerase expression, leading to uncontrolled cell proliferation. Science , this issue p. 1416 ; see also p. 1358

Epigenetic landscapes can shape physiologic and disease phenotypes. We used integrative, high resolution multi-omics methods to characterize the oncogenic drivers of esophageal squamous cell carcinoma (ESCC). We found 98% of CpGs are hypomethylated across the ESCC genome and two-thirds occur in long non-coding (lnc)RNA regions. DNA methylation and epigenetic heterogeneity both coincide with chromosomal topological alterations. Gene body methylation, polycomb repressive complex occupancy, and CTCF binding sites associate with cancer-specific gene regulation. Epigenetically-mediated activation of non-canonical WNT signaling and the lncRNAESCCAL-1were validated as potential ESCC driver alterations. Gene-specific cancer driver roles of epigenetic alterations and heterogeneity are identified.

The use of microRNAs as biomarkers has been proposed for many diseases, including the diagnosis of melanoma. Although hundreds of microRNAs have been identified as differentially expressed in melanomas as compared to benign melanocytic lesions, a limited consensus has been achieved across studies, constraining the effective use of these potentially useful markers. In this study, we applied a machine learning-based pipeline to a dataset consisting of genetic features, clinical features, and next-generation microRNA sequencing from micro-dissected formalin-fixed paraffin embedded melanomas and their adjacent benign precursor nevi. We identified patient age and tumor cellularity as variables that frequently confound the measured expression of potentially diagnostic microRNAs. By employing the ratios of microRNAs that were either enriched or depleted in melanoma compared to the nevi as a normalization strategy, we developed a model that classified all the available published cohorts with an area under the receiver operating characteristic curve of 0.98. External validation on an independent cohort classified lesions with 81% sensitivity and 88% specificity and was uninfluenced by the tumor content of the sample or patient age.

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.

The use of microRNAs as biomarkers has been proposed for many diseases including the diagnosis of melanoma. Although hundreds of microRNAs have been identified as differentially expressed in melanomas as compared to benign melanocytic lesions, limited consensus has been achieved across studies, constraining the effective use of these potentially useful markers. In this study we quantified microRNAs by next-generation sequencing from melanomas and their adjacent benign precursor nevi. We applied a machine learning-based pipeline to identify a microRNA signature that separated melanomas from nevi and was unaffected by confounding variables, such as patient age and tumor cell content. By employing the ratios of microRNAs that were either enriched or depleted in melanoma compared to nevi as a normalization strategy, the classifier performed similarly across multiple published microRNA datasets, obtained by microarray, small RNA sequencing, or RT-qPCR. Validation on separate cohorts of melanomas and nevi correctly classified lesions with 83% sensitivity and 71-83% specificity, independent of variation in tumor cell content of the sample or patient age.

Uveal melanoma is a clinically distinct and particularly lethal subtype of melanoma originating from melanocytes in the eye. Here, we performed multi-region DNA sequencing of primary uveal melanomas and their matched metastases from 35 patients. We observed novel driver mutations and established the order in which these and known driver mutations undergo selection. Metastases had genomic alterations distinct from their primary tumors, and metastatic dissemination sometimes occurred early during the development of the primary tumor. Our study offers new insights into the genetics and evolution of this melanoma subtype, providing potential biomarkers for progression and therapy., Editorial summary: Multi-region sequencing of 35 primary uveal melanomas and their matched metastases yields new insights into the genetics and evolution of these tumors and provides potential biomarkers for progression and therapy.

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.

We elucidated genomic and transcriptomic changes that accompany the evolution of melanoma from pre-malignant lesions by sequencing DNA and RNA from primary melanomas and their adjacent precursors, as well as matched primary tumors and regional metastases. In total, we analyzed 230 histopathologically distinct areas of melanocytic neoplasia from 82 patients. Somatic alterations sequentially induced mitogen-activated protein kinase (MAPK) pathway activation, upregulation of telomerase, modulation of the chromatin landscape, G1/S checkpoint override, ramp-up of MAPK signaling, disruption of the p53 pathway, and activation of the PI3K pathway; no mutations were specifically associated with metastatic progression, as these pathways were perturbed during the evolution of primary melanomas. UV radiation-induced point mutations steadily increased until melanoma invasion, at which point copy-number alterations also became prevalent. Refereed/Peer-reviewed

Deep penetrating nevus (DPN) is characterized by enlarged, pigmented melanocytes that extend through the dermis. DPN can be difficult to distinguish from melanoma but rarely displays aggressive biological behavior. Here, we identify a combination of mutations of the β-catenin and mitogen-activated protein kinase pathways as characteristic of DPN. Mutations of the β-catenin pathway change the phenotype of a common nevus with BRAF mutation into that of DPN, with increased pigmentation, cell volume and nuclear cyclin D1 levels. Our results suggest that constitutive β-catenin pathway activation promotes tumorigenesis by overriding dependencies on the microenvironment that constrain proliferation of common nevi. In melanoma that arose from DPN we find additional oncogenic alterations. We identify DPN as an intermediate stage in the step-wise progression from nevus to melanoma. In summary, we delineate specific genetic alterations and their sequential order, information that can assist in the diagnostic classification and grading of these distinctive neoplasms., Deep penetrating nevi (DPN) are unusual melanocytic neoplasms with unknown genetic drivers. Here the authors show that majority of DPN harbor activating mutations in the β-catenin and the MAP-kinase pathways; this characteristic can help in the classification and grading of these distinctive neoplasms.

Melanomas on sun-exposed skin are heterogeneous tumours, which can be subtyped on the basis of their cumulative levels of exposure to ultraviolet (UV) radiation. A melanocytic neoplasm can also be staged by how far it has progressed, ranging from a benign neoplasm, such as a naevus, to a malignant neoplasm, such as a metastatic melanoma. Each subtype of melanoma can evolve through distinct evolutionary trajectories, passing through (or sometimes skipping over) various stages of transformation. This Review delineates several of the more common progression trajectories that occur in the patient setting and proposes models for tumour evolution that integrate genetic, histopathological, clinical and biological insights from the melanoma literature.

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.

The pathogenic mutations in melanoma have been largely catalogued; however, the order of their occurrence is not known.We sequenced 293 cancer-relevant genes in 150 areas of 37 primary melanomas and their adjacent precursor lesions. The histopathological spectrum of these areas included unequivocally benign lesions, intermediate lesions, and intraepidermal or invasive melanomas.Precursor lesions were initiated by mutations of genes that are known to activate the mitogen-activated protein kinase pathway. Unequivocally benign lesions harbored BRAF V600E mutations exclusively, whereas those categorized as intermediate were enriched for NRAS mutations and additional driver mutations. A total of 77% of areas of intermediate lesions and melanomas in situ harbored TERT promoter mutations, a finding that indicates that these mutations are selected at an unexpectedly early stage of the neoplastic progression. Biallelic inactivation of CDKN2A emerged exclusively in invasive melanomas. PTEN and TP53 mutations were found only in advanced primary melanomas. The point-mutation burden increased from benign through intermediate lesions to melanoma, with a strong signature of the effects of ultraviolet radiation detectable at all evolutionary stages. Copy-number alterations became prevalent only in invasive melanomas. Tumor heterogeneity became apparent in the form of genetically distinct subpopulations as melanomas progressed.Our study defined the succession of genetic alterations during melanoma progression, showing distinct evolutionary trajectories for different melanoma subtypes. It identified an intermediate category of melanocytic neoplasia, characterized by the presence of more than one pathogenic genetic alteration and distinctive histopathological features. Finally, our study implicated ultraviolet radiation as a major factor in both the initiation and progression of melanoma. (Funded by the National Institutes of Health and others.).

Melanoma is a neoplastic disorder that frequently arises on the sun-exposed skin of light-skinned individuals. Melanomas can also occur on anatomic sites that are partially or completely shielded from sunlight, and these subtypes affect world populations at similar incidences, independent of skin color. These sun-shielded melanomas also differ in their clinical and histologic presentations as well as in their patterns of driver mutations and mutational mechanisms. These notable differences have led to classifications distinguishing melanoma subtypes, originally based on clinical and histopathological features1, and later integrating these features with genetic alterations2. This article is protected by copyright. All rights reserved.

Germline copy number variants (CNVs) and somatic copy number alterations (SCNAs) are of significant importance in syndromic conditions and cancer. Massive 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 compatibility with other software. Availability: http://github.com/etal/cnvkit

Background Pancreatic cancer is a deadly disease with a five-year survival of less than 5%. A better understanding of the underlying biology may suggest novel therapeutic targets. Recent surveys of the pancreatic cancer genome have uncovered numerous new alterations; yet systematic functional characterization of candidate cancer genes has lagged behind. To address this challenge, here we have devised a highly-parallel RNA interference-based functional screen to evaluate many genomically-nominated candidate pancreatic cancer genes simultaneously. Results For 185 candidate pancreatic cancer genes, selected from recurrently altered genomic loci, we performed a pooled shRNA library screen of cell growth/viability across 10 different cell lines. Knockdown-associated effects on cell growth were assessed by enrichment or depletion of shRNA hairpins, by hybridization to barcode microarrays. A novel analytical approach (COrrelated Phenotypes for On-Target Effects; COPOTE) was used to discern probable on-target knockdown, based on identifying different shRNAs targeting the same gene and displaying concordant phenotypes across cell lines. Knockdown data were integrated with genomic architecture and gene-expression profiles, and selected findings validated using individual shRNAs and/or independent siRNAs. The pooled shRNA library design delivered reproducible data. In all, COPOTE analysis identified 52 probable on-target gene-knockdowns. Knockdown of known oncogenes (KRAS, MYC, SMURF1 and CCNE1) and a tumor suppressor (CDKN2A) showed the expected contrasting effects on cell growth. In addition, the screen corroborated purported roles of PLEKHG2 and MED29 as 19q13 amplicon drivers. Most notably, the analysis also revealed novel possible oncogenic functions of nucleoporin NUP153 (ostensibly by modulating TGFβ signaling) and Kruppel-like transcription factor KLF5 in pancreatic cancer. Conclusions By integrating physical and functional genomic data, we were able to simultaneously evaluate many candidate pancreatic cancer genes. Our findings uncover new facets of pancreatic cancer biology, with possible therapeutic implications. More broadly, our study provides a general strategy for the efficient characterization of candidate genes emerging from cancer genome studies.

Gene fusions, like BCR/ABL1 in chronic myelogenous leukemia, have long been recognized in hematologic and mesenchymal malignancies. The recent finding of gene fusions in prostate and lung cancers has motivated the search for pathogenic gene fusions in other malignancies. Here, we developed a “breakpoint analysis” pipeline to discover candidate gene fusions by tell-tale transcript level or genomic DNA copy number transitions occurring within genes. Mining data from 974 diverse cancer samples, we identified 198 candidate fusions involving annotated cancer genes. From these, we validated and further characterized novel gene fusions involving ROS1 tyrosine kinase in angiosarcoma (CEP85L/ROS1), SLC1A2 glutamate transporter in colon cancer (APIP/SLC1A2), RAF1 kinase in pancreatic cancer (ATG7/RAF1) and anaplastic astrocytoma (BCL6/RAF1), EWSR1 in melanoma (EWSR1/CREM), CDK6 kinase in T-cell acute lymphoblastic leukemia (FAM133B/CDK6), and CLTC in breast cancer (CLTC/VMP1). Notably, while these fusions involved known cancer genes, all occurred with novel fusion partners and in previously unreported cancer types. Moreover, several constituted druggable targets (including kinases), with therapeutic implications for their respective malignancies. Lastly, breakpoint analysis identified new cell line models for known rearrangements, including EGFRvIII and FIP1L1/PDGFRA. Taken together, we provide a robust approach for gene fusion discovery, and our results highlight a more widespread role of fusion genes in cancer pathogenesis., Author Summary Gene fusions represent an important class of cancer genes, created by rearrangements of the genome that bring together two different genes. Because they are unique to cancer cells, gene fusions are ideal diagnostic markers and therapeutic targets. While gene fusions were once thought restricted mainly to blood cancers, recent discoveries suggest they are more widespread. Here, we have developed an approach for mining DNA microarray data to detect the tell-tale signatures of gene fusions, as “breakpoints” occurring within the encoding DNA or expressed transcripts. We apply this approach to a large collection of nearly 1,000 human cancer specimens. From this analysis, we discover and verify twelve new gene fusions occurring in diverse cancer types. We verify that some of these rearrangements recur in other samples of the same cancer type (supporting a causal role) and that the cancers show dependency on the fusion for cancer cell growth. Notably, some of these fusions (e.g. CEP85L/ROS1 in angiosarcoma) represent the first for that cancer type and thus provide important new biological insight. Some are also good drug targets (including rearrangements of ROS1, RAF1, and CDK6 kinases), with clear implications for therapy.

Pancreatic cancer is a deadly disease, and new therapeutic targets are urgently needed. We previously identified DNA amplification at 7q21-q22 in pancreatic cancer cell lines. Now, by high-resolution genomic profiling of human pancreatic cancer cell lines and human tumors (engrafted in immunodeficient mice to enrich the cancer epithelial fraction), we define a 325 Kb minimal amplicon spanning SMURF1, an E3 ubiquitin ligase and known negative regulator of transforming growth factor β (TGFβ) growth inhibitory signaling. SMURF1 amplification was confirmed in primary human pancreatic cancers by fluorescence in situ hybridization (FISH), where 4 of 95 cases (4.2%) exhibited amplification. By RNA interference (RNAi), knockdown of SMURF1 in a human pancreatic cancer line with focal amplification (AsPC-1) did not alter cell growth, but led to reduced cell invasion and anchorage-independent growth. Interestingly, this effect was not mediated through altered TGFβ signaling, assayed by transcriptional reporter. Finally, overexpression of SMURF1 (but not a catalytic mutant) led to loss of contact inhibition in NIH-3T3 mouse embryo fibroblast cells. Together, these findings identify SMURF1 as an amplified oncogene driving multiple tumorigenic phenotypes in pancreatic cancer, and provide a new druggable target for molecularly directed therapy.

SWI/SNF is a multi-subunit chromatin remodeling complex that uses the energy of ATP hydrolysis to reposition nucleosomes, thereby modulating gene expression. Accumulating evidence suggests that SWI/SNF functions as a tumor suppressor in some cancers. However, the spectrum of SWI/SNF mutations across human cancers has not been systematically investigated. Here, we mined whole-exome sequencing data from 24 published studies representing 669 cases from 18 neoplastic diagnoses. SWI/SNF mutations were widespread across diverse human cancers, with an excess of deleterious mutations, and an overall frequency approaching TP53 mutation. Mutations occurred most commonly in the SMARCA4 enzymatic subunit, and in subunits thought to confer functional specificity (ARID1A, ARID1B, PBRM1, and ARID2). SWI/SNF mutations were not mutually-exclusive of other mutated cancer genes, including TP53 and EZH2 (both previously linked to SWI/SNF). Our findings implicate SWI/SNF as an important but under-recognized tumor suppressor in diverse human cancers, and provide a key resource to guide future investigations.