Efficient identification of neoantigen-specific T-cell responses in advanced human ovarian cancer
Journal for ImmunoTherapy of Cancer volume 7, Article number: 156 (2019)
Efficient identification of neoantigen-specific T-cell responses in epithelial ovarian cancer (EOC) remains a challenge. Existing investigations of spontaneous T-cell response to tumor neoepitope in EOC have taken the approach of comprehensive screening all neoantigen candidates, with a validation rate of 0.5–2%.
Whole-exome and transcriptome sequencing analysis of treatment-naive EOC patients were performed to identify neoantigen candidates, and the immunogenicity of prioritized neoantigens was evaluated by analyzing spontaneous neoantigen-specfic CD4+ and CD8+ T-cell responses in the tumor and/or peripheral blood. The biological relevance of neoantigen-specific T-cell lines and clones were analyzed by evaluating the capacity of autologous ovarian tumor recognition. Genetic transfer of T-cell receptor (TCR) from these neoantigen-specific T-cell clones into peripheral blood T-cells was conducted to generate neoepitope-specific T-cells. The molecular signature associated with positive neoantigen T-cell responses was investigated, and the impacts of expression level and lymphocyte source on neoantigen identification were explored.
Using a small subset of prioritized neoantigen candidates, we were able to detect spontaneous CD4+ and/or CD8+ T-cell responses against neoepitopes from autologous lymphocytes in half of treatment-naïve EOC patients, with a significantly improved validation rate of 19%. Tumors from patients exhibiting neoantigen-specific T-cell responses exhibited a signature of upregulated antigen processing and presentation machinery, which was also associated with favorable patient survival in the TCGA ovarian cohort. T-cells specific against two mutated cancer-associated genes, NUP214 and JAK1, recognized autologous tumors. Gene-engineering with TCR from these neoantigen-specific T-cell clones conferred neoantigen-reactivity to peripheral T-cells.
Our study demonstrated the feasibility of efficiently identifying both CD4+ and CD8+ neoantigen-specific T-cells in EOC. Autologous lymphocytes genetically engineered with tumor antigen-specific TCR can be used to generate cells for use in the personalized adoptive T-cell transfer immunotherapy.
Epithelial ovarian cancer (EOC) is the deadliest gynecological cancer, with approximately 22,240 new cases and 14,070 deaths in the United States for 2018. Approximately 80% of EOCs are diagnosed at an advanced stage, for which the standard treatment is surgery followed by platinum-taxane chemotherapy. Despite the initial efficacy of these standard of care approaches, the overall five-year survival probability is only 28% , and there is an enormous unmet need for the development of alternative therapies. Association between improved clinical outcome and increased intraepithelial CD3+ and/or CD8+ tumor-infiltrating lymphocytes (TILs) suggested antitumor roles of T-cells in EOC [2,3,4,5]. High-affinity neoantigen-reactive T-cells can escape negative selection in the thymus and might have a greater potential to evoke a multi-pronged anti-tumor immune response due to the lack of central tolerance against them . Indeed, recent correlative clinical studies indicate that T-cell reactivity to neoantigens is an important determinant of response to immune checkpoint inhibitors and other immunotherapies , suggesting that efforts to precisely define immunogenic neoantigens either for vaccination [8, 9] or adoptive T-cell therapy (ACT) [10, 11] could potentially provide clinical benefit .
A number of studies have reported T-cells specific to neoantigens in highly mutated tumors such as melanoma and lung cancer [8, 13,14,15,16,17,18,19]. Results from the existing studies of spontaneous T-cell response to tumor neoepitopes in EOC have been mixed [20,21,22,23]. Two earlier studies suggested that EOCs are rarely/unlikely to elicit neoepitope specific spontaneous T-cell response due to the relatively low somatic mutation burden [20, 21]. However, these studies were either limited by the small sample size (1–3 patients) or carried out in a murine ovarian tumor model, which has the concern of neoantigen being silenced and/or lost due to immune editing in the immunocompetent mouse. The studies were also restricted to CD8+ T-cells only, leaving the landscape of spontaneous CD4+ T-cell responses to tumor neoepitopes unexplored. Two more recent studies with larger sample size implied that by comprehensively screening all possible neoantigen candidates, identification of neoepitope specific T-cells is achievable in EOCs [22, 23]. Of the 1714 and 776 putative mutated neoantigens screened, less than 0.5 and 2% were shown to be immunogenic in validation experiments, respectively. Moreover, the capacity of autologous ovarian tumor recognition by neoantigen-specific T-cells has not been addressed.
The aim of this study was to investigate whether the validation rate of neoantigen identification in EOCs can be significantly improved through in silico prioritization. Whole-exome and transcriptome sequencing analysis of treatment-naive EOC patients were performed to identify neoantigen candidates, and the immunogenicity of prioritized neoantigens was evaluated by analyzing spontaneous neoantigen-specfic CD4+ and CD8+ T-cell responses in the tumor and/or peripheral blood. The biological relevance of neoantigen-specific T-cell lines and clones were analyzed by evaluating the capacity of autologous ovarian tumor recognition. Genetic transfer of T-cell receptor (TCR) from these neoantigen-specific T-cell clones into peripheral blood T-cells was conducted to generate neoepitope-specific T-cells. The molecular signature associated with positive neoantigen T-cell responses was investigated, and the impacts of expression level and lymphocyte source on neoantigen identification were discussed .
The characteristics of the study cohort are shown in Additional file 12: Table S1. All 20 treatment-naïve patients underwent maximal debulking surgery (85% were optimally debulked, with 45% complete resection) from which the tumor tissue were procured. Patients in this cohort had typical characteristics of advanced EOC cases: median age at diagnosis of 60 (range 44 to 89), high stage (IIIC, IV; 100%), and the majority with high grade serous histology (75%). The median duration of follow-up was 29.7 months. The median progression-free survival was 18.1 months and the median overall survival was 30.9 months.
Whole-exome sequencing (WES) was performed on 22 pre-therapy biopsies and matched normal samples (Peripheral blood mononuclear cells, PBMC) from the 20 EOC patients, as the first step in our workflow for neoantigen discovery and prioritization (Fig. 1a). The specimens consisted of primary tumor in 9 patients, locally invasive tumor in 9 patients, paired primary and locally invasive tumors in 2 patients. Somatic mutations were identified by comparing the tumor with the matched blood DNA as described [24, 25]. We identified a total of 2096 somatic mutations from the 20 patients, including 1368 non-synonymous somatic mutations (median = 62), and the number of genes with altered amino acid sequence ranged from 9 to 183 per patient (Fig. 1b). TP53 was mutated in 16 patients, including 7 truncating mutations predicted to cause loss-of-function (Additional file 1: Figure S1). Nine genes were mutated in 3 out of 20 patients, including two known Cancer Gene Census (CGC) genes : NF1 and STAG2. Of these nine genes, IL27RA appears to be interesting as two of the three mutations were truncating mutations. The two IL27RA loss-of-function mutations were both identified from locally invasive tumor while the third IL27RA missense mutation was found in primary tumor. In addition, there are 70 genes mutated in two patients, including seven CGC genes (Additional file 13: Table S2). PTEN, BRCA1, and BRCA2 were each mutated in 2 patients, with all of them loss-of-function mutations. There was no gene found to be mutated at a significantly different frequency between primary and locally invasive tumors. In the two patients with both primary and locally invasive tumors, we compared every mutation’s variant allele fraction (VAF) between the two tumors and showed they were overall highly consistent (Additional file 2: Figure S2a).
Identification of neoantigens
Candidate neoantigens were identified using the computational pipeline as outlined in Fig. 1a. We identified neoantigen candidates as mutations harboring mutant peptide whose binding affinity to patients’ MHC was not only strong (< 150 nM), but also specific (i.e., higher affinity than that of matched wild-type peptide). The second condition was included to enrich highly immunogenic neoepitopes because some neoepitope-reactive TCRs were considered to cross-react to the corresponding wild-type epitope if both mutated and wild-type epitopes were similarly presented . In this case, T-cell precursors expressing such TCRs are expected to be eliminated in the thymus, which would result in lower precursor frequencies in the periphery. A total of 449 neoantigen candidates were found to have at least one predicted neopeptide with strong and specific MHC I and/or II-binding affinity. The number of predicted neoantigen in each patient ranged from 4 to 75, with a median of 21 (Additional file 14: Table S3). In the two patients with primary/locally-invasive tumor pair, the majority of neoantigens were shared by both primary and locally invasive tumors (100 and 77.8%, respectively) (Additional file 2: Figure S2b).
These 449 neoantigens included 215 which contained MHC I-binding neopeptides, and 349 which contained class II-binding neopeptides (Fig. 1c). One quarter (115 of 449) of neoantigens were predicted to harbor neopeptides which bind to both class I and II. We classified these 449 neoantigens based on the mutant allele’s expression level in RNAseq data (see Methods). About half (209) neoantigens demonstrated robust expression of the mutant allele while the rest (240) did not.
Prioritization of neoantigens
To prioritize neoantigen candidates for peptide synthesis and T-cell assays, we ranked all predicted neoantigens within each patient based on a pre-defined set of criteria: 1) mutations in CGC genes; 2) MHC binding affinity of the mutant peptide; 3) difference in the binding affinities between mutant and the matched wild type peptides; 4) variant allele fraction (VAF) of the mutation; 5) expression level including both the mutant allele and the overall level of the gene; 6) type of MHC binding (class I only, class II only, or both class I and II). To facilitate this process, we designed a specific type of visualization plot (Christmas Light Plot, or CLP) incorporating all these types of information (Fig. 1d). The final selection of neopeptides involved a target selection board that evaluated the target peptides based on the criteria described above, with additional considerations of biochemical properties related to peptide synthesizability. For the ten patients with autologous PBMC, tumor-derived single-cell suspension and tumor biospecimens available, a total of 75 neopeptides were selected for synthesis, with a median of 7 and a range of 3–12 neopeptides per patient (Additional file 15: Table S4). These include 36, 32, and 7 neopeptides that are predicted to bind to class I only, class II only, and both class I and II, respectively. Twenty five of these 75 neopeptides did not demonstrate robust expression of the mutant allele in RNAseq (Methods), and they were included to investigate the relationship between expression level and induction of T-cell response.
Evaluation of neoepitope immunogenicity
The immunogenicity of the neoepitopes was evaluated in the ten patients from whom live T-cells from both PBMC and tumors were available. CD4+ and CD8+ T-cells were isolated from each specimen and stimulated with T-cell-depleted PBMCs as antigen-presenting cells (APCs) that had been pulsed with pooled neoepitopes. T-cell-depleted PBMCs were used to enrich APCs such as dentritic cells, monocytes/macrophages and B cells. A total of 27 IFN-γ-producing T-cell responses were detected in samples from 5 of 10 patients, including 20 responses against 10 individual neopeptide and 7 responses against 4 pooled neopeptides (Fig. 2a). These positive T-cell responses were highly mutant-specific, with the reactivity against mutated peptide at least two-fold greater than the corresponding wild-type peptide (Additional file 3: Figure S3 and Additional file 4: Figure S4). Both CD8+ and CD4+ T-cells showed neoepitope-specific responses, with 13 responses mediated by CD8+ T-cells and 14 by CD4+ T-cells. The median number of positive T-cell responses against individual neopeptides was 4 in the 5 responders, with a median of 2 reactive neoepitopes per patient.
For the 9 neopeptides (or neopeptide pools) that were recognized by CD8+ T-cells, 5 exclusively elicited mutant-specific T-cell response in either TILs or PBMCs. Likewise, 4 of the 9 neopeptides (or neopeptide pools) that were recognized by CD4+ T-cells exclusively elicited mutant-specific T-cell response in either TILs or PBMCs. This is consistent with prior evidence that there exists a discordance of neoepitope recognition between TILs and PBMCs . Examples included the CD4+ T-cell response against mutant epitopes of JAK1 which were detected only in TILs, and the CD8+ T-cell response against mutant epitopes of TRPC4 which was present only in PBMCs (Fig. 2a). On the other hand, 18 of the 27 responses were detected in both TILs and PBMCs, indicating that both types of patient’s specimen are useful for identifying neoantigen-reactive T-cells (Fig. 2b). Examples include CD8+ and CD4+ T-cell responses against IL27RA2 which were detetcted in both PBMC and TILs. Among the four patients with detected T-cell responses in both TILs and PBMCs, the responses were detected higher in TILs than PBMCs in three patients; while in Pt #5, there were more response in PBMCs than TILs (Fig. 2b). In order to explore immunosuppression status of TILs, we analyzed the expression level of a panel of 10 immune inhibitory molecules from the tumor RNA-Seq data of these four patients. These immune inhibitory molecules include PD1, PDL1, CTLA4, CD80, CD86, LAG3, TIM3, LAGLS9, MYC and FOXP3. Remarkably, Pt #5 showed higher expression of all these immune inhibitory genes than those of the other three patients (Additional file 5: Figure S5).
In total, 7/50 (14%) predicted neoepitopes that showed robust expression in RNAseq data induced T-cell responses. Interestingly, T-cell responses were also detected against 3/25 (12%) predicted neopeptides, KCNH1 in Pt #3, TRPC4 in Pt #15 and DNAH8 in Pt #19, that were not robustly expressed in RNAseq (Additional file 15: Table S4). These 3 genes were all weakly expressed (with RPKM less than 1), although it is possible that the sequence depth of our RNAseq might not be sufficient to detect the mutant allele. It has been suggested that high level of expression may not be required for inclusion of a neoantigen candidate , based on the observations that very low levels (e.g., even a single peptide-HLA complex) may be sufficient for a cell to elicit a cytolytic T-cell response .
Signatures of neoepitope-specific T-cell response
As expected, we found that patients with neoantigen-specific T-cell responses have significantly higher mutation burden and neoantigen load than those without (Fig. 2c). The median number of nonsynomous mutation is 84 in patients with responses, compared with 49 in patients without responses (p = 0.026, one-tailed t-test). The corresponding number of predicted neoantigen is 27 and 8, respectively (p = 0.043, one-tailed t-test). Interestingly, patients with positive T-cell responses are significantly enriched for BRCA1/2 somatic mutations (3/5 vs 0/5, p = 0.038, chi-square test). When gene expression profiles for responders and non-responders were compared, the most significantly and differentially enriched pathways in responders were related to antigen processing and presentation machinery (APPM), suggesting that not only the number of neoantigens, but also antigen processing and presentation in tumor regulate generation of T-cell responses against neoantigens (Fig. 2d). On the other hand, the non-responders were characterized with a unique signature of MYC amplification (Fig. 2e, Additional file 6: Figure S6), which was recently shown to promote immune evasion through the modulation of immune regulatory molecules .
By intersecting the list of differentially expressed genes (responders versus non-responders) with the genes of the APPM pathways, an 31-gene signature of APPM was derived (Fig. 3a, Additional file 16: Table S5). Based on median expression value of the APPM signature, we stratified TCGA ovarian cancer patients into groups of high vs low expression of this signature (Fig. 3b). In consistent with prior evidence that antigen presentation pathway is reduced in high-risk ovarian cancer , patients with higher expression of APPM signature have longer overall survival than patients with lower expression of the signature (p = 5 × 10− 4, Cox Proportional-Hazards Model) (Fig. 3c). Furthermore, patients with higher expression of APPM signature exhibited higher levels of CD4+ memory T-cells, Th1 and Th2 (p = 1.2 × 10− 3, 1 × 10− 13, and 3.1 × 10− 3, respectively, t-test), and lower level of Tregs (p = 4 × 10− 6, t-test) in tumors (Fig. 3d). On the other hand, patients with lower expression of APPM signature have a modest but significant increase in MYC expression (p = 1.7 × 10− 4, t-test) (Fig. 3e). Consistent with these findings, the TCGA ovarian cancer patients displayed a negative association between the expression levels of APPM signature and MYC (p = 1.37 × 10− 6, linear regression) (Fig. 3f). Similar trends were observed in the EOC patients from our study cohort, where the small sample size limited their statistical significance.
Characterization of neoepitope-specific T-cells
In order to further characterize neoepitope-specific T-cells in EOC, we established neoantigen-specific T-cell lines by isolating and expanding peptide-reactive T-cells. T-cell lines that specifically recognized mutated peptides were established in 3 out of 9 cases we attempted. Based on the availability of autologous tumor specimens, we focused on NUP214 neoepitope-specific CD4+ T-cells obtained from TILs of Pt #19 and JAK1 neoepitope-specific CD4+ T-cells obtained from TILs of Pt #11.
About 80% of neoepitope-specific CD4+ T-cell line produced IFN-γ against the mutated NUP214 peptide but not the corresponding wild-type peptide (Fig. 4a). Low-resolution TCR Vβ spectratyping identified the CD4+ T-cell response as oligoclonal (Additional file 7: Fig. S7a), composed of 20% Vβ2+ and 45% Vβ13.1+ T-cells, respectively. Approximately 30% of cells were of other Vβ subtype not identified by this antibody panel. Combination of Vβ-staining and intracellular IFN-γ staining demonstrated that all 3 distinct major subsets recognized the same neopeptide (Additional file 7: Fig. S7b). Therefore, TCR Vβ2+, Vβ13.1+, or Vβ2−Vβ13.1− cells were further isolated by flow cytometric cell sorting to obtain clonal populations (Fig. 4b). Avidity for recognition of mutated peptide by Vβ2+, Vβ13.1+, and Vβ2−Vβ13.1− NUP214-specific CD4+ T-cell clones was similar (Fig. 4c). Responses were strictly specific for mutation as there was no recognition of wild-type peptide even at higher concentrations. As NUP214-specific CD4+ T-cell responses was detected in the tumor from Pt #19 (Fig. 2 a), we reasoned that the mutated NUP214 epitope was naturally presented in the tumor microenvironment. Therefore, we tested whether NUP214-specific CD4+ T-cells are activated by autologous tumor cells. Indeed, NUP214 neoepitope-specific CD4+ T-cells produced IFN-γ specifically against tumor cells, but not against autologous PBMCs (Fig. 4d). These results strongly support that NUP214-specific CD4+ T-cells were activated in the tumor microenvironment. In Pt #19 tumor, both CD45+ hematopoietic cells and EpCAM+ EOC expressed MHC class II (HLA-DR) (Fig. 4e). Therefore, both direct presentation by cancer cells and indirect cross-presentation of tumor-derived NUP214 by hematopoietic antigen-presenting cells are possible as a mechanism for activation of mutated NUP214-specific CD4+ T-cells in the tumor microenvironment.
The JAK1 neoepitope-specific CD4+ T-cells were isolated and expanded (Additional file 8: Figure S8a) and were analyzed for TCR Vβ usage (Additional file 8: Figure S8b). The majority (75%) of JAK1-specific CD4+ T-cell line was Vβ13.6+, indicating a monoclonal response. TCR Vβ13.6+ CD4+ T-cells were further isolated and expanded for functional analyses (Additional file 8: Fig. S8c). Vβ13.6+ CD4+ T-cells specifically recognized JAK1 mutant peptide over the corresponding wild-type peptide (Additional file 8: Figure S8d). Because autologous tumor mononuclear cells (TMCs) were not available for this patient, we tested the reactivity against tumor ascites mononuclear cells (AMCs) and found that CD4+ T-cells produced IFN-γ when co-cultured with the AMCs but not with the autologous PBMCs (Additional file 8: Figure S8e).
Generation of neoepitope-specific T-cells by TCR gene-engineering
To test whether neoantigen-reactivity is solely mediated by TCR and whether neoantigen-specificity can be transferred to other T-cells by TCR gene-engineering, we first cloned TCR gene from 5 neoepitope-specific T-cell clones (3 mutated NUP214-specific CD4-TCR from Pt #19, 1 JAK1-specific CD4-TCR from Pt #11, and 1 TRPC4-specific CD8-TCR from Pt #15) into a retroviral plasmid vector (Fig. 5a) . To test the functionality of the cloned TCR, peripheral T-cells from a healthy donor were polyclonally activated and transduced by the TCR-expressing retroviral vectors. In 4/5 cases, TCR gene-engineered T-cells were demonstrated for mutated peptide-specific reactivity. In the case of 3 NUP214-specific TCRs, two TCRs from Vβ13.1+ or Vβ2−Vβ13.1− T-cell clones provided neoepitope-specific reactivity (Fig. 5b-c), while that from Vβ2+ T-cell clone did not despite of similar reactivity by the parental T-cell clones (Fig. 4c). Both CD4+ and CD8+ T-cells transduced with NUP214-specific TCR showed reactivity against the neoepitope (Additional file 9: Figure S9). Similar observations were made for mutated JAK1-specific TCR, where we have established high-titer retrovirus-packaging PG13 cell clone. After 2 transductions, nearly 60% T-cells expressed transduced TCR, as determined by the increase in TCR Vβ13.6+ expression (Fig. 5d). Mutated JAK1-specific TCR-transduced T-cells also displayed strong and specific reactivity against the mutated JAK1 peptide (Fig. 5e-f). In addition to CD4+ T-cells, we also cloned TCR gene from TRPC4 neoepitope-specific CD8+ T-cells, and confirmed the neoantigen reactivity by TCR gene-engineered T-cells (Additional file 10: Figure S10).
Advanced EOC is a highly lethal disease, with a dismal five-year survival rate below 30% . The strong association of TILs with better survival , which is counter-regulated by innate and adaptive immune escape mechanisms such as expression of immunosuppressive molecules and recruitment of immunosuppressive cells, indicates that patients with advanced EOC might benefit from immunotherapy. Emerging immunotherapies targeting patient-specific neoantigens have attracted considerable interest because immune responses to tumor-specific neoantigens are less constrained by central and peripheral immune tolerance mechanisms .
Given the relatively lower mutation burden of EOC than melanoma and lung cancer, previous investigations of spontaneous T-cell response to tumor neoepitope in EOC have taken the approach of comprehensively screening all identified neoantigen candidates [20,21,22,23]. For example, Bobisse et al screened all 776 predicted HLA class I-binding neopeptide individuals as well as 19 neopeptide pools from 19 patients (mean 46 per patient) and validated 15 of them, or 2% . Dinger et al screened all 1714 putative mutated neoantigens from 7 patients (mean 214 per patient) and validated 8 of them, or 0.5% . In contrast, we only screened a subset (mean 9 per patient) of neopeptide candidates prioritized by our in silico approach, and achieved a substantially higher validation rate (14 of the 75 neopeptide individuals and 10 neopeptide pools, or 19%) for practical use, especially when limited patient-derived specimens are available. On the other hand, our approach might have missed some bona fide neoantigen targets, and therefore underestimated the prevalence of neoepitope specific T-cells in patients with EOC.
As expected, patients exhibiting T-cell responses against neoepitopes were found to have significantly higher mutation burden and neoantigen load. In addition, these responders tended to have higher BRCA1/2 mutation rate than those without. Due to the relatively low mutation burden of EOC and widespread interpatient tumour heterogeneity, our prioritization approach will more likely be applicable to a subset of patients of EOC. In addition, since somatic mutations can vary at different times in tumor evolution (e.g., primary tumor versus tumor relapse), it will be desirable to repeat this procudure sequentially to identify new neoepitope-specific T cells.
We found elevated level of APPM signature in patients who spontaneously developed neoepitope-specific T-cell response, while MYC amplification was observed in patients who did not. Among the five responders in our study, Pt #5 has lower mutation burden than three of the five non-responders (45 vs 49–52, Fig. 2c). Interestingly, the expression level of APPM signature is almost two-fold higher in Pt #5 than the three non-responders (Additional file 11: Figure S11). Our results provided independent edvidence that antigen presentation pathway is associated with ovarian cancer prognosis , and lends additional support for the crucial role of this machinery in TIL infiltration and defining CD8+/CD4+ T-cell responses against neoantigens.
The biological relevance of patient-derived T-cells against neoantigens from known cancer-associated genes, NUP214 and JAK1, were confirmed by the demonstration of autologous tumor recognitions. NUP214 encodes a protein localized to the cytoplasmic face of the nuclear pore complex where it is required for proper cell cycle progression and nucleocytoplasmic transport. NUP214 forms a fusion gene with the DEK gene on chromosome 6 in a t [6, 9] translocation associated with acute myeloid leukemia and myelodysplastic syndrome [35, 36]. Somatic mutations in JAK1/2 have been proposed as a mechanism of evasion to immune recognition, leading to both primary and acquired resistance to PD-1 blockade in some cancer patients [37, 38]. Pre-clinical studies have suggested that targeting JAK/STAT3 could effectively suppress ovarian tumor progression and have therapeutic potential for treating advanced EOC [39,40,41]. Our observations suggest that although cancer cells may evade immunological attack by mutating genes such as JAK1, de novo T-cell responses against such mutations might be exploited to inhibit tumor growth in cancer patients. Interestingly, we discovered IL27RA mutations in 3 patients, including two loss-of-function truncating mutations and one missense mutation. The IL27RA missense mutation was shown to elicit concomitant CD4+ and CD8+ T-cell responses. IL27RA is the alpha subunit of Interleukin 27 (IL27) receptor, which has been reported to demonstrate a dual role of anti-tumor activity and immune regulatory function [42, 43]. A recent proteomic analysis on human ovarian cancer cell lines revealed that IL-27 and IFN-γ shared a broad set of activities, such as the HLA class I antigen presentation . Loss of the IFN-γ pathway gene has been shown as a mechanism of resistance to anti-CTLA-4 therapy . It is possible that mutations in IL27RA provide a novel immune-evasion mechanisms in EOC. Further validations and functional studies will be needed to determine the prognostic and therapeutic potential of IL27-IL27R pathway genes in EOC.
Our study reveals that a significant portion of neoantigen candidates can fail to pass the predefined criteria for mutant allele expression in RNAseq, and we showed some of them could elicit strong mutant-specific CD4+ and/or CD8+ T-cell response in autologous PBMCs and/or TILs. Conceptually, the level of expression of the source protein is a surrogate for determining whether it reaches the threshold required for its efficient processing and presentation in HLA, and whether neoantigen-specific T-cells would efficiently reciognize the tumor target. The latter is mainly determined by affinity of TCR to interact with MHC/neoepitope complex. For virus antigen-specific T-cells, only a single MHC/peptide is sufficient to activate T-cells through high-affinity TCRs . Because affinity for neoantigen-specific TCRs is considered to be high due to absence of central tolerance, it is possible that low level of expression is enough to induce neoantigen-specific T-cell responses. There have been multiple strategies to use expression information for neoantigen prediction, including studies using RNAseq from autologous tumor samples , studies using RNAseq data from different patients of the same cancer type (i.e., TCGA) , and studies not using expression information [22, 28, 48]. Even for the strategy of using RNAseq from autologous tumor samples, the criteria used by different studies to determine neoantigen expression can be vastly different. Therefore, it will be necessary in the future to conduct comprehensive assessments of the role of expression information and the approaches of leveraging it for neoantigen prediction and/or prioritization.
The nature of relatively low mutation burden of EOC also calls for culture enrichment approaches to expand neoepitope specific clones. In the study by Bobisse et al , whole tumors were cultured with high concentration of IL-2 in the presence MHC class I neoepitope pools . In our study, we stimulated isolated CD4+ and CD8+ T-cells from tumors with peptide-pulsed CD4−CD8− cells derived from PBMCs, because APCs in ovarian tumor microenvironment have been shown to be dysfunctional or immunosuppressive . A recent study by Yossef et al showed that the detection of neoantigen-reactive TILs could be enhanced by enriching T-cells that express PD-1 and/or T-cell activation markers followed by microwell culturing to avoid overgrowth of nonreactive T-cells . As the availability of patients’ specimens such as TILs is often limited, future studies are warranted to benchmark the prediction accuracy and effectiveness of different culture enrichment approaches for screening of immunogenic neoantigens in EOC.
We were able to isolate and clone TCRs from in vitro expanded CD4+ and CD8+ T-cell clones reactive against neoepitopes, and demonstrated the feasibility of specifically targeting the neoantigens by TCR gene-engineered T-cells. Our current protocoal takes about 2 weeks to prepare TCR-transduced T-cells from the established neoantigen-specific T-cell lines. Adoptive transfer of autologous tumor antigen-specific T-cells is an effective therapeutic treatment for cancer patients [51, 52]. Using tumor antigen-specific TCR genes, tumor antigen-specific T-cells can be rapidly expanded for infusion into patients in order to mediate immediate elimination of tumors and long-term tumor immune-surveillance. The therapeutic effects of ACT using TCR gene-engineered T-cells have been demonstrated in many clinical trials. High-affinity TCR gene is required to engineer effective T-cells for strong recognition and destruction of cancer cells. It has been recently reported that function of neoepitope-specific CD8+ T-cell from ovarian tumors is significantly higher than PBMCs , suggesting neoepitope-specific CD4+ and/or CD8+ T-cells in TILs could be a useful source of TCR for ACT.
The neoepitope-specific TCRs could be introduced into naïve or central memory T-cells to target driver mutations , simultaneously target multiple mutations, and combine ACT with other immunomodulators such as checkpoint inhibitors. Conceptually multiple neoantigens can be targeted but it is difficult to define an optimal number of neoepitopes. The range of neoepitopes utilized in published clinical trials of neoantigen vaccines varies between 7 and 20 [8, 9, 53]. Due to the intratumoral heterogeneity, it would be desirable to target multiple neoantigens for each patient. Although most previous ACT trials using TCR gene-engineered T cells targeted a single tumor antigen epitope, it is feasible to target multiple tumor antigens using a mixture of engineered T-cell products. As assessment of the immunosuppressive microenvironment was not the main focus of our current study, we did not include exogenous checkpoint blockade (e.g. anti-CTLA4 or anti-PD1) into our neoantigen recognition assays. It is well accepted that immunosuppression within the tumor microenvironment represents a critical barrier for anti-tumor T cell function, both neoantigen-reactive and adoptively transferred engineered T cells. A published report from a recent clinical trial evaluating adoptive transfer of TCR-engineered T cells and ipilimumab (anti-CTLA4) demonstrated that while addition of ipilimumab was safe and tolerable, there was no apparent clinical benefit from the combination . Although the small sample size (n = 4) was insufficiently powered, this study and other preclinical/clinical studies suggest immunomodulatory combinations (e.g. anti-PD1, anti-TGF-β, IDO inhibitors) could still potentiate therapeutic effect of ACT. Additional studies are required to empirically determine the optimal combination for eliciting clinical benefit. Our current study focused on analyzing treatment-naïve ovarian cancer patients at the time of primary debulking surgery. Future studies will be warranted to determine intrinsic factors of tumor and microenvironment as well as the nature of TCR affecting the infiltration of adoptively transferred T cells. In future ACT clinical trials using TCR gene-engineered T cells, it is also feasible to quantify engineered T cells in the tumors using quantitative TCR or digital PCR. Taken together, these results warranted future studies of patient-specific neoantigens as potential targets for downstream translational validation towards adoptive T-cell immunotherapy of ovarian cancer.
In summary, our study demonstrated the feasibility of efficiently identifying both CD4+ and CD8+ neoantigen-specific T-cells in ovarian cancer. Further methodology development , including the improvement of patient response prediction, neoantigen prediction and prioritization, T-cell enrichment and TCR engineering, will be warranted to exploit the therapeutic potential of neoantigen-targeting for personalized immunotherapy in treating advanced EOC.
Patient and sample characteristics
Tumor specimens were collected at the time of primary debulking surgery, flash frozen in liquid nitrogen and stored at − 80 °C. Portions of tumor specimens were dissociated by the Gentle MACS system (Miltenyi Biotec) to obtain single-cell suspensions. After a density gradient isolation, live mononuclear cells were stored in liquid nitrogen until use. PBMC were obtained using a density gradient method and stored in liquid nitrogen. DNA was extracted from the frozen tissues and PBMCs using the GenFIND DNA extraction kit (Agencourt, Pasadena, CA) per the manufacturer’s instructions. All pathology specimens were reviewed by experienced gynecologic pathologists and tumors were classified according to the WHO criteria . Prior to surgery, no patients received neo-adjuvant chemotherapy, and subsequent to surgery all patients received adjuvant platinum/taxane-based chemotherapy. No patient had received immunotherapy treatment prior to surgery.
Whole-exome sequencing and somatic mutation calling
Individual exome capture of each DNA sample followed by single-indexed library generation was carried out using the SureSelect XT Target Enrichment System (Human All Exon V5 kit, Agilent Inc). Cluster generation on cBot was followed by 2 × 100 bp paired-end (PE) sequencing on a HiSeq2500 sequencer according to the manufacturer’s recommended protocol (Illumina Inc.). High quality WES paired-end reads passing Illumina RTA filter were aligned to the NCBI human reference genome (GRCh37) using Burrows-Wheeler Alignment (BWA) tool . PCR duplicated reads were marked and removed using Picard tool. All samples had more than 80% of the targeted regions covered by at least 30-fold. Putative mutations were identified by running variation detection module of Bambino  and Strelka , and then further filtered as previously described [25, 60]. All putative single-nucleotide variants (SNVs) were combined and further filtered based on a standard set of criteria to remove the following common types of false calls:  the alternative allele is present in the matching normal sample and the contingency between the tumor and normal samples is not statistically significant;  the mutant alleles are only present in one strand and the strand bias is statistically significant;  the putative mutation occurs at a site with systematically dropped base quality scores;  the reads harboring the mutant allele are associated with poor mapping quality. The identified somatic mutations were compared to the public human germline databases including dbSNP , 1000 Genomes Project , NHLBI’s Exome Sequencing Project to further exclude remaining germline polymorphisms. All mutations were manually reviewed to ensure accuracy and annotated using ANNOVAR  with the latest NCBI RefSeq database. To extract somatic copy number information based on the sequenced exomes of the samples, Varscan2  was employed, and the results were smoothed and segmented with the Bioconductor DNACopy package . GISTIC2.0  was used to identify recurrent somatic chromosomal alterations in the five samples with mutant-specific T-cell response and the five samples without, respectively.
RNAseq and gene expression analysis
RNAs were purified from fresh frozen tumors using the miRNeasy mini kit (Qiagen). The sequencing libraries were prepared with the TruSeq Stranded Total RNA kit (Illumina Inc) and sequenced for 100 cycle paired-end on a HiSeq2500 sequencer, following the manufacturer’s recommended protocol (Illumina Inc.). Raw reads from RNAseq that passed the Illumina RTA quality filter were first pre-processed using Cutadapt to remove adapter sequences, followed by FASTQC for sequencing base quality control. The remaining reads were mapped to the NCBI human reference genome (GRCh37) and RefSeq annotation database using Tophat . A second round of quality control was performed to identify potential RNAseq library preparation problems by examining mapped BAM files using RSeQC . From the mapping results, the number of reads aligning to each gene was calculated using HTSeq  and RPKM (Reads Per Kilobase of transcript per Million mapped reads) values were obtained using RSEM . ssGSEA , an extension of Gene Set Enrichment Analysis , was performed on gene expression data of each sample using the Reactome gene sets obtained from the MSigDB database . The normalized enrichment score of each gene set in each sample was retrieved, and compared between the five samples with mutant-specific T-cell response and the five samples without. Differentially expressed genes (DEGs) between those patients who responded and those who didn’t were determined using LIMMA . By intersecting the list of DEGs (p < 0.05) with the genes involved in the antigen processing and presentation machinery (APPM) pathways, a signature of APPM constsiting of 31 DEGs were obtained. TCGA gene expression dataset were downloaded from cBioportlal  using Ovarian Serous Cystadenocarcinoma (TCGA Provisional, 307 samples with RNAseq). TCGA ovarian cancer patients were ranked based on median expression value of the 31-gene signature and grouped in top vs bottom tertile (groups of high versus low expression level of signature). Survival analysis was done using the survival package in R. Precalculated dataset of tumor-infiltrating subtypes for the TCGA Ovarian cancer patients were downloaded from xCell  and the comparsion of tumor-infiltrating subtypes between groups of high versus low expression level of signature was performed using Student’s t-test. Correlation between the expression of APPM signature and MYC was calculated using Pearson correlation, and a linear model was built to fit the data and test significance and was plotted as trendline with the confidence intervals.
For each missense SNVs, we obtained 8 to 15-mer peptide sequences containing the mutated amino acid as well as corresponding wild-type ones from RefSeq . Genotypes for patients’ class I and II HLA alleles were determined from next-generation sequencing data using Polysolver and HLAminer, respectively, with default parameters [78, 79]. Class I and II-binding affinities for each combination of peptide/HLA type were predicted using NetMHCpan v3.0 [80,81,82] and NetMHCIIpan v3.1 , respectively, with default parameters. For a mutant peptide to be considered as a neoepitope, we required: 1) mutant binding affinity score is less than 150 nM; 2) the ratio of binding affinity between mutant and the matched wild-type peptides is less than 0.9; and 3) the difference in binding affinity between mutant and wild-type epitopes is at least 100 nM, except for peptides from any Cancer Census genes . If one mutation was predicted to generate multiple neoepitopes, it was counted as one neoantigen . For the predicted neoantigen, we assessed the expression of mutant allele in RNAseq using a set of previously published criteria for neoantigen expression : 1) at least two supporting reads in RNAseq; 2) minimum variant allele fraction (VAF) of 4% for mutations with at least three reads or 20% for mutations with exact two supporting reads; 3) no significant strand bias (p < 0.05).
Synthetic peptides for neoepitopes and the corresponding wild-type epitopes were manufactured at > 90% purity (Genscript, Piscataway, NJ). Reverse phase HPLC produced lyophilized peptides were reconstituted in DMSO (Sigma) and used to test T-cell reactivity.
Monitoring of neoepitope-specific T-cell response
Mutated peptide specific T-cell response was investigated using T-cell presensitization method as previously described . Briefly, CD8+ and CD4+ T-cells were sequentially isolated from PBMCs or tumors of EOC patients using Dynabeads CD8 and CD4 positive isolation kits (ThermoFisher Scientific) and incubated overnight in a 37 °C 5% CO2 incubator. Remaining CD4−CD8− cells from PBMCs were used as APCs for pulsing a pool of 2 μM patient-specific mutated peptides overnight, irradiated at 3000 rad, washed and mixed with CD4+ or CD8+ T-cells. The CD4+ and CD8+ T-cells were cultured in the presence of 10 U/ml IL-2 (Roche) and 10 ng/ml IL-7 (R&D systems). A part of CD4+ T-cells were polyclonally activated by phytohemagglutinin (PHA, Remel-ThermoFisher Scientific) and cultured in the presence of IL-2 and IL-7 to prepare T-APC (antigen-presenting T-cells) . At 13–15 days after the culture, these CD4+ and CD8+ T-cells were harvested and tested reactivity against individual mutated or wild-type peptides, or pooled peptides-loaded on autologous T-APC by IFN-γ ELISPOT assay as previously described . Neoepitope-specific response was considered as positive with a minimum of 25 IFN-γ-spot-forming cells against mutated peptide per 5 × 104 cells as well as the number of spots was 2-fold higher than the corresponding wild-type peptide .
Establishment of neoepitope-specific T-cell clones
In order to characterize in details neoantigen recognition by neoepitope-specific T-cells, we established neoepitope-specific T-cell clones. For CD8+ T-cells, presensitized T-cells that showed neoantigen-specific reactivity in ELISPOT assays were restimulated by mutated peptide-pulsed T-APC and IFN-γ-producing T-cells were labeled using IFN-γ-capture reagent (Miltenyi Biotec) and sorted by flow-cytometry as described . For CD4+ T-cells, presensitized T-cells were similarly restimulated in the presence of monensin (Sigma) and phycoerythrin (PE)-labeled anti-CD154 monoclonal antibody (mAb) as described  and CD154-expressing cells were sorted. Isolated T-cells were polyclonally expanded by PHA stimulation in the presence of allogeneic irradiated PBMC, IL-2 and IL-7. Purity and clonality of T-cells were tested by low-resolution TCR spectratyping using Vβ subtype-specific antibodies (Beckman Coulter). For some oligoclonal T-cell cultures containing different Vβ-expressing T-cells, cells were sorted again based on Vβ expression. Reactivity of neoepitope-specific T-cell clones were tested by ELISA and/or intracellular cytokine staining .
Engineering neoepitope-specific T-cells by TCR gene transduction
Retroviral vectors expressing neoepitope-specific TCR genes were constructed as described previously with modifications . Briefly, a part of sorted neoepitope-specific T-cells (2,000–5,000 T-cells) were lysed in TRIReagent. Total RNA was extracted by columns (Zymo Research), and reverse transcribed using SuperScript IV First-Strand Synthesis System (Thermo Fisher Scientific) using oligo dT primers according to the manufacturer’s instruction. Variable regions for TCR α and β chains were independently PCR-amplified using multiplexed primers, and assembled into a retroviral plasmid vector together with constant regions. Plasmids were amplified in NEBStable competent E. coli (New England Biolabs) and extracted using columns (ZymoResearch). To generate retroviral particles, GP2–293 packaging cell line (Clontech) was co-transfected with TCR-expressing transfer plasmid and envelope (pVSV-G: Clontech) plasmid using Lipofectamine 3000 reagent (Thermo Fisher Scientific). Medium was exchanged 7 h after transfection. Retroviral vectors were harvested at around 36 and 60 h after changing medium. PBMCs from healthy individuals were activated by PHA in the presence of IL-2, IL-7 and 10 ng/ml IL-12 (Peprotech). Activated T-cells were harvested at 36–48 h, and transduced by 125 μl freshly harvested retroviral vectors in a 96-well flat-bottom plate which was coated with Retronectin (TaKaRa Bio) and anti-CD3 mAb (OKT3; eBioscience). Transduction was repeated at 24 h after the first transduction on Retronectin without anti-CD3 mAb. In some experiments, we established PG13 clones stably producing high-titer retroviral vectors as described . Briefly, PG13 (ATCC) was transduced with retroviral supernatant from GP2–293 and cloned by limiting dilution. To infect T-cells with PG13-derived retroviral particles, retroviral supernatant was added in a 96-well plate which was coated with Retronectin with or without anti-CD3 mAb and incubated for 6 h. After removing retroviral supernatant, the plate was rinsed by PBS containing 1% bovine serum albumin (BSA). Activated T-cells were transduced as described above. Function of transduced T-cells was investigated 3–7 days after the second transduction.
Statistical analysis was performed using R. Primary statistical tests include Student’s t test for single comparisons of normally distributed data, Chi-square test for comparsion of BRCA1/2 mutation rate, Pearson’s correlation test for comparisons of APPM signature expression versus MYC expression, and Cox Proportional-Hazards Model to compare survival of TCGA patients with high versus low expression of APPM signature. All statistical tests are 2-tailed unless otherwise specified. A p-value less than 0.05 was considered significant.
Availability of data and materials
The datasets generated and/or analyzed during the current study are available from the corresponding author on reasonable request.
Adoptive cell therapy
Ascites mononuclear cells
Antigen presenting cells
Antigen processing and presentation machinery
Bovine serum albumin
Cancer gene census
Christmas light plot
Differentially expressed genes
Enzyme-linked immune absorbent spot
Epithelial ovarian cancer
International Federation of Gynaecology and Obstetrics
Granulocyte macrophage colony stimulating factor
Human leukocyte antigen
High-performance liquid chromatography
Major histocompatibility complex
Peripheral blood lymphocytes
Peripheral blood mononuclear cells
Phosphate buffered saline
Response Evaluation Criteria In Solid Tumors
Single nucleotide variants
The Cancer Genome Atlas
Tumor infiltrating lymphocytes
Tumor mononuclear cells
Variant allele frequency
Whole exome sequencing
Siegel Rebecca L, Miller Kimberly D, Jemal A. Cancer statistics, 2016. CA Cancer J Clin. 2016;66:7–30.
Curiel TJ, et al. Specific recruitment of regulatory T cells in ovarian carcinoma fosters immune privilege and predicts reduced survival. Nat Med. 2004;10:942–9.
Matsuzaki J, et al. Tumor-infiltrating NY-ESO-1-specific CD8+ T cells are negatively regulated by LAG-3 and PD-1 in human ovarian cancer. Proc Natl Acad Sci U S A. 2010;107:7875–80.
Sato E, et al. Intraepithelial CD8+ tumor-infiltrating lymphocytes and a high CD8+/regulatory T cell ratio are associated with favorable prognosis in ovarian cancer. Proc Natl Acad Sci U S A. 2005;102:18538–43.
Zhang L, et al. Intratumoral T cells, recurrence, and survival in epithelial ovarian cancer. N Engl J Med. 2003;348:203–13.
Schumacher TN, Schreiber RD. Neoantigens in cancer immunotherapy. Science. 2015;348:69.
Snyder A, et al. Genetic basis for clinical response to CTLA-4 blockade in melanoma. N Engl J Med. 2014;371:2189–99.
Ott PA, et al. An immunogenic personal neoantigen vaccine for patients with melanoma. Nature. 2017;547:217–21.
Sahin U, et al. Personalized RNA mutanome vaccines mobilize poly-specific therapeutic immunity against cancer. Nature. 2017;547:222–6.
Tran E, et al. Cancer immunotherapy based on mutation-specific CD4+ T cells in a patient with epithelial cancer. Science. 2014;344:641–5.
Tran E, et al. T-cell transfer therapy targeting mutant KRAS in Cancer. N Engl J Med. 2016;375:2255–62.
Zsiros E, Tanyi J, Balint K, Kandalaft LE. Immunotherapy for ovarian cancer: recent advances and perspectives. Curr Opin Oncol. 2014;26:492–500.
van Rooij N, et al. Tumor exome analysis reveals neoantigen-specific T-cell reactivity in an ipilimumab-responsive melanoma. J Clin Oncol. 2013;31:e439–42.
Lu YC, et al. Efficient identification of mutated cancer antigens recognized by T cells associated with durable tumor regressions. Clin Cancer Res. 2014;20:3401–10.
Carreno BM, et al. Cancer immunotherapy A dendritic cell vaccine increases the breadth and diversity of melanoma neoantigen-specific T cells. Science. 2015;348:803–8.
Robbins PF, et al. Mining exomic sequencing data to identify mutated antigens recognized by adoptively transferred tumor-reactive T cells. Nat Med. 2013;19:747–52.
Linnemann C, et al. High-throughput epitope discovery reveals frequent recognition of neo-antigens by CD4+ T cells in human melanoma. Nat Med. 2015;21:81–5.
Kreiter S, et al. Mutant MHC class II epitopes drive therapeutic immune responses to cancer. Nature. 2015;520:692–6.
McGranahan N, et al. Clonal neoantigens elicit T cell immunoreactivity and sensitivity to immune checkpoint blockade. Science. 2016;351:1463–9.
Wick DA, et al. Surveillance of the tumor mutanome by T cells during progression from primary to recurrent ovarian cancer. Clin Cancer Res. 2014;20:1125–34.
Martin SD, et al. Low mutation burden in ovarian Cancer may limit the utility of Neoantigen-targeted vaccines. PLoS One. 2016;11:e0155189.
Bobisse S, et al. Sensitive and frequent identification of high avidity neo-epitope specific CD8 (+) T cells in immunotherapy-naive ovarian cancer. Nat Commun. 2018;9:1092.
Deniger DC, et al. T-cell responses to TP53 "hotspot" mutations and unique Neoantigens expressed by human ovarian cancers. Clin Cancer Res. 2018;24:5562–73.
Wei L, et al. Intratumoral and Intertumoral genomic heterogeneity of multifocal localized prostate Cancer impacts molecular classifications and genomic prognosticators. Eur Urol. 2016.
Wei L, et al. Whole-genome sequencing of a malignant granular cell tumor with metabolic response to pazopanib. Molecular Case Studies. 2015;1.
Futreal PA, et al. A census of human cancer genes. Nat Rev Cancer. 2004;4:177–83.
Duan F, et al. Genomic and bioinformatic profiling of mutational neoepitopes reveals new rules to predict anticancer immunogenicity. J Exp Med. 2014;211:2231–48.
Rajasagi M, et al. Systematic identification of personal tumor-specific neoantigens in chronic lymphocytic leukemia. Blood. 2014;124:453–62.
Sykulev Y, Joo M, Vturina I, Tsomides TJ, Eisen HN. Evidence that a single peptide-MHC complex on a target cell can elicit a cytolytic T cell response. Immunity. 1996;4:565–71.
Casey SC, et al. MYC regulates the antitumor immune response through CD47 and PD-L1. Science. 2016;352:227–31.
Yoshihara K, et al. High-risk ovarian cancer based on 126-gene expression signature is uniquely characterized by downregulation of antigen presentation pathway. Clin Cancer Res. 2012;18:1374–85.
Tsuji T, et al. Rapid construction of antitumor T-cell receptor vectors from frozen tumors for engineered T-cell therapy. Cancer Immunol Res. 2018.
Ayyoub M, Pignon P, Classe JM, Odunsi K, Valmori D. CD4+ T effectors specific for the tumor antigen NY-ESO-1 are highly enriched at ovarian cancer sites and coexist with, but are distinct from, tumor-associated Treg. Cancer Immunol Res. 2013;1:303–8.
Tanyi JL, et al. Personalized cancer vaccine effectively mobilizes antitumor T cell immunity in ovarian cancer. Sci Transl Med. 2018;10.
Grimwade D, et al. Refinement of cytogenetic classification in acute myeloid leukemia: determination of prognostic significance of rare recurring chromosomal abnormalities among 5876 younger adult patients treated in the United Kingdom Medical Research Council trials. Blood. 2010;116:354–65.
Oancea C, Ruster B, Henschler R, Puccetti E, Ruthardt M. The t (6;9) associated DEK/CAN fusion protein targets a population of long-term repopulating hematopoietic stem cells for leukemogenic transformation. Leukemia. 2010;24:1910–9.
Shin, D.S. et al. Primary Resistance to PD-1 Blockade Mediated by <em>JAK1/2</em> Mutations. Cancer Discovery 7, 188 (2017).
Zaretsky JM, et al. Mutations associated with acquired resistance to PD-1 blockade in melanoma. N Engl J Med. 2016;375:819–29.
Wen W, et al. Targeting JAK1/STAT3 signaling suppresses tumor progression and metastasis in a peritoneal model of human ovarian cancer. Mol Cancer Ther. 2014;13:3037–48.
Gritsina G, et al. Targeted blockade of JAK/STAT3 signaling inhibits ovarian carcinoma growth. Mol Cancer Ther. 2015;14:1035.
Wen W, et al. Synergistic anti-tumor effect of combined inhibition of EGFR and JAK/STAT3 pathways in human ovarian cancer. Mol Cancer. 2015;14:100.
Di Carlo E, et al. The antitumor potential of Interleukin-27 in prostate cancer. Oncotarget. 2014;5:10332–41.
Dibra D, et al. Lack of Immunomodulatory Interleukin-27 enhances oncogenic properties of mutant p53 in vivo. Clin Cancer Res. 2016;22:3876–83.
Petretto A, et al. Proteomic analysis uncovers common effects of IFN-γ and IL-27 on the HLA class I antigen presentation machinery in human cancer cells. Oncotarget. 2016;7:72518–36.
Gao J, et al. Loss of IFN-γ Pathway Genes in Tumor Cells as a Mechanism of Resistance to Anti-CTLA-4 Therapy. Cell. 2016;167:397–404.e9.
Irvine DJ, Purbhoo MA, Krogsgaard M, Davis MM. Direct observation of ligand recognition by T cells. Nature. 2002;419:845–9.
Anagnostou V, et al. Evolution of Neoantigen landscape during immune checkpoint blockade in non-small cell lung Cancer. Cancer Discov. 2017;7:264–76.
Riaz N, et al. Tumor and Microenvironment Evolution during Immunotherapy with Nivolumab. Cell. 2017;171:934–949 e15.
Wilke CM, Kryczek I, Zou W. Antigen-presenting cell (APC) subsets in ovarian cancer. Int Rev Immunol. 2011;30:120–6.
Yossef R, et al. Enhanced detection of neoantigen-reactive T cells targeting unique and shared oncogenes for personalized cancer immunotherapy. JCI Insight. 2018;3.
Rosenberg SA. Cell transfer immunotherapy for metastatic solid cancer--what clinicians need to know. Nat Rev Clin Oncol. 2011;8:577–85.
Yang F, et al. Adoptive cellular therapy (ACT) for Cancer treatment. Adv Exp Med Biol. 2016;909:169–239.
Keskin DB, et al. Neoantigen vaccine generates intratumoral T cell responses in phase Ib glioblastoma trial. Nature. 2019;565:234–9.
Nowicki TS, et al. A pilot trial of the combination of transgenic NY-ESO-1-reactive adoptive cellular therapy with dendritic cell vaccination with or without Ipilimumab. Clin Cancer Res. 2019;25:2096–108.
Schoenberger SP. Is It Possible to Develop Cancer Vaccines to Neoantigens, What Are the Major Challenges, and How Can These Be Overcome? Targeting the Right Antigens in the Right Patients. Cold Spring Harb Perspect Biol. 2018;10.
Tavassoéli, F.A., Devilee, P. Pathology and Genetics of Tumours of the Breast and Female Genital Organs, (2003).
Li H, Durbin R. Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009;25:1754–60.
Edmonson MN, et al. Bambino: a variant detector and alignment viewer for next-generation sequencing data in the SAM/BAM format. Bioinformatics. 2011;27:865–6.
Saunders CT, et al. Strelka: accurate somatic small-variant calling from sequenced tumor-normal sample pairs. Bioinformatics. 2012;28:1811–7.
Holmfeldt L, et al. The genomic landscape of hypodiploid acute lymphoblastic leukemia. Nat Genet. 2013;45:242–52.
Sherry ST, et al. dbSNP: the NCBI database of genetic variation. Nucleic Acids Res. 2001;29:308–11.
Genomes Project C, et al. An integrated map of genetic variation from 1,092 human genomes. Nature. 2012;491:56–65.
Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38:e164.
Koboldt DC, et al. VarScan 2: somatic mutation and copy number alteration discovery in cancer by exome sequencing. Genome Res. 2012;22:568–76.
Huber W, et al. Orchestrating high-throughput genomic analysis with bioconductor. Nat Methods. 2015;12:115–21.
Mermel, C.H. et al. GISTIC2.0 facilitates sensitive and confident localization of the targets of focal somatic copy-number alteration in human cancers. Genome Biol 12, R41 (2011).
Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25:1105–11.
Wang L, Wang S, Li W. RSeQC: quality control of RNA-seq experiments. Bioinformatics. 2012;28:2184–5.
Anders S, Pyl PT, Huber W. HTSeq--a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31:166–9.
Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12:323.
Barbie DA, et al. Systematic RNA interference reveals that oncogenic KRAS-driven cancers require TBK1. Nature. 2009;462:108–12.
Subramanian A, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102:15545–50.
Liberzon A, et al. Molecular signatures database (MSigDB) 3.0. Bioinformatics. 2011;27:1739–40.
Law, C.W., Chen, Y., Shi, W. & Smyth, G.K. voom: Precision weights unlock linear model analysis tools for RNA-seq read counts. Genome Biol 15, R29 (2014).
Cerami E, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012;2:401–4.
Aran D, Hu Z, Butte AJ. xCell: digitally portraying the tissue cellular heterogeneity landscape. Genome Biol. 2017;18:220.
Pruitt KD, Tatusova T, Maglott DR. NCBI reference sequence (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2005;33:D501–4.
Shukla SA, et al. Comprehensive analysis of cancer-associated somatic mutations in class I HLA genes. Nat Biotechnol. 2015;33:1152–8.
Warren RL, et al. Derivation of HLA types from shotgun sequence datasets. Genome Med. 2012;4:95.
Andreatta M, et al. Accurate pan-specific prediction of peptide-MHC class II binding affinity with improved binding core identification. Immunogenetics. 2015;67:641–50.
Hoof I, et al. NetMHCpan, a method for MHC class I binding prediction beyond humans. Immunogenetics. 2009;61:1–13.
Nielsen, M. & Andreatta, M. NetMHCpan-3.0; improved prediction of binding to MHC class I molecules integrating information from multiple receptor and peptide length datasets. Genome Med 8, 33 (2016).
Yadav M, et al. Predicting immunogenic tumour mutations by combining mass spectrometry and exome sequencing. Nature. 2014;515:572–6.
Odunsi K, et al. Efficacy of vaccination with recombinant vaccinia and fowlpox vectors expressing NY-ESO-1 antigen in ovarian cancer and melanoma patients. Proc Natl Acad Sci U S A. 2012;109:5797–802.
Atanackovic D, et al. Monitoring CD4+ T cell responses against viral and tumor antigens using T cells as novel target APC. J Immunol Methods. 2003;278:57–66.
Tsuji T, et al. Antibody-targeted NY-ESO-1 to mannose receptor or DEC-205 in vitro elicits dual human CD8+ and CD4+ T cell responses with broad antigen specificity. J Immunol. 2011;186:1218–27.
Tsuji T, Altorki NK, Ritter G, Old LJ, Gnjatic S. Characterization of preexisting MAGE-A3-specific CD4+ T cells in cancer patients and healthy individuals and their activation by protein vaccination. J Immunol. 2009;183:4800–8.
Matsuzaki J, et al. Nonclassical antigen-processing pathways are required for MHC class II-restricted direct tumor recognition by NY-ESO-1-specific CD4(+) T cells. Cancer Immunol Res. 2014;2:341–50.
Matsuzaki J, et al. Direct tumor recognition by a human CD4(+) T-cell subset potently mediates tumor growth inhibition and orchestrates anti-tumor immune responses. Sci Rep. 2015;5:14896.
We thank Kristen Starbuck, MD for material support.
This work was supported by the RPCI-UPCI Ovarian Cancer SPORE P50 CA159981, NCI Cancer Center Support Grant P30 CA016056, NIH U24 CA232979, NIH R01 CA158318, the Empire State Stem Cell Board from the New York State Health Department (NYSTEM) N14C-002, Ovarian Cancer Research Alliance, and the Roswell Park Alliance Foundation.
Ethics approval and consent to participate
Biospecimens of epithelial ovarian cancer patients were obtained with informed consent at Roswell Park Comprehensive Cancer Center in accordance with a protocol approved by the institutional review board.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. The protein view of selected genes with somatic mutations in the 20 EOC patients studied. (a) TP53. (b) BRCA1. (c) BRCA2. (d) IL27RA. Note that splicing site mutations (c.673-2A > G in TP53 and c.4987-1G > T in BRAC1) are not shown. (PPTX 130 kb)
Figure S2. Comparison of somatic mutation and neoantigen between primary and local invasive tumors. (a) Somatic mutation. (b) Neoantigen. Pt #02 on the left and Pt #10 on the right. X-axis: variant allele frequency (VAF) in WES of primary tumor. Y-axis: VAF in WES of local invasive tumor. (PPTX 260 kb)
Figure S3. Spontaneous CD4+ T-cell response against mutant-specific epitopes from autologous tumor-infiltrating lymphocytes (TIL) and/or peripheral blood lymphocytes (PBL). T-cell reactivity was measured by IFN-γ ELISpot assay (Methods). (−): no peptide; Wt: wildtype peptide; Mt.: Mutant peptide. (PPTX 85 kb)
Figure S4. Spontaneous CD8+ T-cell response against mutant-specific epitopes from autologous tumor-infiltrating lymphocytes (TIL) and/or peripheral blood lymphocytes (PBL). T-cell reactivity was measured by IFN-γ ELISpot assay (Methods). (−): no peptide; Wt: wildtype peptide; Mt.: Mutant peptide. (PPTX 84 kb)
Figure S5. The expression level of a panel of 10 immune inhibitory molecules from the tumor RNA-Seq data. These immune inhibitory molecules include PD1, PDL1, CTLA4, CD80, CD86, LAG3, TIM3, LAGLS9, MYC and FOXP3. The patients shown here include Pt #2, Pt #8, Pt #16 and Pt #5. (PPTX 78 kb)
Figure S6. Recurrent somatic copy number alternations by the GISTC2.0 algorithm. GISTIC deletion (left) and amplification (right) plots using data from the five patients with mutant-specific T-cell response (top), and data from the five patients without mutant-specific T-cell response (bottom). The genome is oriented vertically from top to bottom, and GISTIC q-values at each locus are plotted from left to right on a log scale. The green line represents the default significance threshold (q-value = 0.25). For each plot, known or interesting cancer genes are highlighted. (PPTX 278 kb)
Figure S7. TCR Vβ usage of NUP214 neopeptide-specific CD4+ T-cell line. (a) TCR Vβ usage of NUP214 neopeptide-specific CD4+ T-cell line was determined by flow cytometry. (b) NUP214 neopeptide-specific IFN-γ production from Vβ2+, Vβ13.1+, or Vβ2−Vβ13.1− T-cells was examined by intracellular cytokine staining. (PPTX 221 kb)
Figure S8. Characterization of JAK1 neoepitope-specific CD4+ T-cells. (a) Peptide reactivity of a JAK1 neoepitope-specific CD4+ T-cell line. IFN-γ and GM-CSF production on CD4+ T-cells against JAK1 mutated (IEILRNLYHEIIV) or wild-type (IEILRNLYHENIV) peptide-pulsed autologous EBV-B-cells were determined by intracellular cytokine staining. (b) TCR usage of JAK1 neoepitope-specific CD4+ T-cell line. T-cells were stained with TCR Vβ subtype-specific antibodies and analyzed by flow cytometry. (c) Purity of Vβ13.6+ cells after magnetic-beads sorting. (d) Avidity of JAK1 neoepitope-specific T-cell clone. CD4+ T-cell clones were stimulated with autologous EBV-B-cells pulsed with the indicated concentration of mutated or wild-type peptide for 6 h in the presence of Golgi stop. IFN-γ production from Vβ13.6+ cells were determined by flow cytometry. The data represents mean ± s.d. of duplicate wells. (e) Recognition of autologous tumor-derived cells by Vβ13.6+ T-cell clone. PBMC or AMC were co-cultured with Vβ13.6+ JAK1 neoepitope-specific CD4+ T clones or without T-cells (−) for 24 h. AMC: ascites-derived mononuclear cells. IFN-γ production was measured by ELISA. The data represent mean + s.d. of duplicate wells. *p < 0.05 compared to IFN-γ level against PBMC. (PPTX 250 kb)
Figure S9. Reactivity of NUP214 neoepitope-specific TCR-transduced CD8+ and CD4+ T-cells. IFN-γ production (a) and GM-CSF (b) production from Vβ13.1+ or Vβ2−Vβ13.1− TCR- transduced CD8+ and CD8− (CD4+) T-cells against EBV-B-cells pulsed with or without NUP214 mutated or wild-type peptide were determined by intracellular cytokine staining. (PPTX 683 kb)
Figure S10. Analysis of TRPC4 neoepitope-specific CD8+ T-cell clone. (a) IFN-γ and GM-CSF production on CD8+ T-cells against TRPC4 mutated (QSLFWSIFV) or wild-type (QSLFWSIFG) peptide-pulsed autologous EBV-B-cells were determined by intracellular cytokine staining. (b) Avidity of TRPC4 neoepitope-specific T-cell clone. CD8+ T-cell clone was co-cultured with autologous EBV-B-cells pulsed with the different concentration of mutated or wild-type peptide for 24 h. IFN-γ level in the culture supernatant was measured by ELISA. The data represents mean ± s.d. of duplicate wells. (c) IFN-γ production from TRPC4 neoepitope-specific TCR-transduced T-cells against EBV-B pulsed with or without TRPC4 mutated or wild-type peptide was determined by intracellular cytokine staining. (PPTX 321 kb)
Figure S11. The expression level of APPM signature from the tumor RNA-Seq data. The APPM expression is derived from the median expression value (in RPKM) of the 31 genes within the APPM signature (Methods). The patients shown here include Pt #3, Pt #15, Pt #19 and Pt #5. (PPTX 50 kb)
Table S1. Clinical demographics of treatment-naïve (including both chemotherapy and immunotherapy) epithelial ovarian cancer patients at the time of primary debulking surgery. Patient identification number, age at diagnosis, tumor histologic type, FIGO stage, debulking status, residual tumor mass after debulking surgery, number of recurrences after primary debulking surgery, and RECIST to frontline chemotherapy after primary debulking surgery. (XLSX 10 kb)
Table S2. Somatic point mutations identified from whole-exome sequencing. (a) The 18 patients with either primary tumor or locally invasive tumor; (b) The 2 patients with both primary tumor and locally invasive tumor. AA, amino acid; CGC, cancer gene census; VAF, variant allele frequency. (XLSX 244 kb)
Table S3. Somatic mutation burdens and predicted neoantigen load in the 20 patients. The predicted neoantigens are classified as expressed or non-expressed based on the mutant allele’s expression level in RNAseq data (see Method section). (XLSX 10 kb)
Table S4. Description of the 75 neopeptides screened for immunogenicity. The expression status is based on the mutant allele’s expression level in RNAseq data, and affinity score is predicted by NetMHC algorithm with default setting (see Method section). (XLSX 15 kb)
Table S5. The list of 31 genes within the derived APPM signature (see Method section). (XLSX 11 kb)
About this article
Cite this article
Liu, S., Matsuzaki, J., Wei, L. et al. Efficient identification of neoantigen-specific T-cell responses in advanced human ovarian cancer. j. immunotherapy cancer 7, 156 (2019). https://doi.org/10.1186/s40425-019-0629-6