Incorporating mutational heterogeneity to identify genes that are enriched for synonymous mutations in cancer.

BACKGROUND: Synonymous mutations, which change the DNA sequence but not the encoded protein sequence, can affect protein structure and function, mRNA maturation, and mRNA half-lives. The possibility that synonymous mutations might be enriched in cancer has been explored in several recent studies. However, none of these studies control for all three types of mutational heterogeneity (patient, histology, and gene) that are known to affect the accurate identification of non-synonymous cancer-associated genes. Our goal is to adopt the current standard for non-synonymous mutations in an investigation of synonymous mutations. RESULTS: Here, we create an algorithm, MutSigCVsyn, an adaptation of MutSigCV, to identify cancer-associated genes that are enriched for synonymous mutations based on a non-coding background model that takes into account the mutational heterogeneity across these levels. Using MutSigCVsyn, we first analyzed 2572 cancer whole-genome samples from the Pan-cancer Analysis of Whole Genomes (PCAWG) to identify non-synonymous cancer drivers as a quality control. Indicative of the algorithm accuracy we find that 58.6% of these candidate genes were also found in Cancer Census Gene (CGC) list, and 66.2% were found within the PCAWG cancer driver list. We then applied it to identify 30 putative cancer-associated genes that are enriched for synonymous mutations within the same samples. One of the promising gene candidates is the B cell lymphoma 2 (BCL-2) gene. BCL-2 regulates apoptosis by antagonizing the action of proapoptotic BCL-2 family member proteins. The synonymous mutations in BCL2 are enriched in its anti-apoptotic domain and likely play a role in cancer cell proliferation. CONCLUSION: Our study introduces MutSigCVsyn, an algorithm that accounts for mutational heterogeneity at patient, histology, and gene levels, to identify cancer-associated genes that are enriched for synonymous mutations using whole genome sequencing data. We identified 30 putative candidate genes that will benefit from future experimental studies on the role of synonymous mutations in cancer biology.

Functional genomic analysis of adult and pediatric brain tumor isolates.

Background: Adult and pediatric tumors display stark differences in their mutation spectra and chromosome alterations. Here, we attempted to identify common and unique gene dependencies and their associated biomarkers among adult and pediatric tumor isolates using functional genetic lethal screens and computational modeling. Methods: We performed CRISRP-Cas9 lethality screens in two adult glioblastoma (GBM) tumor isolates and five pediatric brain tumor isolates representing atypical teratoid rhabdoid tumors (ATRT), diffuse intrinsic pontine glioma, GBM, and medulloblastoma. We then integrated the screen results with machine learning-based gene-dependency models generated from data from >900 cancer cell lines. Results: We found that >50% of candidate dependencies of 280 identified were shared between adult GBM tumors and individual pediatric tumor isolates. 68% of screen hits were found as nodes in our network models, along with shared and tumor-specific predictors of gene dependencies. We investigated network predictors associated with ADAR, EFR3A, FGFR1 (pediatric-specific), and SMARCC2 (ATRT-specific) gene dependency among our tumor isolates. Conclusions: The results suggest that, despite harboring disparate genomic signatures, adult and pediatric tumor isolates share a preponderance of genetic dependences. Further, combining data from primary brain tumor lethality screens with large cancer cell line datasets produced valuable insights into biomarkers of gene dependency, even for rare cancers. Importance of the Study: Our results demonstrate that large cancer cell lines data sets can be computationally mined to identify known and novel gene dependency relationships in adult and pediatric human brain tumor isolates. Gene dependency networks and lethality screen results represent a key resource for neuro-oncology and cancer research communities. We also highlight some of the challenges and limitations of this approach.

Genomic and experimental evidence that ALKATI does not predict single agent sensitivity to ALK inhibitors.

Genomic data can facilitate personalized treatment decisions by enabling therapeutic hypotheses in individual patients. Mutual exclusivity has been an empirically useful signal for identifying activating mutations that respond to single agent targeted therapies. However, a low mutation frequency can underpower this signal for rare variants. We develop a resampling based method for the direct pairwise comparison of conditional selection between sets of gene pairs. We apply this method to a transcript variant of anaplastic lymphoma kinase (ALK) in melanoma, termed ALKATI that was suggested to predict sensitivity to ALK inhibitors and we find that it is not mutually exclusive with key melanoma oncogenes. Furthermore, we find that ALKATI is not likely to be sufficient for cellular transformation or growth, and it does not predict single agent therapeutic dependency. Our work strongly disfavors the role of ALKATI as a targetable oncogenic driver that might be sensitive to single agent ALK treatment.