Adoptive transfer of T cells with engineered T-cell receptor (TCR) genes that target tumor-specific antigens can mediate cancer regression. Accumulating evidence suggests that the clinical success of many immunotherapies is mediated by T cells targeting mutated neoantigens unique to the patient. We hypothesized that the most frequent TCR clonotypes infiltrating the tumor were reactive against tumor antigens. To test this hypothesis, we developed a multistep strategy that involved TCRB deep sequencing of the CD8+PD-1+ T-cell subset, matching of TCRA–TCRB pairs by pairSEQ and single-cell RT-PCR, followed by testing of the TCRs for tumor-antigen specificity. Analysis of 12 fresh metastatic melanomas revealed that in 11 samples, up to 5 tumor-reactive TCRs were present in the 5 most frequently occurring clonotypes, which included reactivity against neoantigens. These data show the feasibility of developing a rapid, personalized TCR-gene therapy approach that targets the unique set of antigens presented by the autologous tumor without the need to identify their immunologic reactivity. Cancer Immunol Res; 4(9); 734–43. ©2016 AACR.

The presence of lymphocytes infiltrating into the tumor stroma (tumor-infiltrating lymphocytes; TIL) has been associated with a favorable prognosis in melanoma (1) and other cancer types, including ovarian (2), colon (3), and breast cancer (4). In melanoma, in vitro analysis of expanded TILs revealed a broad specificity of antigen recognition, including melanoma/melanocyte shared differentiation antigens (5–7), cancer germline antigens (8, 9), and mutated neoantigens unique to each patient's tumor (10–12).

Adoptive cell therapy using autologous TILs is an immunotherapeutic approach capable of inducing complete durable regression in 20% of patients with metastatic melanoma (13). However, TILs used for treatment undergo extensive in vivo and in vitro expansion, becoming highly differentiated cells with limited additional proliferative potential (13, 14). Control over which T-cell clonotypes expand in vitro is limited, so the TCR clonotypic repertoire present in the tumor can be altered, potentially leading to decreased frequencies of tumor-reactive clonotypes.

To overcome these problems, we focused our attention on the TCR clonotypes present in the tumor before any in vitro expansion. In melanoma, tumor-specific clonotypes are highly enriched in the fresh CD8+PD-1+ TIL subset (15, 16), which we hypothesize could be due to the oligoclonal expansion that occurs when T cells encounter their specific antigen in the tumor microenvironment in vivo (17), leading to the presence of predominant clonotypes within this population. Thus, the frequency of a clonotype within the TIL repertoire may indicate its tumor reactivity. To test this hypothesis, we analyzed the TCR repertoire of TILs in freshly resected tumors from 12 patients with metastatic melanoma and found that many of the most frequent TCR clonotypes present in the CD8+PD-1+ TIL subset recognized the autologous tumor and either mutated or nonmutated tumor antigens. Thus, it may be possible to efficiently identify tumor-reactive TCRs based solely on their frequency and PD-1 expression in the tumor. This strategy could provide an efficient means to obtain tumor-reactive TCRs that could be genetically engineered into autologous cells with high proliferative potential for use in cell therapy.

Tumor samples

Twelve metastatic melanoma samples were obtained from patients who were not undergoing therapy at the time of sample collection. Patients had undergone a wide range of prior therapies, including surgery, chemotherapy, radiotherapy, immunotherapy, or none of the above. Peripheral blood lymphocytes (PBL) were obtained by either leukapheresis or venipuncture, prepared over Ficoll–Hypaque gradient (LSM; ICN Biomedicals Inc.), and cryopreserved until analysis. After surgical resection, tumor specimens were processed as previously described (18). Briefly, tumor specimens were minced, enzymatically digested overnight at room temperature or for several hours at 37°C [RPMI-1640 with l-glutamine (Lonza), 1 mg/mL collagenase IV (Sigma-Aldrich), 30 U/mL DNAse (Genentech), and antibiotics], and the tissue was separated mechanically using gentleMACS (Miltenyi Biotech). Tumor single-cell suspensions were cryopreserved.

Whole-exome sequencing and RNA sequencing

Genomic DNA purification, library construction, exome capture of approximately 20,000 coding genes, next-generation sequencing of fresh tumor embedded in optimal cutting temperature (Sakura Finetek), and a matched normal pheresis sample were performed as previously described (19). An mRNA sequencing library was prepared from fresh tumors using the Illumina TruSeq RNA library prep kit, as previously described (20). Putative nonsynonymous mutations are defined by ≥3 exome variant reads, ≥8% variant allele fraction (VAF) in the exome, and ≥10 reads in the matched normal sample. Putative mutations with a variant allele frequency of >10% in the tumor exome, as well as mutations that were identified in both transcriptome and exome analysis are initially selected for screening. For some samples (i.e., 3903), the mutations selected based on exome were prioritized only by selecting those with >10 variant reads to increase the confidence of mutation calling.

Antibodies, flow cytometry, and cell sorting

Fluorescently conjugated antibodies were purchased from eBioscience [MIH-4, Anti-Human CD279 conjugated to allophycocyanin (APC) and anti-mouse TCRβ-fluorescein isothiocyanate (FITC)], Miltenyi (4B4-1, anti-human CD137-PE or -APC), BioLegend [anti-human CD8-phycoerythrin (PE)-Cy7, anti-human CD3-APC-Cy7]. For phenotypic characterization and cell sorting of CD8+/−, CD8+PD-1+/− T-cell tumor samples were thawed and rested overnight without cytokines (15). The T cells were sorted by flow cytometry with a modified FACSAria instrument or a BD Jazz instrument (BD Biosciences), gates were set according to isotype and fluorescence minus one controls. The sorting strategy is shown for two representative fresh melanoma samples (3903 and 3998) in Supplementary Fig. S1.

Sample preparation for ImmunoSEQ TCRB deep sequencing and pairSEQ

The T cells were sorted by flow cytometry in comparable numbers for each subset: 100,000 cells for the tumor single-cell suspension bulk TILs, 10,000 cells for the CD8+ and CD8 subsets, and 1,000 to 3,000 cells for the CD8+PD-1+ and CD8+PD-1 subsets. The cells were pelleted and snap frozen. The samples were sent to Adaptive Technologies for genomic DNA extraction and ImmunoSEQ TCRB survey sequencing. Tumor samples were sent to Adaptive Technologies for pairSEQ (21), 1 × 106 total cells from tumor single-cell suspension were pelleted in a table top centrifuge at 6,000 rpm for 30 minutes, resuspended in 200 μL of RNAlater (Invitrogen) and snap frozen.

Single-cell sorting and single-cell RT-PCR

Single-cell sorting was performed using a modified FACSAria instrument or BD Jazz instrument (BD Biosciences) on CD8+PD-1+ TILs; for samples 1913, 2650, 3713, and 3784 CD8+-expanded TILs were used due to limited availability of tumor samples. TCR sequences from the sorted single cells were obtained by a series of two nested PCR reactions. Cells were sorted into RT-PCR buffer. For the first, reverse transcription and amplification reaction were performed with a One-Step RT-PCR kit (Qiagen) using multiplex PCR with multiple Vα and Vβ region primers and one primer for Cα and Cβ regions each (final concentration of each primer is 0.6 μmol/L). The RT-PCR reaction was performed according to the manufacturer's instructions using the following cycling conditions: 50°C for 15 minutes; 95°C for 2 minutes; 95°C for 15 seconds, 60°C for 4 minutes × 18 cycles; 4°C. For the second amplification reaction, 4 μL from the first RT-PCR was used as a template in total 25 μL PCR mix with HotStarTaq DNA polymerase (Qiagen) and multiple internally nested Vα and Vβ region primers and one internally nested primer for Cα and Cβ region each (final concentration of each primer was 0.6 μmol/L). The cycling conditions were 95°C for 15 minutes; 94°C for 30 seconds, 50°C for 30 seconds, 72°C for 1 minute × 50 cycles; 72°C for 10 minutes; 4°C. The PCR products were purified and sequenced by the Sanger method with internally nested Cα and Cβ region primers by Beckmann Coulter. All primers are listed in Supplementary Table S1.

TCR pairs’ reconstruction, cloning into expression vectors, and TCR expression evaluation

In both pairing methods (single-cell RT-PCR and pairSEQ), cDNA is used as template for multiplex PCR using TCRA and TCRB gene–specific primers. The resulting PCR product contains the 3′ end of the variable region and the full CDR3 region of matching TCRA and TCRB genes. These partial TCR sequences were analyzed with IMGT/V-Quest tool (http://www.imgt.org/) which identified the TRAV and TRBV families with the highest likelihood to contain the segment found with our pairing methods. Utilizing the IMGT database, we reconstructed the full-length TRAV and TRBV regions for each pair. In regard to the constant regions, we used modified murine TRAC and TRBC sequences to improve stability and avoid mismatches with the endogenous human TCR after transduction into human T cells (22). Full TCR genes were synthetized, and a 2A peptide (23) was introduced between the TCRB and TCRA chains to ensure a comparable expression efficiency of the two chains. The resulting TCRB–TCRA gene blocks were cloned into either a gamma-retroviral expression vector (24, 25) or for the following TCR pairs: 2650-1, 2650-3, 2650-4, 2650-5, 2650-6, 2650-7, 2650-9, 3903-3A1, 3903-3A2, 3992-1, 3992-2, 3992-3, 3992-4, 3992-5, and 3998-1 into a nonviral Sleeping Beauty transposon system (26, 27). The expression of the TCRB was evaluated with an anti-murine TCRB Ab.

Target cell preparation

Melanoma tumor cell lines (TC 1913, TC 2630, TC 2650, TC 3678, TC 3713, TC 3759, TC 3784, TC 3903, TC 3922, TC 3926, TC 3977, TC 3992, and TC 3998) were established from tumor fragments or from mechanically or enzymatically separated tumor cells and cultured in RPMI 1640 plus 10% FBS (Sigma-Aldrich) supplemented with 100 U/mL penicillin and 100 μg/mL streptomycin at 37°C in 5% CO2. COS-7 cells and COS-7 cells stably transduced with HLA molecules were maintained in DMEM containing 10% FBS (Sigma-Aldrich) supplemented with 100 U/mL penicillin and 100 μg/mL streptomycin at 37°C in 5% CO2. TC 1913 and 1913 tumor–specific neoantigen recognition was previously reported (28). Generation of tandem minigene (TMG) constructs and autologous antigen presenting cells (dendritic cells and CD40L stimulated B cells) was done as previously described (11, 29). Briefly, up to 14 nonsynonymous mutations identified by whole-exome sequencing and RNAseq, each flanked by 12 amino acids of nonmutated protein, were genetically fused together to generate a TMG construct. These constructs were codon optimized, synthesized, and cloned into pcDNA3.1/V5 by Genescript. Autologous antigen-presenting cells were peptide pulsed for 2 hours with 1 μg/mL short peptides (9–10 mers) and overnight with 10 μg/mL of long peptides (25 mers) before coculture.

Target cell recognition functional assay

CD137 upregulation was used to measure target recognition by transduced T cells. CD137 is upregulated transiently in response to TCR stimulation, regardless of the effector cytokines produced or the differentiation state of the cell (30). We used the coexpression of murine TCR constant chain (identified as mTCR) and CD137 to identify the population of transduced antigen-reactive T cells (to be considered reactive the CD137 upregulation had to be greater than 1%, 3 times the background, and inhibited at least 50% by pan–MHC I blocking antibody, clone W6-32). Cells were stained with anti-CD3, anti-CD8, anti-CD137, and anti-murine TCRB antibodies after coculture and acquired by Fortessa (BD Biosciences). Data were analyzed with FlowJo software (Treestar).

Statistical analysis

Wilcoxon signed-rank test was used to determine the statistical significance of the data. P values of 0.05 or less were considered statistically significant. Statistical calculations were performed with Prism program 6.0 (GraphPad Software Inc.).

Study approval

All patient samples were obtained in the course of a National Cancer Institute Institutional Review Board–approved clinical trial. Patients provided informed consent.

CD8+PD-1+ TIL clonotypes are oligoclonal compared with CD8+PD-1 TILs

To characterize the TIL TCR clonotypic repertoire and identify tumor-reactive TCRs, we developed a multistep strategy (Fig. 1A). We first assessed the composition of TILs from 12 fresh metastatic melanoma lesions by flow cytometry (Table 1). The samples varied considerably in the frequency of CD8+ and CD8 lymphocytes (P = 0.15), although the frequency of PD-1 expression, which is a marker for T-cell activation (16, 31), was usually higher on the CD8+ TILs (P = 0.003). TCRB deep sequencing is a robust method for quantifying the frequency of each T-cell clonotype present in different sample types (32–34), which we used to determine whether the CD8+PD-1+ TIL subset displayed evidence of clonal expansion. Genomic DNA extracted from bulk melanoma TILs and from sorted subsets (CD8+, CD8, CD8+PD-1+, and CD8+PD-1) was deep sequenced to determine the number of unique productive (35) TCRB CDR3 sequences that do not contain stop codons or frame shifts (Fig. 1B) in the 10 patients from whom all samples were available. These unique sequences represent a single unique clonotype independent of its frequency in the samples. Different samples can thus have comparable numbers of total reads (Fig. 1C), but it is the different number of unique sequences that determines the level of clonality of each sample. We found a significantly lower number of unique productive sequences present in the CD8+ compared with the CD8 subsets (P = 0.002, Fig. 1B). Within the CD8+ lymphocytes, the PD-1+ population contained a lower number of unique productive sequences compared with PD-1 cells (P = 0.002, Fig. 1B). We also compared the levels of clonal diversity (measured by Shannon entropy; ref. 36) for the samples studied and we found more diversity in the CD8 subset compared with the CD8+ (P = 0.002) and within the CD8+ lymphocytes more diversity in the PD-1 cells compared with the PD-1+ subset (P = 0.002, Supplementary Fig. S2A). We also compared the number of nonsynonymous mutations with the percentage of CD8+PD-1+ TILs and TCR clonal diversity in different TIL subsets (Supplementary Fig. S2B–S2E), but found no significant correlation. Each highly expressed individual clonotype in the CD8+PD-1+ subset was much less frequent in the PD-1 group (P = 0.0003; Supplementary Fig. S3A), confirming that PD-1 is a marker that separates TILs into two separate subsets with different TCR repertoires (15). The same clonotypes that were highly expressed in the CD8+PD-1+ subset were found at low frequency in the total CD8+ subset (Supplementary Fig. S3B, P = 0.001). The highest expressers in the PD-1+ subset though were often also high ranked in the bulk TILs (Supplementary Table S2).

Figure 1.

Strategy overview and TIL characterization. A, schematic representation of the multistep process used to identify tumor-reactive TCRs. (1) TCRB deep sequencing on bulk TIL and sorted CD8+/− TIL and CD8+PD-1+/− TIL populations is used to determine the subset with evidence of clonal expansion and to identify the TCRB sequences of the most dominant clonotypes within that subset. (2) The most dominant TCRB clonotypes in CD8+PD-1+ TILs are paired with TCRA chains identified by single-cell RT-PCR and pairSEQ. (3) TCRA–TCRB pairs are cloned into expression vectors and engineered into T cells. (4) Engineered T cells are tested for tumor reactivity against tumor cell lines, shared tumor antigens and mutated neoantigens. B, unique productive TCRB clonotypes sequences are plotted for bulk TIL and sorted TIL subsets. Productive sequences do not contain stop codons or out of frameshifts so they are likely to be functional. These unique sequences represent a single unique clonotype independent of its frequency in the samples. Wilcoxon matched-pairs signed rank test was applied (n = 10). For samples 1913 and 3922, the CD8 subset was not available. C, total reads of TCRB clonotype sequences are plotted for melanoma bulk TIL and sorted TIL subsets. Wilcoxon matched-pairs signed rank test was applied (n = 10). For samples 1913 and 3922, the CD8 subset was not available.

Figure 1.

Strategy overview and TIL characterization. A, schematic representation of the multistep process used to identify tumor-reactive TCRs. (1) TCRB deep sequencing on bulk TIL and sorted CD8+/− TIL and CD8+PD-1+/− TIL populations is used to determine the subset with evidence of clonal expansion and to identify the TCRB sequences of the most dominant clonotypes within that subset. (2) The most dominant TCRB clonotypes in CD8+PD-1+ TILs are paired with TCRA chains identified by single-cell RT-PCR and pairSEQ. (3) TCRA–TCRB pairs are cloned into expression vectors and engineered into T cells. (4) Engineered T cells are tested for tumor reactivity against tumor cell lines, shared tumor antigens and mutated neoantigens. B, unique productive TCRB clonotypes sequences are plotted for bulk TIL and sorted TIL subsets. Productive sequences do not contain stop codons or out of frameshifts so they are likely to be functional. These unique sequences represent a single unique clonotype independent of its frequency in the samples. Wilcoxon matched-pairs signed rank test was applied (n = 10). For samples 1913 and 3922, the CD8 subset was not available. C, total reads of TCRB clonotype sequences are plotted for melanoma bulk TIL and sorted TIL subsets. Wilcoxon matched-pairs signed rank test was applied (n = 10). For samples 1913 and 3922, the CD8 subset was not available.

Close modal
Table 1.

Characteristics of infiltrating cells in fresh tumors

Patient IDAgeSexTumor locationLymphocytes (% of viable cells)CD3+ (% of lymphocytes)CD3+CD8+ (% of CD3+)CD3+CD8 (% of CD3+)CD3+CD8+ PD-1+ (% of CD3+CD8+)CD3+CD8 PD-1+ (% of CD3+CD)Tumor cell line (TC) availableNonsynonymous mutationsa
1913 40 Subcutaneous 42.9 38.4 67.2 20.1 17.2 24.4 Yes 3,280 
2650 56 Lymph node 9.21 8.52 36.5 62.5 12 2.38 Yes 431 
3678 55 Lymph node 20.7 16.1 51.0 46.5 11.3 12.7 No 504 
3713 53 Lung 28.1 19.7 32.7 33.9 55.2 45.7 Yes 3,976 
3759 21 Subcutaneous 2.06 1.42 53.4 18.4 81.0 12.3 Yes 1,378 
3784 45 Lymph node 18.1 12 35.2 22.2 31.6 15.9 Yes 662 
3903 56 Liver 31.4 26.2 53.3 30.8 74.8 33.7 Yes 385 
3922 65 Lymph node 15.9 7.2 24.5 73.3 11.4 2.84 Yes 449 
3926 37 Lymph node 22.8 15.7 42.3 28.5 34.4 18.5 Yes 340 
3977 51 Lymph node 17.7 15.4 50.3 47 35.7 16.6 Yes N/A 
3992 58 Lymph node 79 16.8 36.4 28.2 15.4 9.78 Yes 159 
3998 30 Liver 17.7 15.8 70.2 15 61.9 11.3 Yes 345 
Patient IDAgeSexTumor locationLymphocytes (% of viable cells)CD3+ (% of lymphocytes)CD3+CD8+ (% of CD3+)CD3+CD8 (% of CD3+)CD3+CD8+ PD-1+ (% of CD3+CD8+)CD3+CD8 PD-1+ (% of CD3+CD)Tumor cell line (TC) availableNonsynonymous mutationsa
1913 40 Subcutaneous 42.9 38.4 67.2 20.1 17.2 24.4 Yes 3,280 
2650 56 Lymph node 9.21 8.52 36.5 62.5 12 2.38 Yes 431 
3678 55 Lymph node 20.7 16.1 51.0 46.5 11.3 12.7 No 504 
3713 53 Lung 28.1 19.7 32.7 33.9 55.2 45.7 Yes 3,976 
3759 21 Subcutaneous 2.06 1.42 53.4 18.4 81.0 12.3 Yes 1,378 
3784 45 Lymph node 18.1 12 35.2 22.2 31.6 15.9 Yes 662 
3903 56 Liver 31.4 26.2 53.3 30.8 74.8 33.7 Yes 385 
3922 65 Lymph node 15.9 7.2 24.5 73.3 11.4 2.84 Yes 449 
3926 37 Lymph node 22.8 15.7 42.3 28.5 34.4 18.5 Yes 340 
3977 51 Lymph node 17.7 15.4 50.3 47 35.7 16.6 Yes N/A 
3992 58 Lymph node 79 16.8 36.4 28.2 15.4 9.78 Yes 159 
3998 30 Liver 17.7 15.8 70.2 15 61.9 11.3 Yes 345 

aPutative nonsynonymous mutations defined by ≥3 exome variant reads, ≥8% variant allele fraction (VAF) in the exome, ≥10 reads in the matched normal sample.

Identification and reconstruction of TCR pairs for the most frequent CD8+PD-1+ TIL clonotypes

To reconstruct a functional TCR, the most frequent TCRB chains present in the CD8+PD-1+ TIL must be paired with the appropriate TCRA chains. One possible approach was to match the most frequent TCRB clonotype with the most frequent TCRA clonotype. However, in the 6 patient samples for which we ultimately matched the correct TCRA with the most frequent TCRB and determined their functionality, the most frequent TCR chains paired together in 3 cases, and the first TCRB in the other 3 cases paired with the 7th, 18th, and 47th TCRA (Supplementary Table S3). Tumor recognition by the first TCRB with the first TCRA occurred in the 3 cases in which the first TCRB clonotype was present in >20% of the PD-1+ population. Discordance of pairing based on frequency of the TCRA was likely due to the presence of more than one α-chain in some cells (37) as well as the variable efficiency of primers used in the TCRA sequencing. Thus, we decided to identify the productive TCR pairs with two different approaches. After identifying the most frequent TCRBs in the CD8+PD-1+ population, we identified the corresponding TCRA with single-cell RT-PCR on CD8+PD-1+ FACS-sorted TILs or CD8+ TILs expanded in vitro. Alternatively, we used the pairSEQ approach (21) on single-cell suspensions from unsorted fresh tumors. The efficiency of the single-cell RT-PCR was between 26% and 90%, depending on the sample. Using this method, we identified a median value of 29 (range, 9–43) unique TCRA–TCRB pairs in each of the CD8+ or CD8+PD-1+ samples. Using pairSEQ on unsorted fresh tumors we identified a median value of 217 (range, 11–883) unique pairs for each sample. A total of 93 (median value 6, range 0–21) TCRA–TCRB pairs from 12 metastatic melanoma patients were identified using both methods (congruent pairs in Table 2; method of identification for each pair in Supplementary Table S4). We generated expression vector constructs encoding the 83 of these pairs ranked within the top 10 CD8+PD-1+ clonotypes and linked them to murine constant chain sequences to improve stability and avoid mismatches with endogenous human TCRs (22), and then introduced them into fresh PBLs. The frequency of T cells that expressed the recombinant TCRs after either retroviral transduction or transfection with a Sleeping Beauty transposon construct ranged between 24.4% and 97.6% (Supplementary Table S5).

Table 2.

Unique TCR pairs identified

Sample IDNumber of unique TCR pairs identified by single-cell RT-PCRaNumber of unique TCR pairs identified by pairSEQbNumber of unique congruent TCR pairsNumber of reconstructed TCR pairs evaluated within the top 10 CD8+PD-1+ clonotypes
1913 29 136 
2650 30 21 
3678 32 11 
3713 21 829 21 
3759 34 133 15 
3784 15 883 
3903 14 156 10 
3922 351 
3926 33 737 
3977 29 21 
3992 20 278 
3998 43 349 19 
Sample IDNumber of unique TCR pairs identified by single-cell RT-PCRaNumber of unique TCR pairs identified by pairSEQbNumber of unique congruent TCR pairsNumber of reconstructed TCR pairs evaluated within the top 10 CD8+PD-1+ clonotypes
1913 29 136 
2650 30 21 
3678 32 11 
3713 21 829 21 
3759 34 133 15 
3784 15 883 
3903 14 156 10 
3922 351 
3926 33 737 
3977 29 21 
3992 20 278 
3998 43 349 19 

aSingle-cell PCR was performed on sorted CD8+PD-1+ TILs and for 1913, 2650, 3713, and 3784 on sorted CD8+-expanded TILs due to the limited availability of tumor samples.

bPairSEQ was performed on bulk TILs.

High-frequency CD8+PD-1+ clonotypes display tumor and mutation reactivity

We then evaluated the antitumor activity of T cells expressing those 83 TCRA–TCRB pairs (Table 2). The TCRs obtained from 10 of the 12 patients were evaluated for response to candidate neoepitopes identified by whole-exome sequencing of autologous tumor (TCRs from samples 2650 and 3977 were evaluated only against the TC line; Table 1). All the TCR pairs were also tested against autologous or HLA-matched antigen-presenting cells transfected with full-length RNA encoding the melanoma/melanocyte shared differentiation antigens MART-1, gp100, and tyrosinase (TYR) and the cancer-germline antigens NY-ESO-1, MAGEA3, and SSX2.

Evaluation of response against the corresponding autologous TC and/or autologous antigen-presenting cells that had either been pulsed with mutated tumor-specific neoantigen minimal epitopes or transfected with TMG (10–12) constructs provided evidence for tumor antigen reactivity in 11 of the 12 patients who were evaluated.

For example, for patient 3998 we initially evaluated the reactivity of some of the top eight most frequent TCR pairs based on the frequency of TCRB (3998-1, 3998-2, 3998-3A1, 3998-3A2, 3998-4, 3998-6, 3998-7, and 3998-8) against the autologous TC (Fig. 2A and C). Six of the TCR pairs tested (3998-1, 3998-2, 3998-4, 3998-6, 3998-7, and 3998-8) showed MHC-restricted recognition of the autologous tumor. The TCRB clonotype ranking third in frequency in the CD8+PD-1+ TILs was associated with two productive TCRA chains but neither of the two combinations (3998-3A1 and 3998-3A2) were tumor reactive. For tumor sample 3998, 345 nonsynonymous mutations were identified (Table 1). We next evaluated the reactivity of the 6 TCR pairs (3998-1, 3998-2, 3998-4, 3998-6, 3998-7, and 3998-8) against 115 mutated antigens encoded by 7 TMGs (Fig. 2B) and six shared melanoma/melanocyte differentiation antigens and cancer-germline antigens (MART-1, gp100, SSX2, TYR, NY-ESO-1, and MAGEA3; Supplementary Fig. S4A and S4B). The 115 mutated antigens were selected for screening from the 345 nonsynonymous mutations based on RNAseq data of their expression. Two TCR pairs (3998-7 and 3998-8) were reactive to TMG-1 (Fig. 2B). Further testing identified MAGEA6E168K as the specific mutation recognized within the antigens encoded by TMG-1 (Fig. 2D and Supplementary Fig. S4C). Reactivity against one shared antigen (NY-ESO-1) was found for TCR pair 3998-5 (Table 3 and Supplementary Fig. S4B).

Figure 2.

Tumor-specific target recognition assay for reconstructed TCR pairs for patient 3998. A, CD137 upregulation on CD8+mTCRB+ cells is shown after coculture with an autologous tumor cell line for 8 of the TCR pairs reconstructed within the top 10 CD8+PD-1+ TILs for this tumor sample: the 1st, 2nd, 3rd (in combination with 2 TCRAs), 4th, 6th, 7th, and 8th most frequent TCRBs. Values are reported as mean ± SEM; the assay was done in duplicate. B, CD137 upregulation on CD8+mTCRB+ cells is shown after coculture with autologous B cells transfected with tandem minigenes (TMG-1 to -7) encoding for 115 nonsynonymous mutations for 6 of the TCR pairs reconstructed (3998-1, 3998-2, 3998-4, 3998-6, 3998-7, and 3998-8). Values are reported as mean ± SEM; the assay was done in duplicate. C, CD137 upregulation is inhibited by pan–MHC I antibody. MHC I–restricted DMF5 TCR is reported as a positive control and MHC II–restricted TCR MAGE-A3 is reported as a negative control. *, greater than 50% inhibition. Values are reported as mean ± SEM; the assay was done in duplicate. D, murine TCRB expression and CD137 upregulation are shown for reconstructed TCR pair 3998-8 after coculture with unpulsed autologous B cells, autologous B cells pulsed with 1 μg/mL, 100 ng/mL, 10 ng/mL, 1 ng/mL, and 0.1 ng/mL of the mutated MAGEA6 peptide (KVDPIGHVY) and wild-type MAGEA6 peptide (EVDPIGHVY), respectively.

Figure 2.

Tumor-specific target recognition assay for reconstructed TCR pairs for patient 3998. A, CD137 upregulation on CD8+mTCRB+ cells is shown after coculture with an autologous tumor cell line for 8 of the TCR pairs reconstructed within the top 10 CD8+PD-1+ TILs for this tumor sample: the 1st, 2nd, 3rd (in combination with 2 TCRAs), 4th, 6th, 7th, and 8th most frequent TCRBs. Values are reported as mean ± SEM; the assay was done in duplicate. B, CD137 upregulation on CD8+mTCRB+ cells is shown after coculture with autologous B cells transfected with tandem minigenes (TMG-1 to -7) encoding for 115 nonsynonymous mutations for 6 of the TCR pairs reconstructed (3998-1, 3998-2, 3998-4, 3998-6, 3998-7, and 3998-8). Values are reported as mean ± SEM; the assay was done in duplicate. C, CD137 upregulation is inhibited by pan–MHC I antibody. MHC I–restricted DMF5 TCR is reported as a positive control and MHC II–restricted TCR MAGE-A3 is reported as a negative control. *, greater than 50% inhibition. Values are reported as mean ± SEM; the assay was done in duplicate. D, murine TCRB expression and CD137 upregulation are shown for reconstructed TCR pair 3998-8 after coculture with unpulsed autologous B cells, autologous B cells pulsed with 1 μg/mL, 100 ng/mL, 10 ng/mL, 1 ng/mL, and 0.1 ng/mL of the mutated MAGEA6 peptide (KVDPIGHVY) and wild-type MAGEA6 peptide (EVDPIGHVY), respectively.

Close modal
Table 3.

Antigen specificity identified

TCRB rank in TIL CD8+PD-1+191326503678371337593784390339223926397739923998
a TC N/A — TC — KIAA-isoM1L MART-1 TC TC — TC 
HLA-A11F33S N/A N/A HELZ2D614N TC N/A — — N/A N/A — TC 
TC — N/A SRPXP55L TC TC — — TC — — — 
HLA-A11F33S — — N/A TC — — N/A N/A N/A — TC 
— TC FBXO21S250Y TC TC — — N/A TC N/A — NY-ESO-1 
— — — WDR46T300I — N/A — N/A TC TC N/A TC 
— — FBXO21S250Y N/A N/A — — N/A TC N/A N/A MAGEA6E168K 
— N/A N/A TC N/A TC — — N/A — N/A MAGEA6 E168K 
— — N/A N/A N/A TC N/A — N/A — N/A N/A 
10 N/A N/A N/A WDR46T300I N/A TC N/A N/A N/A — N/A N/A 
TCRB rank in TIL CD8+PD-1+191326503678371337593784390339223926397739923998
a TC N/A — TC — KIAA-isoM1L MART-1 TC TC — TC 
HLA-A11F33S N/A N/A HELZ2D614N TC N/A — — N/A N/A — TC 
TC — N/A SRPXP55L TC TC — — TC — — — 
HLA-A11F33S — — N/A TC — — N/A N/A N/A — TC 
— TC FBXO21S250Y TC TC — — N/A TC N/A — NY-ESO-1 
— — — WDR46T300I — N/A — N/A TC TC N/A TC 
— — FBXO21S250Y N/A N/A — — N/A TC N/A N/A MAGEA6E168K 
— N/A N/A TC N/A TC — — N/A — N/A MAGEA6 E168K 
— — N/A N/A N/A TC N/A — N/A — N/A N/A 
10 N/A N/A N/A WDR46T300I N/A TC N/A N/A N/A — N/A N/A 

a—, No reactivity was found against autologous tumor cell lines (TC), tumor-specific mutations, melanoma/melanocyte, and cancer-germline antigens tested.

Figure 3 summarizes the findings of all 12 samples. Representative cocultures for all of the tumor samples are shown in Supplementary Figs. S4–S15. For example, in sample 1913 the TCRB rankings 2, 3, and 4 were specific for the autologous TC line, and the clonotypes ranking 2 and 4, as previously found (15, 28), also recognized a mutation in the HLA-11 gene (Table 3). Moreover, the most frequent TCR clonotype was found to be tumor reactive for seven samples (2650, 3759, 3903, 3922, 3926, 3977, and 3998). For all but patient 3992, up to five tumor-reactive TCRs were found among the five most frequently expressed TCRs in the CD8+PD-1+ TILs. Reactivity against autologous neoantigens was found in 5 of the 10 patients whose TCRs were screened against putative autologous mutations. In summary, we found that 36 TCR pairs were reactive against autologous tumor and 11 were directed against mutated tumor-specific neoantigens. This finding indicates that it is possible to identify tumor-reactive TCR pairs in the majority of melanoma samples simply based on their frequency in the CD8+PD-1+ TIL compartment.

Figure 3.

Summary of tumor and mutation reactivity for reconstructed TCR pairs. For every sample analyzed, the graph represents the TCRB frequency of the top 10 CD8+PD-1+ clonotypes and color codes their reactivity against autologous TC lines, shared melanoma/melanocyte and cancer-germline antigens, and tumor-specific mutations. All patients, except 3678, had a corresponding autologous TC line used for testing the TCR pairs. In 11 of 12 patients, up to 5 tumor-reactive TCRs were found in the 5 most frequently expressed TCRs, and this included recognition of mutated neoantigens in 5 of the patients. In 2 patients, reactivity against MART-1 (3922-1) and NY-ESO-1 (3998-5) was also found. The most frequent TCR clonotype was found to be tumor reactive for 7 patients.

Figure 3.

Summary of tumor and mutation reactivity for reconstructed TCR pairs. For every sample analyzed, the graph represents the TCRB frequency of the top 10 CD8+PD-1+ clonotypes and color codes their reactivity against autologous TC lines, shared melanoma/melanocyte and cancer-germline antigens, and tumor-specific mutations. All patients, except 3678, had a corresponding autologous TC line used for testing the TCR pairs. In 11 of 12 patients, up to 5 tumor-reactive TCRs were found in the 5 most frequently expressed TCRs, and this included recognition of mutated neoantigens in 5 of the patients. In 2 patients, reactivity against MART-1 (3922-1) and NY-ESO-1 (3998-5) was also found. The most frequent TCR clonotype was found to be tumor reactive for 7 patients.

Close modal

We have shown that adoptive cell therapy with TILs that appeared to predominantly recognize patient-specific tumor neoantigens (11, 12) or T cells genetically engineered to express TCRs targeting cancer-germline antigens (8) can mediate complete response in patients with metastatic melanoma (10, 12) and in a patient with metastatic cholangiocarcinoma (29). Those studies used a labor-intensive screening approach with tandem minigenes or long peptides representing all known mutations and could identify T cells with reactivity against mutated neoantigens. Tumor-reactive TCRs expressed by mutation-reactive T cells can be isolated, cloned into expression vectors, and can potentially be transferred into autologous cells with high proliferative capacity (13, 38) for use in cell transfer therapy. Here, we demonstrate the identification of tumor and mutation-reactive TCRs from fresh melanoma samples, based on PD-1 expression and on TCRB frequencies as a guide to tumor reactivity.

The development of next-generation TCRB sequencing has allowed a study of the total TCR repertoire in different T-cell compartments in healthy individuals (32) as well as in tumor samples from colorectal (33) and ovarian carcinomas (34). In the present study, we have analyzed the TCR diversity in different subsets of TILs from freshly resected human metastatic melanomas and attempted to determine whether the rank frequency of TCRs was related to their ability to recognize the autologous cancer. This strategy can be used to prospectively identify tumor-reactive TCRs without the prior need to know their antigen specificity.

The major obstacle encountered in the evaluation of the antitumor reactivity of individual TCRs in the fresh tumor prior to in vitro expansion was the technical difficulty in pairing each high-frequency TCRB chain with the correct TCRA. In the present study, we utilized two independent approaches to identify the TCRA–TCRB pairs, single-cell RT-PCR on a specific subset (CD8+PD-1+ TILs or CD8+-expanded TILs) and pairSEQ on unsorted tumors. The single-cell RT-PCR allowed us to directly extract cDNA from single cells and sequence the TCR gene after several rounds of PCR with specific primers using the Sanger method (39), followed by TOPO TA cloning in the event the T cell expressed two different TCRA chains at the same time. This technique, based on a 96-well platform, has also been successful when applied to MHC-multimer sorting on antigen-specific cells (40, 41), although MHC-multimer sorting is feasible only when the HLA restriction element and the minimal epitope of interest are known. In this study, we used a different approach that does not require the knowledge of antigen specificity by sorting specific TIL subsets. Other higher- throughput approaches have been proposed, such as emulsion PCR (42) and “TCR gene capture” that utilizes an RNA-bait library to specifically target the genomic sequence encoding TCR genes (43). Recently, pairSEQ (21), a new high-throughput technology for pairing TCRA and TCRB sequences, has become available, and we tested its feasibility in fresh, unsorted melanoma samples. PairSEQ uses a statistical model for pairing TCRA and TCRB chains and, because it is based on next-generation sequencing, no extra steps are required (such as TOPO TA cloning) to identify two different TCRA genes expressed by the same cell. Combining the single-cell RT-PCR on specific TIL subsets and pairSEQ on unsorted tumors, we successfully identified the majority, but not all, of the top 10 most frequent TIL CD8+PD-1+ TCR pairs in 12 patients and tested them against autologous tumor tissue culture lines and autologous antigen-presenting cells expressing tandem minigenes encoding shared cancer antigens, mutated tumor neoantigens, and/or pulsed with the corresponding mutated peptides (11, 29). Multiple reactive pairs could be identified in the top-ranking TCRs in 11 of 12 patients with metastatic melanoma, including 11 TCRs specific for mutated neoepitopes (Fig. 3). We included the six most common melanoma/melanocyte and cancer-germline antigens commonly recognized by patient TILs in our screening. With this limited screening panel, we identified two TCR pairs that were specific for MART-1 (3922-1) and NY-ESO-1 (3998-5). For the other 25 TCR pairs, we were able to demonstrate MHC-restricted reactivity against the autologous tumor cell line, although their specificity remains undefined.

Nonreactive pairs need to be considered with caution. Our approach is based on PCR; therefore, it is subjected to potential errors that could have altered the original sequence of the TCRA–TCRB pairs. Incorrect pairing is also a possible explanation.

Despite these limitations, we identified tumor-reactive TCRs based on their TCRB frequency in the TIL CD8+PD-1+ population. The major advantage of this approach is the rapidity in finding functional TCRs without the need for further screening. However, once the high-frequency TCRs are identified, they can be tested in vitro in an overnight assay versus the fresh tumor suspension as well as normal autologous PBL to further identify their specificity. This finding opens the possibility for a highly personalized cell transfer cancer therapy in which patients can be treated with autologous, genetically engineered T cells with high proliferative potential. This new approach can potentially be applied to malignancies other than melanoma and is currently under study.

B. Howie has an ownership interest in patents pertaining to pairSEQ Technology. H. Robins is CSO at Adaptive Biotechnologies for which he has ownership interest, including patents. No potential conflicts of interest were disclosed by the other authors.

Conception and design: A. Pasetto, A. Gros, S.A. Rosenberg

Development of methodology: A. Pasetto, A. Gros, D.C. Deniger, R. Matus-Nicodemos, D.C. Douek

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): A. Pasetto, P.F. Robbins, T.D. Prickett, M.R. Parkhurst, J. Gartner, K. Trebska-McGowan, J.S. Crystal

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): A. Pasetto, A. Gros, P.F. Robbins, D.C. Deniger, T.D. Prickett, B. Howie, H. Robins, J. Gartner, J.S. Crystal, S.A. Rosenberg

Writing, review, and/or revision of the manuscript: A. Pasetto, A. Gros, P.F. Robbins, D.C. Deniger, T.D. Prickett, D.C. Douek, H. Robins, J. Gartner, K. Trebska-McGowan, J.S. Crystal, S.A. Rosenberg

Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): J.S. Crystal, S.A. Rosenberg

Study supervision: A. Gros

Other (performed multiple cell assays and analyzed their results): K. Trebska-McGowan

We thank Sanja Stevanović, Eric Tran, William Lu, Mojgan Ahmadzadeh, and Cyril Cohen for helpful discussions; Yang-Li, Mona El-Gamil, and Lien Ngo for technical advice and support; and Arnold Mixon and Shown Farid for flow cytometry technical support and sorting. We also thank the Adelson Medical Research Foundation for their generous support for this study.

This research was supported by the Intramural Research Program of the NIH at the National Cancer Institute.

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1.
Clemente
CG
,
Mihm
MC
 Jr
,
Bufalino
R
,
Zurrida
S
,
Collini
P
,
Cascinelli
N
. 
Prognostic value of tumor infiltrating lymphocytes in the vertical growth phase of primary cutaneous melanoma
.
Cancer
1996
;
77
:
1303
10
.
2.
Sato
E
,
Olson
SH
,
Ahn
J
,
Bundy
B
,
Nishikawa
H
,
Qian
F
, 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
.
3.
Galon
J
,
Costes
A
,
Sanchez-Cabo
F
,
Kirilovsky
A
,
Mlecnik
B
,
Lagorce-Pages
C
, et al
Type, density, and location of immune cells within human colorectal tumors predict clinical outcome
.
Science
2006
;
313
:
1960
4
.
4.
Loi
S
. 
Tumor-infiltrating lymphocytes, breast cancer subtypes and therapeutic efficacy
.
Oncoimmunology
2013
;
2
:
e24720
.
5.
Kawakami
Y
,
Eliyahu
S
,
Delgado
CH
,
Robbins
PF
,
Sakaguchi
K
,
Appella
E
, et al
Identification of a human melanoma antigen recognized by tumor-infiltrating lymphocytes associated with in vivo tumor rejection
.
Proc Natl Acad Sci U S A
1994
;
91
:
6458
62
.
6.
Kawakami
Y
,
Eliyahu
S
,
Delgado
CH
,
Robbins
PF
,
Rivoltini
L
,
Topalian
SL
, et al
Cloning of the gene coding for a shared human melanoma antigen recognized by autologous T cells infiltrating into tumor
.
Proc Natl Acad Sci U S A
1994
;
91
:
3515
9
.
7.
Andersen
RS
,
Thrue
CA
,
Junker
N
,
Lyngaa
R
,
Donia
M
,
Ellebaek
E
, et al
Dissection of T-cell antigen specificity in human melanoma
.
Cancer Res
2012
;
72
:
1642
50
.
8.
Robbins
PF
,
Morgan
RA
,
Feldman
SA
,
Yang
JC
,
Sherry
RM
,
Dudley
ME
, et al
Tumor regression in patients with metastatic synovial cell sarcoma and melanoma using genetically engineered lymphocytes reactive with NY-ESO-1
.
J Clin Oncol
2011
;
29
:
917
24
.
9.
Abate-Daga
D
,
Speiser
DE
,
Chinnasamy
N
,
Zheng
Z
,
Xu
H
,
Feldman
SA
, et al
Development of a T cell receptor targeting an HLA-A*0201 restricted epitope from the cancer-testis antigen SSX2 for adoptive immunotherapy of cancer
.
PLoS One
2014
;
9
:
e93321
.
10.
Robbins
PF
,
Lu
YC
,
El-Gamil
M
,
Li
YF
,
Gross
C
,
Gartner
J
, et al
Mining exomic sequencing data to identify mutated antigens recognized by adoptively transferred tumor-reactive T cells
.
Nat Med
2013
;
19
:
747
52
.
11.
Lu
YC
,
Yao
X
,
Crystal
JS
,
Li
YF
,
El-Gamil
M
,
Gross
C
, et al
Efficient identification of mutated cancer antigens recognized by T cells associated with durable tumor regressions
.
Clin Cancer Res
2014
;
20
:
3401
10
.
12.
Lu
YC
,
Yao
X
,
Li
YF
,
El-Gamil
M
,
Dudley
ME
,
Yang
JC
, et al
Mutated PPP1R3B is recognized by T cells used to treat a melanoma patient who experienced a durable complete tumor regression
.
J Immunol
2013
;
190
:
6034
42
.
13.
Rosenberg
SA
,
Yang
JC
,
Sherry
RM
,
Kammula
US
,
Hughes
MS
,
Phan
GQ
, et al
Durable complete responses in heavily pretreated patients with metastatic melanoma using T-cell transfer immunotherapy
.
Clin Cancer Res
2011
;
17
:
4550
7
.
14.
Robbins
PF
,
Dudley
ME
,
Wunderlich
J
,
El-Gamil
M
,
Li
YF
,
Zhou
J
, et al
Cutting edge: persistence of transferred lymphocyte clonotypes correlates with cancer regression in patients receiving cell transfer therapy
.
J Immunol
2004
;
173
:
7125
30
.
15.
Gros
A
,
Robbins
PF
,
Yao
X
,
Li
YF
,
Turcotte
S
,
Tran
E
, et al
PD-1 identifies the patient-specific CD8(+) tumor-reactive repertoire infiltrating human tumors
.
J Clin Invest
2014
;
124
:
2246
59
.
16.
Inozume
T
,
Hanada
K
,
Wang
QJ
,
Ahmadzadeh
M
,
Wunderlich
JR
,
Rosenberg
SA
, et al
Selection of CD8+PD-1+ lymphocytes in fresh human melanomas enriches for tumor-reactive T cells
.
J Immunother
2010
;
33
:
956
64
.
17.
Nikolich-Zugich
J
,
Slifka
MK
,
Messaoudi
I
. 
The many important facets of T-cell repertoire diversity
.
Nat Rev Immunol
2004
;
4
:
123
32
.
18.
Gros
A
,
Turcotte
S
,
Wunderlich
JR
,
Ahmadzadeh
M
,
Dudley
ME
,
Rosenberg
SA
. 
Myeloid cells obtained from the blood but not from the tumor can suppress T-cell proliferation in patients with melanoma
.
Clin Cancer Res
2012
;
18
:
5212
23
.
19.
Jones
S
,
Wang
TL
,
Shih Ie
M
,
Mao
TL
,
Nakayama
K
,
Roden
R
, et al
Frequent mutations of chromatin remodeling gene ARID1A in ovarian clear cell carcinoma
.
Science
2010
;
330
:
228
31
.
20.
Prickett
TD
,
Crystal
JS
,
Cohen
CJ
,
Pasetto
A
,
Parkhurst
MR
,
Gartner
JJ
, et al
Durable complete response from metastatic melanoma after transfer of autologous T cells recognizing 10 mutated tumor antigens
.
Cancer Immunol Res
2016
Jun 16. [Epub ahead of print]
.
21.
Howie
B
,
Sherwood
AM
,
Berkebile
AD
,
Berka
J
,
Emerson
RO
,
Williamson
DW
, et al
High-throughput pairing of T cell receptor alpha and beta sequences
.
Sci Transl Med
2015
;
7
:
301ra131
.
22.
Cohen
CJ
,
Zhao
Y
,
Zheng
Z
,
Rosenberg
SA
,
Morgan
RA
. 
Enhanced antitumor activity of murine-human hybrid T-cell receptor (TCR) in human lymphocytes is associated with improved pairing and TCR/CD3 stability
.
Cancer Res
2006
;
66
:
8878
86
.
23.
Szymczak
AL
,
Workman
CJ
,
Wang
Y
,
Vignali
KM
,
Dilioglou
S
,
Vanin
EF
, et al
Correction of multi-gene deficiency in vivo using a single ‘self-cleaving’ 2A peptide-based retroviral vector
.
Nat Biotechnol
2004
;
22
:
589
94
.
24.
Jones
S
,
Peng
PD
,
Yang
S
,
Hsu
C
,
Cohen
CJ
,
Zhao
Y
, et al
Lentiviral vector design for optimal T cell receptor gene expression in the transduction of peripheral blood lymphocytes and tumor-infiltrating lymphocytes
.
Hum Gene Ther
2009
;
20
:
630
40
.
25.
Yang
S
,
Rosenberg
SA
,
Morgan
RA
. 
Clinical-scale lentiviral vector transduction of PBL for TCR gene therapy and potential for expression in less-differentiated cells
.
J Immunother
2008
;
31
:
830
9
.
26.
Deniger
DC
,
Yu
J
,
Huls
MH
,
Figliola
MJ
,
Mi
T
,
Maiti
SN
, et al
Sleeping beauty transposition of chimeric antigen receptors targeting receptor tyrosine kinase-like orphan receptor-1 (ROR1) into diverse memory T-cell populations
.
PLoS One
2015
;
10
:
e0128151
.
27.
Deniger
DC
,
Pasetto
A
,
Tran
E
,
Parkhurst
MR
,
Cohen
CJ
,
Robbins
PF
, et al
Stable, nonviral expression of mutated tumor neoantigen-specific T-cell receptors using the sleeping beauty transposon/transposase system
.
Mol Ther
2016
;
24
:
1078
89
.
28.
Huang
J
,
El-Gamil
M
,
Dudley
ME
,
Li
YF
,
Rosenberg
SA
,
Robbins
PF
. 
T cells associated with tumor regression recognize frameshifted products of the CDKN2A tumor suppressor gene locus and a mutated HLA class I gene product
.
J Immunol
2004
;
172
:
6057
64
.
29.
Tran
E
,
Turcotte
S
,
Gros
A
,
Robbins
PF
,
Lu
YC
,
Dudley
ME
, et al
Cancer immunotherapy based on mutation-specific CD4+ T cells in a patient with epithelial cancer
.
Science
2014
;
344
:
641
5
.
30.
Wolfl
M
,
Kuball
J
,
Ho
WY
,
Nguyen
H
,
Manley
TJ
,
Bleakley
M
, et al
Activation-induced expression of CD137 permits detection, isolation, and expansion of the full repertoire of CD8+ T cells responding to antigen without requiring knowledge of epitope specificities
.
Blood
2007
;
110
:
201
10
.
31.
Ahmadzadeh
M
,
Johnson
LA
,
Heemskerk
B
,
Wunderlich
JR
,
Dudley
ME
,
White
DE
, et al
Tumor antigen-specific CD8 T cells infiltrating the tumor express high levels of PD-1 and are functionally impaired
.
Blood
2009
;
114
:
1537
44
.
32.
Robins
HS
,
Campregher
PV
,
Srivastava
SK
,
Wacher
A
,
Turtle
CJ
,
Kahsai
O
, et al
Comprehensive assessment of T-cell receptor beta-chain diversity in alphabeta T cells
.
Blood
2009
;
114
:
4099
107
.
33.
Sherwood
AM
,
Emerson
RO
,
Scherer
D
,
Habermann
N
,
Buck
K
,
Staffa
J
, et al
Tumor-infiltrating lymphocytes in colorectal tumors display a diversity of T cell receptor sequences that differ from the T cells in adjacent mucosal tissue
.
Cancer Immunol Immunother
2013
;
62
:
1453
61
.
34.
Emerson
RO
,
Sherwood
AM
,
Rieder
MJ
,
Guenthoer
J
,
Williamson
DW
,
Carlson
CS
, et al
High-throughput sequencing of T-cell receptors reveals a homogeneous repertoire of tumour-infiltrating lymphocytes in ovarian cancer
.
J Pathol
2013
;
231
:
433
40
.
35.
Srivastava
SK
,
Robins
HS
. 
Palindromic nucleotide analysis in human T cell receptor rearrangements
.
PLoS One
2012
;
7
:
e52250
.
36.
Stewart
JJ
,
Lee
CY
,
Ibrahim
S
,
Watts
P
,
Shlomchik
M
,
Weigert
M
, et al
A Shannon entropy analysis of immunoglobulin and T cell receptor
.
Mol Immunol
1997
;
34
:
1067
82
.
37.
Padovan
E
,
Casorati
G
,
Dellabona
P
,
Meyer
S
,
Brockhaus
M
,
Lanzavecchia
A
. 
Expression of two T cell receptor alpha chains: dual receptor T cells
.
Science
1993
;
262
:
422
4
.
38.
Gattinoni
L
,
Lugli
E
,
Ji
Y
,
Pos
Z
,
Paulos
CM
,
Quigley
MF
, et al
A human memory T cell subset with stem cell-like properties
.
Nat Med
2011
;
17
:
1290
7
.
39.
Wang
GC
,
Dash
P
,
McCullers
JA
,
Doherty
PC
,
Thomas
PG
. 
T cell receptor alphabeta diversity inversely correlates with pathogen-specific antibody levels in human cytomegalovirus infection
.
Sci Transl Med
2012
;
4
:
128ra42
.
40.
Dossinger
G
,
Bunse
M
,
Bet
J
,
Albrecht
J
,
Paszkiewicz
PJ
,
Weissbrich
B
, et al
MHC multimer-guided and cell culture-independent isolation of functional T cell receptors from single cells facilitates TCR identification for immunotherapy
.
PLoS One
2013
;
8
:
e61384
.
41.
Kobayashi
E
,
Mizukoshi
E
,
Kishi
H
,
Ozawa
T
,
Hamana
H
,
Nagai
T
, et al
A new cloning and expression system yields and validates TCRs from blood lymphocytes of patients with cancer within 10 days
.
Nat Med
2013
;
19
:
1542
6
.
42.
Turchaninova
MA
,
Britanova
OV
,
Bolotin
DA
,
Shugay
M
,
Putintseva
EV
,
Staroverov
DB
, et al
Pairing of T-cell receptor chains via emulsion PCR
.
Eur J Immunol
2013
;
43
:
2507
15
.
43.
Linnemann
C
,
Heemskerk
B
,
Kvistborg
P
,
Kluin
RJ
,
Bolotin
DA
,
Chen
X
, et al
High-throughput identification of antigen-specific TCRs by TCR gene capture
.
Nat Med
2013
;
19
:
1534
41
.