Non-viral precision T cell receptor replacement for personalized cell therapy.

T cell receptors (TCRs) enable T cells to specifically recognize mutations in cancer cells1-3. Here we developed a clinical-grade approach based on CRISPR-Cas9 non-viral precision genome-editing to simultaneously knockout the two endogenous TCR genes TRAC (which encodes TCRα) and TRBC (which encodes TCRß). We also inserted into the TRAC locus two chains of a neoantigen-specific TCR (neoTCR) isolated from circulating T cells of patients. The neoTCRs were isolated using a personalized library of soluble predicted neoantigen-HLA capture reagents. Sixteen patients with different refractory solid cancers received up to three distinct neoTCR transgenic cell products. Each product expressed a patient-specific neoTCR and was administered in a cell-dose-escalation, first-in-human phase I clinical trial ( NCT03970382 ). One patient had grade 1 cytokine release syndrome and one patient had grade 3 encephalitis. All participants had the expected side effects from the lymphodepleting chemotherapy. Five patients had stable disease and the other eleven had disease progression as the best response on the therapy. neoTCR transgenic T cells were detected in tumour biopsy samples after infusion at frequencies higher than the native TCRs before infusion. This study demonstrates the feasibility of isolating and cloning multiple TCRs that recognize mutational neoantigens. Moreover, simultaneous knockout of the endogenous TCR and knock-in of neoTCRs using single-step, non-viral precision genome-editing are achieved. The manufacture of neoTCR engineered T cells at clinical grade, the safety of infusing up to three gene-edited neoTCR T cell products and the ability of the transgenic T cells to traffic to the tumours of patients are also demonstrated.

Genome-Wide Polygenic Score Predicts Large Number of High Risk Individuals in Monogenic Undiagnosed Young Onset Parkinson's Disease Patients from India.

Parkinson's disease (PD) is a genetically heterogeneous neurodegenerative disease with poorly defined environmental influences. Genomic studies of PD patients have identified disease-relevant monogenic genes, rare variants of significance, and polygenic risk-associated variants. In this study, whole genome sequencing data from 90 young onset Parkinson's disease (YOPD) individuals are analyzed for both monogenic and polygenic risk. The genetic variant analysis identifies pathogenic/likely pathogenic variants in eight of the 90 individuals (8.8%). It includes large homozygous coding exon deletions in PRKN and SNV/InDels in VPS13C, PLA2G6, PINK1, SYNJ1, and GCH1. Eleven rare heterozygous GBA coding variants are also identified in 13 (14.4%) individuals. In 34 (56.6%) individuals, one or more variants of uncertain significance (VUS) in PD/PD-relevant genes are observed. Though YOPD patients with a prioritized pathogenic variant show a low polygenic risk score (PRS), patients with prioritized VUS or no significant rare variants show an increased PRS odds ratio for PD. This study suggests that both significant rare variants and polygenic risk from common variants together may contribute to the genesis of PD. Further validation using a larger cohort of patients will confirm the interplay between monogenic and polygenic variants and their use in routine genetic PD diagnosis and risk assessment.

Distinct resistance mechanisms arise to allosteric vs. ATP-competitive AKT inhibitors.

The AKT kinases have emerged as promising therapeutic targets in oncology and both allosteric and ATP-competitive AKT inhibitors have entered clinical investigation. However, long-term efficacy of such inhibitors will likely be challenged by the development of resistance. We have established prostate cancer models of acquired resistance to the allosteric inhibitor MK-2206 or the ATP-competitive inhibitor ipatasertib following prolonged exposure. While alterations in AKT are associated with acquired resistance to MK-2206, ipatasertib resistance is driven by rewired compensatory activity of parallel signaling pathways. Importantly, MK-2206 resistance can be overcome by treatment with ipatasertib, while ipatasertib resistance can be reversed by co-treatment with inhibitors of pathways including PIM signaling. These findings demonstrate that distinct resistance mechanisms arise to the two classes of AKT inhibitors and that combination approaches may reverse resistance to ATP-competitive inhibition.

RTK-Dependent Inducible Degradation of Mutant PI3Kα Drives GDC-0077 (Inavolisib) Efficacy.

PIK3CA is one of the most frequently mutated oncogenes; the p110a protein it encodes plays a central role in tumor cell proliferation. Small-molecule inhibitors targeting the PI3K p110a catalytic subunit have entered clinical trials, with early-phase GDC-0077 studies showing antitumor activity and a manageable safety profile in patients with PIK3CA-mutant breast cancer. However, preclinical studies have shown that PI3K pathway inhibition releases negative feedback and activates receptor tyrosine kinase signaling, reengaging the pathway and attenuating drug activity. Here we discover that GDC-0077 and taselisib more potently inhibit mutant PI3K pathway signaling and cell viability through unique HER2-dependent mutant p110a degradation. Both are more effective than other PI3K inhibitors at maintaining prolonged pathway suppression. This study establishes a new strategy for identifying inhibitors that specifically target mutant tumors by selective degradation of the mutant oncoprotein and provide a strong rationale for pursuing PI3Kα degraders in patients with HER2-positive breast cancer. SIGNIFICANCE: The PI3K inhibitors GDC-0077 and taselisib have a unique mechanism of action; both inhibitors lead to degradation of mutant p110a protein. The inhibitors that have the ability to trigger specific degradation of mutant p110a without significant change in wild-type p110a protein may result in improved therapeutic index in PIK3CA-mutant tumors.See related commentary by Vanhaesebroeck et al., p. 20.This article is highlighted in the In This Issue feature, p. 1.

Human ACE2 receptor polymorphisms and altered susceptibility to SARS-CoV-2.

COVID-19 is a respiratory illness caused by a novel coronavirus called SARS-CoV-2. The viral spike (S) protein engages the human angiotensin-converting enzyme 2 (ACE2) receptor to invade host cells with ~10-15-fold higher affinity compared to SARS-CoV S-protein, making it highly infectious. Here, we assessed if ACE2 polymorphisms can alter host susceptibility to SARS-CoV-2 by affecting this interaction. We analyzed over 290,000 samples representing >400 population groups from public genomic datasets and identified multiple ACE2 protein-altering variants. Using reported structural data, we identified natural ACE2 variants that could potentially affect virus-host interaction and thereby alter host susceptibility. These include variants S19P, I21V, E23K, K26R, T27A, N64K, T92I, Q102P and H378R that were predicted to increase susceptibility, while variants K31R, N33I, H34R, E35K, E37K, D38V, Y50F, N51S, M62V, K68E, F72V, Y83H, G326E, G352V, D355N, Q388L and D509Y were predicted to be protective variants that show decreased binding to S-protein. Using biochemical assays, we confirmed that K31R and E37K had decreased affinity, and K26R and T92I variants showed increased affinity for S-protein when compared to wildtype ACE2. Consistent with this, soluble ACE2 K26R and T92I were more effective in blocking entry of S-protein pseudotyped virus suggesting that ACE2 variants can modulate susceptibility to SARS-CoV-2.

Elevated Serum Amino Acids Induce a Subpopulation of Alpha Cells to Initiate Pancreatic Neuroendocrine Tumor Formation.

The cellular origin of sporadic pancreatic neuroendocrine tumors (PNETs) is obscure. Hormone expression suggests that these tumors arise from glucagon-producing alpha cells or insulin-producing ß cells, but instability in hormone expression prevents linage determination. We utilize loss of hepatic glucagon receptor (GCGR) signaling to drive alpha cell hyperproliferation and tumor formation to identify a cell of origin and dissect mechanisms that drive progression. Using a combination of genetically engineered Gcgr knockout mice and GCGR-inhibiting antibodies, we show that elevated plasma amino acids drive the appearance of a proliferative population of SLC38A5+ embryonic progenitor-like alpha cells in mice. Further, we characterize tumors from patients with rare bi-allelic germline GCGR loss-of-function variants and find prominent tumor-cell-associated expression of the SLC38A5 paralog SLC7A8 as well as markers of active mTOR signaling. Thus, progenitor cells arise from adult alpha cells in response to metabolic signals and, when inductive signals are chronically present, drive tumor initiation.

Phenotypic expression and clinical outcomes in a South Asian PRKAG2 cardiomyopathy cohort.

The PRKAG2 syndrome is a rare autosomal dominant phenocopy of sarcomeric hypertrophic cardiomyopathy (HCM), characterized by ventricular pre-excitation, progressive conduction system disease and left ventricular hypertrophy. This study describes the phenotype, genotype and clinical outcomes of a South-Asian PRKAG2 cardiomyopathy cohort over a 7-year period. Clinical, electrocardiographic, echocardiographic, and cardiac MRI data from 22 individuals with PRKAG2 variants (68% men; mean age 39.5 ± 18.1 years), identified at our HCM centre were studied prospectively. At initial evaluation, all of the patients were in NYHA functional class I or II. The maximum left ventricular wall thickness was 22.9 ± 8.7 mm and left ventricular ejection fraction was 53.4 ± 6.6%. Left ventricular hypertrophy was present in 19 individuals (86%) at baseline. 17 patients had an WPW pattern (77%). After a mean follow-up period of 7 years, 2 patients had undergone accessory pathway ablation, 8 patients (36%) underwent permanent pacemaker implantation (atrio-ventricular blocks-5; sinus node disease-2), 3 patients developed atrial fibrillation, 11 patients (50%) developed progressive worsening in NYHA functional class, and 6 patients (27%) experienced sudden cardiac death or equivalent. PRKAG2 cardiomyopathy must be considered in patients with HCM and progressive conduction system disease.

Integrated genomic analysis reveals mutated ELF3 as a potential gallbladder cancer vaccine candidate.

Gallbladder cancer (GBC) is an aggressive gastrointestinal malignancy with no approved targeted therapy. Here, we analyze exomes (n = 160), transcriptomes (n = 115), and low pass whole genomes (n = 146) from 167 gallbladder cancers (GBCs) from patients in Korea, India and Chile. In addition, we also sequence samples from 39 GBC high-risk patients and detect evidence of early cancer-related genomic lesions. Among the several significantly mutated genes not previously linked to GBC are ETS domain genes ELF3 and EHF, CTNNB1, APC, NSD1, KAT8, STK11 and NFE2L2. A majority of ELF3 alterations are frame-shift mutations that result in several cancer-specific neoantigens that activate T-cells indicating that they are cancer vaccine candidates. In addition, we identify recurrent alterations in KEAP1/NFE2L2 and WNT pathway in GBC. Taken together, these define multiple targetable therapeutic interventions opportunities for GBC treatment and management.

Proteogenomics analysis unveils a TFG-RET gene fusion and druggable targets in papillary thyroid carcinomas.

Papillary thyroid cancer (PTC) is the most common type of endocrine malignancy. By RNA-seq analysis, we identify a RET rearrangement in the tumour material of a patient who does not harbour any known RAS or BRAF mutations. This new gene fusion involves exons 1-4 from the 5' end of the Trk fused Gene (TFG) fused to the 3' end of RET tyrosine kinase leading to a TFG-RET fusion which transforms immortalized human thyroid cells in a kinase-dependent manner. TFG-RET oligomerises in a PB1 domain-dependent manner and oligomerisation of TFG-RET is required for oncogenic transformation. Quantitative proteomic analysis reveals the upregulation of E3 Ubiquitin ligase HUWE1 and DUBs like USP9X and UBP7 in both tumor and metastatic lesions, which is further confirmed in additional patients. Expression of TFG-RET leads to the upregulation of HUWE1 and inhibition of HUWE1 significantly reduces RET-mediated oncogenesis.

The Indian cobra reference genome and transcriptome enables comprehensive identification of venom toxins.

Snakebite envenoming is a serious and neglected tropical disease that kills ~100,000 people annually. High-quality, genome-enabled comprehensive characterization of toxin genes will facilitate development of effective humanized recombinant antivenom. We report a de novo near-chromosomal genome assembly of Naja naja, the Indian cobra, a highly venomous, medically important snake. Our assembly has a scaffold N50 of 223.35 Mb, with 19 scaffolds containing 95% of the genome. Of the 23,248 predicted protein-coding genes, 12,346 venom-gland-expressed genes constitute the 'venom-ome' and this included 139 genes from 33 toxin families. Among the 139 toxin genes were 19 'venom-ome-specific toxins' (VSTs) that showed venom-gland-specific expression, and these probably encode the minimal core venom effector proteins. Synthetic venom reconstituted through recombinant VST expression will aid in the rapid development of safe and effective synthetic antivenom. Additionally, our genome could serve as a reference for snake genomes, support evolutionary studies and enable venom-driven drug discovery.

Identification of African-Specific Admixture between Modern and Archaic Humans.

Recent work has demonstrated that two archaic human groups (Neanderthals and Denisovans) interbred with modern humans and contributed to the contemporary human gene pool. These findings relied on the availability of high-coverage genomes from both Neanderthals and Denisovans. Here we search for evidence of archaic admixture from a worldwide panel of 1,667 individuals using an approach that does not require the presence of an archaic human reference genome. We find no evidence for archaic admixture in the Andaman Islands, as previously claimed, or on the island of Flores, where Homo floresiensis fossils have been found. However, we do find evidence for at least one archaic admixture event in sub-Saharan Africa, with the strongest signal in Khoesan and Pygmy individuals from Southern and Central Africa. The locations of these putative archaic admixture tracts are weighted against functional regions of the genome, consistent with the long-term effects of purifying selection against introgressed genetic material.

S100a4 upregulation in Pik3caH1047R;Trp53R270H;MMTV-Cre-driven mammary tumors promotes metastasis.

BACKGROUND: PIK3CA mutations are frequent in human breast cancer. Pik3caH1047R mutant expression in mouse mammary gland promotes tumorigenesis. TP53 mutations co-occur with PIK3CA mutations in human breast cancers. We previously generated a conditionally activatable Pik3caH1047R;MMTV-Cre mouse model and found a few malignant sarcomatoid (spindle cell) carcinomas that had acquired spontaneous dominant-negative Trp53 mutations. METHODS: A Pik3caH1047R;Trp53R270H;MMTV-Cre double mutant mouse breast cancer model was generated. Tumors were characterized by histology, marker analysis, transcriptional profiling, single-cell RNA-seq, and bioinformatics. Cell lines were developed from mutant tumors and used to identify and confirm genes involved in metastasis. RESULTS: We found Pik3caH1047R and Trp53R270H cooperate in driving oncogenesis in mammary glands leading to a shorter latency than either alone. Double mutant mice develop multiple histologically distinct mammary tumors, including adenocarcinoma and sarcomatoid (spindle cell) carcinoma. We found some tumors to be invasive and a few metastasized to the lung and/or the lymph node. Single-cell RNA-seq analysis of the tumors identified epithelial, stromal, myeloid, and T cell groups. Expression analysis of the metastatic tumors identified S100a4 as a top candidate gene associated with metastasis. Metastatic tumors contained a much higher percentage of epithelial-mesenchymal transition (EMT)-signature positive and S100a4-expressing cells. CRISPR/CAS9-mediated knockout of S100a4 in a metastatic tumor-derived cell line disrupted its metastatic potential indicating a role for S100a4 in metastasis. CONCLUSIONS: Pik3caH1047R;Trp53R270H;MMTV-Cre mouse provides a preclinical model to mimic a subtype of human breast cancers that carry both PIK3CA and TP53 mutations. It also allows for understanding the cooperation between the two mutant genes in tumorigenesis. Our model also provides a system to study metastasis and develop therapeutic strategies for PIK3CA/TP53 double-positive cancers. S100a4 found involved in metastasis in this model can be a potential diagnostic and therapeutic target.

SCINA: A Semi-Supervised Subtyping Algorithm of Single Cells and Bulk Samples.

Advances in single-cell RNA sequencing (scRNA-Seq) have allowed for comprehensive analyses of single cell data. However, current analyses of scRNA-Seq data usually start from unsupervised clustering or visualization. These methods ignore prior knowledge of transcriptomes and the probable structures of the data. Moreover, cell identification heavily relies on subjective and possibly inaccurate human inspection afterwards. To address these analytical challenges, we developed SCINA (Semi-supervised Category Identification and Assignment), a semi-supervised model that exploits previously established gene signatures using an expectation-maximization (EM) algorithm. SCINA is applicable to scRNA-Seq and flow cytometry/CyTOF data, as well as other data of similar format. We applied SCINA to a wide range of datasets, and showed its accuracy, stability and efficiency, which exceeded most popular unsupervised approaches. SCINA discovered an intermediate stage of oligodendrocytes from mouse brain scRNA-Seq data. SCINA also detected immune cell population changes in cytometry data in a genetically-engineered mouse model. Furthermore, SCINA performed well with bulk gene expression data. Specifically, we identified a new kidney tumor clade with similarity to FH-deficient tumors (FHD), which we refer to as FHD-like tumors (FHDL). Overall, SCINA provides both methodological advances and biological insights from perspectives different from traditional analytical methods.

The tumor suppressor BAP1 cooperates with BRAFV600E to promote tumor formation in cutaneous 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.

Actionable Activating Oncogenic ERBB2/HER2 Transmembrane and Juxtamembrane Domain Mutations.

Deregulated HER2 is a target of many approved cancer drugs. We analyzed 111,176 patient tumors and identified recurrent mutations in HER2 transmembrane domain (TMD) and juxtamembrane domain (JMD) that include G660D, R678Q, E693K, and Q709L. Using a saturation mutagenesis screen and testing of patient-derived mutations we found several activating TMD and JMD mutations. Structural modeling and analysis showed that the TMD/JMD mutations function by improving the active dimer interface or stabilizing an activating conformation. Further, we found that HER2 G660D employed asymmetric kinase dimerization for activation and signaling. Importantly, anti-HER2 antibodies and small-molecule kinase inhibitors blocked the activity of TMD/JMD mutants. Consistent with this, a G660D germline mutant lung cancer patient showed remarkable clinical response to HER2 blockade.

An Empirical Approach Leveraging Tumorgrafts to Dissect the Tumor Microenvironment in Renal Cell Carcinoma Identifies Missing Link to Prognostic Inflammatory Factors.

By leveraging tumorgraft (patient-derived xenograft) RNA-sequencing data, we developed an empirical approach, DisHet, to dissect the tumor microenvironment (eTME). We found that 65% of previously defined immune signature genes are not abundantly expressed in renal cell carcinoma (RCC) and identified 610 novel immune/stromal transcripts. Using eTME, genomics, pathology, and medical record data involving >1,000 patients, we established an inflamed pan-RCC subtype (IS) enriched for regulatory T cells, natural killer cells, TH1 cells, neutrophils, macrophages, B cells, and CD8+ T cells. IS is enriched for aggressive RCCs, including BAP1-deficient clear-cell and type 2 papillary tumors. The IS subtype correlated with systemic manifestations of inflammation such as thrombocytosis and anemia, which are enigmatic predictors of poor prognosis. Furthermore, IS was a strong predictor of poor survival. Our analyses suggest that tumor cells drive the stromal immune response. These data provide a missing link between tumor cells, the TME, and systemic factors.Significance: We undertook a novel empirical approach to dissect the renal cell carcinoma TME by leveraging tumorgrafts. The dissection and downstream analyses uncovered missing links between tumor cells, the TME, systemic manifestations of inflammation, and poor prognosis. Cancer Discov; 8(9); 1142-55. ©2018 AACR.This article is highlighted in the In This Issue feature, p. 1047.

ERK Mutations and Amplification Confer Resistance to ERK-Inhibitor Therapy.

Purpose: MAPK pathway inhibitors targeting BRAF and MEK have shown clinical efficacy in patients with RAF- and/or RAS-mutated tumors. However, acquired resistance to these agents has been an impediment to improved long-term survival in the clinic. In such cases, targeting ERK downstream of BRAF/MEK has been proposed as a potential strategy for overcoming acquired resistance. Preclinical studies suggest that ERK inhibitors are effective at inhibiting BRAF/RAS-mutated tumor growth and overcome BRAF or/and MEK inhibitor resistance. However, as observed with other MAPK pathway inhibitors, treatment with ERK inhibitors is likely to cause resistance in the clinic. Here, we aimed to model the mechanism of resistance to ERK inhibitors.Experimental Design: We tested five structurally different ATP-competitive ERK inhibitors representing three different scaffolds on BRAF/RAS-mutant cancer cell lines of different tissue types to generate resistant lines. We have used in vitro modeling, structural biology, and genomic analysis to understand the development of resistance to ERK inhibitors and the mechanisms leading to it.Results: We have identified mutations in ERK1/2, amplification and overexpression of ERK2, and overexpression of EGFR/ERBB2 as mechanisms of acquired resistance. Structural analysis of ERK showed that specific compounds that induced on-target ERK mutations were impaired in their ability to bind mutant ERK. We show that in addition to MEK inhibitors, ERBB receptor and PI3K/mTOR pathway inhibitors are effective in overcoming ERK-inhibitor resistance.Conclusions: These findings suggest that combination therapy with MEK or ERBB receptor or PI3K/mTOR and ERK inhibitors may be an effective strategy for managing the emergence of resistance in the clinic. Clin Cancer Res; 24(16); 4044-55. ©2018 AACR.

Comprehensive genomic analysis identifies pathogenic variants in maturity-onset diabetes of the young (MODY) patients in South India.

BACKGROUND: Maturity-onset diabetes of the young (MODY) is an early-onset, autosomal dominant form of non-insulin dependent diabetes. Genetic diagnosis of MODY can transform patient management. Earlier data on the genetic predisposition to MODY have come primarily from familial studies in populations of European origin. METHODS: In this study, we carried out a comprehensive genomic analysis of 289 individuals from India that included 152 clinically diagnosed MODY cases to identify variants in known MODY genes. Further, we have analyzed exome data to identify putative MODY relevant variants in genes previously not implicated in MODY. Functional validation of MODY relevant variants was also performed. RESULTS: We found MODY 3 (HNF1A; 7.2%) to be most frequently mutated followed by MODY 12 (ABCC8; 3.3%). They together account for ~ 11% of the cases. In addition to known MODY genes, we report the identification of variants in RFX6, WFS1, AKT2, NKX6-1 that may contribute to development of MODY. Functional assessment of the NKX6-1 variants showed that they are functionally impaired. CONCLUSIONS: Our findings showed HNF1A and ABCC8 to be the most frequently mutated MODY genes in south India. Further we provide evidence for additional MODY relevant genes, such as NKX6-1, and these require further validation.

Next-Generation Sequencing Reveals Novel Mutations in X-linked Intellectual Disability.

Robust diagnostics for many human genetic disorders are much needed in the pursuit of global personalized medicine. Next-generation sequencing now offers new promise for biomarker and diagnostic discovery, in developed as well as resource-limited countries. In this broader global health context, X-linked intellectual disability (XLID) is an inherited genetic disorder that is associated with a range of phenotypes impacting societies in both developed and developing countries. Although intellectual disability arises due to diverse causes, a substantial proportion is caused by genomic alterations. Studies have identified causal XLID genomic alterations in more than 100 protein-coding genes located on the X-chromosome. However, the causes for a substantial number of intellectual disability and associated phenotypes still remain unknown. Identification of causative genes and novel mutations will help in early diagnosis as well as genetic counseling of families. Advent of next-generation sequencing methods has accelerated the discovery of new genes involved in mental health disorders. In this study, we analyzed the exomes of three families from India with nonsyndromic XLID comprising seven affected individuals. The affected individuals had varying degrees of intellectual disability, microcephaly, and delayed motor and language milestones. We identified potential causal variants in three XLID genes, including PAK3 (V294M), CASK (complex structural variant), and MECP2 (P354T). Our findings reported in this study extend the spectrum of mutations and phenotypes associated with XLID, and calls for further studies of intellectual disability and mental health disorders with use of next-generation sequencing technologies.

Comprehensive genomic analysis of malignant pleural mesothelioma identifies recurrent mutations, gene fusions and splicing alterations.

We analyzed transcriptomes (n = 211), whole exomes (n = 99) and targeted exomes (n = 103) from 216 malignant pleural mesothelioma (MPM) tumors. Using RNA-seq data, we identified four distinct molecular subtypes: sarcomatoid, epithelioid, biphasic-epithelioid (biphasic-E) and biphasic-sarcomatoid (biphasic-S). Through exome analysis, we found BAP1, NF2, TP53, SETD2, DDX3X, ULK2, RYR2, CFAP45, SETDB1 and DDX51 to be significantly mutated (q-score ≥ 0.8) in MPMs. We identified recurrent mutations in several genes, including SF3B1 (∼2%; 4/216) and TRAF7 (∼2%; 5/216). SF3B1-mutant samples showed a splicing profile distinct from that of wild-type tumors. TRAF7 alterations occurred primarily in the WD40 domain and were, except in one case, mutually exclusive with NF2 alterations. We found recurrent gene fusions and splice alterations to be frequent mechanisms for inactivation of NF2, BAP1 and SETD2. Through integrated analyses, we identified alterations in Hippo, mTOR, histone methylation, RNA helicase and p53 signaling pathways in MPMs.