Search

Summary Somatic DNA copy number variations (CNVs) are prevalent in cancer and can drive cancer progression, albeit with often uncharacterized roles in altering cell signaling states. Here, we integrate genomic and proteomic data for 5,598 tumor samples to identify CNVs leading to aberrant signal transduction. The resulting associations recapitulate known kinase-substrate relationships, and further network analysis prioritizes likely causal genes. Of the 303 significant associations we identify from the pan-tumor analysis, 43% are replicated in cancer cell lines, including 44 robust gene-phosphosite associations identified across multiple tumor types. Several predicted regulators of hippo signaling are experimentally validated. Using RNAi, CRISPR, and drug screening data, we find evidence of kinase addiction in cancer cell lines, identifying inhibitors for targeting of kinase-dependent cell lines. We propose copy number status of genes as a useful predictor of differential impact of kinase inhibition, a strategy that may be of use in the future for anticancer therapies., Graphical abstract, Highlights • Alterations in copy number and in signaling pathways are robustly associated in cancer • Differences in essentiality of phosphoproteins show evidence of kinase addiction • Copy number markers of dysregulated signaling suggest treatment stratification options, Memon et al. show that copy number alterations can be used as molecular fingerprints to predict the activity of signaling pathways in cancer. Kinase dependencies of cancer cells can be exploited for targeted anticancer therapy.

Metastatic melanoma carries a poor prognosis despite modern systemic therapies. Understanding the evolution of the disease could help inform patient management. Through whole-genome sequencing of 13 melanoma metastases sampled at autopsy from a treatment naïve patient and by leveraging the analytical power of multi-sample analyses, we reveal evidence of diversification among metastatic lineages. UV-induced mutations dominate the trunk, whereas APOBEC-associated mutations are found in the branches of the evolutionary tree. Multi-sample analyses from a further seven patients confirmed that lineage diversification was pervasive, representing an important mode of melanoma dissemination. Our analyses demonstrate that joint analysis of cancer cell fraction estimates across multiple metastases can uncover previously unrecognised levels of tumour heterogeneity and highlight the limitations of inferring heterogeneity from a single biopsy., Metastatic melanoma is associated with a poor prognosis and understanding the genetic features of metastases may enable better treatment strategies. Here, the authors analyse multiple metastases from individual patients finding high levels of heterogeneity in metastases from different organs.

Myxofibrosarcoma is a common mesenchymal malignancy with complex genomics and heterogeneous clinical outcomes. Through gene-expression profiling of 64 primary high-grade myxofibrosarcomas, we defined an expression signature associated with clinical outcome. The gene most significantly associated with disease-specific death and distant metastasis was ITGA10 (integrin-α10). Functional studies revealed that myxofibrosarcoma cells strongly depended on integrin-α10, whereas normal mesenchymal cells did not. Integrin-α10 transmitted its tumor-specific signal via TRIO and RICTOR, two oncoproteins that are frequently co-overexpressed through gene amplification on chromosome 5p. TRIO and RICTOR activated RAC/PAK and AKT/mTOR to promote sarcoma cell survival. Inhibition of these proteins with EHop-016 (RAC inhibitor) and INK128 (mTOR inhibitor) had antitumor effects in tumor-derived cell lines and mouse xenografts, and combining the drugs enhanced the effects. Our results demonstrate the importance of integrin-α10/TRIO/RICTOR signaling for driving myxofibrosarcoma progression and provide the basis for promising targeted treatment strategies for patients with high-risk disease. Significance: Identifying the molecular pathogenesis for myxofibrosarcoma progression has proven challenging given the highly complex genomic alterations in this tumor type. We found that integrin-α10 promotes tumor cell survival through activation of TRIO–RAC–RICTOR–mTOR signaling, and that inhibitors of RAC and mTOR have antitumor effects in vivo, thus identifying a potential treatment strategy for patients with high-risk myxofibrosarcoma. Cancer Discov; 6(10); 1148–65. ©2016 AACR. This article is highlighted in the In This Issue feature, p. 1069

Simple Summary Clinical responses to the initial treatment of high grade serous ovarian cancer (HGSOC) vary greatly. Widespread intra-site and inter-site genomic heterogeneity presents significant challenges for the development of predictive biomarkers based on pre-treatment sampling of select individual tumors. Non-invasive stratification of patients with HGSOC by risk of outcome could facilitate a higher level of intervention for those with the highest risk of a poor outcome. We developed and validated a machine learning-based integrated marker of HGSOC outcomes to standard chemotherapy that combines a previously developed intra-site and inter-site CT radiomics measure called cluster dissimilarity (cluDiss) with clinical and genomic measures using two retrospective cohorts of internal and external institution datasets. Our approach was more accurate than conventional clinical and average radiomics measures for prognosticating progression-free survival and platinum resistance. Abstract Purpose: Develop an integrated intra-site and inter-site radiomics-clinical-genomic marker of high grade serous ovarian cancer (HGSOC) outcomes and explore the biological basis of radiomics with respect to molecular signaling pathways and the tumor microenvironment (TME). Method: Seventy-five stage III-IV HGSOC patients from internal (N = 40) and external factors via the Cancer Imaging Archive (TCGA) (N = 35) with pre-operative contrast enhanced CT, attempted primary cytoreduction, at least two disease sites, and molecular analysis performed within TCGA were retrospectively analyzed. An intra-site and inter-site radiomics (cluDiss) measure was combined with clinical-genomic variables (iRCG) and compared against conventional (volume and number of sites) and average radiomics (N = 75) for prognosticating progression-free survival (PFS) and platinum resistance. Correlation with molecular signaling and TME derived using a single sample gene set enrichment that was measured. Results: The iRCG model had the best platinum resistance classification accuracy (AUROC of 0.78 [95% CI 0.77 to 0.80]). CluDiss was associated with PFS (HR 1.03 [95% CI: 1.01 to 1.05], p = 0.002), negatively correlated with Wnt signaling, and positively to immune TME. Conclusions: CluDiss and the iRCG prognosticated HGSOC outcomes better than conventional and average radiomic measures and could better stratify patient outcomes if validated on larger multi-center trials.

Somatic mutations can result in the formation of neoantigens, immunogenic peptides that are presented on the tumor cell surface by HLA molecules. These mutations are expected to be under negative selection pressure, but the extent of the resulting neoantigen depletion remains unclear. On the basis of HLA affinity predictions, we annotated the human genome for its translatability to HLA binding peptides and screened for reduced single nucleotide substitution rates in large genomic data sets from untreated cancers. Apparent neoantigen depletion signals become negligible when taking into consideration trinucleotide-based mutational signatures, owing to lack of power or to efficient immune evasion mechanisms that are active early during tumor evolution.

In cancer genomics, recurrence of mutations in independent tumor samples is a strong indicator of functional impact. However, rare functional mutations can escape detection by recurrence analysis owing to lack of statistical power. We enhance statistical power by extending the notion of recurrence of mutations from single genes to gene families that share homologous protein domains. Domain mutation analysis also sharpens the functional interpretation of the impact of mutations, as domains more succinctly embody function than entire genes. By mapping mutations in 22 different tumor types to equivalent positions in multiple sequence alignments of domains, we confirm well-known functional mutation hotspots; identify uncharacterized rare variants in one gene that are equivalent to well-characterized mutations in another gene; detect previously unknown mutation hotspots; and provide hypotheses about molecular mechanisms and downstream effects of domain mutations. With the rapid expansion of cancer genomics projects, protein domain hotspot analysis will likely provide many more leads linking mutations in proteins to the cancer phenotype.

To address limitations of current high-throughput methods for studying cell-cell communication and determining the cell-of-origin of proteins in multicellular environments, we have developed a technique that selectively and continuously labels the proteome of individual cell types in co-culture. Through transgenic expression of exogenous amino acid biosynthesis enzymes, vertebrate cells overcome their dependence on essential amino acids and can be selectively labeled through metabolic incorporation of amino acids produced from heavy isotope-labeled precursors. We have named this method Type specific labeling with Amino acid Precursors (CTAP). Testing CTAP in several human and mouse cell lines, we were able to differentially label the proteome of distinct cell populations in co-culture and determine the relative expression of proteins by quantitative mass spectrometry. In addition, CTAP successfully identified the cell-of-origin of extracellular proteins in co-culture, highlighting its potential use in biomarker discovery for linking secreted factors to their cellular source.

Purpose: Inflammatory bowel disease–associated colorectal cancers (IBD-CRC) are associated with a higher mortality than sporadic colorectal cancers. The poorly defined molecular pathogenesis of IBD-CRCs limits development of effective prevention, detection, and treatment strategies. We aimed to identify biomarkers using whole-exome sequencing of IBD-CRCs to guide individualized management. Experimental Design: Whole-exome sequencing was performed on 34 formalin-fixed paraffin-embedded primary IBD-CRCs and 31 matched normal lymph nodes. Computational methods were used to identify somatic point mutations, small insertions and deletions, mutational signatures, and somatic copy number alterations. Mismatch repair status was examined. Results: Hypermutation was observed in 27% of IBD-CRCs. All hypermutated cancers were from the proximal colon; all but one of the cancers with hypermutation had defective mismatch repair or somatic mutations in the proofreading domain of DNA POLE. Hypermutated IBD-CRCs had increased numbers of predicted neo-epitopes, which could be exploited using immunotherapy. We identified six distinct mutation signatures in IBD-CRCs, three of which corresponded to known mechanisms of mutagenesis. Driver genes were also identified. Conclusions: IBD-CRCs should be evaluated for hypermutation and defective mismatch repair to identify patients with a higher neo-epitope load who may benefit from immunotherapies. Prospective trials are required to determine whether IHC to detect loss of MLH1 expression in dysplastic colonic tissue could identify patients at increased risk of developing IBD-CRC. We identified mutations in genes in IBD-CRCs with hypermutation that might be targeted therapeutically. These approaches would complement and individualize surveillance and treatment programs. Clin Cancer Res; 24(20); 5133–42. ©2018 AACR.

SummaryMembers of the SH2 domain family modulate signal transduction by binding to short peptides containing phosphorylated tyrosines. Each domain displays a distinct preference for the sequence context of the phosphorylated residue. We have developed a high-density peptide chip technology that allows for probing of the affinity of most SH2 domains for a large fraction of the entire complement of tyrosine phosphopeptides in the human proteome. Using this technique, we have experimentally identified thousands of putative SH2-peptide interactions for more than 70 different SH2 domains. By integrating this rich data set with orthogonal context-specific information, we have assembled an SH2-mediated probabilistic interaction network, which we make available as a community resource in the PepspotDB database. A predicted dynamic interaction between the SH2 domains of the tyrosine phosphatase SHP2 and the phosphorylated tyrosine in the extracellular signal-regulated kinase activation loop was validated by experiments in living cells.

Cancer therapy is challenged by the diversity of molecular implementations of oncogenic processes and by the resulting variation in therapeutic responses. Projects such as The Cancer Genome Atlas (TCGA) provide molecular tumor maps in unprecedented detail. The interpretation of these maps remains a major challenge. Here we distilled thousands of genetic and epigenetic features altered in cancers to ∼500 selected functional events (SFEs). Using this simplified description, we derived a hierarchical classification of 3,299 TCGA tumors from 12 cancer types. The top classes are dominated by either mutations (M class) or copy number changes (C class). This distinction is clearest at the extremes of genomic instability, indicating the presence of different oncogenic processes. The full hierarchy shows functional event patterns characteristic of multiple cross-tissue groups of tumors, termed oncogenic signature classes. Targetable functional events in a tumor class are suggestive of class-specific combination therapy. These results may assist in the definition of clinical trials to match actionable oncogenic signatures with personalized therapies.

Transfection of small RNAs (such as small interfering RNAs (siRNAs) and microRNAs (miRNAs)) into cells typically lowers expression of many genes. Unexpectedly, increased expression of genes also occurs. We investigated whether this upregulation results from a saturation effect--that is, competition among the transfected small RNAs and the endogenous pool of miRNAs for the intracellular machinery that processes small RNAs. To test this hypothesis, we analyzed genome-wide transcript responses from 151 published transfection experiments in seven different human cell types. We show that targets of endogenous miRNAs are expressed at significantly higher levels after transfection, consistent with impaired effectiveness of endogenous miRNA repression. This effect exhibited concentration and temporal dependence. Notably, the profile of endogenous miRNAs can be largely inferred by correlating miRNA sites with gene expression changes after transfections. The competition and saturation effects have practical implications for miRNA target prediction, the design of siRNA and short hairpin RNA (shRNA) genomic screens and siRNA therapeutics.

Summary We describe the landscape of genomic alterations in cutaneous melanomas through DNA, RNA, and protein-based analysis of 333 primary and/or metastatic melanomas from 331 patients. We establish a framework for genomic classification into one of four subtypes based on the pattern of the most prevalent significantly mutated genes: mutant BRAF , mutant RAS , mutant NF1 , and Triple-WT (wild-type). Integrative analysis reveals enrichment of KIT mutations and focal amplifications and complex structural rearrangements as a feature of the Triple-WT subtype. We found no significant outcome correlation with genomic classification, but samples assigned a transcriptomic subclass enriched for immune gene expression associated with lymphocyte infiltrate on pathology review and high LCK protein expression, a T cell marker, were associated with improved patient survival. This clinicopathological and multi-dimensional analysis suggests that the prognosis of melanoma patients with regional metastases is influenced by tumor stroma immunobiology, offering insights to further personalize therapeutic decision-making.

Immune checkpoint inhibitors, which unleash a patient’s own T cells to kill tumors, are revolutionizing cancer treatment. To unravel the genomic determinants of response to this therapy, we used whole-exome sequencing of non–small cell lung cancers treated with pembrolizumab, an antibody targeting programmed cell death-1 (PD-1). In two independent cohorts, higher nonsynonymous mutation burden in tumors was associated with improved objective response, durable clinical benefit, and progression-free survival. Efficacy also correlated with the molecular smoking signature, higher neoantigen burden, and DNA repair pathway mutations; each factor was also associated with mutation burden. In one responder, neoantigen-specific CD8+ T cell responses paralleled tumor regression, suggesting that anti–PD-1 therapy enhances neoantigen-specific T cell reactivity. Our results suggest that the genomic landscape of lung cancers shapes response to anti–PD-1 therapy.

Dedifferentiated liposarcoma (DDLS) is a rare but aggressive cancer with high recurrence and low response rates to targeted therapies. Increasing treatment efficacy may require combinations of targeted agents that counteract the effects of multiple abnormalities. To identify a possible multicomponent therapy, we performed a combinatorial drug screen in a DDLS-derived cell line and identified cyclin-dependent kinase 4 (CDK4) and insulin-like growth factor 1 receptor (IGF1R) as synergistic drug targets. We measured the phosphorylation of multiple proteins and cell viability in response to systematic drug combinations and derived computational models of the signaling network. These models predict that the observed synergy in reducing cell viability with CDK4 and IGF1R inhibitors depends on the activity of the AKT pathway. Experiments confirmed that combined inhibition of CDK4 and IGF1R cooperatively suppresses the activation of proteins within the AKT pathway. Consistent with these findings, synergistic reductions in cell viability were also found when combining CDK4 inhibition with inhibition of either AKT or epidermal growth factor receptor (EGFR), another receptor similar to IGF1R that activates AKT. Thus, network models derived from context-specific proteomic measurements of systematically perturbed cancer cells may reveal cancer-specific signaling mechanisms and aid in the design of effective combination therapies.

The Cancer Genome Atlas (TCGA) Research Network has profiled and analyzed large numbers of human tumors to discover molecular aberrations at the DNA, RNA, protein and epigenetic levels. The resulting rich data provide a major opportunity to develop an integrated picture of commonalities, differences and emergent themes across tumor lineages. The Pan-Cancer initiative compares the first 12 tumor types profiled by TCGA. Analysis of the molecular aberrations and their functional roles across tumor types will teach us how to extend therapies effective in one cancer type to others with a similar genomic profile., National Cancer Institute (U.S.), National Human Genome Research Institute (U.S.)

Systematic and quantitative analysis of protein phosphorylation is revealing dynamic regulatory networks underlying cellular responses to environmental cues. However, matching these sites to the kinases that phosphorylate them and the phosphorylation-dependent binding domains that may subsequently bind to them remains a challenge. NetPhorest is an atlas of consensus sequence motifs that covers 179 kinases and 104 phosphorylation-dependent binding domains [Src homology 2 (SH2), phosphotyrosine binding (PTB), BRCA1 C-terminal (BRCT), WW, and 14–3–3]. The atlas reveals new aspects of signaling systems, including the observation that tyrosine kinases mutated in cancer have lower specificity than their non-oncogenic relatives. The resource is maintained by an automated pipe line, which uses phylogenetic trees to structure the currently available in vivo and in vitro data to derive probabilistic sequence models of linear motifs. The atlas is available as a community resource (http://netphorest.info).