pVACview: an interactive visualization tool for efficient neoantigen prioritization and selection.

Background: Neoantigen targeting therapies including personalized vaccines have shown promise in the treatment of cancers, particularly when used in combination with checkpoint blockade therapy. At least 100 clinical trials involving these therapies are underway globally. Accurate identification and prioritization of neoantigens is highly relevant to designing these trials, predicting treatment response, and understanding mechanisms of resistance. With the advent of massively parallel DNA and RNA sequencing technologies, it is now possible to computationally predict neoantigens based on patient-specific variant information. However, numerous factors must be considered when prioritizing neoantigens for use in personalized therapies. Complexities such as alternative transcript annotations, various binding, presentation and immunogenicity prediction algorithms, and variable peptide lengths/registers all potentially impact the neoantigen selection process. There has been a rapid development of computational tools that attempt to account for these complexities. While these tools generate numerous algorithmic predictions for neoantigen characterization, results from these pipelines are difficult to navigate and require extensive knowledge of the underlying tools for accurate interpretation. This often leads to over-simplification of pipeline outputs to make them tractable, for example limiting prediction to a single RNA isoform or only summarizing the top ranked of many possible peptide candidates. In addition to variant detection, gene expression and predicted peptide binding affinities, recent studies have also demonstrated the importance of mutation location, allele-specific anchor locations, and variation of T-cell response to long versus short peptides. Due to the intricate nature and number of salient neoantigen features, presenting all relevant information to facilitate candidate selection for downstream applications is a difficult challenge that current tools fail to address. Results: We have created pVACview, the first interactive tool designed to aid in the prioritization and selection of neoantigen candidates for personalized neoantigen therapies including cancer vaccines. pVACview has a user-friendly and intuitive interface where users can upload, explore, select and export their neoantigen candidates. The tool allows users to visualize candidates across three different levels, including variant, transcript and peptide information. Conclusions: pVACview will allow researchers to analyze and prioritize neoantigen candidates with greater efficiency and accuracy in basic and translational settings The application is available as part of the pVACtools pipeline at pvactools.org and as an online server at pvacview.org.

Neoantigen landscape supports feasibility of personalized cancer vaccine for follicular lymphoma.

ABSTRACT: Personalized cancer vaccines designed to target neoantigens represent a promising new treatment paradigm in oncology. In contrast to classical idiotype vaccines, we hypothesized that "polyvalent" vaccines could be engineered for the personalized treatment of follicular lymphoma (FL) using neoantigen discovery by combined whole-exome sequencing (WES) and RNA sequencing (RNA-seq). Fifty-eight tumor samples from 57 patients with FL underwent WES and RNA-seq. Somatic and B-cell clonotype neoantigens were predicted and filtered to identify high-quality neoantigens. B-cell clonality was determined by the alignment of B-cell receptor (BCR) CDR3 regions from RNA-seq data, grouping at the protein level, and comparison with the BCR repertoire from healthy individuals using RNA-seq data. An average of 52 somatic mutations per patient (range, 2-172) were identified, and ≥2 (median, 15) high-quality neoantigens were predicted for 56 of 58 FL samples. The predicted neoantigen peptides were composed of missense mutations (77%), indels (9%), gene fusions (3%), and BCR sequences (11%). Building off of these preclinical analyses, we initiated a pilot clinical trial using personalized neoantigen vaccination combined with PD-1 blockade in patients with relapsed or refractory FL (#NCT03121677). Synthetic long peptide vaccines targeting predicted high-quality neoantigens were successfully synthesized for and administered to all 4 patients enrolled. Initial results demonstrate feasibility, safety, and potential immunologic and clinical responses. Our study suggests that a genomics-driven personalized cancer vaccine strategy is feasible for patients with FL, and this may overcome prior challenges in the field. This trial was registered at www.ClinicalTrials.gov as #NCT03121677.

Computational prediction of MHC anchor locations guides neoantigen identification and prioritization.

Neoantigens are tumor-specific peptide sequences resulting from sources such as somatic DNA mutations. Upon loading onto major histocompatibility complex (MHC) molecules, they can trigger recognition by T cells. Accurate neoantigen identification is thus critical for both designing cancer vaccines and predicting response to immunotherapies. Neoantigen identification and prioritization relies on correctly predicting whether the presenting peptide sequence can successfully induce an immune response. Because most somatic mutations are single-nucleotide variants, changes between wild-type and mutated peptides are typically subtle and require cautious interpretation. A potentially underappreciated variable in neoantigen prediction pipelines is the mutation position within the peptide relative to its anchor positions for the patient's specific MHC molecules. Whereas a subset of peptide positions are presented to the T cell receptor for recognition, others are responsible for anchoring to the MHC, making these positional considerations critical for predicting T cell responses. We computationally predicted anchor positions for different peptide lengths for 328 common HLA alleles and identified unique anchoring patterns among them. Analysis of 923 tumor samples shows that 6 to 38% of neoantigen candidates are potentially misclassified and can be rescued using allele-specific knowledge of anchor positions. A subset of anchor results were orthogonally validated using protein crystallography structures. Representative anchor trends were experimentally validated using peptide-MHC stability assays and competition binding assays. By incorporating our anchor prediction results into neoantigen prediction pipelines, we hope to formalize, streamline, and improve the identification process for relevant clinical studies.

Neoantigen-targeted CD8+ T cell responses with PD-1 blockade therapy.

Neoantigens are peptides derived from non-synonymous mutations presented by human leukocyte antigens (HLAs), which are recognized by antitumour T cells1-14. The large HLA allele diversity and limiting clinical samples have restricted the study of the landscape of neoantigen-targeted T cell responses in patients over their treatment course. Here we applied recently developed technologies15-17 to capture neoantigen-specific T cells from blood and tumours from patients with metastatic melanoma with or without response to anti-programmed death receptor 1 (PD-1) immunotherapy. We generated personalized libraries of neoantigen-HLA capture reagents to single-cell isolate the T cells and clone their T cell receptors (neoTCRs). Multiple T cells with different neoTCR sequences (T cell clonotypes) recognized a limited number of mutations in samples from seven patients with long-lasting clinical responses. These neoTCR clonotypes were recurrently detected over time in the blood and tumour. Samples from four patients with no response to anti-PD-1 also demonstrated neoantigen-specific T cell responses in the blood and tumour to a restricted number of mutations with lower TCR polyclonality and were not recurrently detected in sequential samples. Reconstitution of the neoTCRs in donor T cells using non-viral CRISPR-Cas9 gene editing demonstrated specific recognition and cytotoxicity to patient-matched melanoma cell lines. Thus, effective anti-PD-1 immunotherapy is associated with the presence of polyclonal CD8+ T cells in the tumour and blood specific for a limited number of immunodominant mutations, which are recurrently recognized over time.

pVACtools: A Computational Toolkit to Identify and Visualize Cancer Neoantigens.

Identification of neoantigens is a critical step in predicting response to checkpoint blockade therapy and design of personalized cancer vaccines. This is a cross-disciplinary challenge, involving genomics, proteomics, immunology, and computational approaches. We have built a computational framework called pVACtools that, when paired with a well-established genomics pipeline, produces an end-to-end solution for neoantigen characterization. pVACtools supports identification of altered peptides from different mechanisms, including point mutations, in-frame and frameshift insertions and deletions, and gene fusions. Prediction of peptide:MHC binding is accomplished by supporting an ensemble of MHC Class I and II binding algorithms within a framework designed to facilitate the incorporation of additional algorithms. Prioritization of predicted peptides occurs by integrating diverse data, including mutant allele expression, peptide binding affinities, and determination whether a mutation is clonal or subclonal. Interactive visualization via a Web interface allows clinical users to efficiently generate, review, and interpret results, selecting candidate peptides for individual patient vaccine designs. Additional modules support design choices needed for competing vaccine delivery approaches. One such module optimizes peptide ordering to minimize junctional epitopes in DNA vector vaccines. Downstream analysis commands for synthetic long peptide vaccines are available to assess candidates for factors that influence peptide synthesis. All of the aforementioned steps are executed via a modular workflow consisting of tools for neoantigen prediction from somatic alterations (pVACseq and pVACfuse), prioritization, and selection using a graphical Web-based interface (pVACviz), and design of DNA vector-based vaccines (pVACvector) and synthetic long peptide vaccines. pVACtools is available at http://www.pvactools.org.

Immunological ignorance is an enabling feature of the oligo-clonal T cell response to melanoma neoantigens.

The impact of intratumoral heterogeneity (ITH) and the resultant neoantigen landscape on T cell immunity are poorly understood. ITH is a widely recognized feature of solid tumors and poses distinct challenges related to the development of effective therapeutic strategies, including cancer neoantigen vaccines. Here, we performed deep targeted DNA sequencing of multiple metastases from melanoma patients and observed ubiquitous sharing of clonal and subclonal single nucleotide variants (SNVs) encoding putative HLA class I-restricted neoantigen epitopes. However, spontaneous antitumor CD8+ T cell immunity in peripheral blood and tumors was restricted to a few clonal neoantigens featuring an oligo-/monoclonal T cell-receptor (TCR) repertoire. Moreover, in various tumors of the 4 patients examined, no neoantigen-specific TCR clonotypes were identified despite clonal neoantigen expression. Mature dendritic cell (mDC) vaccination with tumor-encoded amino acid-substituted (AAS) peptides revealed diverse neoantigen-specific CD8+ T responses, each composed of multiple TCR clonotypes. Isolation of T cell clones by limiting dilution from tumor-infiltrating lymphocytes (TILs) permitted functional validation regarding neoantigen specificity. Gene transfer of TCRαß heterodimers specific for clonal neoantigens confirmed correct TCR clonotype assignments based on high-throughput TCRBV CDR3 sequencing. Our findings implicate immunological ignorance of clonal neoantigens as the basis for ineffective T cell immunity to melanoma and support the concept that therapeutic vaccination, as an adjunct to checkpoint inhibitor treatment, is required to increase the breadth and diversity of neoantigen-specific CD8+ T cells.

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.

Accounting for proximal variants improves neoantigen prediction.

Recent efforts to design personalized cancer immunotherapies use predicted neoantigens, but most neoantigen prediction strategies do not consider proximal (nearby) variants that alter the peptide sequence and may influence neoantigen binding. We evaluated somatic variants from 430 tumors to understand how proximal somatic and germline alterations change the neoantigenic peptide sequence and also affect neoantigen binding predictions. On average, 241 missense somatic variants were analyzed per sample. Of these somatic variants, 5% had one or more in-phase missense proximal variants. Without incorporating proximal variant correction for major histocompatibility complex class I neoantigen peptides, the overall false discovery rate (incorrect neoantigens predicted) and the false negative rate (strong-binding neoantigens missed) across peptides of lengths 8-11 were estimated as 0.069 (6.9%) and 0.026 (2.6%), respectively.

Author Correction: The prognostic effects of somatic mutations in ER-positive breast cancer.

The original version of this Article contained errors in the depiction of confidence intervals in the NF1 BCSS data illustrated in Figure 3b. These have now been corrected in both the PDF and HTML versions of the Article. The incorrect version of Figure 3b is presented in the associated Author Correction.

The prognostic effects of somatic mutations in ER-positive breast cancer.

Here we report targeted sequencing of 83 genes using DNA from primary breast cancer samples from 625 postmenopausal (UBC-TAM series) and 328 premenopausal (MA12 trial) hormone receptor-positive (HR+) patients to determine interactions between somatic mutation and prognosis. Independent validation of prognostic interactions was achieved using data from the METABRIC study. Previously established associations between MAP3K1 and PIK3CA mutations with luminal A status/favorable prognosis and TP53 mutations with Luminal B/non-luminal tumors/poor prognosis were observed, validating the methodological approach. In UBC-TAM, NF1 frame-shift nonsense (FS/NS) mutations were also a poor outcome driver that was validated in METABRIC. For MA12, poor outcome associated with PIK3R1 mutation was also reproducible. DDR1 mutations were strongly associated with poor prognosis in UBC-TAM despite stringent false discovery correction (q = 0.0003). In conclusion, uncommon recurrent somatic mutations should be further explored to create a more complete explanation of the highly variable outcomes that typifies ER+ breast cancer.

Breast Cancer Neoantigens Can Induce CD8+ T-Cell Responses and Antitumor Immunity.

Next-generation sequencing technologies have provided insights into the biology and mutational landscape of cancer. Here, we evaluate the relevance of cancer neoantigens in human breast cancers. Using patient-derived xenografts from three patients with advanced breast cancer (xenografts were designated as WHIM30, WHIM35, and WHIM37), we sequenced exomes of tumor and patient-matched normal cells. We identified 2,091 (WHIM30), 354 (WHIM35), and 235 (WHIM37) nonsynonymous somatic mutations. A computational analysis identified and prioritized HLA class I-restricted candidate neoantigens expressed in the dominant tumor clone. Each candidate neoantigen was evaluated using peptide-binding assays, T-cell cultures that measure the ability of CD8+ T cells to recognize candidate neoantigens, and preclinical models in which we measured antitumor immunity. Our results demonstrate that breast cancer neoantigens can be recognized by the immune system, and that human CD8+ T cells enriched for prioritized breast cancer neoantigens were able to protect mice from tumor challenge with autologous patient-derived xenografts. We conclude that next-generation sequencing and epitope-prediction strategies can identify and prioritize candidate neoantigens for immune targeting in breast cancer. Cancer Immunol Res; 5(7); 516-23. ©2017 AACR.

Temporally Distinct PD-L1 Expression by Tumor and Host Cells Contributes to Immune Escape.

Antibody blockade of programmed death-1 (PD-1) or its ligand, PD-L1, has led to unprecedented therapeutic responses in certain tumor-bearing individuals, but PD-L1 expression's prognostic value in stratifying cancer patients for such treatment remains unclear. Reports conflict on the significance of correlations between PD-L1 on tumor cells and positive clinical outcomes to PD-1/PD-L1 blockade. We investigated this issue using genomically related, clonal subsets from the same methylcholanthrene-induced sarcoma: a highly immunogenic subset that is spontaneously eliminated in vivo by adaptive immunity and a less immunogenic subset that forms tumors in immunocompetent mice, but is sensitive to PD-1/PD-L1 blockade therapy. Using CRISPR/Cas9-induced loss-of-function approaches and overexpression gain-of-function techniques, we confirmed that PD-L1 on tumor cells is key to promoting tumor escape. In addition, the capacity of PD-L1 to suppress antitumor responses was inversely proportional to tumor cell antigenicity. PD-L1 expression on host cells, particularly tumor-associated macrophages (TAM), was also important for tumor immune escape. We demonstrated that induction of PD-L1 on tumor cells was IFNγ-dependent and transient, but PD-L1 induction on TAMs was of greater magnitude, only partially IFNγ dependent, and was stable over time. Thus, PD-L1 expression on either tumor cells or host immune cells could lead to tumor escape from immune control, indicating that total PD-L1 expression in the immediate tumor microenvironment may represent a more accurate biomarker for predicting response to PD-1/PD-L1 blockade therapy, compared with monitoring PD-L1 expression on tumor cells alone. Cancer Immunol Res; 5(2); 106-17. ©2017 AACR.

Truncating Prolactin Receptor Mutations Promote Tumor Growth in Murine Estrogen Receptor-Alpha Mammary Carcinomas.

Estrogen receptor alpha-positive (ERα+) luminal tumors are the most frequent subtype of breast cancer. Stat1(-/-) mice develop mammary tumors that closely recapitulate the biological characteristics of this cancer subtype. To identify transforming events that contribute to tumorigenesis, we performed whole genome sequencing of Stat1(-/-) primary mammary tumors and matched normal tissues. This investigation identified somatic truncating mutations affecting the prolactin receptor (PRLR) in all tumor and no normal samples. Targeted sequencing confirmed the presence of these mutations in precancerous lesions, indicating that this is an early event in tumorigenesis. Functional evaluation of these heterozygous mutations in Stat1(-/-) mouse embryonic fibroblasts showed that co-expression of truncated and wild-type PRLR led to aberrant STAT3 and STAT5 activation downstream of the receptor, cellular transformation in vitro, and tumor formation in vivo. In conclusion, truncating mutations of PRLR promote tumor growth in a model of human ERα+ breast cancer and warrant further investigation.

Comprehensive genomic analysis reveals FLT3 activation and a therapeutic strategy for a patient with relapsed adult B-lymphoblastic leukemia.

The genomic events responsible for the pathogenesis of relapsed adult B-lymphoblastic leukemia (B-ALL) are not yet clear. We performed integrative analysis of whole-genome, whole-exome, custom capture, whole-transcriptome (RNA-seq), and locus-specific genomic assays across nine time points from a patient with primary de novo B-ALL. Comprehensive genome and transcriptome characterization revealed a dramatic tumor evolution during progression, yielding a tumor with complex clonal architecture at second relapse. We observed and validated point mutations in EP300 and NF1, a highly expressed EP300-ZNF384 gene fusion, a microdeletion in IKZF1, a focal deletion affecting SETD2, and large deletions affecting RB1, PAX5, NF1, and ETV6. Although the genome analysis revealed events of potential biological relevance, no clinically actionable treatment options were evident at the time of the second relapse. However, transcriptome analysis identified aberrant overexpression of the targetable protein kinase encoded by the FLT3 gene. Although the patient had refractory disease after salvage therapy for the second relapse, treatment with the FLT3 inhibitor sunitinib rapidly induced a near complete molecular response, permitting the patient to proceed to a matched-unrelated donor stem cell transplantation. The patient remains in complete remission more than 4 years later. Analysis of this patient's relapse genome revealed an unexpected, actionable therapeutic target that led to a specific therapy associated with a rapid clinical response. For some patients with relapsed or refractory cancers, this approach may indicate a novel therapeutic intervention that could alter outcome.

pVAC-Seq: A genome-guided in silico approach to identifying tumor neoantigens.

Cancer immunotherapy has gained significant momentum from recent clinical successes of checkpoint blockade inhibition. Massively parallel sequence analysis suggests a connection between mutational load and response to this class of therapy. Methods to identify which tumor-specific mutant peptides (neoantigens) can elicit anti-tumor T cell immunity are needed to improve predictions of checkpoint therapy response and to identify targets for vaccines and adoptive T cell therapies. Here, we present a flexible, streamlined computational workflow for identification of personalized Variant Antigens by Cancer Sequencing (pVAC-Seq) that integrates tumor mutation and expression data (DNA- and RNA-Seq). pVAC-Seq is available at https://github.com/griffithlab/pVAC-Seq .

Cancer Immunogenomics: Computational Neoantigen Identification and Vaccine Design.

The application of modern high-throughput genomics to the study of cancer genomes has exploded in the past few years, yielding unanticipated insights into the myriad and complex combinations of genomic alterations that lead to the development of cancers. Coincident with these genomic approaches have been computational analyses that are capable of multiplex evaluations of genomic data toward specific therapeutic end points. One such approach is called "immunogenomics" and is now being developed to interpret protein-altering changes in cancer cells in the context of predicted preferential binding of these altered peptides by the patient's immune molecules, specifically human leukocyte antigen (HLA) class I and II proteins. One goal of immunogenomics is to identify those cancer-specific alterations that are likely to elicit an immune response that is highly specific to the patient's cancer cells following stimulation by a personalized vaccine. The elements of such an approach are outlined herein and constitute an emerging therapeutic option for cancer patients.

Association Between Mutation Clearance After Induction Therapy and Outcomes in Acute Myeloid Leukemia.

IMPORTANCE: Tests that predict outcomes for patients with acute myeloid leukemia (AML) are imprecise, especially for those with intermediate risk AML. OBJECTIVES: To determine whether genomic approaches can provide novel prognostic information for adult patients with de novo AML. DESIGN, SETTING, AND PARTICIPANTS: Whole-genome or exome sequencing was performed on samples obtained at disease presentation from 71 patients with AML (mean age, 50.8 years) treated with standard induction chemotherapy at a single site starting in March 2002, with follow-up through January 2015. In addition, deep digital sequencing was performed on paired diagnosis and remission samples from 50 patients (including 32 with intermediate-risk AML), approximately 30 days after successful induction therapy. Twenty-five of the 50 were from the cohort of 71 patients, and 25 were new, additional cases. EXPOSURES: Whole-genome or exome sequencing and targeted deep sequencing. Risk of identification based on genetic data. MAIN OUTCOMES AND MEASURES: Mutation patterns (including clearance of leukemia-associated variants after chemotherapy) and their association with event-free survival and overall survival. RESULTS: Analysis of comprehensive genomic data from the 71 patients did not improve outcome assessment over current standard-of-care metrics. In an analysis of 50 patients with both presentation and documented remission samples, 24 (48%) had persistent leukemia-associated mutations in at least 5% of bone marrow cells at remission. The 24 with persistent mutations had significantly reduced event-free and overall survival vs the 26 who cleared all mutations. Patients with intermediate cytogenetic risk profiles had similar findings. [table: see text]. CONCLUSIONS AND RELEVANCE: The detection of persistent leukemia-associated mutations in at least 5% of bone marrow cells in day 30 remission samples was associated with a significantly increased risk of relapse, and reduced overall survival. These data suggest that this genomic approach may improve risk stratification for patients with AML.

Genome Modeling System: A Knowledge Management Platform for Genomics.

In this work, we present the Genome Modeling System (GMS), an analysis information management system capable of executing automated genome analysis pipelines at a massive scale. The GMS framework provides detailed tracking of samples and data coupled with reliable and repeatable analysis pipelines. The GMS also serves as a platform for bioinformatics development, allowing a large team to collaborate on data analysis, or an individual researcher to leverage the work of others effectively within its data management system. Rather than separating ad-hoc analysis from rigorous, reproducible pipelines, the GMS promotes systematic integration between the two. As a demonstration of the GMS, we performed an integrated analysis of whole genome, exome and transcriptome sequencing data from a breast cancer cell line (HCC1395) and matched lymphoblastoid line (HCC1395BL). These data are available for users to test the software, complete tutorials and develop novel GMS pipeline configurations. The GMS is available at https://github.com/genome/gms.