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.

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 previously unknown driver mutations and established the order in which these and known driver mutations undergo selection. Metastases had genomic alterations distinct from their primary tumors; 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.

Every cell in the human body has a unique set of somatic mutations, but it remains difficult to comprehensively genotype an individual cell1. Here we describe ways to overcome this obstacle in the context of normal human skin, thus offering a glimpse into the genomic landscapes of individual melanocytes from human skin. As expected, sun-shielded melanocytes had fewer mutations than sun-exposed melanocytes. However, melanocytes from chronically sun-exposed skin (for example, the face) had a lower mutation burden than melanocytes from intermittently sun-exposed skin (for example, 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 used to measure cumulative sun damage and risk of skin cancer. Moreover, melanocytes from healthy skin commonly contained pathogenic mutations, although these mutations tended to be weakly oncogenic, probably 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 that are invisible to the naked eye. Overall, our results uncover the genomic landscapes of individual melanocytes, providing key insights into the causes and origins of melanoma.

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.

Desmoplastic melanoma is an uncommon variant of melanoma with sarcomatous histology, distinct clinical behavior and unknown pathogenesis. We performed low-coverage genome and high-coverage exome sequencing of 20 desmoplastic melanomas, followed by targeted sequencing of 293 genes in a validation cohort of 42 cases. A high mutation burden (median of 62 mutations/Mb) ranked desmoplastic melanoma among the most highly mutated cancers. Mutation patterns strongly implicate ultraviolet radiation as the dominant mutagen, indicating a superficially located cell of origin. Newly identified alterations included recurrent promoter mutations of NFKBIE, encoding NF-κB inhibitor ɛ (IκBɛ), in 14.5% of samples. Common oncogenic mutations in melanomas, in particular in BRAF (encoding p.Val600Glu) and NRAS (encoding p.Gln61Lys or p.Gln61Arg), were absent. Instead, other genetic alterations known to activate the MAPK and PI3K signaling cascades were identified in 73% of samples, affecting NF1, CBL, ERBB2, MAP2K1, MAP3K1, BRAF, EGFR, PTPN11, MET, RAC1, SOS2, NRAS and PIK3CA, some of which are candidates for targeted therapies.

Soft-tissue sarcomas are a group of malignant tumours whose clinical management is complicated by morphological heterogeneity, inadequate molecular markers and limited therapeutic options. Receptor tyrosine kinases (RTKs) have been shown to play important roles in cancer, both as therapeutic targets and as prognostic biomarkers. An initial screen of gene expression data for 48 RTKs in 148 sarcomas showed that ROR2 was expressed in a subset of leiomyosarcoma (LMS), gastrointestinal stromal tumour (GIST) and desmoid-type fibromatosis (DTF). This was further confirmed by immunohistochemistry (IHC) on 573 tissue samples from 59 sarcoma tumour types. Here we provide evidence that ROR2 expression plays a role in the invasive abilities of LMS and GIST cells in vitro. We also show that knockdown of ROR2 significantly reduces tumour mass in vivo using a xenotransplantation model of LMS. Lastly, we show that ROR2 expression, as measured by IHC, predicts poor clinical outcome in patients with LMS and GIST, although it was not independent of other clinico-pathological features in a multivariate analysis, and that ROR2 expression is maintained between primary tumours and their metastases. Together, these results show that ROR2 is a useful prognostic indicator in the clinical management of these soft-tissue sarcomas and may represent a novel therapeutic target.

To better understand intrinsic resistance to immune checkpoint blockade (ICB), we established a comprehensive view of the cellular architecture of the treatment-naive melanoma ecosystem and studied its evolution under ICB. Using single-cell, spatial multi-omics, we showed that the tumor microenvironment promotes the emergence of a complex melanoma transcriptomic landscape. Melanoma cells harboring a mesenchymal-like (MES) state, a population known to confer resistance to targeted therapy, were significantly enriched in early on-treatment biopsies from non-responders to ICB. TCF4 serves as the hub of this landscape by being a master regulator of the MES signature and a suppressor of the melanocytic and antigen presentation transcriptional programs. Targeting TCF4 genetically or pharmacologically, using a bromodomain inhibitor, increased immunogenicity and sensitivity of MES cells to ICB and targeted therapy. We thereby uncovered a TCF4-dependent regulatory network that orchestrates multiple transcriptional programs and contributes to resistance to both targeted therapy and ICB in melanoma.

Acral melanoma is an aggressive type of melanoma with unknown origins, arising on the sole, palm, or nail apparatus. It is the most common type of melanoma in individuals with dark skin and is notoriously challenging to treat. Our study examined exome sequencing data from 139 tissue samples, spanning different progression stages, collected from 37 patients. We found that 78.4% of the melanomas displayed one or more clustered copy number transitions with focal amplifications, recurring predominantly on chromosomes 5, 11, 12, and 22. These genomic "hailstorms" were typically shared across all progression stages within individual patients. Genetic alterations known to activate TERT also arose early. By contrast, mutations in the MAP-kinase pathway appeared later during progression, often leading to different tumor areas harboring non-overlapping driver mutations. We conclude that the evolutionary trajectories of acral melanomas substantially diverge from those of melanomas on sun-exposed skin, where MAP-kinase pathway activation initiates the neoplastic cascade followed by immortalization later. The punctuated formation of hailstorms, paired with early TERT activation, suggests a unique mutational mechanism underlying the origins of acral melanoma. Our findings highlight an essential role for telomerase, likely in re-stabilizing tumor genomes after hailstorms have initiated the tumors. The marked genetic heterogeneity, in particular of MAP-kinase pathway drivers, may partly explain the limited success of targeted and other therapies in treating this melanoma subtype.

Oncogenic gene fusions have been identified in many cancers and many serve as biomarkers or targets for therapy. Here we identify six different melanocytic tumours with genomic rearrangements of MET fusing the kinase domain of MET in-frame to six different N-terminal partners. These tumours lack activating mutations in other established melanoma oncogenes. We functionally characterize two of the identified fusion proteins (TRIM4-MET and ZKSCAN1-MET) and find that they constitutively activate the mitogen-activated protein kinase (MAPK), phosphoinositol-3 kinase (PI3K) and phospholipase C gamma 1 (PLCγ1) pathways. The MET inhibitors cabozantinib (FDA-approved for progressive medullary thyroid cancer) and PF-04217903 block their activity at nanomolar concentrations. MET fusion kinases thus provide a potential therapeutic target for a rare subset of melanoma for which currently no targeted therapeutic options currently exist.

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.

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.

A major function of long non-coding RNAs (lncRNAs) is regulating gene expression through changes in chromatin state. Experimental evidence suggests that in cancer, they can influence Polycomb Repressive Complexes (PRC) to retarget to an occupancy pattern resembling that of the embryonic state. We have previously demonstrated that the expression level of lncRNA in the HOX locus, including HOTAIR, is a predictor of breast cancer metastasis. In this current project, RNA in situ hybridization of probes to three different lncRNAs (HOTAIR, ncHoxA1, and ncHoxD4), as well a immunohistochemical staining of EZH2, is undertaken in formalin-fixed paraffin-embedded breast cancer tissues in a high throughput tissue microarray format to correlate expression with clinicopathologic features. Though overall EZH2 and HOTAIR expression levels were highly correlated, the subset of cases with strong HOTAIR expression correlated with ER and PR positivity, while the subset of cases with strong EZH2 expression correlated with an increased proliferation rate, ER and PR negativity, HER2 underexpression, and triple negativity. Co-expression of HOTAIR and EZH2 trended with a worse outcome. In matched primary and metastatic cancers, both HOTAIR and EZH2 had increased expression in the metastatic carcinomas. This is the first study to show that RNA in situ hybridization of formalin fixed paraffin-embedded clinical material can be used to measure levels of long non-coding RNAs. This approach offers a method to make observations on lncRNAs that may influence the cancer epigenome in a tissue-based technique.

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.

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.

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.

Benign melanocytic nevi frequently emerge when an acquired BRAFV600E mutation triggers unchecked proliferation and subsequent arrest in melanocytes. Recent observations have challenged the role of oncogene-induced senescence in melanocytic nevus formation, necessitating investigations into alternative mechanisms for the establishment and maintenance of proliferation arrest in nevi. We compared the transcriptomes of melanocytes from healthy human skin, nevi, and melanomas arising from nevi and identified a set of microRNAs as highly expressed nevus-enriched transcripts. Two of these microRNAs-MIR211-5p and MIR328-3p-induced mitotic failure, genome duplication, and proliferation arrest in human melanocytes through convergent targeting of AURKB. We demonstrate that BRAFV600E induces a similar proliferation arrest in primary human melanocytes that is both reversible and conditional. Specifically, BRAFV600E expression stimulates either arrest or proliferation depending on the differentiation state of the melanocyte. We report genome duplication in human melanocytic nevi, reciprocal expression of AURKB and microRNAs in nevi and melanomas, and rescue of arrested human nevus cells with AURKB expression. Taken together, our data describe an alternative molecular mechanism for melanocytic nevus formation that is congruent with both experimental and clinical observations.

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 p16INK4A 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.

Pancreatic cancer, composed principally of pancreatic adenocarcinoma (PaC), is the fourth leading cause of cancer death in the United States. PaC-associated diabetes may be a marker of early disease. We sought to identify molecules associated with PaC and PaC with diabetes (PaC-DM) using a novel translational bioinformatics approach. We identified fatty acid binding protein-1 (FABP-1) as one of several candidates. The primary aim of this pilot study was to experimentally validate the predicted association between FABP-1 with PaC and PaC with diabetes.

Previously, we showed that the presence of high numbers of macrophages correlates with poor prognosis in nongynecological leiomyosarcoma (LMS). In gynecological LMS, a similar trend was noted but did not reach statistical significance. Colony-stimulating factor-1 (CSF1) is a major chemoattractant for macrophages. Here we show that in a subset of LMS cases, CSF1 is expressed by the malignant cells. Previously, we found that CSF1 is translocated and highly expressed in tenosynovial giant cell tumors (TGCTs), and this observation allowed us to identify genes that showed a coordinate expression with CSF1. Here, we evaluated the expression of CSF1 and TGCT-associated proteins in 149 cases of LMS. The coordinate expression of CSF1 and three TGCT-associated proteins (CD163, FCGR3a, and CTSL1) identified cases with poor prognosis in both gynecological LMS (P = 0.00006) and nongynecological LMS (P = 0.03). In gynecological LMS, the coordinate expression of these four markers was the only independent prognosticator in multivariate analysis (hazard ratio, 4.2; 95% CI, 1.12 to 16; P = 0.03). Our findings indicate that CSF1 may play an important role in the clinical behavior of LMS that may open a window for new therapeutic reagents.

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 catalog of driver genes in cutaneous squamous cell carcinoma, offering points of therapeutic intervention and insights into the biology of cutaneous squamous cell carcinoma.