High-throughput sequencing of the T-cell receptor repertoire: pitfalls and opportunities.

T-cell specificity is determined by the T-cell receptor, a heterodimeric protein coded for by an extremely diverse set of genes produced by imprecise somatic gene recombination. Massively parallel high-throughput sequencing allows millions of different T-cell receptor genes to be characterized from a single sample of blood or tissue. However, the extraordinary heterogeneity of the immune repertoire poses significant challenges for subsequent analysis of the data. We outline the major steps in processing of repertoire data, considering low-level processing of raw sequence files and high-level algorithms, which seek to extract biological or pathological information. The latest generation of bioinformatics tools allows millions of DNA sequences to be accurately and rapidly assigned to their respective variable V and J gene segments, and to reconstruct an almost error-free representation of the non-templated additions and deletions that occur. High-level processing can measure the diversity of the repertoire in different samples, quantify V and J usage and identify private and public T-cell receptors. Finally, we discuss the major challenge of linking T-cell receptor sequence to function, and specifically to antigen recognition. Sophisticated machine learning algorithms are being developed that can combine the paradoxical degeneracy and cross-reactivity of individual T-cell receptors with the specificity of the overall T-cell immune response. Computational analysis will provide the key to unlock the potential of the T-cell receptor repertoire to give insight into the fundamental biology of the adaptive immune system and to provide powerful biomarkers of disease.

Copy number architectures define treatment-mediated selection of lethal prostate cancer clones.

Despite initial responses to hormone treatment, metastatic prostate cancer invariably evolves to a lethal state. To characterize the intra-patient evolutionary relationships of metastases that evade treatment, we perform genome-wide copy number profiling and bespoke approaches targeting the androgen receptor (AR) on 167 metastatic regions from 11 organs harvested post-mortem from 10 men who died from prostate cancer. We identify diverse and patient-unique alterations clustering around the AR in metastases from every patient with evidence of independent acquisition of related genomic changes within an individual and, in some patients, the co-existence of AR-neutral clones. Using the genomic boundaries of pan-autosome copy number changes, we confirm a common clone of origin across metastases and diagnostic biopsies, and identified in individual patients, clusters of metastases occupied by dominant clones with diverged autosomal copy number alterations. These autosome-defined clusters are characterized by cluster-specific AR gene architectures, and in two index cases are topologically more congruent than by chance (p-values 3.07 × 10-8 and 6.4 × 10-4). Integration with anatomical sites suggests patterns of spread and points of genomic divergence. Here, we show that copy number boundaries identify treatment-selected clones with putatively distinct lethal trajectories.

Clinical testing of transcriptome-wide expression profiles in high-risk localized and metastatic prostate cancer starting androgen deprivation therapy: an ancillary study of the STAMPEDE abiraterone Phase 3 trial.

Metastatic and high-risk localized prostate cancer respond to hormone therapy but outcomes vary. Following a pre-specified statistical plan, we used Cox models adjusted for clinical variables to test associations with survival of multi-gene expression-based classifiers from 781 patients randomized to androgen deprivation with or without abiraterone in the STAMPEDE trial. Decipher score was strongly prognostic (p<2×10-5) and identified clinically-relevant differences in absolute benefit, especially for localized cancers. In metastatic disease, classifiers of proliferation, PTEN or TP53 loss and treatment-persistent cells were prognostic. In localized disease, androgen receptor activity was protective whilst interferon signaling (that strongly associated with tumor lymphocyte infiltration) was detrimental. Post-Operative Radiation-Therapy Outcomes Score was prognostic in localized but not metastatic disease (interaction p=0.0001) suggesting the impact of tumor biology on clinical outcome is context-dependent on metastatic state. Transcriptome-wide testing has clinical utility for advanced prostate cancer and identified worse outcomes for localized cancers with tumor-promoting inflammation.

Benchmarking of T cell receptor repertoire profiling methods reveals large systematic biases.

Monitoring the T cell receptor (TCR) repertoire in health and disease can provide key insights into adaptive immune responses, but the accuracy of current TCR sequencing (TCRseq) methods is unclear. In this study, we systematically compared the results of nine commercial and academic TCRseq methods, including six rapid amplification of complementary DNA ends (RACE)-polymerase chain reaction (PCR) and three multiplex-PCR approaches, when applied to the same T cell sample. We found marked differences in accuracy and intra- and inter-method reproducibility for T cell receptor α (TRA) and T cell receptor ß (TRB) TCR chains. Most methods showed a lower ability to capture TRA than TRB diversity. Low RNA input generated non-representative repertoires. Results from the 5' RACE-PCR methods were consistent among themselves but differed from the RNA-based multiplex-PCR results. Using an in silico meta-repertoire generated from 108 replicates, we found that one genomic DNA-based method and two non-unique molecular identifier (UMI) RNA-based methods were more sensitive than UMI methods in detecting rare clonotypes, despite the better clonotype quantification accuracy of the latter.

The naive T-cell receptor repertoire has an extremely broad distribution of clone sizes.

The clone size distribution of the human naive T-cell receptor (TCR) repertoire is an important determinant of adaptive immunity. We estimated the abundance of TCR sequences in samples of naive T cells from blood using an accurate quantitative sequencing protocol. We observe most TCR sequences only once, consistent with the enormous diversity of the repertoire. However, a substantial number of sequences were observed multiple times. We detect abundant TCR sequences even after exclusion of methodological confounders such as sort contamination, and multiple mRNA sampling from the same cell. By combining experimental data with predictions from models we describe two mechanisms contributing to TCR sequence abundance. TCRα abundant sequences can be primarily attributed to many identical recombination events in different cells, while abundant TCRß sequences are primarily derived from large clones, which make up a small percentage of the naive repertoire, and could be established early in the development of the T-cell repertoire.

The T cell differentiation landscape is shaped by tumour mutations in lung cancer.

Tumour mutational burden (TMB) predicts immunotherapy outcome in non-small cell lung cancer (NSCLC), consistent with immune recognition of tumour neoantigens. However, persistent antigen exposure is detrimental for T cell function. How TMB affects CD4 and CD8 T cell differentiation in untreated tumours, and whether this affects patient outcomes is unknown. Here we paired high-dimensional flow cytometry, exome, single-cell and bulk RNA sequencing from patients with resected, untreated NSCLC to examine these relationships. TMB was associated with compartment-wide T cell differentiation skewing, characterized by loss of TCF7-expressing progenitor-like CD4 T cells, and an increased abundance of dysfunctional CD8 and CD4 T cell subsets, with significant phenotypic and transcriptional similarity to neoantigen-reactive CD8 T cells. A gene signature of redistribution from progenitor-like to dysfunctional states associated with poor survival in lung and other cancer cohorts. Single-cell characterization of these populations informs potential strategies for therapeutic manipulation in NSCLC.

Publisher Correction: Spatial heterogeneity of the T cell receptor repertoire reflects the mutational landscape in lung cancer.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Quantitative analysis of the T cell receptor repertoire.

The T cell receptor repertoire provides a window into the cellular adaptive immune response. In the context of cancer, determining the repertoire within a tumor can give important insights into the evolution of the T cell anti-cancer response, and has the potential to identify specific personalized biomarkers for tracking host responses during cancer therapy, including immunotherapy. We describe a protocol for amplifying, sequencing and analyzing T cell receptors which is economical, robust, sensitive and versatile. The key experimental step is the ligation of a single stranded oligonucleotide to the 3' end of the T cell receptor cDNA, which allows easy amplification of all possible rearrangements using only a single set of primers per locus, while simultaneously introducing a unique molecular identifier to label each starting cDNA molecule. After sequencing, this molecular identifier can be used to correct both sequence errors and the effects of differential PCR amplification efficiency, thus producing a more accurate measure of the true T cell receptor frequency within the sample. We describe a detailed protocol describing this method to create libraries of T cell receptors from in vitro T cell cultures, blood or tissue samples. We combine this with a computational pipeline, which incorporates sample multiplexing, T cell receptor annotation and error correction to provide accurate counts of individual T cell receptor sequences within samples. The integrated experimental and computational pipeline should be of value to researchers interested in documenting and understanding the T cell immune response to cancer, and in manipulating it for therapeutic purposes.

Spatial heterogeneity of the T cell receptor repertoire reflects the mutational landscape in lung cancer.

Somatic mutations together with immunoediting drive extensive heterogeneity within non-small-cell lung cancer (NSCLC). Herein we examine heterogeneity of the T cell antigen receptor (TCR) repertoire. The number of TCR sequences selectively expanded in tumors varies within and between tumors and correlates with the number of nonsynonymous mutations. Expanded TCRs can be subdivided into TCRs found in all tumor regions (ubiquitous) and those present in a subset of regions (regional). The number of ubiquitous and regional TCRs correlates with the number of ubiquitous and regional nonsynonymous mutations, respectively. Expanded TCRs form part of clusters of TCRs of similar sequence, suggestive of a spatially constrained antigen-driven process. CD8+ tumor-infiltrating lymphocytes harboring ubiquitous TCRs display a dysfunctional tissue-resident phenotype. Ubiquitous TCRs are preferentially detected in the blood at the time of tumor resection as compared to routine follow-up. These findings highlight a noninvasive method to identify and track relevant tumor-reactive TCRs for use in adoptive T cell immunotherapy.

Urine-derived lymphocytes as a non-invasive measure of the bladder tumor immune microenvironment.

Despite the advances in cancer immunotherapy, only a fraction of patients with bladder cancer exhibit responses to checkpoint blockade, highlighting a need to better understand drug resistance and identify rational immunotherapy combinations. However, accessibility to the tumor prior and during therapy is a major limitation in understanding the immune tumor microenvironment (TME). Herein, we identified urine-derived lymphocytes (UDLs) as a readily accessible source of T cells in 32 patients with muscle invasive bladder cancer (MIBC). We observed that effector CD8+ and CD4+ cells and regulatory T cells within the urine accurately map the immune checkpoint landscape and T cell receptor repertoire of the TME. Finally, an increased UDL count, specifically high expression of PD-1 (PD-1hi) on CD8+ at the time of cystectomy, was associated with a shorter recurrence-free survival. UDL analysis represents a dynamic liquid biopsy that is representative of the bladder immune TME that may be used to identify actionable immuno-oncology (IO) targets with potential prognostic value in MIBC.

The T Cell Response to the Contact Sensitizer Paraphenylenediamine Is Characterized by a Polyclonal Diverse Repertoire of Antigen-Specific Receptors.

Paraphenylenediamine (PPD) is a common component of hair dyes and black henna tattoos and can cause skin sensitization and allergic contact dermatitis (ACD). The cutaneous inflammatory reaction associated with ACD is driven by both CD4+ and CD8+ T cells. However, the characteristics of such responses with respect to clonal breadth and magnitude are poorly defined. In this study, we have characterized the in vitro recall response of peripheral blood T cells prepared from PPD-allergic individuals to a PPD-human serum albumin (HSA) conjugate (PPD-HSA). Quantitative high throughput sequencing was used to characterize the changes in the repertoire of T cell receptor (TCR) α and ß genes after exposure to antigen in vitro. The PPD conjugate induced expansion of T cells carrying selected TCRs, with around 800 sequences (around 1%) being 8 or more times as abundant after culture than before. The expanded sequences showed strong skewing of V and J usage, consistent with an antigen-driven clonal expansion. The complementarity-determining region 3 sequences of the expanded TCRs could be grouped into several families of related amino acid sequence, but the overall diversity of the expanded sample was not much less than that of a random sample of the same size. The results suggest a model in which PPD-HSA conjugate stimulates a broad diversity of TCRs, with a wide range of stimulation strengths, which manifest as different degrees of in vitro expansion.

Quantitative Characterization of the T Cell Receptor Repertoire of Naïve and Memory Subsets Using an Integrated Experimental and Computational Pipeline Which Is Robust, Economical, and Versatile.

The T cell receptor (TCR) repertoire can provide a personalized biomarker for infectious and non-infectious diseases. We describe a protocol for amplifying, sequencing, and analyzing TCRs which is robust, sensitive, and versatile. The key experimental step is ligation of a single-stranded oligonucleotide to the 3' end of the TCR cDNA. This allows amplification of all possible rearrangements using a single set of primers per locus. It also introduces a unique molecular identifier to label each starting cDNA molecule. This molecular identifier is used to correct for sequence errors and for effects of differential PCR amplification efficiency, thus producing more accurate measures of the true TCR frequency within the sample. This integrated experimental and computational pipeline is applied to the analysis of human memory and naive subpopulations, and results in consistent measures of diversity and inequality. After error correction, the distribution of TCR sequence abundance in all subpopulations followed a power law over a wide range of values. The power law exponent differed between naïve and memory populations, but was consistent between individuals. The integrated experimental and analysis pipeline we describe is appropriate to studies of T cell responses in a broad range of physiological and pathological contexts.