Zachary Sethna

Abstract Pancreatic ductal adenocarcinoma (PDAC) is lethal in 88% of patients 1 , yet harbours mutation-derived T cell neoantigens that are suitable for vaccines 2,3 . Here in a phase I trial of adjuvant autogene cevumeran, an individualized neoantigen vaccine based on uridine mRNA–lipoplex nanoparticles, we synthesized mRNA neoantigen vaccines in real time from surgically resected PDAC tumours. After surgery, we sequentially administered atezolizumab (an anti-PD-L1 immunotherapy), autogene cevumeran (a maximum of 20 neoantigens per patient) and a modified version of a four-drug chemotherapy regimen (mFOLFIRINOX, comprising folinic acid, fluorouracil, irinotecan and oxaliplatin). The end points included vaccine-induced neoantigen-specific T cells by high-threshold assays, 18-month recurrence-free survival and oncologic feasibility. We treated 16 patients with atezolizumab and autogene cevumeran, then 15 patients with mFOLFIRINOX. Autogene cevumeran was administered within 3 days of benchmarked times, was tolerable and induced de novo high-magnitude neoantigen-specific T cells in 8 out of 16 patients, with half targeting more than one vaccine neoantigen. Using a new mathematical strategy to track T cell clones (CloneTrack) and functional assays, we found that vaccine-expanded T cells comprised up to 10% of all blood T cells, re-expanded with a vaccine booster and included long-lived polyfunctional neoantigen-specific effector CD8 + T cells. At 18-month median follow-up, patients with vaccine-expanded T cells (responders) had a longer median recurrence-free survival (not reached) compared with patients without vaccine-expanded T cells (non-responders; 13.4 months, P = 0.003). Differences in the immune fitness of the patients did not confound this correlation, as responders and non-responders mounted equivalent immunity to a concurrent unrelated mRNA vaccine against SARS-CoV-2. Thus, adjuvant atezolizumab, autogene cevumeran and mFOLFIRINOX induces substantial T cell activity that may correlate with delayed PDAC recurrence.

MOTIVATION: High-throughput sequencing of large immune repertoires has enabled the development of methods to predict the probability of generation by V(D)J recombination of T- and B-cell receptors of any specific nucleotide sequence. These generation probabilities are very non-homogeneous, ranging over 20 orders of magnitude in real repertoires. Since the function of a receptor really depends on its protein sequence, it is important to be able to predict this probability of generation at the amino acid level. However, brute-force summation over all the nucleotide sequences with the correct amino acid translation is computationally intractable. The purpose of this paper is to present a solution to this problem. RESULTS: We use dynamic programming to construct an efficient and flexible algorithm, called OLGA (Optimized Likelihood estimate of immunoGlobulin Amino-acid sequences), for calculating the probability of generating a given CDR3 amino acid sequence or motif, with or without V/J restriction, as a result of V(D)J recombination in B or T cells. We apply it to databases of epitope-specific T-cell receptors to evaluate the probability that a typical human subject will possess T cells responsive to specific disease-associated epitopes. The model prediction shows an excellent agreement with published data. We suggest that OLGA may be a useful tool to guide vaccine design. AVAILABILITY AND IMPLEMENTATION: Source code is available at https://github.com/zsethna/OLGA. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A fundamental challenge for cancer vaccines is to generate long-lived functional T cells that are specific for tumour antigens. Here we find that mRNA–lipoplex vaccines against somatic mutation-derived neoantigens may solve this challenge in pancreatic ductal adenocarcinoma (PDAC), a lethal cancer with few mutations. At an extended 3.2-year median follow-up from a phase 1 trial of surgery, atezolizumab (PD-L1 inhibitory antibody), autogene cevumeran1 (individualized neoantigen vaccine with backbone-optimized uridine mRNA–lipoplex nanoparticles) and modified (m) FOLFIRINOX (chemotherapy) in patients with PDAC, we find that responders with vaccine-induced T cells (n = 8) have prolonged recurrence-free survival (RFS; median not reached) compared with non-responders without vaccine-induced T cells (n = 8; median RFS 13.4 months; P = 0.007). In responders, autogene cevumeran induces CD8+ T cell clones with an average estimated lifespan of 7.7 years (range 1.5 to roughly 100 years), with approximately 20% of clones having latent multi-decade lifespans that may outlive hosts. Eighty-six percent of clones per patient persist at substantial frequencies approximately 3 years post-vaccination, including clones with high avidity to PDAC neoepitopes. Using PhenoTrack, a novel computational strategy to trace single T cell phenotypes, we uncover that vaccine-induced clones are undetectable in pre-vaccination tissues, and assume a cytotoxic, tissue-resident memory-like T cell state up to three years post-vaccination with preserved neoantigen-specific effector function. Two responders recurred and evidenced fewer vaccine-induced T cells. Furthermore, recurrent PDACs were pruned of vaccine-targeted cancer clones. Thus, in PDAC, autogene cevumeran induces de novo CD8+ T cells with multiyear longevity, substantial magnitude and durable effector functions that may delay PDAC recurrence. Adjuvant mRNA–lipoplex neoantigen vaccines may thus solve a pivotal obstacle for cancer vaccination. In a phase 1 trial, patients with pancreatic ductal adenocarcinoma who were treated with surgery and bespoke neoantigen mRNA vaccines combined with anti-PD-L1 and chemotherapy exhibited marked long-lived persistence of neoantigen-specific CD8+ T cell clones, which correlated with prolonged recurrence-free survival at a 3.2-year follow-up.

We quantify the VDJ recombination and somatic hypermutation processes in human B cells using probabilistic inference methods on high-throughput DNA sequence repertoires of human B-cell receptor heavy chains. Our analysis captures the statistical properties of the naive repertoire, first after its initial generation via VDJ recombination and then after selection for functionality. We also infer statistical properties of the somatic hypermutation machinery (exclusive of subsequent effects of selection). Our main results are the following: the B-cell repertoire is substantially more diverse than T-cell repertoires, owing to longer junctional insertions; sequences that pass initial selection are distinguished by having a higher probability of being generated in a VDJ recombination event; somatic hypermutations have a non-uniform distribution along the V gene that is well explained by an independent site model for the sequence context around the hypermutation site.