Cancer genomes contain large numbers of somatic mutations but few of these mutations drive tumor development. Current approaches either identify driver genes on the basis of mutational recurrence or ...approximate the functional consequences of nonsynonymous mutations by using bioinformatic scores. Passenger mutations are enriched in characteristic nucleotide contexts, whereas driver mutations occur in functional positions, which are not necessarily surrounded by a particular nucleotide context. We observed that mutations in contexts that deviate from the characteristic contexts around passenger mutations provide a signal in favor of driver genes. We therefore developed a method that combines this feature with the signals traditionally used for driver-gene identification. We applied our method to whole-exome sequencing data from 11,873 tumor-normal pairs and identified 460 driver genes that clustered into 21 cancer-related pathways. Our study provides a resource of driver genes across 28 tumor types with additional driver genes identified according to mutations in unusual nucleotide contexts.
Small cell lung carcinoma (SCLC) is highly mutated, yet durable response to immune checkpoint blockade (ICB) is rare. SCLC also exhibits cellular plasticity, which could influence its immunobiology. ...Here we discover that a distinct subset of SCLC uniquely upregulates MHC I, enriching for durable ICB benefit.
modeling confirms epigenetic recovery of MHC I in SCLC following loss of neuroendocrine differentiation, which tracks with derepression of STING. Transient EZH2 inhibition expands these nonneuroendocrine cells, which display intrinsic innate immune signaling and basally restored antigen presentation. Consistent with these findings, murine nonneuroendocrine SCLC tumors are rejected in a syngeneic model, with clonal expansion of immunodominant effector CD8 T cells. Therapeutically, EZH2 inhibition followed by STING agonism enhances T-cell recognition and rejection of SCLC in mice. Together, these data identify MHC I as a novel biomarker of SCLC immune responsiveness and suggest novel immunotherapeutic approaches to co-opt SCLC's intrinsic immunogenicity. SIGNIFICANCE: SCLC is poorly immunogenic, displaying modest ICB responsiveness with rare durable activity. In profiling its plasticity, we uncover intrinsically immunogenic MHC I
subpopulations of nonneuroendocrine SCLC associated with durable ICB benefit. We also find that combined EZH2 inhibition and STING agonism uncovers this cell state, priming cells for immune rejection.
.
Despite continued widespread use, the genomic effects of cisplatin-based chemotherapy and implications for subsequent treatment are incompletely characterized. Here, we analyze whole exome sequencing ...of matched pre- and post-neoadjuvant cisplatin-based chemotherapy primary bladder tumor samples from 30 muscle-invasive bladder cancer patients. We observe no overall increase in tumor mutational burden post-chemotherapy, though a significant proportion of subclonal mutations are unique to the matched pre- or post-treatment tumor, suggesting chemotherapy-induced and/or spatial heterogeneity. We subsequently identify and validate a novel mutational signature in post-treatment tumors consistent with known characteristics of cisplatin damage and repair. We find that post-treatment tumor heterogeneity predicts worse overall survival, and further observe alterations in cell-cycle and immune checkpoint regulation genes in post-treatment tumors. These results provide insight into the clinical and genomic dynamics of tumor evolution with cisplatin-based chemotherapy, suggest mechanisms of clinical resistance, and inform development of clinically relevant biomarkers and trials of combination therapies.
The diversity of clinical tumor profiling approaches (small panels to whole exomes with matched or unmatched germline analysis) may engender uncertainty about their benefits and liabilities, ...particularly in light of reported germline false positives in tumor-only profiling and use of global mutational and/or neoantigen data. The goal of this study was to determine the impact of genomic analysis strategies on error rates and data interpretation across contexts and ancestries.
We modeled common tumor profiling modalities-large (n = 300 genes), medium (n = 48 genes), and small (n = 15 genes) panels-using clinical whole exomes (WES) from 157 patients with lung or colon adenocarcinoma. We created a tumor-only analysis algorithm to assess germline false positive rates, the impact of patient ancestry on tumor-only results, and neoantigen detection.
After optimizing a germline filtering strategy, the germline false positive rate with tumor-only large panel sequencing was 14 % (144/1012 variants). For patients whose tumor-only results underwent molecular pathologist review (n = 91), 50/54 (93 %) false positives were correctly interpreted as uncertain variants. Increased germline false positives were observed in tumor-only sequencing of non-European compared with European ancestry patients (p < 0.001; Fisher's exact) when basic germline filtering approaches were used; however, the ExAC database (60,706 germline exomes) mitigated this disparity (p = 0.53). Matched and unmatched large panel mutational load correlated with WES mutational load (r(2) = 0.99 and 0.93, respectively; p < 0.001). Neoantigen load also correlated (r(2) = 0.80; p < 0.001), though WES identified a broader spectrum of neoantigens. Small panels did not predict mutational or neoantigen load.
Large tumor-only targeted panels are sufficient for most somatic variant identification and mutational load prediction if paired with expanded germline analysis strategies and molecular pathologist review. Paired germline sequencing reduced overall false positive mutation calls and WES provided the most neoantigens. Without patient-matched germline data, large germline databases are needed to minimize false positive mutation calling and mitigate ethnic disparities.
Immune checkpoint inhibitors (ICIs) have minimal therapeutic effect in hormone receptor-positive (HR+ ) breast cancer. We present final overall survival (OS) results (n = 88) from a randomized phase ...2 trial of eribulin ± pembrolizumab for patients with metastatic HR+ breast cancer, computationally dissect genomic and/or transcriptomic data from pre-treatment tumors (n = 52) for molecular associations with efficacy, and identify cytokine changes differentiating response and ICI-related toxicity (n = 58). Despite no improvement in OS with combination therapy (hazard ratio 0.95, 95% CI 0.59-1.55, p = 0.84), immune infiltration and antigen presentation distinguished responding tumors, while tumor heterogeneity and estrogen signaling independently associated with resistance. Moreover, patients with ICI-related toxicity had lower levels of immunoregulatory cytokines. Broadly, we establish a framework for ICI response in HR+ breast cancer that warrants diagnostic and therapeutic validation. ClinicalTrials.gov Registration: NCT03051659.
Squamous cell carcinoma of the anal canal (ASCC) accounts for 2% to 4% of gastrointestinal malignancies in the United States and is increasing in incidence; however, genomic features of ASCC are ...incompletely characterized. Primary treatment of ASCC involves concurrent chemotherapy and radiation (CRT), but the mutational landscape of resistance to CRT is unknown. Here, we aim to compare mutational features of ASCC in the pre- and post-CRT setting.
We perform whole-exome sequencing of primary (
= 31) and recurrent (
= 30) ASCCs and correlate findings with clinical data. We compare genomic features of matched pre- and post-CRT tumors to identify genomic features of CRT response. Finally, we investigate the mutational underpinnings of an extraordinary ASCC response to immunotherapy.
We find that both primary and recurrent ASCC tumors harbor mutations in genes, such as
and
, that are also mutated in other HPV-associated cancers. Overall mutational burden was not significantly different in pre- versus post-CRT tumors, and several examples of shared clonal driver mutations were identified. In two cases, clonally related pre- and post-CRT tumors harbored distinct oncogenic driver mutations in the same cancer gene (
or
). A patient with recurrent disease achieved an exceptional response to anti-programmed death (PD-1) therapy, and genomic dissection revealed high mutational burden and predicted neoantigen load.
We perform comprehensive mutational analysis of ASCC and characterize mutational features associated with CRT. Although many primary and recurrent tumors share driver events, we identify several unique examples of clonal evolution in response to treatment.
.
This study highlights the significance of genetic ancestry in renal cell carcinoma (RCC) risk assessment, urges carful interpretation of CHEK2’s role in RCC predisposition, and uncovers understudied ...inherited mutations potentially increasing RCC risk.
Previous germline studies on renal cell carcinoma (RCC) have usually pooled clear and non–clear cell RCCs and have not adequately accounted for population stratification, which might have led to an inaccurate estimation of genetic risk. Here, we aim to analyze the major germline drivers of RCC risk and clinically relevant but underexplored germline variant types.
We first characterized germline pathogenic variants (PVs), cryptic splice variants, and copy number variants (CNVs) in 1436 unselected RCC patients. To evaluate the enrichment of PVs in RCC, we conducted a case-control study of 1356 RCC patients ancestry matched with 16512 cancer-free controls using approaches accounting for population stratification and histological subtypes, followed by characterization of secondary somatic events.
Clear cell RCC patients (n = 976) exhibited a significant burden of PVs in VHL compared with controls (odds ratio OR: 39.1, p = 4.95e-05). Non–clear cell RCC patients (n = 380) carried enrichment of PVs in FH (OR: 77.9, p = 1.55e-08) and MET (OR: 1.98e11, p = 2.07e-05). In a CHEK2-focused analysis with European participants, clear cell RCC (n = 906) harbored nominal enrichment of low-penetrance CHEK2 variants—p.Ile157Thr (OR: 1.84, p = 0.049) and p.Ser428Phe (OR: 5.20, p = 0.045), while non–clear cell RCC (n = 295) exhibited nominal enrichment of CHEK2 loss of function PVs (OR: 3.51, p = 0.033). Patients with germline PVs in FH, MET, and VHL exhibited significantly earlier age of cancer onset than patients without germline PVs (mean: 46.0 vs 60.2 yr, p < 0.0001), and more than half had secondary somatic events affecting the same gene (n = 10/15, 66.7%). Conversely, CHEK2 PV carriers exhibited a similar age of onset to patients without germline PVs (mean: 60.1 vs 60.2 yr, p = 0.99), and only 30.4% carried somatic events in CHEK2 (n = 7/23). Finally, pathogenic germline cryptic splice variants were identified in SDHA and TSC1, and pathogenic germline CNVs were found in 18 patients, including CNVs in FH, SDHA, and VHL.
This analysis supports the existing link between several RCC risk genes and RCC risk manifesting in earlier age of onset. It calls for caution when assessing the role of CHEK2 due to the burden of founder variants with varying population frequency. It also broadens the definition of the RCC germline landscape of pathogenicity to incorporate previously understudied types of germline variants.
In this study, we carefully compared the frequency of rare inherited mutations with a focus on patients’ genetic ancestry. We discovered that subtle variations in genetic background may confound a case-control analysis, especially in evaluating the cancer risk associated with specific genes, such as CHEK2. We also identified previously less explored forms of rare inherited mutations, which could potentially increase the risk of kidney cancer.
Abstract
Motivation
Large-scale sequencing studies have created a need to succinctly visualize genomic characteristics of patient cohorts linked to widely variable phenotypic information. This is ...often done by visualizing the co-occurrence of variants with comutation plots. Current tools lack the ability to create highly customizable and publication quality comutation plots from arbitrary user data.
Results
We developed CoMut, a stand-alone, object-oriented Python package that creates comutation plots from arbitrary input data, including categorical data, continuous data, bar graphs, side bar graphs and data that describes relationships between samples.
Availability and implementation
The CoMut package is open source and is available at https://github.com/vanallenlab/comut under the MIT License, along with documentation and examples. A no installation, easy-to-use implementation is available on Google Colab (see GitHub).