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.

Well-differentiated papillary mesothelioma is an uncommon mesothelial neoplasm that most frequently arises in the peritoneal cavity of women of reproductive age. Whereas malignant mesothelioma is an aggressive tumor associated with poor outcome, well-differentiated papillary mesothelioma typically exhibits indolent behavior. However, histologically differentiating between these two entities can be challenging, necessitating the development of distinguishing biomarkers. While the genetic alterations that drive malignant mesothelioma have recently been determined, the molecular pathogenesis of well-differentiated papillary mesothelioma is unknown. Here we performed genomic profiling on a cohort of ten well-differentiated papillary mesothelioma of the peritoneum. We identified that all tumors harbored somatic missense mutations in either the TRAF7 or CDC42 genes, and lacked alterations involving BAP1, NF2, CDKN2A, DDX3X, SETD2, and ALK that are frequent in malignant mesothelioma. We recently identified that another mesothelial neoplasm, adenomatoid tumor of the genital tract, is genetically defined by somatic missense mutations in the TRAF7 gene, indicating a shared molecular pathogenesis between well-differentiated papillary mesothelioma and adenomatoid tumors. To the best of our knowledge, well-differentiated papillary mesothelioma is the first human tumor type found to harbor recurrent mutations in the CDC42 gene, which encodes a Rho family GTPase. Immunohistochemistry demonstrated intact BAP1 expression in all cases of well-differentiated papillary mesothelioma, indicating that this is a reliable marker for distinguishing well-differentiated papillary mesothelioma from malignant mesotheliomas that frequently display loss of expression. Additionally, all well-differentiated papillary mesothelioma demonstrated robust expression of L1 cell adhesion molecule (L1CAM), a marker of NF-kB pathway activation, similar to that observed in adenomatoid tumors. In contrast, we have previously shown that L1CAM staining is not observed in normal mesothelial cells and malignant mesotheliomas of the peritoneum. Together, these studies demonstrate that well-differentiated papillary mesothelioma is genetically defined by mutually exclusive mutations in TRAF7 and CDC42 that molecularly distinguish this entity from malignant mesothelioma.

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.

Metagenomic next-generation sequencing (mNGS) of body fluids is an emerging approach to identify occult pathogens in undiagnosed patients. We hypothesized that metagenomic testing can be simultaneously used to detect malignant neoplasms in addition to infectious pathogens.

Knowledge of key drivers and therapeutic targets in mucosal melanoma is limited due to the paucity of comprehensive mutation data on this rare tumor type. To better understand the genomic landscape of mucosal melanoma, here we describe whole genome sequencing analysis of 67 tumors and validation of driver gene mutations by exome sequencing of 45 tumors. Tumors have a low point mutation burden and high numbers of structural variants, including recurrent structural rearrangements targeting TERT, CDK4 and MDM2. Significantly mutated genes are NRAS, BRAF, NF1, KIT, SF3B1, TP53, SPRED1, ATRX, HLA-A and CHD8. SF3B1 mutations occur more commonly in female genital and anorectal melanomas and CTNNB1 mutations implicate a role for WNT signaling defects in the genesis of some mucosal melanomas. TERT aberrations and ATRX mutations are associated with alterations in telomere length. Mutation profiles of the majority of mucosal melanomas suggest potential susceptibility to CDK4/6 and/or MEK inhibitors.

A subset of Spitz tumors harbor fusions of NTRK3 with ETV6, MYO5A, and MYH9. We evaluated a series of 22 melanocytic tumors in which an NTRK3 fusion was identified as part of the diagnostic workup. Tumors in which NTRK3 was fused to ETV6 occurred in younger patients were predominantly composed of epithelioid melanocytes and were classified by their histopathologic features as Spitz tumors. In contrast, those in which NTRK3 was fused to MYO5A were predominantly composed of spindled melanocytes arrayed in fascicles with neuroid features such as pseudo-Verocay bodies. To further investigate the effects of the fusion kinases ETV6-NTRK3 and MYO5A-NTRK3 in melanocytes, we expressed them in immortalized melanocytes and determined their subcellular localization by immunofluorescence. ETV6-NTRK3 was localized to the nucleus and diffusely within the cytoplasm and caused melanocytes to adopt an epithelioid cytomorphology. In contrast, MYO5A-NTRK3, appeared excluded from the nucleus of melanocytes, was localized to dendrites, and resulted in a highly dendritic cytomorphology. Our findings indicate that ETV6-NTRK3 and MYO5A-NTRK3 have distinct subcellular localizations and effects on cellular morphology.

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.

KRAS mutations drive a quarter of cancer mortality, and most are undruggable. Several inhibitors of the MAPK pathway are FDA approved but poorly tolerated at the doses needed to adequately extinguish RAS/RAF/MAPK signaling in the tumor cell. We found that oncogenic KRAS signaling induced ferrous iron (Fe2+) accumulation early in and throughout mutant KRAS-mediated transformation. We converted an FDA-approved MEK inhibitor into a ferrous iron-activatable drug conjugate (FeADC) and achieved potent MAPK blockade in tumor cells while sparing normal tissues. This innovation allowed sustainable, effective treatment of tumor-bearing animals, with tumor-selective drug activation, producing superior systemic tolerability. Ferrous iron accumulation is an exploitable feature of KRAS transformation, and FeADCs hold promise for improving the treatment of KRAS-driven solid tumors.

Melanomas that have histopathologic features that overlap with those of Spitz nevus are referred to as spitzoid melanomas. However, the diagnostic concept is used inconsistently and genomic analyses suggest it is a heterogeneous category. Spitz tumors, the spectrum of melanocytic neoplasms extending from Spitz nevi to their malignant counterpart Spitz melanoma, are defined in the 2018 WHO classification of skin tumors by the presence of specific genetic alterations, such as kinase fusions or HRAS mutations. It is unclear what fraction of "spitzoid melanomas" defined solely by their histopathologic features belong to the category of Spitz melanoma or to other melanoma subtypes. We assembled a cohort of 25 spitzoid melanomas diagnosed at a single institution over an 8-year period and performed high-coverage DNA sequencing of 480 cancer related genes. Transcriptome wide RNA sequencing was performed for select cases. Only nine cases (36%) had genetic alterations characteristic of Spitz melanoma, including HRAS mutation or fusion involving BRAF, ALK, NTRK1, or MAP3K8. The remaining cases were divided into those with an MAPK activating mutation and those without an MAPK activating mutation. Both Spitz melanoma and spitzoid melanomas in which an MAPK-activating mutation could not be identified tended to occur in younger patients on skin with little solar elastosis, infrequently harbored TERT promoter mutations, and had a lower burden of pathogenic mutations than spitzoid melanomas with non-Spitz MAPK-activating mutations. The MAPK-activating mutations identified affected non-V600 residues of BRAF as well as NRAS, MAP2K1/2, NF1, and KIT, while BRAF V600 mutations, the most common mutations in melanomas of the WHO low-CSD category, were entirely absent. While the "spitzoid melanomas" comprising our cohort were enriched for bona fide Spitz melanomas, the majority of melanomas fell outside of the genetically defined category of Spitz melanomas, indicating that histomorphology is an unreliable predictor of Spitz lineage.

Mucosal melanoma is a rare and poorly characterized subtype of human melanoma. Here we perform a cross-species analysis by sequencing tumor-germline pairs from 46 primary human muscosal, 65 primary canine oral and 28 primary equine melanoma cases from mucosal sites. Analysis of these data reveals recurrently mutated driver genes shared between species such as NRAS, FAT4, PTPRJ, TP53 and PTEN, and pathogenic germline alleles of BRCA1, BRCA2 and TP53. We identify a UV mutation signature in a small number of samples, including human cases from the lip and nasal mucosa. A cross-species comparative analysis of recurrent copy number alterations identifies several candidate drivers including MDM2, B2M, KNSTRN and BUB1B. Comparison of somatic mutations in recurrences and metastases to those in the primary tumor suggests pervasive intra-tumor heterogeneity. Collectively, these studies suggest a convergence of some genetic changes in mucosal melanomas between species but also distinctly different paths to tumorigenesis.

Adenomatoid tumors are the most common neoplasm of the epididymis, and histologically similar adenomatoid tumors also commonly arise in the uterus and fallopian tube. To investigate the molecular pathogenesis of these tumors, we performed genomic profiling on a cohort of 31 adenomatoid tumors of the male and female genital tracts. We identified that all tumors harbored somatic missense mutations in the TRAF7 gene, which encodes an E3 ubiquitin ligase belonging to the family of tumor necrosis factor receptor-associated factors (TRAFs). These mutations all clustered into one of five recurrent hotspots within the WD40 repeat domains at the C-terminus of the protein. Functional studies in vitro revealed that expression of mutant but not wild-type TRAF7 led to increased phosphorylation of nuclear factor-kappa B (NF-kB) and increased expression of L1 cell adhesion molecule (L1CAM), a marker of NF-kB pathway activation. Immunohistochemistry demonstrated robust L1CAM expression in adenomatoid tumors that was absent in normal mesothelial cells, malignant peritoneal mesotheliomas and multilocular peritoneal inclusion cysts. Together, these studies demonstrate that adenomatoid tumors of the male and female genital tract are genetically defined by TRAF7 mutation that drives aberrant NF-kB pathway activation.

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.