Large-scale cancer sequencing studies of patient cohorts have statistically implicated many genes driving cancer growth and progression, and their identification has yielded substantial translational ...impact. However, a remaining challenge is to increase the resolution of driver prediction from the level of genes to mutations because mutation-level predictions are more closely aligned with the goal of precision cancer medicine. Here, we present CHASMplus, a computational method that is uniquely capable of identifying driver missense mutations, including those specific to a cancer type, as evidenced by significantly superior performance on diverse benchmarks. Applied to 8,657 tumor samples across 32 cancer types in The Cancer Genome Atlas (TCGA), CHASMplus identifies over 4,000 unique driver missense mutations in 240 genes, supporting a prominent role for rare driver mutations. We show which TCGA cancer types are likely to yield discovery of new driver missense mutations by additional sequencing, which has important implications for public policy.
Display omitted
•CHASMplus identifies driver missense mutations in a cancer-type-specific manner•Rare driver mutations are common in cancer when considered as a group•In some cancers, further sequencing will yield insight into rare driver mutations•Comprehensive resource of driver propensity for all possible missense mutations
Tokheim et al. introduce a computational approach to accurately separate driver from passenger mutations in cancer. Their analysis revealed that most driver mutations occur only in a few patients, presenting a challenge for precision medicine, and several cancer types will benefit from additional sequencing to identify these rare driver mutations.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ, ZRSKP
Immune checkpoint blockade (ICB) has demonstrated efficacy by reinvigorating immune cytotoxicity against tumors. However, the mechanisms underlying how ICB induces responses in a subset of patients ...remain unclear. Using bulk and single-cell transcriptomic cohorts of melanoma patients receiving ICB, we proposed a clustering model based on the expression of an antigen-presenting machinery (APM) signature consisting of 23 genes in a forward-selection manner. We characterized four APM clusters associated with distinct immune characteristics, cancer hallmarks, and patient prognosis in melanoma. The model predicts differential regulation of APM genes during ICB, which shaped ICB responsiveness. Surprisingly, while immunogenically hot tumors with high baseline APM expression prior to treatment are correlated with a better response to ICB than cold tumors with low APM expression, a subset of hot tumors with the highest pre-ICB APM expression fail to upregulate APM expression during treatment. In addition, they undergo immunoediting and display infiltration of exhausted T cells. In comparison, tumors associated with the best patient prognosis demonstrate significant APM upregulation and immune infiltration following ICB. They also show infiltration of tissue-resident memory T cells, shaping prolonged antitumor immunity. Using only pre-treatment transcriptomic data, our model predicts the dynamic APM-mediated tumor-immune interactions in response to ICB and provides insights into the immune escape mechanisms in hot tumors that compromise the ICB efficacy. We highlight the prognostic value of APM expression in predicting immune response in chronic diseases.
Full text
Available for:
IZUM, KILJ, NUK, PILJ, PNG, SAZU, UL, UM, UPUK
3.
Evaluating the evaluation of cancer driver genes Tokheim, Collin J.; Papadopoulos, Nickolas; Kinzler, Kenneth W. ...
Proceedings of the National Academy of Sciences - PNAS,
12/2016, Volume:
113, Issue:
50
Journal Article
Peer reviewed
Open access
Sequencing has identified millions of somatic mutations in human cancers, but distinguishing cancer driver genes remains a major challenge. Numerous methods have been developed to identify driver ...genes, but evaluation of the performance of these methods is hindered by the lack of a gold standard, that is, bona fide driver gene mutations. Here, we establish an evaluation framework that can be applied to driver gene prediction methods. We used this framework to compare the performance of eight such methods. One of these methods, described here, incorporated a machine-learning–based ratiometric approach. We show that the driver genes predicted by each of the eight methods vary widely. Moreover, the P values reported by several of the methods were inconsistent with the uniform values expected, thus calling into question the assumptions that were used to generate them. Finally, we evaluated the potential effects of unexplained variability in mutation rates on false-positive driver gene predictions. Our analysis points to the strengths and weaknesses of each of the currently available methods and offers guidance for improving them in the future.
Full text
Available for:
BFBNIB, NMLJ, NUK, PNG, SAZU, UL, UM, UPUK
Accurately predicting regulatory sequences and enhancers in entire genomes is an important but difficult problem, especially in large vertebrate genomes. With the advent of ChIP-seq technology, ...experimental detection of genome-wide EP300/CREBBP bound regions provides a powerful platform to develop predictive tools for regulatory sequences and to study their sequence properties. Here, we develop a support vector machine (SVM) framework which can accurately identify EP300-bound enhancers using only genomic sequence and an unbiased set of general sequence features. Moreover, we find that the predictive sequence features identified by the SVM classifier reveal biologically relevant sequence elements enriched in the enhancers, but we also identify other features that are significantly depleted in enhancers. The predictive sequence features are evolutionarily conserved and spatially clustered, providing further support of their functional significance. Although our SVM is trained on experimental data, we also predict novel enhancers and show that these putative enhancers are significantly enriched in both ChIP-seq signal and DNase I hypersensitivity signal in the mouse brain and are located near relevant genes. Finally, we present results of comparisons between other EP300/CREBBP data sets using our SVM and uncover sequence elements enriched and/or depleted in the different classes of enhancers. Many of these sequence features play a role in specifying tissue-specific or developmental-stage-specific enhancer activity, but our results indicate that some features operate in a general or tissue-independent manner. In addition to providing a high confidence list of enhancer targets for subsequent experimental investigation, these results contribute to our understanding of the general sequence structure of vertebrate enhancers.
Cooperative dysregulation of gene sequence and expression may contribute to cancer formation and progression. The Cancer Genome Atlas (TCGA) Network recently catalogued gene sequence and expression ...data for a collection of glioblastoma multiforme (GBM) tumors. We developed an automated, model-free method to rapidly and exhaustively examine the correlation among somatic mutation and gene expression and interrogated 149 GBM tumor samples from the TCGA. The method identified 41 genes whose mutation status is highly correlated with drastic changes in the expression (z-score ± 2.0), across tumor samples, of other genes. Some of the 41 genes have been previously implicated in GBM pathogenesis (e.g., NF1, TP53, RB1, and IDH1) and others, while implicated in cancer, had not previously been highlighted in studies using TCGA data (e.g., SYNE1, KLF6, FGFR4, and EPHB4). The method also predicted that known oncogenes and tumor suppressors participate in GBM via drastic over- and underexpression, respectively. In addition, the method identified a known synthetic lethal interaction between TP53 and PLK1, other potential synthetic lethal interactions with TP53, and correlations between IDH1 mutation status and the overexpression of known GBM survival genes.
High-grade serous ovarian carcinoma (HGSOC) is the most frequent type of ovarian cancer and has a poor outcome. It has been proposed that fallopian tube cancers may be precursors of HGSOC but ...evolutionary evidence for this hypothesis has been limited. Here, we perform whole-exome sequence and copy number analyses of laser capture microdissected fallopian tube lesions (p53 signatures, serous tubal intraepithelial carcinomas (STICs), and fallopian tube carcinomas), ovarian cancers, and metastases from nine patients. The majority of tumor-specific alterations in ovarian cancers were present in STICs, including those affecting TP53, BRCA1, BRCA2 or PTEN. Evolutionary analyses reveal that p53 signatures and STICs are precursors of ovarian carcinoma and identify a window of 7 years between development of a STIC and initiation of ovarian carcinoma, with metastases following rapidly thereafter. Our results provide insights into the etiology of ovarian cancer and have implications for prevention, early detection and therapeutic intervention of this disease.
Immune checkpoint inhibitors have shown significant therapeutic responses against tumors containing increased mutation-associated neoantigen load. We have examined the evolving landscape of tumor ...neoantigens during the emergence of acquired resistance in patients with non-small cell lung cancer after initial response to immune checkpoint blockade with anti-PD-1 or anti-PD-1/anti-CTLA-4 antibodies. Analyses of matched pretreatment and resistant tumors identified genomic changes resulting in loss of 7 to 18 putative mutation-associated neoantigens in resistant clones. Peptides generated from the eliminated neoantigens elicited clonal T-cell expansion in autologous T-cell cultures, suggesting that they generated functional immune responses. Neoantigen loss occurred through elimination of tumor subclones or through deletion of chromosomal regions containing truncal alterations, and was associated with changes in T-cell receptor clonality. These analyses provide insight into the dynamics of mutational landscapes during immune checkpoint blockade and have implications for the development of immune therapies that target tumor neoantigens.
Acquired resistance to immune checkpoint therapy is being recognized more commonly. This work demonstrates for the first time that acquired resistance to immune checkpoint blockade can arise in association with the evolving landscape of mutations, some of which encode tumor neoantigens recognizable by T cells. These observations imply that widening the breadth of neoantigen reactivity may mitigate the development of acquired resistance.
.
Despite the initial successes of immunotherapy, there is an urgent clinical need for molecular assays that identify patients more likely to respond. Here, we report that ultrasensitive measures of ...circulating tumor DNA (ctDNA) and T-cell expansion can be used to assess responses to immune checkpoint blockade in metastatic lung cancer patients (
= 24). Patients with clinical response to therapy had a complete reduction in ctDNA levels after initiation of therapy, whereas nonresponders had no significant changes or an increase in ctDNA levels. Patients with initial response followed by acquired resistance to therapy had an initial drop followed by recrudescence in ctDNA levels. Patients without a molecular response had shorter progression-free and overall survival compared with molecular responders 5.2 vs. 14.5 and 8.4 vs. 18.7 months; HR 5.36; 95% confidence interval (CI), 1.57-18.35;
= 0.007 and HR 6.91; 95% CI, 1.37-34.97;
= 0.02, respectively, which was detected on average 8.7 weeks earlier and was more predictive of clinical benefit than CT imaging. Expansion of T cells, measured through increases of T-cell receptor productive frequencies, mirrored ctDNA reduction in response to therapy. We validated this approach in an independent cohort of patients with early-stage non-small cell lung cancer (
= 14), where the therapeutic effect was measured by pathologic assessment of residual tumor after anti-PD1 therapy. Consistent with our initial findings, early ctDNA dynamics predicted pathologic response to immune checkpoint blockade. These analyses provide an approach for rapid determination of therapeutic outcomes for patients treated with immune checkpoint inhibitors and have important implications for the development of personalized immune targeted strategies.
Rapid and sensitive detection of circulating tumor DNA dynamic changes and T-cell expansion can be used to guide immune targeted therapy for patients with lung cancer.
.
Major efforts to sequence cancer genomes are now occurring throughout the world. Though the emerging data from these studies are illuminating, their reconciliation with epidemiologic and clinical ...observations poses a major challenge. In the current study, we provide a mathematical model that begins to address this challenge. We model tumors as a discrete time branching process that starts with a single driver mutation and proceeds as each new driver mutation leads to a slightly increased rate of clonal expansion. Using the model, we observe tremendous variation in the rate of tumor development—providing an understanding of the heterogeneity in tumor sizes and development times that have been observed by epidemiologists and clinicians. Furthermore, the model provides a simple formula for the number of driver mutations as a function of the total number of mutations in the tumor. Finally, when applied to recent experimental data, the model allows us to calculate the actual selective advantage provided by typical somatic mutations in human tumors in situ. This selective advantage is surprisingly small—0.004 ± 0.0004—and has major implications for experimental cancer research.
Full text
Available for:
BFBNIB, NMLJ, NUK, PNG, SAZU, UL, UM, UPUK
Metastases are responsible for the majority of cancer-related deaths. Although genomic heterogeneity within primary tumors is associated with relapse, heterogeneity among treatment-naïve metastases ...has not been comprehensively assessed. We analyzed sequencing data for 76 untreated metastases from 20 patients and inferred cancer phylogenies for breast, colorectal, endometrial, gastric, lung, melanoma, pancreatic, and prostate cancers. We found that within individual patients, a large majority of driver gene mutations are common to all metastases. Further analysis revealed that the driver gene mutations that were not shared by all metastases are unlikely to have functional consequences. A mathematical model of tumor evolution and metastasis formation provides an explanation for the observed driver gene homogeneity. Thus, single biopsies capture most of the functionally important mutations in metastases and therefore provide essential information for therapeutic decision-making.