Pathogenic Germline Variants in 10,389 Adult Cancers.

We conducted the largest investigation of predisposition variants in cancer to date, discovering 853 pathogenic or likely pathogenic variants in 8% of 10,389 cases from 33 cancer types. Twenty-one genes showed single or cross-cancer associations, including novel associations of SDHA in melanoma and PALB2 in stomach adenocarcinoma. The 659 predisposition variants and 18 additional large deletions in tumor suppressors, including ATM, BRCA1, and NF1, showed low gene expression and frequent (43%) loss of heterozygosity or biallelic two-hit events. We also discovered 33 such variants in oncogenes, including missenses in MET, RET, and PTPN11 associated with high gene expression. We nominated 47 additional predisposition variants from prioritized VUSs supported by multiple evidences involving case-control frequency, loss of heterozygosity, expression effect, and co-localization with mutations and modified residues. Our integrative approach links rare predisposition variants to functional consequences, informing future guidelines of variant classification and germline genetic testing in cancer.

Functional Genomic Landscape of Human Breast Cancer Drivers, Vulnerabilities, and Resistance.

Large-scale genomic studies have identified multiple somatic aberrations in breast cancer, including copy number alterations and point mutations. Still, identifying causal variants and emergent vulnerabilities that arise as a consequence of genetic alterations remain major challenges. We performed whole-genome small hairpin RNA (shRNA) "dropout screens" on 77 breast cancer cell lines. Using a hierarchical linear regression algorithm to score our screen results and integrate them with accompanying detailed genetic and proteomic information, we identify vulnerabilities in breast cancer, including candidate "drivers," and reveal general functional genomic properties of cancer cells. Comparisons of gene essentiality with drug sensitivity data suggest potential resistance mechanisms, effects of existing anti-cancer drugs, and opportunities for combination therapy. Finally, we demonstrate the utility of this large dataset by identifying BRD4 as a potential target in luminal breast cancer and PIK3CA mutations as a resistance determinant for BET-inhibitors.

New Brain Tumor Entities Emerge from Molecular Classification of CNS-PNETs.

Primitive neuroectodermal tumors of the central nervous system (CNS-PNETs) are highly aggressive, poorly differentiated embryonal tumors occurring predominantly in young children but also affecting adolescents and adults. Herein, we demonstrate that a significant proportion of institutionally diagnosed CNS-PNETs display molecular profiles indistinguishable from those of various other well-defined CNS tumor entities, facilitating diagnosis and appropriate therapy for patients with these tumors. From the remaining fraction of CNS-PNETs, we identify four new CNS tumor entities, each associated with a recurrent genetic alteration and distinct histopathological and clinical features. These new molecular entities, designated "CNS neuroblastoma with FOXR2 activation (CNS NB-FOXR2)," "CNS Ewing sarcoma family tumor with CIC alteration (CNS EFT-CIC)," "CNS high-grade neuroepithelial tumor with MN1 alteration (CNS HGNET-MN1)," and "CNS high-grade neuroepithelial tumor with BCOR alteration (CNS HGNET-BCOR)," will enable meaningful clinical trials and the development of therapeutic strategies for patients affected by poorly differentiated CNS tumors.

Pan-cancer ion transport signature reveals functional regulators of glioblastoma aggression.

Ion channels, transporters, and other ion-flux controlling proteins, collectively comprising the "ion permeome", are common drug targets, however, their roles in cancer remain understudied. Our integrative pan-cancer transcriptome analysis shows that genes encoding the ion permeome are significantly more often highly expressed in specific subsets of cancer samples, compared to pan-transcriptome expectations. To enable target selection, we identified 410 survival-associated IP genes in 33 cancer types using a machine-learning approach. Notably, GJB2 and SCN9A show prominent expression in neoplastic cells and are associated with poor prognosis in glioblastoma, the most common and aggressive brain cancer. GJB2 or SCN9A knockdown in patient-derived glioblastoma cells induces transcriptome-wide changes involving neuron projection and proliferation pathways, impairs cell viability and tumor sphere formation in vitro, perturbs tunneling nanotube dynamics, and extends the survival of glioblastoma-bearing mice. Thus, aberrant activation of genes encoding ion transport proteins appears as a pan-cancer feature defining tumor heterogeneity, which can be exploited for mechanistic insights and therapy development.

Candidate Cancer Driver Mutations in Distal Regulatory Elements and Long-Range Chromatin Interaction Networks.

A comprehensive catalog of cancer driver mutations is essential for understanding tumorigenesis and developing therapies. Exome-sequencing studies have mapped many protein-coding drivers, yet few non-coding drivers are known because genome-wide discovery is challenging. We developed a driver discovery method, ActiveDriverWGS, and analyzed 120,788 cis-regulatory modules (CRMs) across 1,844 whole tumor genomes from the ICGC-TCGA PCAWG project. We found 30 CRMs with enriched SNVs and indels (FDR < 0.05). These frequently mutated regulatory elements (FMREs) were ubiquitously active in human tissues, showed long-range chromatin interactions and mRNA abundance associations with target genes, and were enriched in motif-rewiring mutations and structural variants. Genomic deletion of one FMRE in human cells caused proliferative deficiencies and transcriptional deregulation of cancer genes CCNB1IP1, CDH1, and CDKN2B, validating observations in FMRE-mutated tumors. Pathway analysis revealed further sub-significant FMREs at cancer genes and processes, indicating an unexplored landscape of infrequent driver mutations in the non-coding genome.

Therapeutic targeting of ependymoma as informed by oncogenic enhancer profiling.

Genomic sequencing has driven precision-based oncology therapy; however, the genetic drivers of many malignancies remain unknown or non-targetable, so alternative approaches to the identification of therapeutic leads are necessary. Ependymomas are chemotherapy-resistant brain tumours, which, despite genomic sequencing, lack effective molecular targets. Intracranial ependymomas are segregated on the basis of anatomical location (supratentorial region or posterior fossa) and further divided into distinct molecular subgroups that reflect differences in the age of onset, gender predominance and response to therapy. The most common and aggressive subgroup, posterior fossa ependymoma group A (PF-EPN-A), occurs in young children and appears to lack recurrent somatic mutations. Conversely, posterior fossa ependymoma group B (PF-EPN-B) tumours display frequent large-scale copy number gains and losses but have favourable clinical outcomes. More than 70% of supratentorial ependymomas are defined by highly recurrent gene fusions in the NF-κB subunit gene RELA (ST-EPN-RELA), and a smaller number involve fusion of the gene encoding the transcriptional activator YAP1 (ST-EPN-YAP1). Subependymomas, a distinct histologic variant, can also be found within the supratetorial and posterior fossa compartments, and account for the majority of tumours in the molecular subgroups ST-EPN-SE and PF-EPN-SE. Here we describe mapping of active chromatin landscapes in 42 primary ependymomas in two non-overlapping primary ependymoma cohorts, with the goal of identifying essential super-enhancer-associated genes on which tumour cells depend. Enhancer regions revealed putative oncogenes, molecular targets and pathways; inhibition of these targets with small molecule inhibitors or short hairpin RNA diminished the proliferation of patient-derived neurospheres and increased survival in mouse models of ependymomas. Through profiling of transcriptional enhancers, our study provides a framework for target and drug discovery in other cancers that lack known genetic drivers and are therefore difficult to treat.

Human phospho-signaling networks of SARS-CoV-2 infection are rewired by population genetic variants.

SARS-CoV-2 infection hijacks signaling pathways and induces protein-protein interactions between human and viral proteins. Human genetic variation may impact SARS-CoV-2 infection and COVID-19 pathology; however, the genetic variation in these signaling networks remains uncharacterized. Here, we studied human missense single nucleotide variants (SNVs) altering phosphorylation sites modulated by SARS-CoV-2 infection, using machine learning to identify amino acid substitutions altering kinase-bound sequence motifs. We found 2,033 infrequent phosphorylation-associated SNVs (pSNVs) that are enriched in sequence motif alterations, potentially reflecting the evolution of signaling networks regulating host defenses. Proteins with pSNVs are involved in viral life cycle and host responses, including RNA splicing, interferon response (TRIM28), and glucose homeostasis (TBC1D4) with potential associations with COVID-19 comorbidities. pSNVs disrupt CDK and MAPK substrate motifs and replace these with motifs of Tank Binding Kinase 1 (TBK1) involved in innate immune responses, indicating consistent rewiring of signaling networks. Several pSNVs associate with severe COVID-19 and hospitalization (STARD13, ARFGEF2). Our analysis highlights potential genetic factors contributing to inter-individual variation of SARS-CoV-2 infection and COVID-19 and suggests leads for mechanistic and translational studies.

Chromatin accessibility of primary human cancers ties regional mutational processes and signatures with tissues of origin.

Somatic mutations in cancer genomes are associated with DNA replication timing (RT) and chromatin accessibility (CA), however these observations are based on normal tissues and cell lines while primary cancer epigenomes remain uncharacterised. Here we use machine learning to model megabase-scale mutation burden in 2,500 whole cancer genomes and 17 cancer types via a compendium of 900 CA and RT profiles covering primary cancers, normal tissues, and cell lines. CA profiles of primary cancers, rather than those of normal tissues, are most predictive of regional mutagenesis in most cancer types. Feature prioritisation shows that the epigenomes of matching cancer types and organ systems are often the strongest predictors of regional mutation burden, highlighting disease-specific associations of mutational processes. The genomic distributions of mutational signatures are also shaped by the epigenomes of matched cancer and tissue types, with SBS5/40, carcinogenic and unknown signatures most accurately predicted by our models. In contrast, fewer associations of RT and regional mutagenesis are found. Lastly, the models highlight genomic regions with overrepresented mutations that dramatically exceed epigenome-derived expectations and show a pan-cancer convergence to genes and pathways involved in development and oncogenesis, indicating the potential of this approach for coding and non-coding driver discovery. The association of regional mutational processes with the epigenomes of primary cancers suggests that the landscape of passenger mutations is predominantly shaped by the epigenomes of cancer cells after oncogenic transformation.

Divergent clonal selection dominates medulloblastoma at recurrence.

The development of targeted anti-cancer therapies through the study of cancer genomes is intended to increase survival rates and decrease treatment-related toxicity. We treated a transposon-driven, functional genomic mouse model of medulloblastoma with 'humanized' in vivo therapy (microneurosurgical tumour resection followed by multi-fractionated, image-guided radiotherapy). Genetic events in recurrent murine medulloblastoma exhibit a very poor overlap with those in matched murine diagnostic samples (<5%). Whole-genome sequencing of 33 pairs of human diagnostic and post-therapy medulloblastomas demonstrated substantial genetic divergence of the dominant clone after therapy (<12% diagnostic events were retained at recurrence). In both mice and humans, the dominant clone at recurrence arose through clonal selection of a pre-existing minor clone present at diagnosis. Targeted therapy is unlikely to be effective in the absence of the target, therefore our results offer a simple, proximal, and remediable explanation for the failure of prior clinical trials of targeted therapy.

A transcriptome-based signature of pathological angiogenesis predicts breast cancer patient survival.

The specific genes and molecules that drive physiological angiogenesis differ from those involved in pathological angiogenesis, suggesting distinct mechanisms for these seemingly related processes. Unveiling genes and pathways preferentially associated with pathologic angiogenesis is key to understanding its mechanisms, thereby facilitating development of novel approaches to managing angiogenesis-dependent diseases. To better understand these different processes, we elucidated the transcriptome of the mouse retina in the well-accepted oxygen-induced retinopathy (OIR) model of pathological angiogenesis. We identified 153 genes changed between normal and OIR retinas, which represent a molecular signature relevant to other angiogenesis-dependent processes such as cancer. These genes robustly predict the survival of breast cancer patients, which was validated in an independent 1,000-patient test cohort (40% difference in 15-year survival; p = 2.56 x 10-21). These results suggest that the OIR model reveals key genes involved in pathological angiogenesis, and these may find important applications in stratifying tumors for treatment intensification or for angiogenesis-targeted therapies.

Phosphoproteome and drug-response effects mediated by the three protein phosphatase 2A inhibitor proteins CIP2A, SET, and PME-1.

Protein phosphatase 2A (PP2A) critically regulates cell signaling and is a human tumor suppressor. PP2A complexes are modulated by proteins such as cancerous inhibitor of protein phosphatase 2A (CIP2A), protein phosphatase methylesterase 1 (PME-1), and SET nuclear proto-oncogene (SET) that often are deregulated in cancers. However, how they impact cellular phosphorylation and how redundant they are in cellular regulation is poorly understood. Here, we conducted a systematic phosphoproteomics screen for phosphotargets modulated by siRNA-mediated depletion of CIP2A, PME-1, and SET (to reactivate PP2A) or the scaffolding A-subunit of PP2A (PPP2R1A) (to inhibit PP2A) in HeLa cells. We identified PP2A-modulated targets in diverse cellular pathways, including kinase signaling, cytoskeleton, RNA splicing, DNA repair, and nuclear lamina. The results indicate nonredundancy among CIP2A, PME-1, and SET in phosphotarget regulation. Notably, PP2A inhibition or reactivation affected largely distinct phosphopeptides, introducing a concept of nonoverlapping phosphatase inhibition- and activation-responsive sites (PIRS and PARS, respectively). This phenomenon is explained by the PPP2R1A inhibition impacting primarily dephosphorylated threonines, whereas PP2A reactivation results in dephosphorylation of clustered and acidophilic sites. Using comprehensive drug-sensitivity screening in PP2A-modulated cells to evaluate the functional impact of PP2A across diverse cellular pathways targeted by these drugs, we found that consistent with global phosphoproteome effects, PP2A modulations broadly affect responses to more than 200 drugs inhibiting a broad spectrum of cancer-relevant targets. These findings advance our understanding of the phosphoproteins, pharmacological responses, and cellular processes regulated by PP2A modulation and may enable the development of combination therapies.

Phosphoproteomics Analysis Identifies Novel Candidate Substrates of the Nonreceptor Tyrosine Kinase, Src-related Kinase Lacking C-terminal Regulatory Tyrosine and N-terminal Myristoylation Sites (SRMS).

SRMS (Src-related kinase lacking C-terminal regulatory tyrosine and N-terminal myristoylation sites), also known as PTK 70 (Protein tyrosine kinase 70), is a non-receptor tyrosine kinase that belongs to the BRK family of kinases (BFKs). To date less is known about the cellular role of SRMS primarily because of the unidentified substrates or signaling intermediates regulated by the kinase. In this study, we used phosphotyrosine antibody-based immunoaffinity purification in large-scale label-free quantitative phosphoproteomics to identify novel candidate substrates of SRMS. Our analyses led to the identification of 1258 tyrosine-phosphorylated peptides which mapped to 663 phosphoproteins, exclusively from SRMS-expressing cells. DOK1, a previously characterized SRMS substrate, was also identified in our analyses. Functional enrichment analyses revealed that the candidate SRMS substrates were enriched in various biological processes including protein ubiquitination, mitotic cell cycle, energy metabolism and RNA processing, as well as Wnt and TNF signaling. Analyses of the sequence surrounding the phospho-sites in these proteins revealed novel candidate SRMS consensus substrate motifs. We utilized customized high-throughput peptide arrays to validate a subset of the candidate SRMS substrates identified in our MS-based analyses. Finally, we independently validated Vimentin and Sam68, as bona fide SRMS substrates through in vitro and in vivo assays. Overall, our study identified a number of novel and biologically relevant SRMS candidate substrates, which suggests the involvement of the kinase in a vast array of unexplored cellular functions.

ActiveDriverDB: human disease mutations and genome variation in post-translational modification sites of proteins.

Interpretation of genetic variation is needed for deciphering genotype-phenotype associations, mechanisms of inherited disease, and cancer driver mutations. Millions of single nucleotide variants (SNVs) in human genomes are known and thousands are associated with disease. An estimated 21% of disease-associated amino acid substitutions corresponding to missense SNVs are located in protein sites of post-translational modifications (PTMs), chemical modifications of amino acids that extend protein function. ActiveDriverDB is a comprehensive human proteo-genomics database that annotates disease mutations and population variants through the lens of PTMs. We integrated >385,000 published PTM sites with ∼3.6 million substitutions from The Cancer Genome Atlas (TCGA), the ClinVar database of disease genes, and human genome sequencing projects. The database includes site-specific interaction networks of proteins, upstream enzymes such as kinases, and drugs targeting these enzymes. We also predicted network-rewiring impact of mutations by analyzing gains and losses of kinase-bound sequence motifs. ActiveDriverDB provides detailed visualization, filtering, browsing and searching options for studying PTM-associated mutations. Users can upload mutation datasets interactively and use our application programming interface in pipelines. Integrative analysis of mutations and PTMs may help decipher molecular mechanisms of phenotypes and disease, as exemplified by case studies of TP53, BRCA2 and VHL. The open-source database is available at https://www.ActiveDriverDB.org.

MIMP: predicting the impact of mutations on kinase-substrate phosphorylation.

Protein phosphorylation is important in cellular pathways and altered in disease. We developed MIMP (http://mimp.baderlab.org/), a machine learning method to predict the impact of missense single-nucleotide variants (SNVs) on kinase-substrate interactions. MIMP analyzes kinase sequence specificities and predicts whether SNVs disrupt existing phosphorylation sites or create new sites. This helps discover mutations that modify protein function by altering kinase networks and provides insight into disease biology and therapy development.

Pathway and network analysis of cancer genomes.

Genomic information on tumors from 50 cancer types cataloged by the International Cancer Genome Consortium (ICGC) shows that only a few well-studied driver genes are frequently mutated, in contrast to many infrequently mutated genes that may also contribute to tumor biology. Hence there has been large interest in developing pathway and network analysis methods that group genes and illuminate the processes involved. We provide an overview of these analysis techniques and show where they guide mechanistic and translational investigations.

Global phosphoproteomic analysis identifies SRMS-regulated secondary signaling intermediates.

BACKGROUND: The non-receptor tyrosine kinase, SRMS (Src-related kinase lacking C-terminal regulatory tyrosine and N-terminal myristoylation sites) is a member of the BRK family kinases (BFKs) which represents an evolutionarily conserved relative of the Src family kinases (SFKs). Tyrosine kinases are known to regulate a number of cellular processes and pathways via phosphorylating substrate proteins directly and/or by partaking in signaling cross-talks leading to the indirect modulation of various signaling intermediates. In a previous study, we profiled the tyrosine-phosphoproteome of SRMS and identified multiple candidate substrates of the kinase. The broader cellular signaling intermediates of SRMS are unknown. METHODS: In order to uncover the broader SRMS-regulated phosphoproteome and identify the SRMS-regulated indirect signaling intermediates, we performed label-free global phosphoproteomics analysis on cells expressing wild-type SRMS. Using computational database searching and bioinformatics analyses we characterized the dataset. RESULTS: Our analyses identified 60 hyperphosphorylated (phosphoserine/phosphothreonine) proteins mapped from 140 hyperphosphorylated peptides. Bioinfomatics analyses identified a number of significantly enriched biological and cellular processes among which DNA repair pathways were found to be upregulated while apoptotic pathways were found to be downregulated. Analyses of motifs derived from the upregulated phosphosites identified Casein kinase 2 alpha (CK2α) as one of the major potential kinases contributing to the SRMS-dependent indirect regulation of signaling intermediates. CONCLUSIONS: Overall, our phosphoproteomics analyses identified serine/threonine phosphorylation dynamics as important secondary events of the SRMS-regulated phosphoproteome with implications in the regulation of cellular and biological processes.

Subgroup-specific structural variation across 1,000 medulloblastoma genomes.

Medulloblastoma, the most common malignant paediatric brain tumour, is currently treated with nonspecific cytotoxic therapies including surgery, whole-brain radiation, and aggressive chemotherapy. As medulloblastoma exhibits marked intertumoural heterogeneity, with at least four distinct molecular variants, previous attempts to identify targets for therapy have been underpowered because of small samples sizes. Here we report somatic copy number aberrations (SCNAs) in 1,087 unique medulloblastomas. SCNAs are common in medulloblastoma, and are predominantly subgroup-enriched. The most common region of focal copy number gain is a tandem duplication of SNCAIP, a gene associated with Parkinson's disease, which is exquisitely restricted to Group 4α. Recurrent translocations of PVT1, including PVT1-MYC and PVT1-NDRG1, that arise through chromothripsis are restricted to Group 3. Numerous targetable SCNAs, including recurrent events targeting TGF-ß signalling in Group 3, and NF-κB signalling in Group 4, suggest future avenues for rational, targeted therapy.

g:Profiler-a web server for functional interpretation of gene lists (2016 update).

Functional enrichment analysis is a key step in interpreting gene lists discovered in diverse high-throughput experiments. g:Profiler studies flat and ranked gene lists and finds statistically significant Gene Ontology terms, pathways and other gene function related terms. Translation of hundreds of gene identifiers is another core feature of g:Profiler. Since its first publication in 2007, our web server has become a popular tool of choice among basic and translational researchers. Timeliness is a major advantage of g:Profiler as genome and pathway information is synchronized with the Ensembl database in quarterly updates. g:Profiler supports 213 species including mammals and other vertebrates, plants, insects and fungi. The 2016 update of g:Profiler introduces several novel features. We have added further functional datasets to interpret gene lists, including transcription factor binding site predictions, Mendelian disease annotations, information about protein expression and complexes and gene mappings of human genetic polymorphisms. Besides the interactive web interface, g:Profiler can be accessed in computational pipelines using our R package, Python interface and BioJS component. g:Profiler is freely available at http://biit.cs.ut.ee/gprofiler/.

Evolutionary constraint and disease associations of post-translational modification sites in human genomes.

Interpreting the impact of human genome variation on phenotype is challenging. The functional effect of protein-coding variants is often predicted using sequence conservation and population frequency data, however other factors are likely relevant. We hypothesized that variants in protein post-translational modification (PTM) sites contribute to phenotype variation and disease. We analyzed fraction of rare variants and non-synonymous to synonymous variant ratio (Ka/Ks) in 7,500 human genomes and found a significant negative selection signal in PTM regions independent of six factors, including conservation, codon usage, and GC-content, that is widely distributed across tissue-specific genes and function classes. PTM regions are also enriched in known disease mutations, suggesting that PTM variation is more likely deleterious. PTM constraint also affects flanking sequence around modified residues and increases around clustered sites, indicating presence of functionally important short linear motifs. Using target site motifs of 124 kinases, we predict that at least ∼180,000 motif-breaker amino acid residues that disrupt PTM sites when substituted, and highlight kinase motifs that show specific negative selection and enrichment of disease mutations. We provide this dataset with corresponding hypothesized mechanisms as a community resource. As an example of our integrative approach, we propose that PTPN11 variants in Noonan syndrome aberrantly activate the protein by disrupting an uncharacterized cluster of phosphorylation sites. Further, as PTMs are molecular switches that are modulated by drugs, we study mutated binding sites of PTM enzymes in disease genes and define a drug-disease network containing 413 novel predicted disease-gene links.

Single cell-derived clonal analysis of human glioblastoma links functional and genomic heterogeneity.

Glioblastoma (GBM) is a cancer comprised of morphologically, genetically, and phenotypically diverse cells. However, an understanding of the functional significance of intratumoral heterogeneity is lacking. We devised a method to isolate and functionally profile tumorigenic clones from patient glioblastoma samples. Individual clones demonstrated unique proliferation and differentiation abilities. Importantly, naïve patient tumors included clones that were temozolomide resistant, indicating that resistance to conventional GBM therapy can preexist in untreated tumors at a clonal level. Further, candidate therapies for resistant clones were detected with clone-specific drug screening. Genomic analyses revealed genes and pathways that associate with specific functional behavior of single clones. Our results suggest that functional clonal profiling used to identify tumorigenic and drug-resistant tumor clones will lead to the discovery of new GBM clone-specific treatment strategies.