Article Text
Abstract
T-cell receptor sequencing (TCRseq) enables tracking of T-cell clonotypes recognizing the same antigen over time and across biological compartments. TCRseq has been used to test if cross-reactive antitumor T cells are responsible for development of immune-related adverse events (irAEs) following immune checkpoint blockade. Prior studies have interpreted T-cell clones shared among the tumor and irAE as evidence supporting this, but interpretations of these findings are challenging, given the constraints of TCRseq. Here we capitalize on a rare opportunity to understand the impact of potential confounders, such as sample size, tissue compartment, and collection batch/timepoint, on the relative proportion of shared T-cell clones between an irAE and tumor specimens. TCRseq was performed on tumor-involved and -uninvolved tissues, including an irAE, that were obtained throughout disease progression and at the time of rapid autopsy from a patient with renal cell carcinoma treated with programmed death-1 (PD-1) blockade. Our analyses show significant effects of these confounders on our ability to understand T-cell receptor overlap, and we present mitigation strategies and study design recommendations to reduce these errors. Implementation of these strategies will enable more rigorous TCRseq-based studies of immune responses in human tissues, particularly as they relate to antitumor T-cell cross-reactivity in irAEs following checkpoint blockade.
- immunotherapy
- receptors
- immunologic
- t-lymphocytes
- immunologic techniques
- biostatistics
This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See http://creativecommons.org/licenses/by-nc/4.0/.
Statistics from Altmetric.com
Introduction
PD(L)-1 checkpoint blockade is complicated by the development of immune-related adverse events (irAEs) in 5%–20% of treated patients.1 Severe irAEs have been reported in up to 10% of patients and can result in hospitalization, interruption or discontinuation of therapy, and rarely, death.1 Notably, irAEs increase in prevalence and severity with combination immunotherapy, an approach likely required to improve disappointingly low response rates.2 Prediction, prevention, and treatment of irAEs will require delineation of etiological mechanisms. One hypothesis is that irAEs are the result of cross-reactivity of an antigen-specific antitumor immune response. Supporting this hypothesis, improvements in response rates and survival in patients who develop irAEs have been observed in studies across tumor types.3–6 T-cell repertoire profiling with T-cell receptor sequencing (TCRseq) enables tracking of individual T-cell clones recognizing the same antigen.7 Multiple case series have identified shared T-cell clones in tumor and irAE tissues, thereby providing a foundation for shared antigen specificity between the tumor and irAE.8–10 Unfortunately, these types of analyses may be subject to several sources of confounding that are rarely considered and are often difficult to address in the context of human immune-oncology.
The current study evaluates notable sources of confounding in the analysis of T-cell receptor (TCR) repertoire overlap between tumor-involved and irAE specimens in a patient who developed a refractory irAE (dermatitis) while receiving PD-1 blockade for metastatic renal cell carcinoma (RCC). We demonstrate how these analytical pitfalls could lead to erroneous interpretation of TCRseq data obtained from distinct biological compartments and timepoints within the same patient. These factors have previously been considered as potentiators of confoundment in large-scale genomic datasets11 but have not been evaluated or implemented in TCR studies.
Case report
A woman in her early 70s underwent radical nephrectomy for clear cell RCC followed by systemic therapy, including anti-PD-1. The patient’s clinical course, including development of an irAE in the form of a lichenoid dermatitis (LD),12–16 and biospecimen collection are shown in figure 1. The immune-related LD (irAE) persisted despite two treatment breaks and systemic prednisone (figure 1). Anti-PD-1 therapy was discontinued because of the severity of the irAE. Progressing metastases in the small bowel and a new brain metastasis were confirmed by biopsy. The patient died approximately 2 months after cessation of anti-PD-1 therapy and a rapid autopsy was performed. Findings included metastatic RCC involving the brain (specimen TM1), jejunum (specimen TM2), and mesentery (specimen TM3). Of the three mediastinal lymph nodes (LNs) sampled, one was histologically unremarkable (LN1); one showed multiple large fibrotic nodules (LN2); and one showed multiple small subcapsular fibrotic nodules (LN3). An inflamed seborrheic keratosis (benign skin lesion) was sampled from the skin (SK) as well as uninvolved normal tissues from the left kidney and normal small bowel (NSB).
Supplemental material
To test for possible T-cell cross-reactivity among the tumor and irAE (figure 1D), TCR Vβ CDR3 sequencing was performed on all specimens (online supplemental table S1).17 18 The irAE shared 147 unique clonotypes (4.7%) with the pretreatment primary tumor (Tp) and 118 unique clonotypes (3.7%) with jejunal metastasis (TM2) (online supplemental figure S1A,B). In total, 127 unique T-cell clones present in the irAE (4.0%) were also found in at least one tumor specimen and absent in all healthy, non-lymphoid specimens (online supplemental figure S1C). We next tested if library size (the total number of productive sequencing reads) influences T-cell repertoire overlap among specimens. The number of clones shared with a given specimen was highly correlated with library size, illustrated for the irAE (Spearman’s rho, R=0.7, p=0.031; online supplemental figure S1D, left) and LN2 (R=0.78, p=0.012; online supplemental figure S1D, right). Random subsampling weighted by clonal abundance within each specimen was used to equalize library sizes to eliminate this confounding (online supplemental figure S1E).7 11 Not surprisingly, a strong correlation was observed between library size and the number of unique clonotypes (R=0.93, p<2.2e−16; online supplemental figure S2). The degree of clonal sharing was also correlated with the number of unique clonotypes in each specimen (R=0.85 and p=0.0035 for the irAE, and R=0.92 and p=0.00047 for LN2; online supplemental figure S3A). Using weighted downsampling,7 11 we normalized specimens to the same library size, which eliminated the correlation between the number of unique clonotypes in a specimen and clonal sharing (online supplemental figure S3B). Therefore, we used weighted downsampling for the remainder of our analyses.
Tumor-specific T cells can be detected in paired uninvolved tissue, even when the normal tissue is collected 10–15 cm from the tumor itself.19–21 Likewise, clonotype sharing analyses could be confounded by bypassing viral-specific T cells. We performed a reanalysis of a previously published functional assay20 21 and found that viral-specific T-cell clones showed notable clonotype sharing across multiple tissue compartments, including the tumor, in a patient with non-small cell lung cancer. A similar pattern was previously observed with neoantigen-specific T-cell clones,20 21 suggesting that T cells can traffic across tumor involved/uninvolved compartments regardless of the presence of antigen (online supplemental figure S4A,B). Indeed, after mapping the TCRs from our present study to a public TCR database with annotated antigen specificity (vdjdb, https://vdjdb.cdr3.net/), we found an Epstein-Barr Virus (EBV)-specific clone that was detected in tumor-involved tissues from our study participant (online supplemental figure S4C–D). This emphasizes that clonotype sharing alone is not necessarily associated with biological relevance. Additional abundance measurement and antigen specificity analyses are warranted to assist further interpretation. The SK and LN outliers in online supplemental figure S1D (red arrows) also indicate increased T-cell repertoire sharing between specimens from the same tissue compartment, even when collected at different locations and timepoints (ie, the irAE and SK). Tissue compartment confounding—a greater degree of T-cell repertoire overlap between specimens collected from the same tissue site—is illustrated with pairwise comparisons in online supplemental figure S5A,B. While T-cell repertoire sharing is reduced in samples from different tissue compartments, shared clones are still detected between seemingly unrelated specimens. These could reflect circulating clones at the time of tissue collection and are not necessarily reflective of biologically meaningful clonal sharing, (ie, batch effect confounding). We used the Morisita Overlap Index (MI), which incorporates relative clonal abundance and is not influenced by library size in our dataset (online supplemental figure S6A), to calculate the overlap between the irAE and all other specimens. The relative clonal sharing among all specimens is illustrated in a chord diagram (online supplemental figure S6B), in which the width of the bands is proportional to the MI values. The highest MI was observed between specimens collected from the same batch and from the same tissue compartment (online supplemental figure S6C). These population-level comparisons likely capture a combination of biological and batch effects, which cannot be distinguished in this dataset.
We next tested for evidence of cross-reactive T cells in the tumor and irAE while accounting for the sources of confounding identified previously (figure 1E). Within the limitations of specimen availability, tissue compartment and batch effect were considered in selecting comparator specimens for meaningful analyses. First, the T-cell repertoire overlap between the primary tumor and the metastases was quantified using MI, which demonstrated a significantly higher overlap among the progressing metastases relative to the mediastinal LNs and normal tissues (p=0.044; figure 1F,G). A chord diagram highlights population-level sharing among all tumor specimens (figure 1G), with the greatest TCR repertoire overlap observed between the metastases from the brain (TM1) and mesentery (TM3, MI 0.47). Although batch effect is a potential confounder, this degree of overlap is not observed with the small bowel metastasis (TM2, MI 0.03 and 0.04 with TM1 and TM3, respectively), also collected at autopsy. Since the primary tumor and metastases were collected at different time points and from different tissue sites, batch and tissue compartment effects could not be confounders in this analysis. We recognize that a small sample size limits interpretation of the aforementioned findings, given that TM1 and TM3 have the smallest library sizes (online supplemental table S1), even though they satisfied our criteria for inclusion. Consequently, we focused on the metastasis with the largest library size (TM2, 6945 reads) for additional analyses.
Following library size normalization, 9.1% of unique primary tumor clonotypes are shared with TM2 relative to 5.8% (95% CI 5.2% to 6.8%) and 3.7% (95% CI 3.3% to 4.3%) of Tp clonotypes shared with LN2 and NSB, respectively (figure 1H). Clonal expansion of shared clones in the primary tumor was greatest for those shared with TM2, with shared clones representing 27.1% of total primary tumor reads. The TCR repertoire overlap between the primary tumor and LN2 may suggest an antitumor signature in the mediastinal LNs, a site of radiographical tumor regression. The same approach was used to evaluate TCR repertoire overlap between the irAE and tumor specimens. The irAE repertoire was most similar to TM3, TM1, and the regression site LNs (figure 1I), with intermediate overlap with Tp and the least overlap with the normal control specimens. Although the relative degrees of overlap with the irAE potentially suggest a biologically relevant pattern, the magnitude of population-level TCR repertoire sharing is quite small relative to values observed among other specimens (figure 1J).
Specimen libraries were then normalized to allow direct pairwise comparison of clonal sharing with the irAE. Sharing between the irAE and Tp (largest library size of 12,588 reads) was compared with LN2 and NSB. A similar degree of irAE clonotype sharing was observed in Tp and LN2, which was greater than that observed for NSB (figure 1K). There was no evidence of clonal proliferation in the irAE or Tp (figure 1K). Finally, we evaluated the most abundant clonotypes in each specimen for overlap with the irAE, given that prior studies have implicated the highest-frequency intratumor clonotypes in mediating antitumor immunity.22 There was no enrichment of irAE-shared clones in the tumor relative to the non-tumor specimens (online supplemental figure S7).
Based on observations that antigen specificity may be determined by limited contact sites in the TCR CDR3, we applied the grouping lymphocyte interactions by paratope hotspots 2 (GLIPH2) algorithm23 24 to identify and cluster TCR sequences into possible antigen specificity groups. In order to be included in downstream analyses, clusters had to contain ≥3 unique CDR3s, ≥10 reads for each CDR3, a variable gene beta (vb) score of <0.05, and a length score <0.05. One cluster with significant enrichment in the irAE was identified. Notably, the primary tumor (3.66%) had the highest abundance of T-cell clones in the ‘SSQD’ CDR3 motif cluster (figure 1L), followed by the irAE (2.42%), which were both higher than representation of this motif in non-diseased TCR repertoires from four healthy donors (range: 0.03%–1.67%, online supplemental figure S8).25 Though the human leukocyte antigen (HLA) information is unknown, three of the four healthy donors had common clonotypes shared with the patient in our study, indicating that at least one HLA allele was shared among them. By querying additional published skin/tumor-reactive TCR data, the specific motif SSQD was reported in a T-cell clone recognizing an epitope derived from Maspin, which functions as a tumor suppressor gene in epithelial cells.26 Collectively, this indicates that, though analyses of the total TCR repertoire and a subset of high abundance clones do not show a signature of enriched sharing between the irAE and tumor specimens relative to non-tumor specimens at the clonotype level, more ‘antigen-driven’ approaches may be useful to identify potential specificity clusters, especially when coupled with functional assays to confirm antigen specificity and cross-reactivity between irAEs and tumors.
Discussion
As immune checkpoint blocking agents become first-line and second-line therapies for a growing number of tumor types, we are faced with an increasing number of diverse irAEs that may develop during or after treatment. The association of cutaneous irAEs with clinical benefit in some patients suggests that there may be a common antigen that may underlie both durable antitumor responses and clinically significant irAEs. It is conceivable that T cells with a common TCR could mediate both tumor regression and irAE development and progression, as has been evidenced by prior studies evaluating clonal overlap of TCR clonotypes between tumor and irAE tissues8–10 and that expansion of peripheral blood T-cell clones prior to irAE onset positively correlates with irAE severity during checkpoint blockade treatment.27
The large number and circulating nature of T cells predispose these studies to detecting false positive signals, that is, detection of differential or statistically significant clonal overlap that is not necessarily of pathogenic relevance. Biological differences exacerbate this issue, including variation in T-cell numbers and clonality in different tissue types. In addition, due to differences in sampling, clonotype detection can be limited, particularly for rare/low-frequency clonotypes. The analysis pitfalls and mitigation strategies identified in this study are summarized in table 1, and we present considerations for prospective specimen collection in online supplemental figure S9. Many of these factors are already considered as a standard part of large-scale genomic analyses, but they are not yet routinely applied to immune receptor sequencing datasets and, to date, no studies have demonstrated the differential outcomes when these important sources of confounding are not acknowledged. Strengths of this study include the rare opportunity to analyze the TCR repertoire in the same patient across time, tissue compartments, and disease states, and the ability to compare with published tumor-reactive/skin-reactive TCRs and non-irAE skin TCRs. We recognize that we are limited in our ability to comprehensively dissect all potential sources of confounding owing to limited sample availability. Lastly, the data-driven recommendations made in this study highlight the scientific value of rapid autopsy to answer complex questions using human tissue specimens.
Methods
Case selection
Specimens from the underlying primary tumor and/or metastatic site and from skin affected by the cutaneous irAE were collected from the Johns Hopkins Hospital surgical pathology archives and the Rapid Autopsy program and Franklin Square Hospital. Overall patient response to anti-PD-1 therapy was classified according to Response Evaluation Criteria in Solid Tumors V.1.1.
TCRseq and bioinformatic analysis
DNA extraction from formalin-fixed paraffin-embedded (FFPE)-preserved tumor and skin biopsy specimens was performed using the DNeasy Blood and Tissue Kit (Qiagen). The TCR-B locus was amplified and sequenced using the ImmunoSEQ assay (Adaptive Biotechnologies). Non-productive TCR CDR3 sequences (premature stop or frameshift), sequences with amino acid length less than 7, and sequences not starting with ‘C’ or ending with ‘F/W’ were excluded from the final analyses. Specimens with at least 1000 reads were included in the final analysis. To focus on T cells recognizing the same antigen, we analyzed amino acid clonotypes exclusively.
The degree of clonality for each specimen was assessed by the productive clonality matrix, which is defined as 1-Pielou’s evenness.28 Values near one represent samples with one or a few predominant clones (monoclonal or oligoclonal samples), whereas values near 0 represent a polyclonal population.
A random subsampling approach weighted by clonal abundance was used to equalize library sizes for relative comparisons of TCR repertoire overlap. For subsampling, each clonotype at amino acid level was treated as a sample and specimens were randomly sampled with replacement and weighted by clonal abundance (or frequency) until the total read count equaled that of the comparator library. To account for subsampling variation, the procedure was repeated 100 times and the 95%CIs for all subsampled comparisons are reported.
The degree of T cell clone overlap at the species level was evaluated using the Morisita overlap index.29 30 This measurement accounts for differences in library size and diversity per specimen, values near one the species occur in the same proportion in both samples, whereas values near 0 implies the two samples do not overlap in terms of species. Clonotypic sharing at the individual clone level was assessed in pairwise biological compartments before and after normalization to the same library size. Clones that were copresented in any of the compartment pairs are defined as shared clones. Based on the clonal frequency distributions (online supplemental figure 4), we assessed the top 40 clones in each specimen for overlap with the irAE repertoire. GLIPH223 24 was used for antigen-specific clustering. The motif significantly-enriched in the irAE was queried in a published dataset of tumor-reactive/skin-reactive TCRs in lung cancer26 and in a dataset of dermal/epidermal TCRs from healthy donors.25
Statistical analysis
Statistical analysis was performed using R software. The Mann Whitney U test was used for comparison of 2-group data. For analysis of >2 group data, Kruskal-Wallis was used. Spearman’s rho correlation was used to determine correlation significance. TCR preprocessing was performed using tcR package. Chord diagram was performed using the circlize package.31 32 p<0.05 was considered significant.
Data and code availability
Bulk TCR Vβ sequencing data generated by Adaptive Biotechnologies are available in the Adaptive Biotechnologies ImmuneACCESS repository at DOI: 10.21417/TRCJZ2021JITC. The code to perform downsampling of the TCR repertoire to the same library size and relevant figures are available online (https://github.com/BKI-immuno/dermatitis/).
Ethics statements
Ethics approval
This study was approved by the institutional review board (IRB) at Johns Hopkins University (JHU) and was conducted in accordance with the Declaration of Helsinki and the International Conference on Harmonization Good Clinical Practice guidelines. The patient described in this study provided written informed consent as approved by the IRB of JHU.
Acknowledgments
We thank the patient and the patient’s family for participation in this study, members of our research and administrative teams who contributed to this study, and also Fiamma Berner and Lukas Flatz for generous and prompt tumor-reactive/skin-reactive T-cell receptor sequencing data sharing.
References
Supplementary materials
Supplementary Data
This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.
Footnotes
Twitter @jihk99, @SmithImmunology
TC and JZ contributed equally.
Contributors TC, GJK, and H-YC conceived of and conducted the experiments. KNS and JT oversaw the study design, data interpretation, and manuscript preparation. JZ, BZ, PB, and HJ led the bioinformatic analyses. FV, JEH, HH, and MEA oversaw the clinical care of the patient and led the specimen acquisition. All authors contributed to and edited the manuscript.
Funding KNS was supported by the Lung Cancer Foundation of America, the IASLC Foundation, Swim Across America, and The Commonwealth Foundation. KNS, JT, JZ, and BZ were supported by the Mark Foundation for Cancer Research. HJ was partially supported by the National Institutes of Health (NIH)/National Human Genome Research Institute (grant R01HG009518). TC was supported by NIH (T32 CA193145). KNS and HJ were supported by R37 CA251447. This research was funded in part through the Bloomberg-Kimmel Institute for Cancer Immunotherapy, Bloomberg Philanthropies, and P30CA006973.
Competing interests HH has received clinical research funding from Bristol-Myers Squibb and Merck and serves in an advisory role for Pfizer, Merck, and Bristol-Myers Squibb. JT receives research funding from Bristol-Myers Squibb and serves a consulting/advisory role for Bristol-Myers Squibb, Merck, and Astra Zeneca. KNS has received travel support/honoraria from Illumina, Inc., receives research funding from Bristol-Myers Squibb, Enara Bio, and Astra Zeneca, and owns founder’s equity in manaT Bio. The terms of all these arrangements are being managed by the investigators’ respective institutions in accordance with their conflict of interest policies.
Provenance and peer review Not commissioned; externally peer reviewed.
Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.