
==== Front
eLife
Elife
eLife
eLife
2050-084X
eLife Sciences Publications, Ltd

39312285
91705
10.7554/eLife.91705
version of record
Research Article
Cancer Biology
Genetics and Genomics
Patient-derived xenografts and single-cell sequencing identifies three subtypes of tumor-reactive lymphocytes in uveal melanoma metastases
Karlsson Joakim W 12†
Sah Vasu R 2†
Olofsson Bagge Roger https://orcid.org/0000-0001-5795-0355
234
Kuznetsova Irina 1
Iqba Munir 5
Alsen Samuel 2
Stenqvist Sofia 2
Saxena Alka https://orcid.org/0000-0001-5683-0618
5
Ny Lars 26
Nilsson Lisa M 12
Nilsson Jonas A https://orcid.org/0000-0003-0346-6837
jonas.a.nilsson@surgery.gu.se
12
1 https://ror.org/02xz7d723 Harry Perkins Institute of Medical Research and University of Western Australia Perth Australia
2 https://ror.org/01tm6cn81 Sahlgrenska Center for Cancer Research, Institute of Clinical Sciences, Sahlgrenska Academy, University of Gothenburg Gothenburg Sweden
3 https://ror.org/04vgqjj36 Department of Surgery, Sahlgrenska University Hospital Gothenburg Sweden
4 https://ror.org/01tm6cn81 Wallenberg Centre for Molecular and Translational Medicine, University of Gothenburg Gothenburg Sweden
5 https://ror.org/01dbmzx78 Genomics WA, Telethon Kids Institute, Harry Perkins Institute of Medical Research and University of Western Australia Nedlands Australia
6 https://ror.org/04vgqjj36 Department of Oncology, Sahlgrenska University Hospital Gothenburg Sweden
Yang Yongliang Reviewing Editor https://ror.org/023hj5876 Dalian University of Technology China

Liu Caigang Senior Editor https://ror.org/032d4f246 Shengjing Hospital of China Medical University China

† These authors contributed equally to this work.

23 9 2024
2024
12 RP9170522 8 2023
This manuscript was published as a preprint.18 5 2023

This manuscript was published as a reviewed preprint.16 11 2023

The reviewed preprint was revised.11 7 2024

© 2023, Karlsson, Sah et al
2023
Karlsson, Sah et al
https://creativecommons.org/licenses/by/4.0/ This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Uveal melanoma (UM) is a rare melanoma originating in the eye’s uvea, with 50% of patients experiencing metastasis predominantly in the liver. In contrast to cutaneous melanoma, there is only a limited effectiveness of combined immune checkpoint therapies, and half of patients with uveal melanoma metastases succumb to disease within 2 years. This study aimed to provide a path toward enhancing immunotherapy efficacy by identifying and functionally validating tumor-reactive T cells in liver metastases of patients with UM. We employed single-cell RNA-seq of biopsies and tumor-infiltrating lymphocytes (TILs) to identify potential tumor-reactive T cells. Patient-derived xenograft (PDX) models of UM metastases were created from patients, and tumor sphere cultures were generated from these models for co-culture with autologous or MART1-specific HLA-matched allogenic TILs. Activated T cells were subjected to TCR-seq, and the TCRs were matched to those found in single-cell sequencing data from biopsies, expanded TILs, and in livers or spleens of PDX models injected with TILs. Our findings revealed that tumor-reactive T cells resided not only among activated and exhausted subsets of T cells, but also in a subset of cytotoxic effector cells. In conclusion, combining single-cell sequencing and functional analysis provides valuable insights into which T cells in UM may be useful for cell therapy amplification and marker selection.

patient-derived xenograft
tumor-infiltrating lymphocytes
cellular immunotherapy
Research organism

Human
Mouse
http://dx.doi.org/10.13039/501100002794 Cancerfonden Olofsson Bagge Roger Ny Lars Nilsson Jonas A http://dx.doi.org/10.13039/100007436 Erling-Persson Family Foundation Nilsson Jonas A http://dx.doi.org/10.13039/501100004359 Vetenskapsrådet Olofsson Bagge Roger Nilsson Jonas A http://dx.doi.org/10.13039/100016669 Sjöbergstiftelsen Nilsson Jonas A Harry Perkins Institute of Medical Research start-up grants Nilsson Jonas A Västragötaland Region ALF grant Nilsson Jonas A Olofsson Bagge Roger Ny Lars http://dx.doi.org/10.13039/501100004063 Knut and Alice Wallenberg Foundation Nilsson Jonas A Olofsson Bagge Roger BioCARE strategic grants Nilsson Jonas A Lion’s Cancerfond Väst Alsen Samuel Ny Lars http://dx.doi.org/10.13039/501100005009 Assar Gabrielssons Foundation Sah Vasu R Gustaf V Jubileumsklinikens forskningsfond Ny Lars Wilhelm och Martina Lundgrens Vetenskapsfond Karlsson Joakim W The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.Author impact statementSingle-cell sequencing and functional analysis identify T cells that can be useful for marker selection and cell therapy in uveal melanoma, a disease largely unresponsive to conventional immune checkpoint therapies.
publishing-routeprc
==== Body
pmcIntroduction

Uveal melanoma (UM) is a rare form of melanoma and the most common primary malignancy of the eye (Jager et al., 2020). It develops in the uvea of the eye, most often in the choroid and ciliary body, and more infrequently in the iris. Primary UM can be cured by brachytherapy or enucleation of the eye, but 50% of patients will develop metastasis (Buder et al., 2013). The most common route for metastatic spread is the liver for reasons that are poorly understood. Patients with metastatic disease have a median survival of less than a year (Khoja et al., 2019) however, ongoing clinical studies suggest that some treatments could prolong survival (Carvajal et al., 2017). Monotherapy immune checkpoint inhibitors (ICI) are markedly less effective in patients with UM (Johnson et al., 2019; Algazi et al., 2016) than in those with metastatic cutaneous melanomas. However, combination treatments with PD-1/CTLA-4 inhibitors (Piulats et al., 2021; Pelster et al., 2021; Najjar et al., 2020) or PD-1/HDAC inhibitors (Jespersen et al., 2019; Ny et al., 2021) have demonstrated longer overall survival than historic benchmark data (Khoja et al., 2019). Tebentafusp, a T cell engager suitable for patients with the HLA-A2 genotype (Damato et al., 2019; Middleton et al., 2020), is the first therapy to show a prolonged overall survival in a phase 3 randomized trial, increasing median survival from 16.0 months to 21.7 months (Nathan et al., 2021). Another development is locoregional treatments, where a recent phase 3 trial demonstrated that isolated hepatic perfusion (IHP) triples hepatic progression-free survival (Olofsson Bagge et al., 2023) compared with the best alternative care and historic benchmark data. Although retrospective data suggest an overall survival benefit with liver-directed therapy (Ben-Shabat et al., 2016), in the SCANDIUM phase 3 trial, overall survival of patients treated with IHP was not statistically significantly better that for patients treated with best alternative care (Olofsson Bagge et al., 2024). Moreover, invariably, all patients progressed on both tebentafusp and IHP; therefore, more research is needed.

Adoptive cell therapy (ACT) with TILs or CAR-T cells has not been extensively studied in UM. Pilot data from a clinical trial (Chandran et al., 2017) showed that TILs can cause responses in patients, and we have previously showed that HER2 CAR-T cells can eradicate UM in xenografts (Forsberg et al., 2019). There is, however, a lack of robust ex vivo screening models and very few patient-derived xenograft (PDX) mouse models from metastases (Carita et al., 2015) to use in studies to improve ACT. Part of the lack of success for ACT may be attributable to a lack of tumor-reactive lymphocytes, possibly owing to lower mutation burden in UM compared to cutaneous melanoma (van der Kooij et al., 2019). However, both the clinical trial and previous analyses of TILs in UM liver metastases have suggested that at least some TILs are tumor reactive (Chandran et al., 2017; Durante et al., 2020; Karlsson et al., 2020; Rothermel et al., 2016).

Defining tumor-reactive T lymphocytes (TRLs) among TILs would enable the identification of candidate biomarkers for selective expansion or TCR transgenics in cell therapy experiments in mouse models and cell therapy trials. To decipher how immunotherapy can be made more effective, the aim of this study was to identify and functionally validate TRLs in metastases of patients with UM. To this end, we used paired single-cell RNA (scRNA)- and TCR-seq, as well as functional experiments using PDX models and autologous TIL cultures.

Results

We previously collected metastases of UM from patients in routine clinical care and two clinical trials. These were the SCANDIUM phase 3 randomized trial comparing IHP to best alternative care (Olofsson et al., 2014) and the PEMDAC phase 2 single-arm trial where patients received a combination of the PD-1 inhibitor pembrolizumab and the HDAC inhibitor entinostat (Ny et al., 2021). Owing to the limited amount available of these biopsies, we prioritized one part of the biopsy for DNA/RNA preparation to sequence (Karlsson et al., 2020), one part for transplantation into immunocompromised mice to generate PDX mouse models, one part for generation of TILs, one part for cryopreservation of finely minced tumor (flow cytometry and single-cell sequencing), and one part for formalin-fixed paraffin-embedded (FFPE) blocks (Figure 1a).

Figure 1. Single-cell RNA sequencing of uveal melanoma metastases.

(a) Schematic showing the study workflow. Tumors extracted from patient liver metastases were subjected to single-cell RNA-seq (scRNA-seq), formalin fixation for immunohistochemical (IHC) analysis, and patient-derived xenograft (PDX) generation. (b) Uniform Manifold Approximation and Projection (UMAP) dimensionality-reduced expression profiles of single cells in biopsies from 14 patients with uveal melanoma (UM). scRNA-seq data from all samples were integrated using FastMNN, Louvain clustering performed, and clusters annotated using cell-type marker genes compiled from the literature. Some cell types were annotated after additional subclustering (Figure 1—figure supplement 1a–f). (c) Cluster proportions within each sample, for all cell types. Subcutaneous biopsies are indicated with asterisk, remaining samples being liver biopsies. Panel a generated with BioRender.com.

Figure 1—figure supplement 1. Identification of cell types in single-cell RNA sequencing data.

(a, b) Subclustering of cells that initially grouped together with endothelial cells (a). Shown per tissue of origin (b). (c, d) Subclustering of monocyte-like cells. Shown per tissue of origin (d). (e, f) Marker genes used to label clusters in (a) and (c), respectively. (g) Sample contribution to each cluster from Figure 1b. (h) Copy number profiles for uveal melanoma (UM) cells identified in each sample, as determined by inferCNV (Tickle et al., 2019) analysis. Blue represents copy number loss and red copy number gain. (i) Sample contribution to each cluster from Figure 2a. (j, k) Cell contribution to each CD8+ T cell cluster from biopsies of liver and subcutaneous tumors, respectively (Haghverdi et al., 2018).

Single-cell sequencing defines an atlas of T cell states in UM metastases

We obtained biopsies suitable for sequencing from 14 patients, two of which were subcutaneous and 12 were liver metastases. To identify and profile the T cells present in these samples, we performed paired scRNA- and TCR-seq, integrated the samples, and annotated the cell types. This has resulted in the largest single-cell atlas of metastatic UM to date, complementing prior efforts, including our own study of TILs from UM metastases (Durante et al., 2020; Karlsson et al., 2020; Lin et al., 2021; Pandiani et al., 2021).

Uniform Manifold Approximation and Projection (UMAP) dimensionality-reduced expression profiles revealed clusters containing melanoma, CD8+ or CD4+ T cells, NK cells, monocytes, and small clusters of B cells, endothelial cells, and other minor cell populations (Figure 1b and c, Figure 1—figure supplement 1a–f). Clusters were generally well mixed with respect to contributions from different samples, suggesting few patient-specific or batch effects after data integration (Figure 1—figure supplement 1g). The identity of cells labeled as UM was further confirmed by inferring copy number changes. This revealed characteristic aberrations seen in metastatic UM, including monosomy of chromosome 3 and gain of 8q, as well as subclonal variation within patients (Figure 1—figure supplement 1h).

Further subclustering of CD8+ T cells showed the presence of 12 different groups (Figure 2a and b and Figure 1—figure supplement 1i–k), all characterized by the differential expression of immunological marker genes (Figure 2c and d and Supplementary file 1). The clusters can be divided into two axes. One of these is characterized by expression of genes associated with late activation and exhaustion, such as PDCD1, HAVCR2, TNFRSF9, and HLA-DRA. Additional expression of ITGAE (CD103) suggests that these clusters contain tissue-resident memory cells (Ganesan et al., 2017). Multiplex immunofluorescence showed that PD-1- and ICOS-expressing T cells were in close proximity to the tumor cells (Figure 2e). The other axis expressed genes associated with naïve/memory-like or early activated phenotypes, such as IL7R, TCF7, CCR7, and NR4A1 (Figure 2d).

Figure 2. CD8 T cells in uveal melanoma metastases.

(a) Uniform Manifold Approximation and Projection (UMAP) dimensionality reduction and more detailed re-annotation of the CD8+ T cell subsets from Figure 1b, using the marker genes shown in (c). (b) CD8+ T cell cluster proportions within each sample. (c) Expression levels of markers used for CD8+ T cell subset annotation. (d) Expression of marker genes signifying early (NR4A1) activation, naïve/memory-like phenotypes (TCF7, IL7R), late activation (HLA-DRA, PDCD1, ICOS), and progression toward exhaustion/dysfunction (LAG3, TIGIT, HAVCR2). Additional statistically identified genes differing between clusters can be found in Supplementary file 1. (e) Multiplex immunohistochemical (IHC) staining for T cell (CD3, ICOS, PD-1) and cancer marker genes (β-catenin) in a patient biopsy (UM24).

Individual clusters within these axes mostly correspond to intermediate phenotypes that progress from activation to exhaustion. However, a few seem to represent distinct states. Two of these were within the IL7+ axis. Of these, one expressed KLRB1 and NCR3 and the other expressed NK-associated markers, such as FCGR3A, FCRL6, and NKG7 (Figure 2c and d). Although both clusters expressed genes related to cytotoxicity, only the latter expressed GZMB and GZMH. Other markers associated with this group are CX3CR1 and FGFBP2. Previous studies labeled these cells as cytotoxic NK-like T cells (Li et al., 2019; Huuhtanen, 2023).

In summary, this map reveals the cell types present in liver and subcutaneous metastases of UM and divides infiltrating T cells into two major phenotypic groups, characterized by the mutually exclusive expression of HAVCR2 and IL7R, but also reveals distinct cytotoxic subsets expressing NK-associated markers.

TILs from UM metastases can kill their cognate tumors in PDX models

Given the presence of activated T cells in tumors, we investigated whether these could eradicate metastatic UM in ACT experiments using PDX models. These models are useful personalized models for studying tumor biology and more recently, tumor immunology (Patton et al., 2021). However, the biopsies we received were generally very small, and the initial take rate using the same single-cell dispersion method used for cutaneous melanoma metastases (Patton et al., 2021) was unsatisfactory. Therefore, we transplanted small tumor pieces or used splenic injections, which led to successful PDX generation for some patients. Ten of these are presented in this study. All the PDX tumors exhibited confirmatory markers of melanoma, such as SOX10, MELAN-A (MART1), HMB45, or S100B. Sanger sequencing further revealed mutually exclusive mutations in the UM oncogenes GNAQ, GNA11, CYSLTR2, and PLCB4, as found in biopsies (Karlsson et al., 2020; Figure 3a and Figure 3—figure supplement 1). The UM tumors grew at a slower rate compared to our experience with cutaneous melanoma, which has a take rate of 95% and where a model is usually established within 3–12 weeks (Figure 3b; Einarsdottir et al., 2014).

Figure 3. PDX models and 3D tumorsphere cultures of uveal melanoma to find TRLs.

(a) Immunohistochemical (IHC) staining of HMB45, SOX10, S100, MELAN-A (MART1), and PRAME expression in one patient-derived xenograft (PDX) (see also Figure 1—figure supplement 1 for additional PDX models), as well as Sanger sequencing to verify mutation of GNA11 at Q209L. Color reactions were with magenta. (b) Tumor growth curves of established PDX models. (c–d) Tumor-infiltrating lymphocytes (TIL) therapy experiments in PDX models. UM1 and UM22 were transplanted subcutaneously into NOG or hIL2-NOG mice. When tumors were palpable, 20 million autologous TILs were injected in the tail vein. Tumor growth was monitored by a caliper. Tumor sizes are average size ± SEM. (e) Images of tumor spheroids created (scale bar represents 100 µM). UM23 was unable to form spheres. (f–g) Activation (41BB+) and degranulation (CD107a+) status for spheres and TIL co-cultures, (f) and (g) respectively, using n=2 or 3 replicates.

Figure 3—figure supplement 1. Immunohistochemical (IHC) showing expression of HMB45, SOX10, S100, MELAN-A, PRAME expression in patient-derived xenograft (PDX) models.

Staining color used was either diaminobenzidine (DAB) (top four samples) or magenta.

Figure 3—figure supplement 2. Characterization of TILs used in sphere culture experiments.

(a) Gating strategy for Far Red Dye tagged tumor cells (APC) and Jade dye tagged tumor-infiltrating lymphocytes (TILs) (FITC). (b) Far Red fluorescent dye incorporated spheres co-cultured with Jade dye labeled TILs (green) co-culture with a co-localized population (yellow) at 24 hr time-point.

We have previously developed a method to study ACT in PDX models, termed PDXv2 (Forsberg et al., 2019; Jespersen et al., 2017; Webb et al., 2010). Survival, expansion, and activity of TILs and CAR-T cells were improved in this model using human IL2-expressing transgenic NOD SCID IL2 receptor gamma knockout mice (hIL2-NOG) (Forsberg et al., 2019; Jespersen et al., 2017; Ny et al., 2020; Huuhtanen et al., 2022; Nilsson et al., 2022). Tumor eradication by TILs correlates with immunotherapy response in patients with melanoma (Jespersen et al., 2017; Ny et al., 2020). We conducted a PDXv2 experiment using PDX models from two patients with UM. The injection of autologous TILs into hIL2-NOG mice carrying UM1 resulted in tumor rejection, whereas UM22 tumors could not be eradicated (Figure 3c and d). This constitutes proof of concept that TILs from UM metastases can kill autologous tumors in a PDXv2 model.

TILs contain subsets capable of ex vivo expansion and killing in 3D tumor sphere cultures

To gain further insight into the capacity of T cells in UM biopsies to kill tumors, we next investigated the tumor reactivity of ex vivo expanded TILs on autologous UM cells. Since most UM PDX models did not grow fast enough to routinely conduct ACT experiments, we used the tumors grown in the mice to generate 3D tumor sphere cultures (Ma et al., 2021; Weiswald et al., 2015) of high viability for TIL experiments (Gopal et al., 2021; Al Hity et al., 2021). We generated tumor sphere cultures from 10 PDX models and expanded their autologous TILs (Figure 3e). When tumor spheres and TILs were co-cultured, TILs from patients UM1 and UM9 were activated the most (surface 4-1BB+) and also degranulated (surface CD107a+; Figure 3f and g). These TILs also bound to the tumor spheres, as assessed by labeling tumor spheres and TILs with different vital dyes and then imaging with a real-time cell culture imager (Figure 3—figure supplement 2). The UM22 TILs were hardly able to get activated by their autologous tumor spheres (Figure 3b), which likely explains their lack of activity in the PDX ACT experiment (Figure 3d). However, since UM22 is capable of presenting MART1 on HLA-A2, and since we previously identified a minor population of MART1-reactive TILs in this biopsy (Karlsson et al., 2020), we decided to expand and enrich this subset using MART1 peptide HLA-A2 dextramers. When using these cells, we were able to achieve potent activation and degranulation (Figure 3f and g). This confirms that T cells extracted from biopsies can contain subsets capable of recognizing and killing tumors, which can be expanded ex vivo.

Determining the phenotypes of TRLs in biopsies using scRNA- and TCR-seq

To gain clues on how to identify reactive T cells from biopsies through their transcriptomic phenotypes, we co-cultured UM1 and UM9 TILs with their autologous tumors and sorted 4-1BB activated T cells by cell sorting. We prepared RNA from these cells, and as a control, from MART1-specific allogenic UM46 TILs and sequenced their TCRs. Then, using the single-cell atlas built from biopsies, we mapped those TCRs back to T cell subsets from their original tumors. This showed that those TRLs predominantly emanated from exhausted and late-activated cells (Figure 4a and b). UM9 also contained a sizeable amount of TRLs in the NK-like/cytolytic cluster (Figure 4b). The lack of potential TRL clones identified in the TCF7+ axis of the CD8+ T cell atlas from biopsies is also compatible with a previous report showing that TCF7+ T cells in melanoma tend to be non-reactive bystanders (Li et al., 2019).

Figure 4. TCR sequencing of CD8 T cells.

(a) Proportions of cells from each sample that match clonotypes found in experimentally identified 4-1BB+ and MART1-reactive T cells and which are present in a given biopsy CD8+ T cell cluster. The difference between matching and non-matching cells in each cluster is shown on the right, highlighting subsets that are enriched among the reactive T cells. (b) Biopsy CD8+ T cells with clonotypes matching identified reactive T cells, highlighted in the Uniform Manifold Approximation and Projection (UMAP) representation. (c) As in (b), but highlighting cells with clonotypes shared between biopsies and tumor-infiltrating lymphocytes (TILs). (d) Proportions of biopsy CD8+ T cells in each cluster matching either clonotypes from TIL cultures or experimentally identified reactive T cells. Statistical differences in frequency were determined using binomial tests between the frequency of the latter in each cluster relative to the background frequency of all TIL clonotypes present in the same cluster. Frequencies were calculated as number of cells from a given category in a given cluster divided by the total number of cells from that category, where category refers to either TILs or 4-1BB+/MART1-reactive cells. p-Values were adjusted for multiple testing using Bonferroni correction. (e) Distributions of cells matching the two different categories of experimentally identified reactive cells among biopsy CD8+ T cells clusters. (f) Cells from single-cell RNA-seq (scRNA-seq) of TIL cultures that have clonotypes matching experimentally identified reactive cells. (g) Shared clonotypes among all biopsy and TIL samples. (h) GLIPH2 clusters of significantly similar and HLA-restricted clonotypes (Huang et al., 2020), mapped to known antigens in public databases using TCRMatch (Chronister et al., 2021), based on CDR3β sequences. (i) Biopsy CD8+ T cells with clonotypes in either high- or low-confidence GLIPH2 clusters that match MART1 motifs in public databases.

Figure 4—figure supplement 1. Mapping TRLs in single-cell data and databases.

(a) Biopsy CD8+ T cells with clonotypes that match single-cell RNA (scRNA)-sequenced tumor-infiltrating lymphocytes (TILs), as in Figure 4c but shown separately for each sample. (b) TIL cells with clonotypes that match any corresponding biopsy. (c) Proportions within each CD8+ T cell cluster that match either 4-1BB+ or MART1-selected TILs from UM1, UM9, and UM46, respectively. (d) Overlap among TCRβ chains from the reactive cells sorted out from these three samples. (e) Shared TCRβ chains among the subset of cells from biopsies and TIL cultures that also have a match to any TCRβ chain among experimentally identified reactive TILs. (f, g) Matches to any CDR3β-antigen pair in public databases, as inferred by TCRMatch, for CD8+ T cells from biopsies and TILs. Samples from liver and subcutaneous metastases shown in (f) and (g), respectively. (h) As in (f–g), for experimentally identified reactive TILs. (i) Cells from biopsy CD8+ T cells that were included in any high-confidence GLIPH2 cluster.

To evaluate which T cells in biopsies were reflective of those capable of generating TIL cultures, we further mapped our previously published TCR-seq data of UM TILs to the biopsy data (Figure 4c) to compare them to the clones matching experimentally identified TRLs in the biopsy (Figure 4b–d, Figure 4—figure supplement 1a and b, and Supplementary file 2). Cells from the TIL dataset tended to correspond to biopsy T cells with exhausted and late-activated profiles, similar to the experimentally identified TRLs. This might be expected since both have undergone expansion, which might favor T cells in certain phenotypic states. However, TRLs were significantly over- and underrepresented in some phenotypic clusters compared to TILs in general. Overrepresented clusters included proliferative exhausted and IL7R+ exhausted T cells. The underrepresented ones were activated (non-exhausted), naïve/memory-like, and FCGR3A+ NK-like T cells. The distribution of TRLs across the clusters was also different for cells sorted by 4-1BB positivity versus MART1 recognition. MART1-reactive cells were more likely to be found in exhausted and late-activated clusters than 4-1BB+ cells, which tended to map to earlier activation states as well to the NK-like cluster (Figure 4e). The latter may be partly due to their prevalence in UM9 (Figure 4—figure supplement 1c).

Mapping TCRs from enriched 4-1BB+ cells to TILs in corresponding patients from previously published TCR-seq data revealed that a large fraction of TILs in these cultures was among the identified TRL clones (Figure 4f). In six of the eight previously published TIL cultures, we were able to recover and sequence the corresponding cryopreserved biopsy of that patient. Five of the six cultures had many TCRs that overlapped between TILs and biopsies (Figure 4g and Figure 4—figure supplement 1a and b). UM3 did not have any matching TCRs between the biopsies and expanded TIL cultures. No overlap was observed between samples from different patients with the 4-1BB+/MART1+-selected TCRs (Figure 4—figure supplement 1d). This is reasonable given the different HLA genotypes of these patients. In general, very few clonotypes were shared between samples from different patients (Figure 4—figure supplement 1e). In summary, experimentally identified TRLs mapped back to biopsies were overrepresented in proliferative exhausted and IL7R+ exhausted T cells, and these T cells could also be identified in sequenced cultures of expanded TILs.

Identifying potential antigens recognized by TRLs

To understand potential antigens recognized by a given TCR, we searched for matches in public databases of known TCRβ chains, based on CDR3β motifs, using TCRMatch (Chronister et al., 2021). The results confirmed that a number of TCRs matched known MART1-reactive clonotypes but also revealed matching TCRs that might recognize other melanoma antigens, such as PMEL and MAGEA1 (Figure 4—figure supplement 1f, g, and h). Overall, a large number of TCRs also matched clonotypes recognizing viral antigens such as influenza, EBV, and human coronaviruses. However, an even greater number lacked known matches, suggesting they might recognize tumor-specific antigens (neoantigens).

However, public databases are currently limited. As an alternative approach to nominate tumor-reactive clonotypes, we used an algorithm designed to predict clusters of TCRs that are sufficiently similar and HLA-restricted to potentially recognize the same antigens. GLIPH2 (Huang et al., 2020) nominated a few high-confidence clusters of similar TCR sequences based on the biopsy scRNA-seq data (Figure 4h and Supplementary file 3). To further test the hypothesis that each cluster recognizes the same antigen, we mapped the CDR3β motifs of these sequences to the same public databases. Interestingly, the highest-ranked cluster was overrepresented for sequences matching MART1, and this cluster also exhibited significant restriction toward HLA-A*02:01, on which MART1 is known to be presented. The remaining high-confidence clusters mainly matched antigens of viral origin or lacked matches to known antigens (Figure 4h).

Similar to the experimentally determined MART1-reactive T cells, the bioinformatically determined cluster of potential MART1-recognizing clones mostly belonged to the exhausted CD8+ T cell cluster (Figure 4i, Figure 4—figure supplement 1i). This supports the finding that TRLs found in biopsies may have predominantly exhausted profiles. In addition, a high-confidence group of clonotypes co-localized to the NK-like cluster was nominated by GLIPH2, although no sequences from this cluster matched the public TCR-antigen pairs (Figure 4h, Figure 4—figure supplement 1i). Given that the cells from this cluster were experimentally identified as 4-1BB+ upon tumor co-culture (Figure 4b), one might speculate that they could recognize unknown tumor antigens. Taken together, these findings suggest that both bioinformatically inferred TRLs from biopsies and ex vivo expanded TRLs tend to have exhausted phenotypes. Among these, T cells recognize known melanoma antigens such as MART1, PMEL, and MAGEA1, as well as unknown sequences that might be neoantigens.

TRLs home to the tumor of PDX mice bearing liver metastases of UM

TIL cultures contain both TRLs, virus-specific T cells, and bystander T cells that grow during IL2-containing culture conditions. We had previously found that TILs can home to subcutaneously growing skin melanoma metastases (Forsberg et al., 2019; Jespersen et al., 2017) but we were interested in determining whether they could identify and home to UM growing in the liver. The PDX models used to generate sphere cultures were prepared by subcutaneous transplantation. To instead generate a liver metastasis model, we performed splenic injection of UM22 into NOG mice (Sugase et al., 2020). Tumor cells forming liver metastases were further transplanted via tail vein injection into hIL2-NOG mice. After verifying tumors by ultrasound, autologous bulk or MART1-selected UM22 TILs were intravenously injected for TIL therapy studies (Figure 5a). The mice were harvested 3 weeks later, and tumors recovered from the liver were analyzed by flow cytometry, IHC, scRNA, and TCR-seq. As expected from Figure 3f and g, MART1-selected TILs had a higher degree of activation than unselected bulk UM22 TILs (Figure 5b), despite similar numbers of TILs present in the livers. There were no visible differences in the size of the liver tumors between mice receiving unselected or MART1-selected UM22 TILs. This may be due to T cell exclusion, since most TILs colonized the rim of the tumor (Figure 5c). More MART1-specific TILs than unselected TILs were able to get activated, as assessed by the activation markers 4-1BB and CD69 (Figure 5b, d, and e). This shows that T cells recognizing a tumor antigen can successfully migrate to liver metastases of UM in mice.

Figure 5. CD8 T cells get activated by tumor cells in vivo.

(a) Bulk unsorted tumor-infiltrating lymphocytes (TILs) or MART1-selected TILs were injected into hIL2-NOG mice carrying liver tumors from patient UM22. (b) Flow cytometry analysis of single-cell suspension liver metastasis, comparing treatment of UM22 TILs and UM22 MART1-specific TILs for CD3+, CD3+CD8+CD69+, or CD3+CD8+CD137+. (c) Immunohistochemical (IHC) with diaminobenzidine (DAB) showing tumor (SOX10) and TILs (CD3) within a liver metastasis, (d, e) corresponding analysis of 4-1BB+ TILs in a section, (d) and in image analysis comparing both treatments (e). Statistical tests in b and e were unpaired two-tailed t-tests, assuming equal variance. *: p<0.05; **: p<0.01. (f) Samples of tumor and TILs from the liver and spleen, respectively, were sequenced with single-cell RNA-seq (scRNA-seq). n=3 biological replicates were performed for each group of liver samples, and n=2 for spleen samples (out of which one spleen sample for MART1-selected TILs was the pooled material of two independent mice). Sequencing reads mapping to human and mouse were separated with XenoCell (Cheloni et al., 2021), after which cells were clustered and annotated as described in Methods. (g) Subclustering of CD8+ T cells identified three overall clusters, one of which represented a mixed profile of the other two, but with marked cell cycle activity. (h) Contributions of the different experimental conditions to each CD8+ T cell cluster. (i) Markers distinguishing CD8+ T cell clusters, identified using the FindAllMarkers function of Seurat (the union of the top 25 genes per condition, ranked by log2 fold change). Expression per experimental condition is shown below. (j) All TCRβ chains identified in each experimental condition. Subsets found in TIL culture scRNA-seq data are highlighted, as are any matches to antigens in public databases. (k) Differentially expressed genes between bulk TIL mixtures or MART1-selected TILs present in the livers of mice. A pseudo-bulk approach was used, summing read counts across all cells within a given replicate, and statistical testing performed with DESeq2 (Love et al., 2014). Genes with q<0.05 after Benjamini-Hochberg correction were considered significant. (l) Expression of KRT86, DUSP4, and LAYN in biopsy CD8+ T cells. (m) Mapping indiciated genes to phenotypic clusters from the single-cell RNAseq UMAP in Figure 2a. Figure generated with BioRender.com.

Figure 5—figure supplement 1. Characterization of TRLs in the UM22 PDX model.

(a) All cells from single-cell RNA-seq (scRNA-seq) analysis of patient-derived xenograft (PDX) samples, shown separately per experimental condition. (b) Copy number profiles of uveal melanoma (UM) cells in PDX samples, as determined by inferCNV (Tickle et al., 2019) analysis. Blue represents copy number loss and red copy number gain. (c) Markers discriminating CD8+ T cell clusters in biopsies, contrasted with CD8+ T cell clusters identified in PDX samples. (d) TCRβ chains identified in the PDX samples that match clonotypes in a separate scRNA-seq experiment of UM22 tumor-infiltrating lymphocytes (TILs). (e) Shared unique TCRβ chains between PDX samples from each experimental condition. (f) As in (e), shown relative to PDX CD8+ T cell clusters. (g) As in (f), but only including cells that also match any clonotype found in scRNA-seq from UM22 TIL cultures. (h) Arithmetic difference between (g) and (f), showing cells with TCRs uniquely detected PDX samples compared to the sequenced UM22 TIL culture. (i) All TCRβ chains identified in each PDX CD8+ T cell cluster. Clonotypes found in TIL cultures are highlighted, as well as any matches to public antigen databases. (j) Differentially expressed genes between bulk TIL mixtures residing in liver and spleen of PDX models. A pseudo-bulk approach (Squair et al., 2021) was used, summing read counts across all cells within a given replicate, and statistical testing performed with DESeq2 (Love et al., 2014). Genes with q<0.05 were considered significant.

Figure 5—figure supplement 2. Marker genes on T cells in the UM22 PDX model.

(a) Expression of marker genes related to a previously described tumor-reactive set of tissue-resident CD8+ T cells (Clarke et al., 2019) , among CD8+ T cells found in patient-derived xenograft (PDX) models. (b) Expression of KRT86, DUSP4, LAYN, HAVCR2, and IL7R in tumor-infiltrating lymphocyte (TIL) culture from UM22. Whether a given cell has a clonotype identified in PDX samples is indicated.

Determining the phenotypes of TRLs in PDX liver metastases using scRNA- and TCR-seq

Next, we profiled the livers and spleens of mice injected with unselected bulk or MART1-selected UM22 TILs by scRNA- and TCR-seq (Figure 5f, Figure 5—figure supplement 1a). After bioinformatically separating human- from mouse-derived cells and integrating the data, clearly distinguishable clusters emerged containing UM cells, all characterized by the expression of the differentiation genes TYR, MLANA, and PMEL. Their designation as UM cells was further confirmed by inference of their copy number profiles (Figure 5—figure supplement 1b). As expected, all UM cells resided in the liver (Figure 5—figure supplement 1a), despite the mice being injected with UM cells via the tail vein. The TILs discovered in these samples constituted three clusters: C1, C2, and a proliferative cluster with a mixed set of cells from the two main clusters (Figure 5g). Compared with CD8+ T cell clusters from the biopsies, the clusters found in mice did not separate with respect to the same set of markers. Rather, all clusters expressed markers compatible with late-activated or exhausted cytotoxic cells to some extent (Figure 5—figure supplement 1c). These TILs were detectable in both the liver and the spleen, although there was a notable correspondence between, on one hand, liver TILs and C1, and on the other, spleen TILs and C2 (Figure 5g–i).

Investigating the TCR clonotypes represented in the four experimental conditions revealed several clones unique to or overrepresented in the liver biopsies of mice with MART1-selected T cells, as compared to unselected TILs or spleen samples (Figure 5j). Given that these clones have been sorted to enrich for reactivity toward an antigen known to be presented by UM22 cells, and that they have successfully migrated to the liver tumor, this subset is likely to be strongly enriched for TRLs. Among these, some were also found in previous scRNA-seq data from expanded UM22 TILs (Figure 5—figure supplement 1d–i).

Differential expression analysis between this TRL subset and unselected TILs found in the livers of PDX mice revealed statistically significant genes that may be potential markers of the tumor-reactive phenotype (Figure 5k, Figure 5—figure supplement 1j, Supplementary file 4). This included the activation marker TNFRFS9 (4-1BB) as well as KIR2DL4, which has been found to be expressed in certain exhausted subsets of CD8+ T cells (Figure 2c; Li et al., 2019). Interestingly, this subset also bears a striking resemblance to a recently described group of TIM-3+IL7R- tumor-reactive tissue-resident memory cells in lung cancer, further characterized by the specific expression of KRT86, DUSP4, and LAYN (Figure 5—figure supplement 2; Clarke et al., 2019). This set was enriched in patients responding to anti-PD1 therapy (Clarke et al., 2019). Cells expressing this marker combination were also present in biopsy samples, predominantly in exhausted clusters (Figure 5l and m), as well as in TIL cultures from UM22 (Figure 5—figure supplement 2b). This might suggest avenues for enriching TRLs from biopsies for use with ex vivo expansion and ACT based on their gene expression profiles.

CD8 cells diminish after tumor eradication in PDX models

Although the UM22 PDX model demonstrated signs of TIL homing and activation, TILs were unable to cure the mice from the tumor, neither when grown subcutaneously (Figure 3d) nor in the liver (Figure 5c). We therefore developed a liver metastasis model of UM1, which responds to TILs in a subcutaneous PDX model (Figure 3d). Injection of TILs into hIL2-NOG mice bearing UM1 liver metastases resulted in a complete clearance or tumors as assessed macroscopically and by ultrasound imaging. To characterize the cellular composition of T cells in the livers and spleens of the PDX mice treated with TILs, we performed scRNA- and TCR-seq. Clustering analysis identified 15 clusters of human cells that were annotated by known marker genes defined by differential expression analysis and supported by literature (Figure 6b and c, Figure 5—figure supplement 1c, Supplementary files 5 and 6). Notably, and at variance with the UM22 model, none of the clusters contained any melanoma cells, confirming that the TILs had eradicated the tumor.

Figure 6. Characterization of CD8 T cells in the UM1 PDX model following tumor eradication.

(a) Establishment of another uveal melanoma (UM) liver metastasis model. UM1 patient-derived xenograft (PDX) cells were injected in the tail vein after having been serially transplanted in spleen followed by harvesting from liver. Ultrasound confirmed growth in liver before injection of autologous UM1 tumor-infiltrating lymphocytes (TILs) or HER2 CAR-T cells as controls. (b) Response to TILs as assessed by ultrasound monitoring. (c) Uniform Manifold Approximation and Projection (UMAP) of T cells in the UM1 liver metastatic model, showing 15 different cell populations. A list of marker genes was used to annotate clusters. The marker genes list was compiled from differential expression analysis and literature. (d) Dot plot showing an average expression of marker genes and detection rate of cells in which the marker gene is detected across 15 cell populations.

Figure 6—figure supplement 1. Characterization of TILs in vivo in the UM1 PDX.

(a) Uniform Manifold Approximation and Projection (UMAP) split on spleen and liver. Liver has 20,418 cells and spleen has 18,208 cells. (b) Bar plot showing number of cells across each cell type comparing between spleen and liver. (c) Violin plot and UMAP showing expression distribution for CD4 and CD8A genes across T cell clusters in liver and spleen. (d) Heatmap showing gene expression (z-score) values of pseudo-bulk (‘summed counts’) for each cell population. (e) Feature plot showing gene expression of two common marker genes for each cell population.

We compared number of cells coming from liver and spleen samples. A total of 20,418 cells were found in liver and 18,208 from spleen (Figure 6—figure supplement 1a and b). Strikingly, even though the TILs injected were primarily CD8 positive, the T cells remaining in liver and spleen after tumor eradication were primarily CD4+ positive cells (Figure 6—figure supplement 1b and c). Most CD8 cells resided in cluster 11 with fewer cells also in clusters 8 and 12 due to expression of cell cycle genes. Re-clustering of cluster 11 resulted in two subpopulations (11.0 and 11.1) annotated as CD8+ cytotoxic expressing killer cell lectin like receptor G1 (KLRG1) and CD8+ cytotoxic expressing higher levels of CD8A and CD8B (Figure 6—figure supplement 1d). These cells also expressed higher levels of CTLA4, TIGIT, and TNFRSF9 (encoding activation marker CD137/4-1BB) and the γδTcell markers TRGC2/TRGV9. Cluster 8 expressed high levels of Ki-67 (MKI67) assembly factor for spindle microtubules (ASPM), stathmin 1 (STMN1), DNA topoisomerase II alpha (TOP2A), TPX2 microtubule nucleation factor (TPX2) corresponding to the proliferative cellular state. Cluster 12 expressed minichromosome maintenance complex components 2, 6, 7 (MCM2, MCM6, MCM5, MCM7), GINS complex subunit 2 (GINS2), helicase, lymphoid specific (HELLS) that are associated with the proliferation and cell cycle state.

Next, we investigated the CD4+ cells and found that clusters 1, 3, 5, 9 expressed programmed cell death 1 (PD-1, also known as PDCD1), lymphocyte activation gene 3 (LAG3), cytotoxic T lymphocyte-associated protein 4 (CTLA-4), T cell immunoreceptor with immunoglobulin and ITIM domains (TIGIT), thymocyte selection associated high mobility group box (TOX), and eomesodermin (EOMES) marker genes that are associated with a dysfunctional or exhausted state and was annotated as a CD4+ exhaustion TOX+, CD4+ exhaustion PTCD1+, CD4+ exhaustion LAG3+, and CD4+ exhaustion TIGIT+, respectively. Three adjacent clusters 0, 2, 10 expressed interleukin 7 receptor (IL7R), rhotekin 2 (RTKN2), C-C motif chemokine receptor 3 (CCR3), and protein tyrosine phosphatase receptor type C (PTPRC, also known as CD45) marker genes corresponding to a memory T cells phenotype and were annotated as a CD4+ memory RTKN2+, CD4+ memory CCR3+, and CD4+ memory PTPRC+, respectively. More liver cells contributed to the CD4+ exhaustion cluster whereas the CD4+ memory cluster was primarily made up of spleen cells (Figure 6—figure supplement 1b). Finally, three adjacent clusters 7, 6, 4 were associated with a CD4+ Th17 phenotype that play important role in the inflammatory response expressing the transcription factor retinoic acid orphan receptor (ROR)γt (RORC), chemokine CC receptor 6 (CCR6), a killer cell lectin-like receptor B1 (KLRB1, also known as CD161), interleukin 4-induced 1 (IL4I1), cyclin-dependent kinase inhibitor 1A (CDKN1A), aquaporin 3 (AQP3) genes (Figure 6—figure supplement 1e) and annotated as a CD4+ Th17 CD161+ (KLRB1), CD4+ Th17 PTPN13+, and CD4+ Th17 PTPRM+ respectively.

Finally, we investigated the TCR clonotype abundance and visualized top 10 the most abundant clonotypes (Figure 7a). Clonotypes 302, 11155, and 21797 were the most abundant, and as expected, they were present in CD4+ cells. Clonotypes 302 mapped to the CD4+ exhaustion site, 111555 to the CD4+ memory, and 21797 to the CD4+ Th17 cluster (Figure 7b). To investigate if there were any differences in gene expression between these clonally expanded T cells in liver versus spleen, we performed pseudo-bulk analysis. Genes such as FASLG, ICOS, TIGIT, IL2RA that were higher in T cells residing in liver are known to be expressed during activation and exhaustion (Figure 7c). We also mapped where TRLs are residing by extracting the top 10 T cell receptor beta sequences obtained from the CD137-sorted UM1 TILs (Figure 3f–g and Figure 4). We identified a few of these in the PDX data and mapped them on the UMAP (Figure 7d). TRLs were residing on the CD4+ Th17 and proliferation/cell cycle clusters.

Figure 7. Mapping clonotypes in single-cell data from the UM1 PDX model.

(a) Bar plot showing top 10 the most abundant clonotypes for liver and spleen samples detected by the TCR-seq. (b) Uniform Manifold Approximation and Projection (UMAP) showing where clonotypes 302, 11155, and 21797 are residing. (c) Volcano plot highlighting differentially expressed genes between liver and spleen detected for the clonotype 302 using pseudo-bulk approach. log2FC positive means liver is upregulated relative to the control (spleen). (d) UMAP showing where tumor-reactive T lymphocytes (TRLs) (from Figure 3f and g, Figure 4) are residing.

Discussion

UM continues to be a condition with unmet medical need, despite recent advances in locoregional therapy and immunotherapy (Fu et al., 2022). However, the fact that a T cell engager is prolonging survival in patients with UM having an HLA-A2 genotype (Nathan et al., 2021) and that immune checkpoint inhibitors in various combination therapies on occasion can yield clinical benefit (Piulats et al., 2021; Pelster et al., 2021; Ny et al., 2021) suggests that if we learn more about the tumor immunology of UM, then we can devise new immunotherapies. Here, we generated datasets and models that can be used not only to characterize tumor immunity, but also to functionally test hypotheses and targets in autologous culture and in vivo systems.

ACT with TILs has also been tested in patients with metastatic UM (Chandran et al., 2017). Although there were responses in a few patients, the data clearly show room for improvement. The HAITILS trial (NCT04812470) addresses whether hepatic arterial infusion results in local accumulation of TILs in liver metastases. However, this will not circumvent the low mutational burden of UM (other than iris melanoma; Karlsson et al., 2020) and the potentially low number of TRLs in the cell product (that we see in this study). A greater understanding of TRLs in patients’ metastases is therefore of interest to learn how to enhance the activity of TIL products and how to expand, for instance, neoantigen- or melanoma-associated antigen-specific TRLs (van den Berg et al., 2020; Seiter et al., 2002). Here, we aimed to identify tumor-reactive lymphocytes in terms of marker characteristics and functional assessments of TILs. These analyses suggest that some of the TILs in our cultures that became activated (4-1BB positive) when binding tumor cells were already TRLs (4-1BB or PD-1 positive) when they resided in the tumor. However, they were never able to clear the tumor, and most likely became exhausted or dysfunctional in the microenvironment.

Presence of TRL-specific markers adds validity to a strategy whereby enrichment of TRLs could be done directly from the tumor by sorting for those markers. Indeed, previous studies have employed activation, exhaustion, or tissue-resident markers such as 4-1BB (Ye et al., 2014), PD-1 (Inozume et al., 2010; Gros et al., 2014; Zheng et al., 2022), CD39 (Zheng et al., 2022; Simoni et al., 2018), and CD103 (Webb et al., 2010; Duhen et al., 2018), or the specificity of a TCR for a melanoma-associated or melanoma-specific antigen (Yee et al., 2002; Dudley et al., 2002). However, a risk with this strategy is that some TRLs are missed. Most notably, we revealed TRLs in clusters of T cells negative for activation and exhaustion markers. These cells are interesting since they express high levels of perforin and granzyme B, suggesting that they are cytotoxic. It is tempting to speculate that these cells are effector or effector memory cytotoxic T cells since they express RNA encoding CD16, FCRL6, and CX3CR1 (Clémenceau et al., 2008; Schreeder et al., 2008; Gerlach et al., 2016). In the UM9 sample, most of the T cells in the biopsy that mapped to TRLs identified in the tumor sphere assay resided in this cluster.

PDX models have been developed for UM, but to our knowledge, this is the first time they have been used to study immunotherapy using autologous TILs and single-cell sequencing. Although one model responded to TIL therapy, the other model did not. In the UM22 liver metastasis model, we found evidence that T cells were pushed out of the tumor. This may partly explain the failure of immune checkpoint inhibitors in many patients (Algazi et al., 2016) and why ACT with TILs does not work for the majority of patients with metastatic UM (Chandran et al., 2017). Nevertheless, we profiled TRLs that preferentially migrated to liver metastases in mice, as opposed to non-reactive T cells and those that migrate to the spleen. We identified marker genes and phenotypic states that might be used to identify relevant subsets present in patient biopsies and TIL cultures. These T cells appeared to resemble a recently described IL7R-HAVCR2+KRT86+DUSP4+LAYN+ subset, which has been associated with response to anti-PD1 therapy in lung cancer (Clarke et al., 2019). They also tended to express ITGAE (CD103), suggesting that they may be tissue-resident memory cells. This subset was identified in both PDX liver metastases, cultured TILs, and patient biopsies. The phenotypes of cells with this marker combination were predominantly exhaustion-like in the PDX samples as well as in biopsies. This is compatible with the results from the profiling of biopsy CD8+ T cells, where mapping of experimentally identified MART1-reactive T cells back to biopsies and computational prediction of tumor antigen-recognizing T cells suggested an overrepresentation of TRL candidates in exhausted and late-activated gene expression-based clusters. Whether cells expressing these canonical exhaustion-related markers are truly dysfunctional or capable of eliciting an antitumor response remains to be explored (Clarke et al., 2019).

In the UM1 liver metastasis model, the tumor was eradicated, demonstrating that TILs can cause a complete response also of liver metastases. The single-cell data obtained confirmed the loss of melanoma cells but it also resulted in a loss of CD8+ T cells. It is well known in immunology that the fate of CD8 T cells after an acute immunological reaction is either death of the effector cells or development of memory. A strength of the PDX experiment is that we could model. But since we harvested the livers after tumor eradication, a limitation is that we could not investigate the expression of genes in T cells that are actively combating the tumor. However, since the CD4+ T cells were present in high numbers, we could compute differences between the T cells of the same clonotype residing either in liver or in spleen. This analysis did indeed suggest that some level of activation/exhaustion remained even after tumor clearance. Moreover, some of the TCRs from the co-culture experiment of tumorspheres and TILs designed to reveal TRLs were indeed identified in vivo among the CD8 and the proliferative cells. We are therefore confident that PDX models such as UM22 and UM1 will be useful in future studies where T cells are sorted for potential markers prior to expansion or infusion, where the T cells are genetically engineered or grown in cocktails of factors to enhance potency and stemness. We also propose the model to be appropriate for further humanization to unravel the interplay between tumor cells, liver cells and stroma cells, myeloid cells, and other immune cells.

When mapping TCRs to TCR databases, we noted not only similarities to known melanoma-associated antigens, but also to viruses. Importantly, TCRs that recognize different antigens can share CDR3β chains. The TCR alpha chain is even more promiscuous but not always reported in these databases, which makes it difficult to find accurate TCR pair matches. Therefore, until current databases have been expanded and improved, we cannot be certain about the exact antigen matches of some of these receptors. Nevertheless, an increasing body of literature suggests that bystander cells in tumors are often antigen-experienced, can be directed against bacteria or viruses, and even participate in tumor control, for instance during immune checkpoint inhibitor therapy (Rosato et al., 2019; Çuburu et al., 2022; Hu et al., 2022). Whether the presence of foreign antigen-specific T cells is beneficial for cellular immunotherapy remains to be evaluated.

In summary, we profiled the tumor microenvironment of metastatic UM and identified the phenotypes of tumor-reactive T cell subsets. The results of this study provide new avenues for targeted expansion of T cells to improve cellular and immune checkpoint immunotherapy. Importantly, however, to potentiate TIL therapies in UM, immune cold tumors likely need to be turned hot, i.e., to create an inflammatory environment promoting T cell infiltration. Simultaneously, lymphodepletion is needed to avoid existing T cells sequestering IL2 needed for the TIL infusion product. Tumor-targeting treatments such as IHP or percutaneous hepatic perfusion may fulfill these needs and could act as adjuvant therapy for both immune checkpoint inhibitor treatment and TIL therapy. As melphalan triggers immunogenic death, influx of new T cells into cold tumors may occur, partly explaining the very high response rate of these therapies (Olofsson Bagge et al., 2023). Ongoing clinical trials are combining locoregional therapy with ICI (NCT04463368) to capitalize on this concept. However, data from these trials have not yet been published. Combination of liver-directed therapy and TIL therapy is also warranted in future animal experiments and clinical trials. We propose that PDX models are optimal to capture the intra-patient heterogeneity of cancer and that multiple mechanisms in tumor immunology and treatment can be modeled at single-cell resolution.

Methods

Patient samples

According to the ethical approval of the Regional Ethical Review Board (#289-12, #44-18, and #144-13), the patients provided oral and written information with signed informed consent agreements. Biopsies were extracted from either subcutaneous or liver metastases (Supplementary file 7). Tumor biopsies were divided into pieces that were snap-frozen, minced, embedded in formalin blocks, or used for TIL cultures and cryopreservation. Snap-frozen tumor fragments were homogenized using a Bullet Blender (Next Advance, NY, USA). DNA and RNA were extracted using the AllPrep DNA/RNA kit (QIAGEN, Germany).

Immunostainings of sections

FFPE patient samples were evaluated using H&E staining and IHC. IHC was performed using an autostainer (Autostainer Link 48, Dako) with primary anti-human antibodies against Melan-A (clone A103; Dako), S100 (IR504; Dako), PRAME (E7I1B, CST), SOX10 (E6B6I, CST), HMB45 (M0634; Dako), CD3 (D4V8L, CST), and CD137 (E6Z7F, CST). HRP DAB or HRP Magenta (Dako) was used to stain the proteins of interest, and counterstaining was performed using hematoxylin.

We also performed multiplex immunofluorescence using The PhenoCycler-Fusion system (Akoya Biosciences) at the Spatial Proteomics Facility at SciLifeLab, Stockholm. The data for CD3, ICOS, and PD-1 primary antibodies are shown here; however, the full dataset will be published elsewhere.

PDX models

Animal experiments were conducted in conformity with EU directive 2010/63 (Regional Animal Ethics Committee of Gothenburg approvals no. 2014-36 and 2018-1183). PDX models were generated by transplanting small pieces into the flank (subcutaneous) of 6- to 8-week-old female immunocompromised, non-obese, severe combined immune deficient interleukin-2 chain receptor γ knockout mice (NOG mice) or hIL2-NOG mice (Taconic, Denmark). Cryopreserved single cells from biopsies were thawed and transplanted by splenic injection to create liver metastases in NOG mice or hIL2-NOG mice (Taconic, Denmark). All PDX tumors were analyzed using IHC, and their identity was verified by staining patient biopsies using clinically graded antibodies. For subcutaneous implants, tumor growth was monitored twice a week following TILS treatment. After establishment, the liver metastasis models formed liver metastases equally well by splenic injection and tail vein injection. The mice were evaluated using ultrasound scans before TIL treatment and grossly assessed at the time of sacrifice.

Generation of TILs and MART1-specific T cells

The metastatic fragments were placed in culture medium (90% RPMI 1640 [Invitrogen], 10% heat inactivated human AB serum [HS, Sigma-Aldrich], 6000 IU/ml recombinant human IL2 [Peprotech], penicillin and streptomycin [Invitrogen]) to generate young TILs (yTILs) as previously described (Jespersen et al., 2017). MART1-specific yTILs were identified as previously described (Karlsson et al., 2020) and sorted using FACSAria III (BD Biosciences). yTILs were expanded using a standard REP (Rapid Expansion Protocol), constituting irradiated (40 Gy) allogeneic feeder cells, CD3 antibody (clone: OKT3) (30 ng/ml) (Miltenyi), medium (AIM-V, Invitrogen) supplemented with 10% HS and 6000 IU/ml IL2 and monitored for 14 days with media changes (Jespersen et al., 2017). After completion of the REP, MART1 specificity was confirmed using MART1-specific dextramers (Immudex, Copenhagen, Denmark). REP and REP MART1 TILs were used for downstream, cell line, spheroid, and in vivo experiments.

Spheroid assay

PDX models were harvested, and tumor tissue was digested using a human tumor dissociation kit (Miltenyi) according to the manufacturer’s instructions in conjunction with the gentleMACS Octo tissue dissociator (Miltenyi). This single-cell suspension was used to prepare spheroids using an ultra-low attachment, 96-well, clear, round-bottomed plates (Corning, NY, USA), and GBM-MG serum-free media (CLS Cell Lines Service, Germany) for 72–96 hr in a 37°C incubator. To monitor treatment responses, cells were labeled with CellVue Jade Cell Labeling Kit (Thermo Fisher Scientific) and/or CellVue Claret Far Red Fluorescent Cell Linker Mini Kit (Sigma) and used for co-culture experiments. The 3D images were captured using IncuCyte (Essen BioScience)±TILs treatment, followed by imaging and analysis at 72 hr with IncuCyte Live Cell Analysis. Tumor spheres stained with a far-red dye were measured for the red object area/µm (Buder et al., 2013) mask with >100 μm. Flow cytometry was used to confirm the positively stained immune cells (Jade-FITC channel) and cancer cells (far-red APC channel).

Co-culture and TCR-seq

A human UM cell line derived from patient UM22 (Karlsson et al., 2020) was grown in complete medium (RPMI-1640 supplemented with 10% FBS, glutamine, and gentamycin) and cultured at 37°C with 5% CO2. The cell line or spheroids were analyzed with or without TILs, and the co-culture was analyzed using human CD3 (FITC and HIT3a, BioLegend), CD8 (APC and SK1, BioLegend), CD107a (APC and H4A3, BD Biosciences), CD137 (PE and 4B4-1, BioLegend), active caspase-3 (clone C92-605; BD Biosciences), and granzyme B (clone GB11; BD Biosciences) antibodies for 30 min at 4°C. Flow cytometry data were acquired using the BD Accuri C6 and BD Accuri C6 Plus (BD Biosciences). The supernatants of these samples were subjected to Granzyme B secretion analysis using an ELISA kit (R&D Systems). Following detection of the activated cell populations in certain samples, a microfluidics-based flow cell sorting instrument (WOLF NanoCellect, San Diego, CA, USA) was used to sort CD3+CD137+ cells from the co-culture of tumor spheres and TILs. The RNA of the sorted cells was extracted using a Nucleospin RNA XS kit (Macherey-Nagel), and its concentration was measured using QuBit. Libraries were prepared for TCR profiling using a SMARTer Human TCR a/b Profiling Kit v2 (catalog number 63478) and subsequent NGS library preparation and bead purification; Tapestation 4200 (Agilent) was used to confirm sample quality control with automated electrophoresis.

Statistical analysis for spheroid assays

Statistical tests for the spheroid assays were done using two-tailed unpaired t-tests, with p-values represented as * for p<0.05 and ** for p<0.01; all error bars represent standard error (SEM), unless otherwise stated.

scRNA-seq data analysis

Library preparation

5’ GEX and TCR-seq: Single-cell suspensions prepared from cryopreserved UM biopsies were subjected to dead-cell removal using beads (Miltenyi Biotech) before quality and viability validation and loading on Chromium instrument (10x Genomics). Single-cell 5’ gene expression and TCR libraries were prepared using vendor protocol with the Next GEM Single Cell 5' Kit v2 (10x Genomics, catalog number 1000263), Library Construction Kit (catalog number 1000190), Single Cell Human TCR Amplification Kit (catalog number 1000252), and Dual Index kit TT Set A (catalog number 1000215). Libraries were sequenced for quality check on iSeq100 (Illumina) before deep sequencing on NovaSeq 6000 (Illumina) in the required format to yield at least 20,000 reads per cell for gene expression and 5000 reads per cell for TCR libraries.

Alignment and gene expression quantification

Paired scRNA-seq and TCR-seq data from biopsies were aligned to the 10x Genomics-provided GRCh38 reference genome (v. 2020A) and V(D)J reference (v. 7.0.0) using the Cell Ranger pipeline (v. 7.0.1) and the ‘multi’ command, with both ‘Gene Expression’ and ‘VDJ-T’ analyses activated, as well as with the setting ‘check-library-compatibility,false’.

Removal of ambient RNA

To remove ambient RNA, CellBender (v. 0.2.0) (Fleming et al., 2022) was used on the output from Cell Ranger, with the parameters ‘--fpr 0.01 --epocs 150’ together with expected and total cell numbers estimated from Cell Ranger quality control plots, as described in the instructions at https://cellbender.readthedocs.io/en/latest/usage/index.html (accessed September 15, 2022). The values for the parameters ‘--low-count-threshold’ and ‘--learning-rate’ were changed from their defaults for some problematic samples and iteratively adapted as needed after inspecting CellBender output quality control plots, as described in the documentation.

Separation of human and mouse reads in PDX scRNA-seq data

To separate reads of mouse and human origin in the scRNA-seq data from PDX samples, XenoCell (v. 1.0) (Cheloni et al., 2021; Conway et al., 2012) was used. As input, compressed fastq files, merged across lanes, were provided, together with reference genomes for human and mouse (the same genomes as for Cell Ranger). An index was created with the command ‘xenocell.py generate_index’, with the parameters ‘--threads 1 --memory 24’. Reads were classified according to host (mouse) or graft (human) with the command ‘xenocell.py classify_reads’ and the parameters ‘--barcode_start 1 --barcode_length 16 --threads 1 -memory 24 --compression_level 1’. Reads belonging to human and mouse, respectively, were then extracted into new fastq files using the command ‘xenocell.py extract_cellular_barcodes’ with the same values for barcode_start, barcode_length, threads, memory, and compression_level parameters as in the previous step, but with ‘--lower_threshold 0 --upper_threshold 0.1’ for graft and ‘--lower_threshold 0.9 --upper_threshold 1.0’ for host. The workflow was implemented in a Nextflow script available on GitHub, (copy archived at Karlsson, 2022). The new fastq files were then used for separate Cell Ranger runs for human and mouse reads for a given sample, as described above for human biopsy samples, except with the corresponding 10x Genomics-provided mouse reference genome instead for reads determined to be of mouse origin.

Quality control and cell filtering

Gene counts and assembled TCR sequences were imported to R with Seurat (Satija et al., 2015). An initial object was created and normalized with the NormalizeData function of the same package. Samples were integrated with fastMNN (Haghverdi et al., 2018) (batchelor package v. 1.10.0) via the convenience function RunFastMNN from the SeuratWrappers (v. 0.3.0) package (parameters: ‘features = 4000, auto.merge=T, d = 150, k = 15’). Assembled and annotated TCR sequences from Cell Ranger were further added to this object using djvdj (Sheridan, 2020). Cells with <200 reads were then removed.

A first set of potential doublets were identified using scDblFinder (v. 1.8.0), using default parameters. In addition to this, cells expressing conflicting lineage-associated markers were also marked as doublets. Incompatible marker sets were defined as any among CD3D, CD4, CD8A, CD8B or with assembled TCR receptor expressed together with any of MLANA, PMEL, TYR, HBA1, HBA2, or HBB; any cell with assembled TCR receptor expressed together with any of CD14, CD19, MS4A1, or JCHAIN; or any among MLANA, PMEL or TYR expressed together with CD14, CD19, MS4A1, JCHAIN, or NCR1. For the purpose of further refining doublet detection, an initial clustering was performed using the Seurat functions FindNeighbors (parameters: ‘reduction = ‘mnn’, dims = 1:150, k.param=5’) and FindClusters (‘resolution = 5’). For each cluster, a two-sided Fisher’s exact test was carried out, assessing whether that cluster was associated with doublets. p-Values were adjusted for multiple testing with Benjamini-Hochberg correction. In addition, each cluster was also assessed for higher than expected UMI counts per cell. Doublet clusters to be removed were nominated from any cluster with Fisher’s test odds ratio <1 and q<0.05, percent doublets greater than or equal to the 95% percentile of doublet percentages, or UMI counts above the 95% percentile. These thresholds were chosen in conjunction with visual inspection of plots regarding how these candidate clusters associated with other metrics of poor quality, such as proximity in UMAPs to cells expressing incompatible lineage markers and higher/lower than usual ribosomal or mitochondrial read percentage.

Global thresholds to keep cells fulfilling the following criteria were set: number of detected genes >200, UMI counts >500, ribosomal read percentage <40%, and percentage of genes in the top 100 expressed genes <80%. However, further specific filtering was made for two cell types, erythrocytes and plasma cells, known to have generally lower RNA content. Since common global thresholds on RNA and number of detected genes may filter out these cell types entirely, lower thresholds were motivated specifically for these. To establish suitable thresholds, candidate members of these cell types were first identified based on the distinctly expressed marker genes HLA1 (erythrocytes) and JCHAIN (plasma cells, with IL3RA and LILRA4 additionally demanded to be zero), respectively. Cells in each respective group (without having removed candidate doublets) were clustered based on Spearman correlation coefficients with the pheatmap function in the R package of the same name (default parameters) and visualized together. In both cases, cells were divided into clear clusters, some of which had predicted doublets were overrepresented. The latter clusters also tended to have higher than expected RNA content. Cluster without overrepresented candidate doublets were retained and the rest discarded from the dataset. Thresholds to keep cells with UMI counts and number of detected genes greater than the first quartile were defined for each of the two cell types based on cells in clusters not associating with doublets.

Further to this, additional quality filtering thresholds for remaining cells were determined adaptively using miQC (v. 1.2.0) (Hippen et al., 2021), which first creates mixture models to assess the distribution of number of genes detected versus percentage of mitochondrial reads in each sample and then filters cells based on a posterior probability of being compromised derived from this distribution. For a few samples that failed to converge with the default model, a one-dimensional mixture model was used instead. For all samples, the parameters ‘posterior_cutoff = 0.9’ and ‘keep_all_below_boundary = TRUE’ were used to filter out low-quality cells.

Sample integration, unsupervised clustering, and cell-type annotation

After removing low-quality cells and doublets, all samples were reintegrated using fastMNN, with the same parameters as described above. A UMAP embedding was created using the Seurat function RunUMAP (‘parameters: reduction = ‘mnn’, dims = 1:150’). Unsupervised clustering was performed with FindNeighbors (‘reduction = ‘mnn’, dims = 1:150’) and FindClusters (‘resolution = 2’). Genes differing between clusters were assessed using the FindAllMarkers function (default parameters). Clusters were annotated based on a combination of generally known cell-type markers (CD3D and CD3G for T cells, CD8A for CD8+ T cells, CD4 for CD4+ T cells, FOXP3 for regulatory T cells, NCAM1 and NCR1 for NK cells, CD19 and MS4A1 for B cells, HBA1, HBA2, and HBB for erythrocytes, CD14 for general monocytic cells, PMEL, MLANA, and TYR for melanocytic cells, as well as MKI67 for actively cycling cells) and markers compiled from literature (Li et al., 2019; Durante et al., 2020; Malone, 2021; Collin and Bigley, 2018; Buonomo et al., 2022). CD8+ T cells were determined and studied in greater detail to identify further subgroups. After subsetting the main Seurat object to only contain CD8+ T cells, these were subclustered by rerunning FindNeighbors (‘reduction = ‘mnn’, k.param = 15, dims = 1:150’) and FindClusters (‘algorithm = 1, resolution = 1.3’). A new UMAP was created using RunUMAP (‘reduction = ‘mnn’, dims = 1:150, min.dist = 0.05, n.neighbors = 15’). 16 clusters were obtained, among which some were merged based on similar expression of key marker genes. The new CD8+ T cell clusters were annotated based on a combination of literature-derived markers (Li et al., 2019; Zhang et al., 2021; Zhao et al., 2020; Crespo et al., 2013) for CD8+ T cell phenotypes as well as genes found to distinguish the clusters that could be related to specific phenotypes.

Processing of scRNA-seq from TIL cultures

scRNA-seq data of TIL cultures from a previous study (Karlsson, 2022) were obtained, loaded into a Seurat object in R, and normalized with the NormalizeData function (default parameters). Variable features were determined using FindVariableFeatures (‘selection.method = ‘vst’, nfeatures = 4000’). Cell cycle scores were then calculated using CellCycleScoring (default provided gene sets) and data scaled using ScaleData (‘vars.to.regress=c(‘S.Score’, ‘G2M.Score’)). Principal components were calculated using the RunPCA function (npcs = 30) and a UMAP embedding calculated using the RunUMAP function (dims = 1:30). These steps were performed on each sample separately, without sample integration.

GLIPH2 analysis

TCR sequences were clustered into inferred specificity groups using GLIPH2 (Huang et al., 2020). As input, unique TCR sequences and their counts from annotated CD8+ T cells were given, together with previously published HLA genotypes from bulk RNA and WGS (Karlsson, 2022). HLA genotypes for two additional samples were inferred from previously published exome data (Ny et al., 2021) using nf-core/hlatyping (v2.0.0) (Ewels, 2022), with the parameter ‘--genome GRCh37’. Clonotypes with more than one TCR alpha chain were excluded from analysis, as these are not supported by GLIPH2. Parameters provided were the option to use the v2 human reference for CD8 T cells and with ‘all_aa_interchangeable’ set to Yes. The tool was run from the web interface and according to instructions at http://50.255.35.37:8080/ (accessed December 21, 2022). High-confidence TCR clusters were considered those with V-gene bias score <0.05, at least three unique TCRs and at least three unique individuals contributing TCRs.

Copy number analysis

To determine copy number changes in UM cells from scRNA-seq data, the inferCNV R package (v. 1.10.1) (Tickle et al., 2019) was used according to the instructions provided at https://github.com/broadinstitute/infercnv/wiki (accessed November 2, 2022). A genomic position file was created from the reference annotation provided with Cell Ranger (v. 2020A) using the inferCNV helper script gtf_to_position_file.py. All cells annotated as UM were marked as ‘malignant’ and remaining cells were used as a reference diploid background as input to inferCNV, together with the following parameters: ‘cutoff = 0.1, cluster_by_groups = TRUE, denoise = TRUE, HMM = FALSE, num_threads = 1, analysis_mode=’samples’, output_format = ‘pdf’.

Differential expression analysis

Differentially expressed genes between experimental conditions in scRNA-seq data from PDX samples were assessed using a pseudo-bulk approach together. First, all the human CD8+ T cells for each sample were used to create gene subset expression matrix. All cells in a given sample were summed on a per-gene basis to create pseudo-bulk samples, imitating a bulk RNA-seq experiment. These samples were normalized and analyzed for differential expression between the conditions represented by the samples using DESeq2 (Love et al., 2014), with the parameters test = ‘LRT’, reduced = ‘~1’, relative to a full model specification encoding the two experimental conditions to be compared. For each pairwise comparison of conditions, only samples in those two conditions were included in the DESeq2 main object. To extract a table of results, the ‘results’ function from the same R package was used, with the parameter ‘alpha = 0.05’. Significant genes were visualized with the pheatmap function from the R package of the same name (v. 1.0.12), with the parameter ‘scale = ‘row’.

scRNA- and TCR-seq analysis for UM1 sample

FASTQs files lane-1 and lane-2 were concatenated and used as an input for the XenoCell (v. 1.0) tool to split scRNA-seq data on host (mm) and graft (hs). XenoCell classify_reads function was used to calculate the percentage of graft- and host-specific reads for each cellular barcode. Then, extract_cellular_barcodes function was used to extract cellular barcodes for the host (mm), and the same function was used to extract cellular barcodes for the graft (hs). The feature-barcode matrix for graft (hs) was generated using Cell Ranger multi (v. 7.0.1) pipeline. Parameters --localcores and --localmem were set to 32 and 240 respectively and human reference genome GRCh38 (v. 2020-A) and V(D)J reference (v. 7.0.0) were used. Next, CellBender remove-background CellBender (v. 0.2.0) function was used to remove ambient RNA. For each sample CellBender parameters were adjusted as follows --expected-cells (25L: 7573; 25S: 2274; 39L: 7973; 39S: 10527; 45L: 8438; 45S: 10447; 55L: 6318; 55S: 4143), --total-droplets-included (25L: 17000; 25S: 10000; 39L: 12000; 39S: 17000; 45L: 16000; 45S: 19000; 55L: 15000; 55S: 13000), --fpr (25L: 0.01; 25S: 0.01; 39L: 0.01; 39S: 0.01; 45L: 0.01; 45S: 0.01; 55L: 0.01; 55S: 0.01), --epochs (25L: 150; 25S: 150; 39L: 200; 39S: 150; 45L: 150; 45S: 200; 55L: 150; 55S: 150), --low-count-threshold (25L: 3; 25S: 5; 39L: 3; 39S: 3; 45L: 3; 45S: 5; 55L: 5; 55S: 5), --learning-rate (25L: 0.0000125; 25S: 0.00005; 39L: 0.00000625; 39S: 0.0000125; 45L: 0.00000625; 45S: 0.00000625; 55L: 0.00005; 55S: 0.00005).

R (v. 4.3.2) was used for downstream analysis. Data generated by CellBender was parsed using Read_CellBender_h5_Multi_Directory function from scCustomize (v. 2.0.1) package. Seurat object was created using CreateSeuratObject function from Seurat (v. 5.0.1) package with parameters min.cells and min.features set to 3 and 200 respectively. This resulted in the Seurat object with 23,884 genes across 64,091 cells. V(D)J data was added to the Seurat object using import_vdj function from djvdj (v. 0.1.0) library. Doublets were removed using scDblFinder function from scDblFinder (v. 1.16.0) package that resulted in 57,200 cells. Cells with a low coverage were removed before using scDblFinder (nCount<200, Germain, 2024). QCs filtering was performed to remove dying/dead, any remaining empty droplets or multiples cells. Cells with a high gene deletion number could indicate doublets or multiplets, thus, cells that have unique gene counts over 6500 or less than 400 were removed. Cells with a total read count (or UMI) less than 500 and more than 25,000 were filtered out. Finally, cells with mitochondrial counts more than 20% and ribosomal counts more than 15% were removed. Next, the data was normalized, scaled, and dimensionality reduction applied using NormalizeData, FindVariableFeatures (nfeatures = 2000), ScaleData, RunPCA. This step followed standardized Seurat workflow (‘Data analysis workflow, Integrative analysis in Seurat v5’, https://satijalab.org/seurat/articles/seurat5_integration, accessed 2024, Satija, 2024). Integration was performed using IntegrateLayers function from Seurat (v. 5.0.1) package with parameters set to method = CCAIntegration, orig.reduction = ‘pca’, new.reduction = ‘integrated.cca’. FindNeighbors (reduction = ‘integrated.cca’, dims = 1:70) and FindClusters (resolution = c(0.1–2.0), cluster.name = ‘cca_clusters’) functions from Seurat (v. 5.0.1) package were used and resolution 0.6 (15 clusters) was selected. FindAllMarkers (logfc.threshold = 0.3, min.pct = 0.3, only.pos = TRUE, test.use = ‘wilcox’) function from Seurat (v. 5.0.1) package was used to find differentially expressed genes. To visualize results DimPlot, VlnPlot, FeaturePlot, and DotPlot Seurat (v. 5.0.1) functions were used. For plotting the volcano plot we used EnhancedVolcano (Blighe, 2018); for plotting a bar plot ggplot2 (v. 3.5.0) was used. To pseudo-bulk, AggregateExpression function from Seurat (v. 5.0.1) package was used. The gene list was compiled using marker genes coming from differential expression analysis and from the literature. Results were visualized as a heatmap using plotGroupedHeatmap function from the scater (v. 1.30.1) package. To find clonotype abundance calc_frequency function from djvdj (v. 0.1.0) package was used; to visualize results plot_clone_frequency function was used from the same package.

TCR sequence analysis

Bulk TCR-seq data analysis

TCR sequences sequenced with the Takara SMARTer Human TCR a/b Profiling Kit v2 were aligned and assembled using the open source Cogent NGS Immune Profiler software (v. 1.0) provided by the same company, using the parameters ‘-r TCR’ and ‘-t Both’.

Matching of TCR receptors between samples

TCRs from scRNA-seq samples and/or bulk TCR-seq samples were matched based on identical combinations of CDR3β, Vb and Jb. This was the minimal common combination of criteria available on which to consistently compare TCRs across all TCR-sequenced samples, since the bulk TCR-seq protocol does not provide information in pairing of alpha and beta chains.

Matching of TCR receptors to public databases

To match TCRs to public databases with information on possible antigens they might recognize, TCRMatch (v. 1.0.2) (Chronister et al., 2021) was used (parameter: -s 0.97) together with a modified version of the reference database from IEDB (Vita et al., 2019) provided with the tool, which we expanded with additional information from VDJdb (Shugay et al., 2018), McPAS-TCR (Tickotsky et al., 2017), and TCR3d (Gowthaman and Pierce, 2019) (using the datasets for both cancer and virus-targeting TCRs). In addition, this database file was further modified to replace internal commas in antigen names, since the default formatting conflicts with the use of comma as the main delimiter in the database file. This facilitated easier analysis of the results. TCR matches were visualized with Sankey diagrams using the ggsankey R package (Sjoberg, 2021).

Overrepresentation analysis for TCRs in CD8+ T cell clusters

To assess whether the 41BB- or MART1 reactivity experiment would enrich for clones mapping to any phenotype relative to cultured TILs in general, a binomial test was used. The test considers the observed fraction of cells with matches to the reactivity screen versus an expected probability. Given that cells from the reactivity experiment were taken from a TIL culture, we wanted to assess whether these had a meaningfully different distribution compared to cultured TILs in general in regard to cluster representation. Therefore, the expected probability was estimated as the fraction of TILs that match any given cluster. Thus, these two frequencies were compared with one two-sided binomial test (binom.test in R) per cluster, and p-values further adjusted for multiple testing with Bonferroni correction.

Funding Information

This paper was supported by the following grants:

http://dx.doi.org/10.13039/501100002794 Cancerfonden to Roger Olofsson Bagge, Lars Ny, Jonas A Nilsson.

http://dx.doi.org/10.13039/100007436 Erling-Persson Family Foundation to Jonas A Nilsson.

http://dx.doi.org/10.13039/501100004359 Vetenskapsrådet to Roger Olofsson Bagge, Jonas A Nilsson.

http://dx.doi.org/10.13039/100016669 Sjöbergstiftelsen to Jonas A Nilsson.

Harry Perkins Institute of Medical Research start-up grants to Jonas A Nilsson.

Västragötaland Region ALF grant to Jonas A Nilsson, Roger Olofsson Bagge, Lars Ny.

http://dx.doi.org/10.13039/501100004063 Knut and Alice Wallenberg Foundation to Jonas A Nilsson, Roger Olofsson Bagge.

BioCARE strategic grants to Jonas A Nilsson.

Lion’s Cancerfond Väst to Samuel Alsen, Lars Ny.

http://dx.doi.org/10.13039/501100005009 Assar Gabrielssons Foundation to Vasu R Sah.

Gustaf V Jubileumsklinikens forskningsfond to Lars Ny.

Wilhelm och Martina Lundgrens Vetenskapsfond to Joakim W Karlsson.

Acknowledgements

We thank Carina Karlsson for providing technical support. Grant support came from Cancerfonden (to JAN and ROB), Familjen Erling Persson (to JAN), Knut and Alice Wallenberg Foundation (to JAN and ROB), Vetenskapsrådet (to JAN and ROB), Sjöbergstiftelsen (to JAN), BioCARE Strategic grants (to JAN), Lion’s Cancerfond Väst (to SA and LN), Västra Götaland Regionen ALF grant (to JAN, ROB, and LN), Assar Gabrielsson fond (to VS), Gustaf V Jubileumsklinikens forskningsfond (to LN) and Wilhelm & Martina Lundgrens Vetenskapsfond (to JK). JAN is the Inaugural Chair of Melanoma Discovery, primarily supported by donations to Perkins from family and friends of Scott Kirkbride, the MACA Cancer 200 ride for cancer research, the HBF Foundation and Perpetual. The Genomics WA facility (AS) is supported by BioPlatforms Australia, State Government Western Australia, Australian Cancer Research Foundation, Cancer Research Trust, Harry Perkins Institute of Medical Research, Telethon Kids Institute and the University of Western Australia. We also thank the Pawsey Supercomputing Centre for providing computational resources for this project.

Additional information

Competing interests

Author contributions

Ethics

Additional files

Supplementary file 1. Statistically identified marker genes for each CD8+ T cell cluster in Figure 2a, using the FindMarkers function in the Seurat R package.

Supplementary file 2. Complete statistics for the analysis of cluster representation among reactivity-screened TCRs matched to biopsy CD8+ T cells, as compared to representation among tumor-infiltrating lymphocytes (TILs) matched to the same clusters (Figure 4d).

Supplementary file 3. Complete statistical output from the GLIPH2 analysis referred to in Figure 4h and i and Figure 4—figure supplement 1i.

Supplementary file 4. Complete statistics for the differential expression analyses in Figure 5k and Figure 5—figure supplement 1j.

Supplementary file 5. Differentially expressed genes in different clusters of the Uniform Manifold Approximation and Projection (UMAP) in Figure 5.

Supplementary file 6. References for marker gene selection and annotation of clusters in Figure 5.

Supplementary file 7. Description of the patient samples analyzed in this study.

MDAR checklist

Data availability

All newly generated single-cell RNA/TCR sequencing data are available at NCBI Sequence Read Archive under the collective accession number PRJNA1135711. Single-cell RNA/TCR sequencing data from TILs were from a previous study and are available via controlled access from European Genome-Phenome Archive (Dataset ID EGAD00001006031).

The following dataset was generated:

Karlsson J Sah VR Bagge RO Kuznetsova I Iqbal M Alsén S Stenqvis S Saxena A Ny L Nilsson LM 2024 Patient-derived xenografts and single-cell sequencing identifies three subtypes of tumor-reactive lymphocytes in uveal melanoma metastases NCBI BioProject PRJNA1135711

The following previously published dataset was used:

Joakim K Jonas AN 2020 Sequencing of uveal melanoma metastases European Genome-phenome Archive EGAD00001006031

10.7554/eLife.91705.3.sa0
eLife assessment
Yang Yongliang Reviewing Editor Dalian University of Technology China

Solid
Valuable
This study presents valuable findings on tumor-reactive T cells in liver metastases of uveal melanoma (UM). The authors conducted single-cell RNA sequencing to identify potential tumor-reactive T cells and used PDX models for functional analysis. The evidence supporting their claims is solid. The work will be of interest to scientists working in the field of uveal melanoma.

10.7554/eLife.91705.3.sa1
Reviewer #1 (Public Review):
Reviewer
This work successfully identified and validated TRLs in hepatic metastatic uveal melanoma, providing new horizons for enhanced immunotherapy. Uveal melanoma is a highly metastatic cancer that, unlike cutaneous melanoma, has a limited effect on immune checkpoint responses, and thus there is a lack of formal clinical treatment for metastatic UM. In this manuscript, the authors described the immune microenvironmental profile of hepatic metastatic uveal melanoma by sc-RNAseq, TCR-seq, and PDX models. Firstly, they identified and defined the phenotypes of tumor-reactive T lymphocytes (TRLs). Moreover, they validated the activity of TILs by in vivo PDX modeling as well as in vitro co-culture of 3D tumorsphere cultures and autologous TILs. Additionally, the authors found that TRLs are mainly derived from depleted and late-activated T cells, which recognize melanoma antigens and tumor-specific antigens. Most importantly, they identified TRLs-associated phenotypes, which provide new avenues for targeting expanded T cells to improve cellular and immune checkpoint immunotherapy.

Comments on revised manuscript

The revised manuscript has addressed all my concerns.

10.7554/eLife.91705.3.sa2
Reviewer #2 (Public Review):
Reviewer
Summary:

The study's goal is to characterize and validate tumor-reactive T cells in liver metastases of uveal melanoma (UM), which could contribute to enhancing immunotherapy for these patients. The authors used single-cell RNA and TCR sequencing to find potential tumor-reactive T cells and then used patient-derived xenograft (PDX) models and tumor sphere cultures for functional analysis. They discovered that tumor-reactive T cells exist in activated/exhausted T cell subsets and in cytotoxic effector cells. Functional experiments with isolated TILs show that they are capable of killing UM cells in vivo and ex vivo.

Strengths:

The study highlights the potential of using single-cell sequencing and functional analysis to identify T cells that can be useful for cell therapy and marker selection in UM treatment. This is important and novel as conventional immune checkpoint therapies are not highly effective in treating UM. Additionally, the study's strength lies in its validation of findings through functional assays, which underscores the clinical relevance of the research.

Weaknesses:

The manuscript may pose challenges for individuals with limited knowledge of single-cell analysis and immunology markers, making it less accessible to a broader audience.

10.7554/eLife.91705.3.sa3
Author response
Karlsson Joakim W Author Harry Perkins Institute of Medical Research Perth Australia

Sah Vasu R Author University of Gothenburg Gothenburg Sweden

Olofsson Bagge Roger Author University of Gothenburg Gothenburg Sweden

Kuznetsova Irina Author Harry Perkins Institute of Medical Research Perth Australia

Iqba Munir Author University of Western Australia Perth Australia

Alsen Samuel Author University of Gothenburg Gothenburg Sweden

Stenqvist Sofia Author University of Gothenburg Gothenburg Sweden

Saxena Alka Author Harry Perkins Institute of Medical Research Perth Australia

Ny Lars Author University of Gothenburg Gothenburg Sweden

Nilsson Lisa Author Harry Perkins Institute of Medical Research Perth Australia

Nilsson Jonas A Author Harry Perkins Institute of Medical Research Perth Australia

The following is the authors’ response to the original reviews.

Reviewer #1 (Public Review):

Summary:

This work successfully identified and validated TRLs in hepatic metastatic uveal melanoma, providing new horizons for enhanced immunotherapy. Uveal melanoma is a highly metastatic cancer that, unlike cutaneous melanoma, has a limited effect on immune checkpoint responses, and thus there is a lack of formal clinical treatment for metastatic UM. In this manuscript, the authors described the immune microenvironmental profile of hepatic metastatic uveal melanoma by sc-RNAseq, TCR-seq, and PDX models. Firstly, they identified and defined the phenotypes of tumor-reactive T lymphocytes (TRLs). Moreover, they validated the activity of TILs by in vivo PDX modelling as well as in vitro coculture of 3D tumorsphere cultures and autologous TILs. Additionally, the authors found that TRLs are mainly derived from depleted and late activated T cells, which recognize melanoma antigens and tumor-specific antigens. Most importantly, they identified TRLs associated phenotypes, which provide new avenues for targeting expanded T cells to improve cellular and immune checkpoint immunotherapy.

Strengths:

Jonas A. Nilsson, et al. has been working on new therapies for melanoma. The team has also previously performed the most comprehensive genome-wide analysis of uveal melanoma available, presenting the latest insights into metastatic disease. In this work, the authors performed paired sc-RNAseq and TCR-seq on 14 patients with metastatic UM, which is the largest single-cell map of metastatic UM available. This provides huge data support for other studies of metastatic UM.

We thank the reviewer for these kind words about our work.

Weaknesses:

Although the paper does have strengths in principle, the weaknesses of the paper are that these strengths are not directly demonstrated. That is, insufficient analyses are performed to fully support the key claims in the manuscript by the data presented. In particular:

The author's description of the overall results of the article should be logical, not just a description of the observed phenomena. For example, the presentation related to the results of TRLs lacked logic. In addition, the title of the article emphasizes the three subtypes of hepatic metastatic UM TRLs, but these three subtypes are not specifically discussed in the results as well as the discussion section. The title of the article is not a very comprehensive generalization and should be carefully considered by the authors.

We thank the reviewer for the critical reading of our work. We have added more data and more discussion.

The authors' claim that they are the first to use autologous TILs and sc-RNAseq to study immunotherapy needs to be supported by the corresponding literature to be more convincing. This can help the reader to understand the innovation and importance of the methodology.

We have gone through the manuscript and found that we only refer to being first in using PDX models and autologous TILs to study immunotherapy responses by single-cell sequencing. While there are data to be deduced from other studies, we still believe this to be an accurate statement.

In addition, the authors argue that TILs from metastatic UM can kill tumor cells. This is the key and bridging point to the main conclusion of the article. Therefore, the credibility of this conclusion should be considered. Metastatic UM1 and UM9 remain responsive to autologous tumors under in vitro conditions with their autologous TILs.

UM1 responds also in vivo in the subcutaneous model in the paper. We have also finished an experiment where we show that this model also responds in a liver metastasis model. These data have been added in this revised version of the paper. We add two main figures and one supplementary figure where we characterize the response in vivo and also by single-cell sequencing of TILs.

In contrast, UM22, also as a metastatic UM, did not respond to TIL treatment. In particular, the presence of MART1-responsive TILs. The reliability of the results obtained by the authors in the model of only one case of UM22 liver metastasis should be considered. The authors should likewise consider whether such a specific cellular taxon might also exist in other patients with metastatic UM, producing an immune response to tumor cells. The results would be more comprehensive if supported by relevant data.

The reviewer has interpreted the results absolutely right, the allogenic and autologous MART1-specific TILs cells while reactive in vitro against UM22, cannot kill this tumor either in a subcutaneous or liver metastases model. We hypothesize this has to do with an immune exclusion phenotype and show weak immunohistochemistry that suggest this. We hope the addition of more UM1 data can be viewed as supportive of tumor-reactivity also in vivo.

In addition, the authors in that study used previously frozen biopsy samples for TCR-seq, which may be associated with low-quality sequencing data, high risk of outcome indicators, and unfriendly access to immune cell information. The existence of these problems and the reliability of the results should be considered. If special processing of TCR-seq data from frozen samples was performed, this should also be accounted for.

We agree with the reviewers and acknowledge we never anticipated the development of single-cell sequencing techniques when we started biobank 2013. We performed dead cell removal before the 10x Genomics experiment. We have also done extensive quality controls and believe that the data from the biopsies should be viewed as a whole and that quantitative intra-patient comparisons cannot be done.

Reviewer #2 (Public Review):

Summary:

The study's goal is to characterize and validate tumor-reactive T cells in liver metastases of uveal melanoma (UM), which could contribute to enhancing immunotherapy for these patients. The authors used single-cell RNA and TCR sequencing to find potential tumor-reactive T cells and then used patientderived xenograft (PDX) models and tumor sphere cultures for functional analysis. They discovered that tumor-reactive T cells exist in activated/exhausted T cell subsets and in cytotoxic effector cells. Functional experiments with isolated TILs show that they are capable of killing UM cells in vivo and ex vivo.

Strengths:

The study highlights the potential of using single-cell sequencing and functional analysis to identify T cells that can be useful for cell therapy and marker selection in UM treatment. This is important and novel as conventional immune checkpoint therapies are not highly effective in treating UM. Additionally, the study's strength lies in its validation of findings through functional assays, which underscores the clinical relevance of the research.

We thank the reviewer for these kind words about our work.

Weaknesses:

The manuscript may pose challenges for individuals with limited knowledge of single-cell analysis and immunology markers, making it less accessible to a broader audience.

The first draft of the manuscript (excluding methods) was written by a person (J.A.N) who is not a bioinformatician. It has been corrected to include the correct nomenclature where applicable but overall it is written with the aim to be understandable. We have made an additional effort in this version.

Reviewer #1 (Recommendations For The Authors):

(1) Firstly, the authors should provide high-resolution pictures to ensure readability for readers.

We have converted to pdf ourselves and that improved resolution. We are happy to provide high-resolution to the office if needed for the printing.

(2) Furthermore, some parts of the article are more colloquial, and the authors should consider the logic and academic nature of the overall writing of the article. For example, authors should double-check whether the relevant expressions in the results are correct. For example, 'TCR' in the fourth part of the results should be 'TRLs'.

We thank the reviewer for the recommendations and have gone through the manuscript.

(3) Moreover, UM22 is described several times in the results as a metastatic UM and should be clearly defined in the methodology.

The UM22 and UM1 samples are described in-depth in Karlsson et al., Nature Communications, 2020, a paper that is cited in the beginning of Results as part of the narrative. The current work can be viewed as an extension of that work.

(4) Finally, it is recommended that authors describe a part of the results in full before citing the corresponding picture, otherwise, it will lead to confusion among readers.

We have made an effort in the revised version to describe the new data in more detail.

Reviewer #2 (Recommendations For The Authors):

The manuscript is very interesting and important to understanding key aspects of uveal melanoma immune profile and functionality. However, in my opinion, there are a few aspects that could be addressed.

- The manuscript lacks comprehensive details about the samples used, such as their disease progression, response to treatment, or any relevant information that could shed light on potential differences between samples. It would be valuable to know whether these samples were collected before any systemic treatment or if any of the patients underwent immunotherapy post-sample collection, along with the outcomes of such treatments. Providing this information would enrich the manuscript and provide a more holistic view of the research.

We thank the reviewer for the recommendation and have included a new Supplementary table 7 with information about the samples. We have also pasted in individual samples’ contribution to the UMAP to add further holistic view.

- The results presented and discussed in the manuscript seem to indicate that there were no significant differences across the various samples, including comparisons between lymph-node and liver metastases. However, this lack of variation or the reasons for not discussing any observed differences should be clarified. If there are distinctions between the samples, it would be beneficial to discuss these findings in the manuscript.

We thank the reviewer for the recommendation. Whereas 14 samples are many for a uveal melanoma study it is not really powered to do intra-patient comparisons.

- The manuscript may pose difficulties for individuals with limited knowledge of single-cell analysis and immunology markers, potentially limiting its accessibility. To make the research more inclusive, the authors might consider presenting the technical aspects of their work in a less descriptive manner and providing explanations for those less familiar with the technology. This would help a broader audience grasp the significance of the study's findings.

The manuscript is from a multidisciplinary team where all have read and commented. The draft was written by a tumor biologist and edited by a bioinformatician for accuracy. We honestly think it is more understandable than most studies in this bioinformatics era. But we have tried to describe the new data in an easier way.

No competing interests declared.

Data curation, Software, Formal analysis, Validation, Investigation, Visualization, Writing – original draft, Writing – review and editing, Conducted all bioinformatic analysis except for the UM1 liver metastasis model data, generated figures and wrote several sections of the paper.

Formal analysis, Validation, Investigation, Methodology, Writing – review and editing, Conducted functional experiments, generated figures and wrote some sections of the paper.

Resources, Writing – review and editing, Coordinated patient consents, surgeries and collection of liver metastases of uveal melanoma.

Data curation, Formal analysis, Validation, Investigation, Writing – review and editing, Conducted bioinformatic analysis for the liver metastatic model of the UM1 sample.

Methodology, Writing – review and editing, Conducted single-cell sequencing.

Methodology, Writing – review and editing, Performed cell sorting and flow analyses.

Investigation, Methodology, Writing – review and editing, Generated PDX models and developed new transplantation and imaging techniques.

Resources, Supervision, Writing – review and editing, Managed single-cell sequencing workflow.

Resources, Writing – review and editing, Performed patient recruitment and supervision.

Conceptualization, Resources, Supervision, Investigation, Methodology, Project administration, Writing – review and editing, Supervised, managed the biobank and sequencing, generated TILs, managed, conceived and visualized data from animal experiments.

Conceptualization, Data curation, Supervision, Funding acquisition, Investigation, Visualization, Writing – original draft, Project administration, Writing – review and editing, Conceived and supervised the study, wrote the paper and generated figures.

Human subjects: According to the ethical approval of the Regional Ethical Review Board in Gothenburg (#289-12, #44-18 and #144-13), the patients provided oral and written information with signed informed consent agreements.

Animal experiments were conducted in conformity with E.U. directive 2010/63 (Regional Animal Ethics Committee of Gothenburg approvals no. 2014-36 and 2018-1183).
==== Refs
References

Algazi AP Tsai KK Shoushtari AN Munhoz RR Eroglu Z Piulats JM Ott PA Johnson DB Hwang J Daud AI Sosman JA Carvajal RD Chmielowski B Postow MA Weber JS Sullivan RJ 2016 Clinical outcomes in metastatic uveal melanoma treated with PD-1 and PD-L1 antibodies Cancer 122 3344 3353 10.1002/cncr.30258 27533448
Al Hity G Yang F Campillo-Funollet E Greenstein AE Hunt H Mampay M Intabli H Falcinelli M Madzvamuse A Venkataraman C Flint MS 2021 An integrated framework for quantifying immune-tumour interactions in a 3D co-culture model Communications Biology 4 781 10.1038/s42003-021-02296-7 34168276
Ben-Shabat I Belgrano V Ny L Nilsson J Lindnér P Olofsson Bagge R 2016 Long-Term follow-up evaluation of 68 Patients with uveal melanoma liver metastases treated with isolated hepatic perfusion Annals of Surgical Oncology 23 1327 1334 10.1245/s10434-015-4982-5 26628434
Blighe K 2018 EnhancedVolcano v.1.20.0 GitHub https://github.com/kevinblighe/EnhancedVolcano
Buder K Gesierich A Gelbrich G Goebeler M 2013 Systemic treatment of metastatic uveal melanoma: Review of literature and future perspectives Cancer Medicine 2 674 686 10.1002/cam4.133 24403233
Buonomo EL Mei S Guinn SR Leo IR Peluso MJ Nolan M-A Schildberg FA Zhao L Lian C Xu S Misdraji J Kharchenko PV Sharpe AH 2022 Liver stromal cells restrict macrophage maturation and stromal IL-6 limits the differentiation of cirrhosis-linked macrophages Journal of Hepatology 76 1127 1137 10.1016/j.jhep.2021.12.036 35074474
Carita G Némati F Decaudin D 2015 Uveal melanoma patient-derived xenografts Ocular Oncology and Pathology 1 161 169 10.1159/000370154 27172261
Carvajal RD Schwartz GK Tezel T Marr B Francis JH Nathan PD 2017 Metastatic disease from uveal melanoma: treatment options and future prospects The British Journal of Ophthalmology 101 38 44 10.1136/bjophthalmol-2016-309034 27574175
Chandran SS Somerville RPT Yang JC Sherry RM Klebanoff CA Goff SL Wunderlich JR Danforth DN Zlott D Paria BC Sabesan AC Srivastava AK Xi L Pham TH Raffeld M White DE Toomey MA Rosenberg SA Kammula US 2017 Treatment of metastatic uveal melanoma with adoptive transfer of tumour-infiltrating lymphocytes: A single-centre, two-stage, single-arm, phase 2 study The Lancet. Oncology 18 792 802 10.1016/S1470-2045(17)30251-6 28395880
Cheloni S Hillje R Luzi L Pelicci PG Gatti E 2021 XenoCell: Classification of cellular barcodes in single cell experiments from xenograft samples BMC Medical Genomics 14 34 10.1186/s12920-021-00872-8 33514375
Chronister WD Crinklaw A Mahajan S Vita R Koşaloğlu-Yalçın Z Yan Z Greenbaum JA Jessen LE Nielsen M Christley S Cowell LG Sette A Peters B 2021 TCRMatch: Predicting T-cell receptor specificity based on sequence similarity to previously characterized receptors Frontiers in Immunology 12 640725 10.3389/fimmu.2021.640725 33777034
Clarke J Panwar B Madrigal A Singh D Gujar R Wood O Chee SJ Eschweiler S King EV Awad AS Hanley CJ McCann KJ Bhattacharyya S Woo E Alzetani A Seumois G Thomas GJ Ganesan A-P Friedmann PS Sanchez-Elsner T Ay F Ottensmeier CH Vijayanand P 2019 Single-cell transcriptomic analysis of tissue-resident memory T cells in human lung cancer The Journal of Experimental Medicine 216 2128 2149 10.1084/jem.20190249 31227543
Clémenceau B Vivien R Berthomé M Robillard N Garand R Gallot G Vollant S Vié H 2008 Effector memory alphabeta T lymphocytes can express FcgammaRIIIa and mediate antibody-dependent cellular cytotoxicity Journal of Immunology 180 5327 5334 10.4049/jimmunol.180.8.5327 18390714
Collin M Bigley V 2018 Human dendritic cell subsets: an update Immunology 154 3 20 10.1111/imm.12888 29313948
Conway T Wazny J Bromage A Tymms M Sooraj D Williams ED Beresford-Smith B 2012 Xenome--a tool for classifying reads from xenograft samples Bioinformatics 28 i172 i178 10.1093/bioinformatics/bts236 22689758
Crespo J Sun H Welling TH Tian Z Zou W 2013 T cell anergy, exhaustion, senescence, and stemness in the tumor microenvironment Current Opinion in Immunology 25 214 221 10.1016/j.coi.2012.12.003 23298609
Çuburu N Bialkowski L Pontejo SM Sethi SK Bell ATF Kim R Thompson CD Lowy DR Schiller JT 2022 Harnessing anti-cytomegalovirus immunity for local immunotherapy against solid tumors PNAS 119 e2116738119 10.1073/pnas.2116738119 35749366
Damato BE Dukes J Goodall H Carvajal RD 2019 Tebentafusp: T cell redirection for the treatment of metastatic uveal melanoma Cancers 11 971 10.3390/cancers11070971 31336704
Dudley ME Wunderlich JR Robbins PF Yang JC Hwu P Schwartzentruber DJ Topalian SL Sherry R Restifo NP Hubicki AM Robinson MR Raffeld M Duray P Seipp CA Rogers-Freezer L Morton KE Mavroukakis SA White DE Rosenberg SA 2002 Cancer regression and autoimmunity in patients after clonal repopulation with antitumor lymphocytes Science 298 850 854 10.1126/science.1076514 12242449
Duhen T Duhen R Montler R Moses J Moudgil T de Miranda NF Goodall CP Blair TC Fox BA McDermott JE Chang S-C Grunkemeier G Leidner R Bell RB Weinberg AD 2018 Co-expression of CD39 and CD103 identifies tumor-reactive CD8 T cells in human solid tumors Nature Communications 9 2724 10.1038/s41467-018-05072-0 30006565
Durante MA Rodriguez DA Kurtenbach S Kuznetsov JN Sanchez MI Decatur CL Snyder H Feun LG Livingstone AS Harbour JW 2020 Single-cell analysis reveals new evolutionary complexity in uveal melanoma Nature Communications 11 496 10.1038/s41467-019-14256-1 31980621
Einarsdottir BO Bagge RO Bhadury J Jespersen H Mattsson J Nilsson LM Truvé K López MD Naredi P Nilsson O Stierner U Ny L Nilsson JA 2014 Melanoma patient-derived xenografts accurately model the disease and develop fast enough to guide treatment decisions Oncotarget 5 9609 9618 10.18632/oncotarget.2445 25228592
Ewels P 2022 The nf-core framework for community-curated bioinformatics pipelines 0.1 Zenodo 10.5281/zenodo.7116979
Fleming SJ Chaffin MD Arduini A Akkad AD Banks E Marioni JC Philippakis AA Ellinor PT Babadi M 2022 Unsupervised removal of systematic background noise from droplet-based single-cell experiments using CellBender bioRxiv 10.1101/791699
Forsberg EMV Lindberg MF Jespersen H Alsén S Bagge RO Donia M Svane IM Nilsson O Ny L Nilsson LM Nilsson JA 2019 HER2 CAR-T Cells eradicate uveal melanoma and T-cell therapy-resistant human melanoma in IL2 transgenic NOD/SCID IL2 receptor knockout mice Cancer Research 79 899 904 10.1158/0008-5472.CAN-18-3158 30622115
Fu Y Xiao W Mao Y 2022 Recent advances and challenges in uveal melanoma immunotherapy Cancers 14 3094 10.3390/cancers14133094 35804863
Ganesan A-P Clarke J Wood O Garrido-Martin EM Chee SJ Mellows T Samaniego-Castruita D Singh D Seumois G Alzetani A Woo E Friedmann PS King EV Thomas GJ Sanchez-Elsner T Vijayanand P Ottensmeier CH 2017 Tissue-resident memory features are linked to the magnitude of cytotoxic T cell responses in human lung cancer Nature Immunology 18 940 950 10.1038/ni.3775 28628092
Gerlach C Moseman EA Loughhead SM Alvarez D Zwijnenburg AJ Waanders L Garg R de la Torre JC von Andrian UH 2016 The chemokine receptor CX3CR1 defines three antigen-experienced CD8 T cell subsets with distinct roles in immune surveillance and homeostasis Immunity 45 1270 1284 10.1016/j.immuni.2016.10.018 27939671
Germain PL 2024 scDblFinder 1.18.0 BioConductor https://bioconductor.org/packages/release/bioc/vignettes/scDblFinder/inst/doc/scDblFinder.html
Gopal S Kwon S-J Ku B Lee DW Kim J Dordick JS 2021 3D tumor spheroid microarray for high-throughput, high-content natural killer cell-mediated cytotoxicity Communications Biology 4 893 10.1038/s42003-021-02417-2 34290356
Gowthaman R Pierce BG 2019 TCR3d: The T cell receptor structural repertoire database Bioinformatics 35 5323 5325 10.1093/bioinformatics/btz517 31240309
Gros A Robbins PF Yao X Li YF Turcotte S Tran E Wunderlich JR Mixon A Farid S Dudley ME Hanada KI Almeida JR Darko S Douek DC Yang JC Rosenberg SA 2014 PD-1 identifies the patient-specific CD8 The Journal of Clinical Investigation 124 2246 2259 10.1172/JCI73639 24667641
Haghverdi L Lun ATL Morgan MD Marioni JC 2018 Batch effects in single-cell RNA-sequencing data are corrected by matching mutual nearest neighbors Nature Biotechnology 36 421 427 10.1038/nbt.4091 29608177
Hippen AA Falco MM Weber LM Erkan EP Zhang K Doherty JA Vähärautio A Greene CS Hicks SC 2021 miQC: An adaptive probabilistic framework for quality control of single-cell RNA-sequencing data PLOS Computational Biology 17 e1009290 10.1371/journal.pcbi.1009290 34428202
Hu ZI Link VM Lima-Junior DS Delaleu J Bouladoux N Han S-J Collins N Belkaid Y 2022 Immune checkpoint inhibitors unleash pathogenic immune responses against the microbiota PNAS 119 e2200348119 10.1073/pnas.2200348119 35727974
Huang H Wang C Rubelt F Scriba TJ Davis MM 2020 Analyzing the Mycobacterium tuberculosis immune response by T-cell receptor clustering with GLIPH2 and genome-wide antigen screening Nature Biotechnology 38 1194 1202 10.1038/s41587-020-0505-4 32341563
Huuhtanen J Chen L Jokinen E Kasanen H Lönnberg T Kreutzman A Peltola K Hernberg M Wang C Yee C Lähdesmäki H Davis MM Mustjoki S 2022 Evolution and modulation of antigen-specific T cell responses in melanoma patients Nature Communications 13 5988 10.1038/s41467-022-33720-z 36220826
Huuhtanen J 2023 Single-cell characterization of anti-LAG3+anti-PD1 treatment in melanoma patients The Journal of Clinical Investigation 1 e164809 10.1172/JCI164809
Inozume T Hanada K-I Wang QJ Ahmadzadeh M Wunderlich JR Rosenberg SA Yang JC 2010 Selection of CD8+PD-1+ lymphocytes in fresh human melanomas enriches for tumor-reactive T cells Journal of Immunotherapy 33 956 964 10.1097/CJI.0b013e3181fad2b0 20948441
Jager MJ Shields CL Cebulla CM Abdel-Rahman MH Grossniklaus HE Stern M-H Carvajal RD Belfort RN Jia R Shields JA Damato BE 2020 Uveal melanoma Nature Reviews. Disease Primers 6 24 10.1038/s41572-020-0158-0 32273508
Jespersen H Lindberg MF Donia M Söderberg EMV Andersen R Keller U Ny L Svane IM Nilsson LM Nilsson JA 2017 Clinical responses to adoptive T-cell transfer can be modeled in an autologous immune-humanized mouse model Nature Communications 8 707 10.1038/s41467-017-00786-z 28955032
Jespersen H Olofsson Bagge R Ullenhag G Carneiro A Helgadottir H Ljuslinder I Levin M All-Eriksson C Andersson B Stierner U Nilsson LM Nilsson JA Ny L 2019 Concomitant use of pembrolizumab and entinostat in adult patients with metastatic uveal melanoma (PEMDAC study): Protocol for a multicenter phase II open label study BMC Cancer 19 415 10.1186/s12885-019-5623-3 31046743
Johnson DB Bao R Ancell KK Daniels AB Wallace D Sosman JA Luke JJ 2019 Response to Anti-PD-1 in uveal melanoma without high-volume liver metastasis Journal of the National Comprehensive Cancer Network 17 114 117 10.6004/jnccn.2018.7070 30787124
Karlsson J Nilsson LM Mitra S Alsén S Shelke GV Sah VR Forsberg EMV Stierner U All-Eriksson C Einarsdottir B Jespersen H Ny L Lindnér P Larsson E Olofsson Bagge R Nilsson JA 2020 Molecular profiling of driver events in metastatic uveal melanoma Nature Communications 11 1894 10.1038/s41467-020-15606-0 32313009
Karlsson J 2022 Xenocell_nextflow swh:1:rev:bdf2a305d8482a7ba5d9ffaf173a6823f78f8eaa Software Heritage https://archive.softwareheritage.org/swh:1:dir:f5b2ae98e976d8491223e3fe9a90319e73d6f5c2;origin=https://github.com/jowkar/xenocell_nextflow;visit=swh:1:snp:cf8e682450fb736fd559e1adbc1e4a9dce518e1b;anchor=swh:1:rev:bdf2a305d8482a7ba5d9ffaf173a6823f78f8eaa
Khoja L Atenafu EG Suciu S Leyvraz S Sato T Marshall E Keilholz U Zimmer L Patel SP Piperno-Neumann S Piulats J Kivelä TT Pfoehler C Bhatia S Huppert P Van Iersel LBJ De Vries IJM Penel N Vogl T Cheng T Fiorentini G Mouriaux F Tarhini A Patel PM Carvajal R Joshua AM 2019 Meta-analysis in metastatic uveal melanoma to determine progression free and overall survival benchmarks: an international rare cancers initiative (IRCI) ocular melanoma study Annals of Oncology 30 1370 1380 10.1093/annonc/mdz176 31150059
Li H van der Leun AM Yofe I Lubling Y Gelbard-Solodkin D van Akkooi ACJ van den Braber M Rozeman EA Haanen JBAG Blank CU Horlings HM David E Baran Y Bercovich A Lifshitz A Schumacher TN Tanay A Amit I 2019 Dysfunctional CD8 T cells form a proliferative, dynamically regulated compartment within human melanoma Cell 176 775 789 10.1016/j.cell.2018.11.043 30595452
Lin W Beasley AB Ardakani NM Denisenko E Calapre L Jones M Wood BA Warburton L Forrest ARR Gray ES 2021 Intra- and intertumoral heterogeneity of liver metastases in a patient with uveal melanoma revealed by single-cell RNA sequencing Cold Spring Harbor Molecular Case Studies 7 a006111 10.1101/mcs.a006111 34470851
Love MI Huber W Anders S 2014 Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 Genome Biology 15 1 34 10.1186/s13059-014-0550-8
Ma D Hernandez GA Lefebvre AEYT Alshetaiwi H Blake K Dave KR Rauf M Williams JW Davis RT Evans KT Longworth A Masoud MYG Lee R Edwards RA Digman MA Kessenbrock K Lawson DA 2021 Patient-derived xenograft culture-transplant system for investigation of human breast cancer metastasis Communications Biology 4 1268 10.1038/s42003-021-02596-y 34741115
Malone AF 2021 Monocytes and Macrophages in kidney transplantation and insights from single cell RNA-seq studies Kidney360 2 1654 1659 10.34067/KID.0003842021 35372970
Middleton MR McAlpine C Woodcock VK Corrie P Infante JR Steven NM Evans TRJ Anthoney A Shoushtari AN Hamid O Gupta A Vardeu A Leach E Naidoo R Stanhope S Lewis S Hurst J O’Kelly I Sznol M 2020 Tebentafusp, A TCR/Anti-CD3 bispecific fusion protein targeting gp100, potently activated antitumor immune responses in patients with metastatic melanoma Clinical Cancer Research 26 5869 5878 10.1158/1078-0432.CCR-20-1247 32816891
Najjar YG Navrazhina K Ding F Bhatia R Tsai K Abbate K Durden B Eroglu Z Bhatia S Park S Chowdhary A Chandra S Kennedy J Puzanov I Ernstoff M Vachhani P Drabick J Singh A Xu T Yang J Carvajal R Manson D Kirkwood JM Cohen J Sullivan R Johnson D Funchain P Shoushtari A 2020 Ipilimumab plus nivolumab for patients with metastatic uveal melanoma: a multicenter, retrospective study Journal for Immunotherapy of Cancer 8 e000331 10.1136/jitc-2019-000331 32581057
Nathan P Hassel JC Rutkowski P Baurain J-F Butler MO Schlaak M Sullivan RJ Ochsenreither S Dummer R Kirkwood JM Joshua AM Sacco JJ Shoushtari AN Orloff M Piulats JM Milhem M Salama AKS Curti B Demidov L Gastaud L Mauch C Yushak M Carvajal RD Hamid O Abdullah SE Holland C Goodall H Piperno-Neumann S IMCgp100-202 Investigators 2021 Overall survival benefit with tebentafusp in metastatic uveal melanoma The New England Journal of Medicine 385 1196 1206 10.1056/NEJMoa2103485 34551229
Nilsson LM Vilhav C Karlsson JW Fagman JB Giglio D Engström CE Naredi P Nilsson JA 2022 Genetics and therapeutic responses to tumor-infiltrating lymphocyte therapy of pancreatic cancer patient-derived xenograft models Gastro Hep Advances 1 1037 1048 10.1016/j.gastha.2022.07.006 39131259
Ny L Rizzo LY Belgrano V Karlsson J Jespersen H Carstam L Bagge RO Nilsson LM Nilsson JA 2020 Supporting clinical decision making in advanced melanoma by preclinical testing in personalized immune-humanized xenograft mouse models Annals of Oncology 31 266 273 10.1016/j.annonc.2019.11.002 31959343
Ny L Jespersen H Karlsson J Alsén S Filges S All-Eriksson C Andersson B Carneiro A Helgadottir H Levin M Ljuslinder I Olofsson Bagge R Sah VR Stierner U Ståhlberg A Ullenhag G Nilsson LM Nilsson JA 2021 The PEMDAC phase 2 study of pembrolizumab and entinostat in patients with metastatic uveal melanoma Nature Communications 12 5155 10.1038/s41467-021-25332-w 34453044
Olofsson R Ny L Eilard MS Rizell M Cahlin C Stierner U Lönn U Hansson J Ljuslinder I Lundgren L Ullenhag G Kiilgaard JF Nilsson J Lindnér P 2014 Isolated hepatic perfusion as a treatment for uveal melanoma liver metastases (the SCANDIUM trial): study protocol for a randomized controlled trial Trials 15 317 10.1186/1745-6215-15-317 25106493
Olofsson Bagge R Nelson A Shafazand A All-Eriksson C Cahlin C Elander N Helgadottir H Kiilgaard JF Kinhult S Ljuslinder I Mattsson J Rizell M Sternby Eilard M Ullenhag GJ Nilsson JA Ny L Lindnér P 2023 Isolated hepatic perfusion with melphalan for patients with isolated uveal melanoma liver metastases: A multicenter, randomized, open-label, phase III trial (the SCANDIUM Trial) Journal of Clinical Oncology 41 3042 3050 10.1200/JCO.22.01705 36940407
Olofsson Bagge R Nelson A Shafazand A All-Eriksson C Cahlin C Elander N Gustavsson A Helgadottir H Kiilgaard JF Kinhult S Ljuslinder I Mattsson J Rizell M Sternby Eilard M Ullenhag GJ Nilsson JA Ny L Lindnér P 2024 Survival and Quality of Life after Isolated Hepatic Perfusion with Melphalan as a Treatment for Uveal Melanoma Liver Metastases - Final Results from the Phase III Randomized Controlled Trial SCANDIUM Annals of Surgery 10.1097/SLA.0000000000006255 38420778
Pandiani C Strub T Nottet N Cheli Y Gambi G Bille K Husser C Dalmasso M Béranger G Lassalle S Magnone V Pédeutour F Irondelle M Maschi C Nahon-Estève S Martel A Caujolle J-P Hofman P LeBrigand K Davidson I Baillif S Barbry P Ballotti R Bertolotto C 2021 Single-cell RNA sequencing reveals intratumoral heterogeneity in primary uveal melanomas and identifies HES6 as a driver of the metastatic disease Cell Death and Differentiation 28 1990 2000 10.1038/s41418-020-00730-7 33462406
Patton EE Mueller KL Adams DJ Anandasabapathy N Aplin AE Bertolotto C Bosenberg M Ceol CJ Burd CE Chi P Herlyn M Holmen SL Karreth FA Kaufman CK Khan S Kobold S Leucci E Levy C Lombard DB Lund AW Marie KL Marine J-C Marais R McMahon M Robles-Espinoza CD Ronai ZA Samuels Y Soengas MS Villanueva J Weeraratna AT White RM Yeh I Zhu J Zon LI Hurlbert MS Merlino G 2021 Melanoma models for the next generation of therapies Cancer Cell 39 610 631 10.1016/j.ccell.2021.01.011 33545064
Pelster MS Gruschkus SK Bassett R Gombos DS Shephard M Posada L Glover MS Simien R Diab A Hwu P Carter BW Patel SP 2021 Nivolumab and ipilimumab in metastatic uveal melanoma: Results From a single-arm phase II study Journal of Clinical Oncology 39 599 607 10.1200/JCO.20.00605 33125309
Piulats JM Espinosa E de la Cruz Merino L Varela M Alonso Carrión L Martín-Algarra S López Castro R Curiel T Rodríguez-Abreu D Redrado M Gomà M Rullán AJ Calvo González A Berrocal-Jaime A 2021 Nivolumab plus ipilimumab for treatment-naïve metastatic uveal melanoma: An open-label, multicenter, phase II trial by the Spanish multidisciplinary melanoma group (GEM-1402) Journal of Clinical Oncology 39 586 598 10.1200/JCO.20.00550 33417511
Rosato PC Wijeyesinghe S Stolley JM Nelson CE Davis RL Manlove LS Pennell CA Blazar BR Chen CC Geller MA Vezys V Masopust D 2019 Virus-specific memory T cells populate tumors and can be repurposed for tumor immunotherapy Nature Communications 10 567 10.1038/s41467-019-08534-1 30718505
Rothermel LD Sabesan AC Stephens DJ Chandran SS Paria BC Srivastava AK Somerville R Wunderlich JR Lee C-CR Xi L Pham TH Raffeld M Jailwala P Kasoji M Kammula US 2016 Identification of an immunogenic subset of metastatic uveal melanoma Clinical Cancer Research 22 2237 2249 10.1158/1078-0432.CCR-15-2294 26712692
Satija R Farrell JA Gennert D Schier AF Regev A 2015 Spatial reconstruction of single-cell gene expression data Nature Biotechnology 33 495 502 10.1038/nbt.3192 25867923
Satija R 2024 Seurat v5 Satija Lab https://satijalab.org/seurat/
Schreeder DM Pan J Li FJ Vivier E Davis RS 2008 FCRL6 distinguishes mature cytotoxic lymphocytes and is upregulated in patients with B-cell chronic lymphocytic leukemia European Journal of Immunology 38 3159 3166 10.1002/eji.200838516 18991291
Seiter S Monsurro V Nielsen M-B Wang E Provenzano M Wunderlich JR Rosenberg SA Marincola FM 2002 Frequency of MART-1/MelanA and gp100/PMel17-specific T cells in tumor metastases and cultured tumor-infiltrating lymphocytes Journal of Immunotherapy 25 252 263 10.1097/00002371-200205000-00008 12000867
Sheridan R 2020 djvdj v0.0.0.9000 GitHub https://github.com/rnabioco/djvdj
Shugay M Bagaev DV Zvyagin IV Vroomans RM Crawford JC Dolton G Komech EA Sycheva AL Koneva AE Egorov ES Eliseev AV Van Dyk E Dash P Attaf M Rius C Ladell K McLaren JE Matthews KK Clemens EB Douek DC Luciani F van Baarle D Kedzierska K Kesmir C Thomas PG Price DA Sewell AK Chudakov DM 2018 VDJdb: a curated database of T-cell receptor sequences with known antigen specificity Nucleic Acids Research 46 D419 D427 10.1093/nar/gkx760 28977646
Simoni Y Becht E Fehlings M Loh CY Koo S-L Teng KWW Yeong JPS Nahar R Zhang T Kared H Duan K Ang N Poidinger M Lee YY Larbi A Khng AJ Tan E Fu C Mathew R Teo M Lim WT Toh CK Ong B-H Koh T Hillmer AM Takano A Lim TKH Tan EH Zhai W Tan DSW Tan IB Newell EW 2018 Bystander CD8+ T cells are abundant and phenotypically distinct in human tumour infiltrates Nature 557 575 579 10.1038/s41586-018-0130-2 29769722
Sjoberg D 2021 ggsankey v.0.0.99999 GitHub https://github.com/davidsjoberg/ggsankey
Squair JW Gautier M Kathe C Anderson MA James ND Hutson TH Hudelle R Qaiser T Matson KJE Barraud Q Levine AJ La Manno G Skinnider MA Courtine G 2021 Confronting false discoveries in single-cell differential expression Nature Communications 12 5692 10.1038/s41467-021-25960-2 34584091
Sugase T Lam BQ Danielson M Terai M Aplin AE Gutkind JS Sato T 2020 Development and optimization of orthotopic liver metastasis xenograft mouse models in uveal melanoma Journal of Translational Medicine 18 208 10.1186/s12967-020-02377-x 32434572
Tickle T Tirosh I Georgescu C Brown M Haas B 2019 InferCNV of the trinity CTAT project. preprint at 124eec0 Github https://github.com/broadinstitute/inferCNV
Tickotsky N Sagiv T Prilusky J Shifrut E Friedman N 2017 McPAS-TCR: A manually curated catalogue of pathology-associated T cell receptor sequences Bioinformatics 33 2924 2929 10.1093/bioinformatics/btx286 28481982
van den Berg JH Heemskerk B van Rooij N Gomez-Eerland R Michels S van Zon M de Boer R Bakker NAM Jorritsma-Smit A van Buuren MM Kvistborg P Spits H Schotte R Mallo H Karger M van der Hage JA Wouters MWJM Pronk LM Geukes Foppen MH Blank CU Beijnen JH Nuijen B Schumacher TN Haanen JBAG 2020 Tumor infiltrating lymphocytes (TIL) therapy in metastatic melanoma: boosting of neoantigen-specific T cell reactivity and long-term follow-up Journal for Immunotherapy of Cancer 8 e000848 10.1136/jitc-2020-000848 32753545
van der Kooij MK Speetjens FM van der Burg SH Kapiteijn E 2019 Uveal versus cutaneous melanoma; same origin, very distinct tumor types Cancers 11 845 10.3390/cancers11060845 31248118
Vita R Mahajan S Overton JA Dhanda SK Martini S Cantrell JR Wheeler DK Sette A Peters B 2019 The immune epitope database (IEDB): 2018 update Nucleic Acids Research 47 D339 D343 10.1093/nar/gky1006 30357391
Webb JR Wick DA Nielsen JS Tran E Milne K McMurtrie E Nelson BH 2010 Profound elevation of CD8+ T cells expressing the intraepithelial lymphocyte marker CD103 (alphaE/beta7 Integrin) in high-grade serous ovarian cancer Gynecologic Oncology 118 228 236 10.1016/j.ygyno.2010.05.016 20541243
Weiswald LB Bellet D Dangles-Marie V 2015 Spherical cancer models in tumor biology Neoplasia 17 1 15 10.1016/j.neo.2014.12.004 25622895
Ye Q Song DG Poussin M Yamamoto T Best A Li C Coukos G Powell Jr DJ 2014 CD137 accurately identifies and enriches for naturally occurring tumor-reactive T cells in tumor Clinical Cancer Research 20 44 55 10.1158/1078-0432.CCR-13-0945 24045181
Yee C Thompson JA Byrd D Riddell SR Roche P Celis E Greenberg PD 2002 Adoptive T cell therapy using antigen-specific CD8+ T cell clones for the treatment of patients with metastatic melanoma: in vivo persistence, migration, and antitumor effect of transferred T cells PNAS 99 16168 16173 10.1073/pnas.242600099 12427970
Zhang J He T Xue L Guo H 2021 Senescent T cells: a potential biomarker and target for cancer therapy EBioMedicine 68 103409 10.1016/j.ebiom.2021.103409 34049248
Zhao Y Shao Q Peng G 2020 Exhaustion and senescence: two crucial dysfunctional states of T cells in the tumor microenvironment Cellular & Molecular Immunology 17 27 35 10.1038/s41423-019-0344-8 31853000
Zheng C Fass JN Shih Y-P Gunderson AJ Sanjuan Silva N Huang H Bernard BM Rajamanickam V Slagel J Bifulco CB Piening B Newell PHA Hansen PD Tran E 2022 Transcriptomic profiles of neoantigen-reactive T cells in human gastrointestinal cancers Cancer Cell 40 410 423 10.1016/j.ccell.2022.03.005 35413272
