Key Parameters of Tumor Epitope Immunogenicity Revealed Through a Consortium Approach Improve Neoantigen Prediction.

Many approaches to identify therapeutically relevant neoantigens couple tumor sequencing with bioinformatic algorithms and inferred rules of tumor epitope immunogenicity. However, there are no reference data to compare these approaches, and the parameters governing tumor epitope immunogenicity remain unclear. Here, we assembled a global consortium wherein each participant predicted immunogenic epitopes from shared tumor sequencing data. 608 epitopes were subsequently assessed for T cell binding in patient-matched samples. By integrating peptide features associated with presentation and recognition, we developed a model of tumor epitope immunogenicity that filtered out 98% of non-immunogenic peptides with a precision above 0.70. Pipelines prioritizing model features had superior performance, and pipeline alterations leveraging them improved prediction performance. These findings were validated in an independent cohort of 310 epitopes prioritized from tumor sequencing data and assessed for T cell binding. This data resource enables identification of parameters underlying effective anti-tumor immunity and is available to the research community.

Non-viral precision T cell receptor replacement for personalized cell therapy.

T cell receptors (TCRs) enable T cells to specifically recognize mutations in cancer cells1-3. Here we developed a clinical-grade approach based on CRISPR-Cas9 non-viral precision genome-editing to simultaneously knockout the two endogenous TCR genes TRAC (which encodes TCRα) and TRBC (which encodes TCRß). We also inserted into the TRAC locus two chains of a neoantigen-specific TCR (neoTCR) isolated from circulating T cells of patients. The neoTCRs were isolated using a personalized library of soluble predicted neoantigen-HLA capture reagents. Sixteen patients with different refractory solid cancers received up to three distinct neoTCR transgenic cell products. Each product expressed a patient-specific neoTCR and was administered in a cell-dose-escalation, first-in-human phase I clinical trial ( NCT03970382 ). One patient had grade 1 cytokine release syndrome and one patient had grade 3 encephalitis. All participants had the expected side effects from the lymphodepleting chemotherapy. Five patients had stable disease and the other eleven had disease progression as the best response on the therapy. neoTCR transgenic T cells were detected in tumour biopsy samples after infusion at frequencies higher than the native TCRs before infusion. This study demonstrates the feasibility of isolating and cloning multiple TCRs that recognize mutational neoantigens. Moreover, simultaneous knockout of the endogenous TCR and knock-in of neoTCRs using single-step, non-viral precision genome-editing are achieved. The manufacture of neoTCR engineered T cells at clinical grade, the safety of infusing up to three gene-edited neoTCR T cell products and the ability of the transgenic T cells to traffic to the tumours of patients are also demonstrated.

CIViCdb 2022: evolution of an open-access cancer variant interpretation knowledgebase.

CIViC (Clinical Interpretation of Variants in Cancer; civicdb.org) is a crowd-sourced, public domain knowledgebase composed of literature-derived evidence characterizing the clinical utility of cancer variants. As clinical sequencing becomes more prevalent in cancer management, the need for cancer variant interpretation has grown beyond the capability of any single institution. CIViC contains peer-reviewed, published literature curated and expertly-moderated into structured data units (Evidence Items) that can be accessed globally and in real time, reducing barriers to clinical variant knowledge sharing. We have extended CIViC's functionality to support emergent variant interpretation guidelines, increase interoperability with other variant resources, and promote widespread dissemination of structured curated data. To support the full breadth of variant interpretation from basic to translational, including integration of somatic and germline variant knowledge and inference of drug response, we have enabled curation of three new Evidence Types (Predisposing, Oncogenic and Functional). The growing CIViC knowledgebase has over 300 contributors and distributes clinically-relevant cancer variant data currently representing >3200 variants in >470 genes from >3100 publications.

ERK mediates interferon gamma-induced melanoma cell death.

BACKGROUND: Interferon-gamma (IFNγ) exerts potent growth inhibitory effects on a wide range of cancer cells through unknown signaling pathways. We pursued complementary screening approaches to characterize the growth inhibition pathway. METHODS: We performed chemical genomics and whole genome targeting CRISPR/Cas9 screens using patient-derived melanoma lines to uncover essential nodes in the IFNγ-mediated growth inhibition pathway. We used transcriptomic profiling to identify cell death pathways activated upon IFNγ exposure. Live imaging experiments coupled with apoptosis assays confirmed the involvement of these pathways in IFNγ-mediated cell death. RESULTS: We show that IFNγ signaling activated ERK. Blocking ERK activation rescued IFNγ-mediated apoptosis in 17 of 23 (~ 74%) cell lines representing BRAF, NRAS, NF1 mutant, and triple wild type subtypes of cutaneous melanoma. ERK signaling induced a stress response, ultimately leading to apoptosis through the activity of DR5 and NOXA proteins. CONCLUSIONS: Our results provide a new understanding of the IFNγ growth inhibition pathway, which will be crucial in defining mechanisms of immunotherapy response and resistance.

Standard operating procedure for somatic variant refinement of sequencing data with paired tumor and normal samples.

PURPOSE: Following automated variant calling, manual review of aligned read sequences is required to identify a high-quality list of somatic variants. Despite widespread use in analyzing sequence data, methods to standardize manual review have not been described, resulting in high inter- and intralab variability. METHODS: This manual review standard operating procedure (SOP) consists of methods to annotate variants with four different calls and 19 tags. The calls indicate a reviewer's confidence in each variant and the tags indicate commonly observed sequencing patterns and artifacts that inform the manual review call. Four individuals were asked to classify variants prior to, and after, reading the SOP and accuracy was assessed by comparing reviewer calls with orthogonal validation sequencing. RESULTS: After reading the SOP, average accuracy in somatic variant identification increased by 16.7% (p value = 0.0298) and average interreviewer agreement increased by 12.7% (p value < 0.001). Manual review conducted after reading the SOP did not significantly increase reviewer time. CONCLUSION: This SOP supports and enhances manual somatic variant detection by improving reviewer accuracy while reducing the interreviewer variability for variant calling and annotation.

DGIdb 2.0: mining clinically relevant drug-gene interactions.

The Drug-Gene Interaction Database (DGIdb, www.dgidb.org) is a web resource that consolidates disparate data sources describing drug-gene interactions and gene druggability. It provides an intuitive graphical user interface and a documented application programming interface (API) for querying these data. DGIdb was assembled through an extensive manual curation effort, reflecting the combined information of twenty-seven sources. For DGIdb 2.0, substantial updates have been made to increase content and improve its usefulness as a resource for mining clinically actionable drug targets. Specifically, nine new sources of drug-gene interactions have been added, including seven resources specifically focused on interactions linked to clinical trials. These additions have more than doubled the overall count of drug-gene interactions. The total number of druggable gene claims has also increased by 30%. Importantly, a majority of the unrestricted, publicly-accessible sources used in DGIdb are now automatically updated on a weekly basis, providing the most current information for these sources. Finally, a new web view and API have been developed to allow searching for interactions by drug identifiers to complement existing gene-based search functionality. With these updates, DGIdb represents a comprehensive and user friendly tool for mining the druggable genome for precision medicine hypothesis generation.

GenVisR: Genomic Visualizations in R.

UNLABELLED: Visualizing and summarizing data from genomic studies continues to be a challenge. Here, we introduce the GenVisR package to addresses this challenge by providing highly customizable, publication-quality graphics focused on cohort level genome analyses. GenVisR provides a rapid and easy-to-use suite of genomic visualization tools, while maintaining a high degree of flexibility by leveraging the abilities of ggplot2 and Bioconductor. AVAILABILITY AND IMPLEMENTATION: GenVisR is an R package available via Bioconductor (https://bioconductor.org/packages/GenVisR) under GPLv3. Support is available via GitHub (https://github.com/griffithlab/GenVisR/issues) and the Bioconductor support website. CONTACTS: obigriffith@wustl.edu or mgriffit@wustl.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

RNAi screen in Drosophila larvae identifies histone deacetylase 3 as a positive regulator of the hsp70 heat shock gene expression during heat shock.

The transcription regulation of the Drosophila hsp70 gene is a complex process that involves the regulation of multiple steps, including the establishment of paused Pol II and release of Pol II into elongation upon heat shock activation. While the major players involved in the regulation of gene expression have been studied in detail, additional factors involved in this process continue to be discovered. To identify factors involved in hsp70 expression, we developed a screen that capitalizes on a visual assessment of heat shock activation using a hsp70-beta galactosidase reporter and publicly available RNAi fly lines to deplete candidate proteins. We validated the screen by showing that the depletion of HSF, CycT, Cdk9, Nurf 301, or ELL prevented the full induction of hsp70 by heat shock. Our screen also identified the histone deacetylase HDAC3 and its associated protein SMRTER as positive regulators of hsp70 activation. Additionally, we show that HDAC3 and SMRTER contribute to hsp70 gene expression at a step subsequent to HSF-mediated activation and release of the paused Pol II that resides at the promoter prior to heat shock induction.

Integrative systems biology reveals NKG2A-biased immune responses correlate with protection in infectious disease, autoimmune disease, and cancer.

Infection, autoimmunity, and cancer are principal human health challenges of the 21st century. Often regarded as distinct ends of the immunological spectrum, recent studies hint at potential overlap between these diseases. For example, inflammation can be pathogenic in infection and autoimmunity. T resident memory (TRM) cells can be beneficial in infection and cancer. However, these findings are limited by size and scope; exact immunological factors shared across diseases remain elusive. Here, we integrate large-scale deeply clinically and biologically phenotyped human cohorts of 526 patients with infection, 162 with lupus, and 11,180 with cancer. We identify an NKG2A+ immune bias as associative with protection against disease severity, mortality, and autoimmune/post-acute chronic disease. We reveal that NKG2A+ CD8+ T cells correlate with reduced inflammation and increased humoral immunity and that they resemble TRM cells. Our results suggest NKG2A+ biases as a cross-disease factor of protection, supporting suggestions of immunological overlap between infection, autoimmunity, and cancer.

Opposing tumor-cell-intrinsic and -extrinsic roles of the IRF1 transcription factor in antitumor immunity.

Type I interferon (IFN-I) and IFN-γ foster antitumor immunity by facilitating T cell responses. Paradoxically, IFNs may promote T cell exhaustion by activating immune checkpoints. The downstream regulators of these disparate responses are incompletely understood. Here, we describe how interferon regulatory factor 1 (IRF1) orchestrates these opposing effects of IFNs. IRF1 expression in tumor cells blocks Toll-like receptor- and IFN-I-dependent host antitumor immunity by preventing interferon-stimulated gene (ISG) and effector programs in immune cells. In contrast, expression of IRF1 in the host is required for antitumor immunity. Mechanistically, IRF1 binds distinctly or together with STAT1 at promoters of immunosuppressive but not immunostimulatory ISGs in tumor cells. Overexpression of programmed cell death ligand 1 (PD-L1) in Irf1-/- tumors only partially restores tumor growth, suggesting multifactorial effects of IRF1 on antitumor immunity. Thus, we identify that IRF1 expression in tumor cells opposes host IFN-I- and IRF1-dependent antitumor immunity to facilitate immune escape and tumor growth.

Antitumor Immune Responses in B2M-Deficient Cancers.

ß2-microglobulin (B2M) is a critical component of the MHC class I molecule and is required to present tumor antigens to T cells. Its loss results in acquired resistance to immune checkpoint blockade (ICB) therapies. However, there have been well-documented cases of B2M-inactivated tumors responding to ICB, justifying investigation of how an antitumor immune response can be generated to tumors without surface MHC class I. We knocked out B2M in three murine models with varying baseline MHC class I expression and sensitivity to anti-programmed death receptor (PD-1) therapy and analyzed the immune responses. MC38 and YUMMER2.1 without B2M responded to anti-PD-1 alone or with an IL2 agonist, and this was mediated by CD4+ T cells and natural killer (NK) cells. The more aggressive B16 without B2M expression only partially responded to the IL2 agonist, and this was dependent on NK cells. When analyzing nearly 300 pretreatment biopsies from patients with melanoma receiving PD-1 blockade-based therapies, we found infrequent B2M mutations or homozygous loss but more frequent LOH or copy-number gains. B2M LOH was enriched in biopsies from patients without response to therapy, and these biopsies were more frequently infiltrated by activated NK cells. We conclude that in the absence of B2M, activation of CD4+ T cells and NK cells can mediate responses to murine models of PD-1 blockade therapy. In addition, in human melanoma, the intratumoral presence of activated NK cells upon partial B2M loss likely selects against tumor escape through low surface MHC class I expression.

A Pilot Study of Neoadjuvant Nivolumab, Ipilimumab, and Intralesional Oncolytic Virotherapy for HER2-negative Breast Cancer.

Purpose: Neoadjuvant combination immune checkpoint blockade and intralesional oncolytic virotherapy have the potential to activate antitumor responses in patients with breast cancer. Experimental Design: Eligibility for this pilot phase I trial included patients with localized HER2-negative breast cancer who received systemic nivolumab and ipilimumab and intratumor talimogene laherparepvec (T-VEC; NCT04185311). The primary objective was to evaluate the safety and adverse event profile of immunotherapy combined with T-VEC in patients with localized, HER2-negative breast cancer. Results: Six patients were enrolled, 4 having relapses after prior neoadjuvant chemotherapy and 2 who were previously untreated. Toxicities included 1 patient having grade 3 hypotension and type 1 diabetes mellitus, 3 patients with hypothyroidism, and all patients having constitutional symptoms known to be associated with the administration of T-VEC. One patient had a pathologic complete response, 3 patients had pathologic partial responses, 1 showed no significant response, and 1 had disease progression. Biopsies demonstrated increased immune cell infiltration in samples from patients who responded to therapy. Conclusions: This triple immunotherapy regimen provided responses in patients with advanced or relapsed HER2-negative breast cancer, at the expense of long-term toxicities. Significance: Systemic immune checkpoint blockade with a programmed death receptor 1 and a CTL antigen-4 blocking antibody, combined with intralesional oncolytic virotherapy, is a chemotherapy-free combination aimed at inducing an antitumor immune response locally and systemic immunity.

Prior anti-CTLA-4 therapy impacts molecular characteristics associated with anti-PD-1 response in advanced melanoma.

Immune checkpoint inhibitors (ICIs), including CTLA-4- and PD-1-blocking antibodies, can have profound effects on tumor immune cell infiltration that have not been consistent in biopsy series reported to date. Here, we analyze seven molecular datasets of samples from patients with advanced melanoma (N = 514) treated with ICI agents to investigate clinical, genomic, and transcriptomic features of anti-PD-1 response in cutaneous melanoma. We find that prior anti-CTLA-4 therapy is associated with differences in genomic, individual gene, and gene signatures in anti-PD-1 responders. Anti-CTLA-4-experienced melanoma tumors that respond to PD-1 blockade exhibit increased tumor mutational burden, inflammatory signatures, and altered cell cycle processes compared with anti-CTLA-4-naive tumors or anti-CTLA-4-experienced, anti-PD-1-nonresponsive melanoma tumors. We report a harmonized, aggregate resource and suggest that prior CTLA-4 blockade therapy is associated with marked differences in the tumor microenvironment that impact the predictive features of PD-1 blockade therapy response.

Integrated analysis of genomic and transcriptomic data for the discovery of splice-associated variants in cancer.

Somatic mutations within non-coding regions and even exons may have unidentified regulatory consequences that are often overlooked in analysis workflows. Here we present RegTools ( www.regtools.org ), a computationally efficient, free, and open-source software package designed to integrate somatic variants from genomic data with splice junctions from bulk or single cell transcriptomic data to identify variants that may cause aberrant splicing. We apply RegTools to over 9000 tumor samples with both tumor DNA and RNA sequence data. RegTools discovers 235,778 events where a splice-associated variant significantly increases the splicing of a particular junction, across 158,200 unique variants and 131,212 unique junctions. To characterize these somatic variants and their associated splice isoforms, we annotate them with the Variant Effect Predictor, SpliceAI, and Genotype-Tissue Expression junction counts and compare our results to other tools that integrate genomic and transcriptomic data. While many events are corroborated by the aforementioned tools, the flexibility of RegTools also allows us to identify splice-associated variants in known cancer drivers, such as TP53, CDKN2A, and B2M, and other genes.

Ipilimumab with or without nivolumab in PD-1 or PD-L1 blockade refractory metastatic melanoma: a randomized phase 2 trial.

In this randomized phase 2 trial, blockade of cytotoxic T-lymphocyte protein 4 (CTLA-4) with continuation of programmed death protein 1 (PD-1) blockade in patients with metastatic melanoma who had received front-line anti-PD-1 or therapy against programmed cell death 1 ligand 1 and whose tumors progressed was tested in comparison with CTLA-4 blockade alone. Ninety-two eligible patients were randomly assigned in a 3:1 ratio to receive the combination of ipilimumab and nivolumab, or ipilimumab alone. The primary endpoint was progression-free survival. Secondary endpoints included the difference in CD8 T cell infiltrate among responding and nonresponding tumors, objective response rate, overall survival and toxicity. The combination of nivolumab and ipilimumab resulted in a statistically significant improvement in progression-free survival over ipilimumab (hazard ratio = 0.63, 90% confidence interval (CI) = 0.41-0.97, one-sided P = 0.04). Objective response rates were 28% (90% CI = 19-38%) and 9% (90% CI = 2-25%), respectively (one-sided P = 0.05). Grade 3 or higher treatment-related adverse events occurred in 57% and 35% of patients, respectively, which is consistent with the known toxicity profile of these regimens. The change in intratumoral CD8 T cell density observed in the present analysis did not reach statistical significance to support the formal hypothesis tested as a secondary endpoint. In conclusion, primary resistance to PD-1 blockade therapy can be reversed in some patients with the combination of CTLA-4 and PD-1 blockade. Clinicaltrials.gov identifier: NCT03033576 .

Distinct clonal identities of B-ALLs arising after lenolidomide therapy for multiple myeloma.

Patients with multiple myeloma (MM) who are treated with lenalidomide rarely develop a secondary B-cell acute lymphoblastic leukemia (B-ALL). The clonal and biological relationship between these sequential malignancies is not yet clear. We identified 17 patients with MM treated with lenalidomide, who subsequently developed B-ALL. Patient samples were evaluated through sequencing, cytogenetics/fluorescence in situ hybridization (FISH), immunohistochemical (IHC) staining, and immunoglobulin heavy chain (IgH) clonality assessment. Samples were assessed for shared mutations and recurrently mutated genes. Through whole exome sequencing and cytogenetics/FISH analysis of 7 paired samples (MM vs matched B-ALL), no mutational overlap between samples was observed. Unique dominant IgH clonotypes between the tumors were observed in 5 paired MM/B-ALL samples. Across all 17 B-ALL samples, 14 (83%) had a TP53 variant detected. Three MM samples with sufficient sequencing depth (>500×) revealed rare cells (average of 0.6% variant allele frequency, or 1.2% of cells) with the same TP53 variant identified in the subsequent B-ALL sample. A lack of mutational overlap between MM and B-ALL samples shows that B-ALL developed as a second malignancy arising from a founding population of cells that likely represented unrelated clonal hematopoiesis caused by a TP53 mutation. The recurrent variants in TP53 in the B-ALL samples suggest a common path for malignant transformation that may be similar to that of TP53-mutant, treatment-related acute myeloid leukemia. The presence of rare cells containing TP53 variants in bone marrow at the initiation of lenalidomide treatment suggests that cellular populations containing TP53 variants expand in the presence of lenalidomide to increase the likelihood of B-ALL development.

Spatial profiling reveals association between WNT pathway activation and T-cell exclusion in acquired resistance of synovial sarcoma to NY-ESO-1 transgenic T-cell therapy.

BACKGROUND: Genetically engineered T-cell immunotherapies for adoptive cell transfer (ACT) have emerged as a promising form of cancer treatment, but many of these patients develop recurrent disease. Furthermore, delineating mechanisms of resistance may be challenging since the analysis of bulk tumor profiling can be complicated by spatial heterogeneity. METHODS: Tumor samples were collected from a patient with synovial sarcoma who developed acquired resistance to ACT targeting NY-ESO-1. Biopsies (primary, progressive metastasis, and recurrence) were subjected to bulk tumor DNA and RNA sequencing, as well as high-dimensional spatial profiling of RNA and protein targets. Untreated and progressive lesions were compared with identified patterns associated with acquired resistance to ACT. RESULTS: Gene expression patterns due to immune activity and infiltration were diluted in bulk tumor sequencing. The metastasis was enriched for tumor regions with increased CTNNB1 (encoding beta-catenin), which were negatively associated with the expression of T-cell surface proteins and antigen presentation machinery. Spatial profiling was most highly concordant with bulk sequencing in the lesions with decreased spatial heterogeneity. CONCLUSIONS: Complementary use of bulk and spatial profiling enables more accurate interrogation of tumor specimens, particularly to address complex questions regarding immunotherapeutic mechanisms. Our study uses this approach to demonstrate a mechanism of T-cell exclusion and resistance to cellular immunotherapy in synovial sarcoma.

Remodeling of the tumor microenvironment through PAK4 inhibition sensitizes tumors to immune checkpoint blockade.

PAK4 inhibition can sensitize tumors to immune checkpoint blockade (ICB) therapy, however, the underlying mechanisms remain unclear. We report that PAK4 inhibition reverses immune cell exclusion by increasing the infiltration of CD8 T cells and CD103+ dendritic cells (DCs), a specific type of DCs that excel at cross-presenting tumor antigens and constitute a source of CXCL10. Interestingly, in melanoma clinical datasets, PAK4 expression levels negatively correlate with the presence of CCL21, the ligand for CCR7 expressed in CD103+ DCs. Furthermore, we extensively characterized the transcriptome of PAK4 knock out (KO) tumors, in vitro and in vivo, and established the importance of PAK4 expression in the regulation of the extracellular matrix, which can facilitate immune cell infiltration. Comparison between PAK4 wild type (WT) and KO anti-PD-1 treated tumors revealed how PAK4 deletion sensitizes tumors to ICB from a transcriptomic perspective. In addition, we validated genetically and pharmacologically that inhibition of PAK4 kinase activity is sufficient to improve anti-tumor efficacy of anti-PD-1 blockade in multiple melanoma mouse models. Therefore, this study provides novel insights into the mechanism of action of PAK4 inhibition and provides the foundation for a new treatment strategy that aims to overcome resistance to PD-1 blockade by combining anti-PD-1 with a small molecule PAK4 kinase inhibitor.

Checkpoint blockade-induced CD8+ T cell differentiation in head and neck cancer responders.

BACKGROUND: Immune checkpoint blockade (ICB) response in recurrent/metastatic head and neck squamous cell carcinoma (HNSCC) is limited to 15%-20% of patients and underpinnings of resistance remain undefined. METHODS: Starting with an anti-PD1 sensitive murine HNSCC cell line, we generated an isogenic anti-PD1 resistant model. Mass cytometry was used to delineate tumor microenvironments of both sensitive parental murine oral carcinoma (MOC1) and resistant MOC1esc1 tumors. To examine heterogeneity and clonal dynamics of tumor infiltrating lymphocytes (TILs), we applied paired single-cell RNA and TCR sequencing in three HNSCC models. RESULTS: Anti-PD1 resistant MOC1esc1 line displayed a conserved cell intrinsic immune evasion signature. Immunoprofiling showed distinct baseline tumor microenvironments of MOC1 and MOC1esc1, as well as the remodeling of immune compartments on ICB in MOC1esc1 tumors. Single cell sequencing analysis identified several CD8 +TIL subsets including Tcf7 +Pd1- (naïve/memory-like), Tcf7 +Pd1+ (progenitor), and Tcf7-Pd1+ (differentiated effector). Mapping TCR shared fractions identified that successful anti-PD1 or anti-CTLA4 therapy-induced higher post-treatment T cell lineage transitions. CONCLUSIONS: These data highlight critical aspects of CD8 +TIL heterogeneity and differentiation and suggest facilitation of CD8 +TIL differentiation as a strategy to improve HNSCC ICB response.