tcR: an R package for T cell receptor repertoire advanced data analysis

T-cell Receptor Repertoire Analysis in the Context of Transarterial Chemoembolization Synergy with Systemic Therapy for Hepatocellular Carcinoma

T-cell receptor (TCR) sequencing provides a novel platform for insight into and characterization of intricate T-cell profiles, advancing the understanding of tumor immune heterogeneity. Recently, transarterial chemoembolization (TACE) combined with systemic therapy has become the recommended regimen for advanced hepatocellular carcinoma. The regulation of the immune microenvironment after TACE and its impact on tumor progression and recurrence has been a focus of research. By examining and tracking fluctuations in the TCR repertoire following combination treatment, novel perspectives on the modulation of the tumor microenvironment post-TACE and the underlying mechanisms governing tumor progression and recurrence can be gained. Clarifying the distinctive metrics and dynamic alterations of the TCR repertoire within the context of combination therapy is imperative for understanding the mechanisms of anti-tumor immunity, assessing efficacy, exploiting novel treatments, and further advancing precision oncology in the treatment of hepatocellular carcinoma. In this review, we initially summarized the fundamental characteristics of TCR repertoire and depicted immune microenvironment remodeling after TACE. Ultimately, we illustrated the prospective applications of TCR repertoires in TACE combined with systemic therapy.

The physiological landscape and specificity of antibody repertoires are consolidated by multiple immunizations

Diverse antibody repertoires spanning multiple lymphoid organs (i.e., bone marrow, spleen, lymph nodes) form the foundation of protective humoral immunity. Changes in their composition across lymphoid organs are a consequence of B-cell selection and migration events leading to a highly dynamic and unique physiological landscape of antibody repertoires upon antigenic challenge (e.g., vaccination). However, to what extent B cells encoding identical or similar antibody sequences (clones) are distributed across multiple lymphoid organs and how this is shaped by the strength of a humoral response remains largely unexplored. Here, we performed an in-depth systems analysis of antibody repertoires across multiple distinct lymphoid organs of immunized mice and discovered that organ-specific antibody repertoire features (i.e., germline V-gene usage and clonal expansion profiles) equilibrated upon a strong humoral response (multiple immunizations and high serum titers). This resulted in a surprisingly high degree of repertoire consolidation, characterized by highly connected and overlapping B-cell clones across multiple lymphoid organs. Finally, we revealed distinct physiological axes indicating clonal migrations and showed that antibody repertoire consolidation directly correlated with antigen specificity. Our study uncovered how a strong humoral response resulted in a more uniform but redundant physiological landscape of antibody repertoires, indicating that increases in antibody serum titers were a result of synergistic contributions from antigen-specific B-cell clones distributed across multiple lymphoid organs. Our findings provide valuable insights for the assessment and design of vaccine strategies.

TCRosetta: An Integrated Analysis and Annotation Platform for T-cell Receptor Sequences

T cells and T-cell receptors (TCRs) are essential components of the adaptive immune system. Characterization of the TCR repertoire offers a promising and highly informative source for understanding the functions of T cells in the immune response and immunotherapy. Although TCR repertoire studies have attracted much attention, there are few online servers available for TCR repertoire analysis, especially for TCR sequence annotation or advanced analyses. Therefore, we developed TCRosetta, a comprehensive online server that integrates analytical methods for TCR repertoire analysis and visualization. TCRosetta combines general feature analysis, large-scale sequence clustering, network construction, peptide-TCR binding prediction, generation probability calculation, and k-mer motif analysis for TCR sequences, making TCR data analysis as simple as possible. The TCRosetta server accepts multiple input data formats and can analyze ∼ 20,000 TCR sequences in less than 3 min. TCRosetta is the most comprehensive web server available for TCR repertoire analysis and is freely available at https://guolab.wchscu.cn/TCRosetta/.

A20 haploinsufficiency disturbs immune homeostasis and drives the transformation of lymphocytes with permissive antigen receptors

Genetic TNFAIP3 (A20) inactivation is a classical somatic lymphoma lesion and the genomic trait in haploinsufficiency of A20 (HA20). In a cohort of 34 patients with HA20, we show that heterozygous TNFAIP3 loss skews immune repertoires toward lymphocytes with classical self-reactive antigen receptors typically found in B and T cell lymphomas. This skewing was mediated by a feed-forward tumor necrosis factor (TNF)/A20/nuclear factor κB (NF-κB) loop that shaped pre-lymphoma transcriptome signatures in clonally expanded B (CD81, BACH2, and NEAT1) or T (GATA3, TOX, and PDCD1) cells. The skewing was reversed by anti-TNF treatment but could also progress to overt lymphoma. Analysis of conditional TNFAIP3 knock-out mice reproduced the wiring of the TNF/A20/NF-κB signaling axis with permissive antigen receptors and suggested a distinct regulation in B and T cells. Together, patients with the genetic disorder HA20 provide an exceptional window into A20/TNF/NF-κB-mediated control of immune homeostasis and early steps of lymphomagenesis that remain clinically unrecognized.

Genomic analysis and identification of a novel superantigen, SargEY, in Staphylococcus argenteus isolated from atopic dermatitis lesions

During surveillance of Staphylococcus aureus in lesions from patients with atopic dermatitis (AD), we isolated Staphylococcus argenteus, a species registered in 2011 as a new member of the genus Staphylococcus and previously considered a lineage of S. aureus. Genome sequence comparisons between S. argenteus isolates and representative S. aureus clinical isolates from various origins revealed that the S. argenteus genome from AD patients closely resembles that of S. aureus causing skin infections. We previously reported that 17%-22% of S. aureus isolated from skin infections produce staphylococcal enterotoxin Y (SEY), which predominantly induces T-cell proliferation via the T-cell receptor (TCR) Vα pathway. Complete genome sequencing of S. argenteus isolates revealed a gene encoding a protein similar to superantigen SEY, designated as SargEY, on its chromosome. Population structure analysis of S. argenteus revealed that these isolates are ST2250 lineage, which was the only lineage positive for the SEY-like gene among S. argenteus. Recombinant SargEY demonstrated immunological cross-reactivity with anti-SEY serum. SargEY could induce proliferation of human CD4+ and CD8+ T cells, as well as production of TNF-α and IFN-γ. SargEY showed emetic activity in a marmoset monkey model. SargEY and SET (a phylogenetically close but uncharacterized SE) revealed their dependency on TCR Vα in inducing human T-cell proliferation. Additionally, TCR sequencing revealed other previously undescribed Vα repertoires induced by SEH. SargEY and SEY may play roles in exacerbating the respective toxin-producing strains in AD.

IMPORTANCE

Staphylococcus aureus is frequently isolated from active lesions of atopic dermatitis (AD) patients. We reported that 17%-22% of S. aureus isolated from AD patients produced a novel superantigen staphylococcal enterotoxin Y (SEY). Unlike many S. aureus superantigens that activate T cells via T-cell receptor (TCR) Vß, SEY activates T cells via TCR Vα and stimulates cytokine secretion. Staphylococcus argenteus was isolated from AD patients during the surveillance for S. aureus. Phylogenetic comparison of the genome indicated that the isolate was very similar to S. aureus causing skin infections. The isolate encoded a SEY-like protein, designated SargEY, which, like SEY, activated T cells via the TCR Vα. ST2250 is the only lineage positive for SargEY gene. ST2250 S. argenteus harboring a superantigen SargEY gene may be a novel staphylococcal clone that infects human skin and is involved in the exacerbation of AD.

PVRIG is Expressed on Stem-Like T Cells in Dendritic Cell-Rich Niches in Tumors and Its Blockade May Induce Immune Infiltration in Non-Inflamed Tumors

Cancers that are poorly immune infiltrated pose a substantial challenge, with current immunotherapies yielding limited clinical success. Stem-like memory T cells (TSCM) have been identified as a subgroup of T cells that possess strong proliferative capacity and that can expand and differentiate following interactions with dendritic cells (DCs). In this study, we explored the pattern of expression of a recently discovered inhibitory receptor poliovirus receptor-related immunoglobulin domain protein (PVRIG) and its ligand, poliovirus receptor-related ligand 2 (PVRL2), in the human tumor microenvironment. Using spatial and single-cell RNA transcriptomics data across diverse cancer indications, we found that among the T-cell checkpoints, PVRIG is uniquely expressed on TSCM and PVRL2 is expressed on DCs in immune aggregate niches in tumors. PVRIG blockade could therefore enhance TSCM-DC interactions and efficiently drive T-cell infiltration to tumors. Consistent with these data, following PVRIG blockade in patients with poorly infiltrated tumors, we observed immune modulation including increased tumor T-cell infiltration, T-cell receptor (TCR) clonality, and intratumoral T-cell expansion, all of which were associated with clinical benefit. These data suggest PVRIG blockade as a promising strategy to induce potent antitumor T-cell responses, providing a novel approach to overcome resistance to immunotherapy in immune-excluded tumors.

Shape of the art: TCR-repertoire after allogeneic hematopoietic cell transplantation

The human adaptive immune repertoire is characterized by specificity and diversity to provide immunity against past and future tasks. Such tasks are mainly infections but also malignant transformations of cells. With its multiple lines of defense, the human immune system contains both, rapid reaction forces and the potential to capture, disassemble and analyze strange structures in order to teach the adaptive immune system and mount a specific immune response. Prevention and mitigation of autoimmunity is of equal importance. In the context of allogeneic hematopoietic cell transplantation (HCT) specific challenges exist with the transfer of cells from the adapted donor immune system to the immunosuppressed recipient. Those challenges are immunogenetic disparity between donor and host, reconstitution of immunity early after HCT by expansion of mature immune effector cells, and impaired thymic function, if the recipient is an adult (as it is the case in most HCTs). The possibility to characterize the adaptive immune repertoire by massively parallel sequencing of T-cell receptor gene rearrangements allows for a much more detailed characterization of the T-cell repertoire. In addition, high-dimensional characterization of immune effector cells based on their immunophenotype and single cell RNA sequencing allow for much deeper insights in adaptive immune responses. We here review, existing - still incomplete - information on immune reconstitution after allogeneic HCT. Building on the technological advances much deeper insights into immune recovery after HCT and adaptive immune responses and can be expected in the coming years.

Characterization of Th17 tissue-resident memory cells in non-inflamed intestinal tissue of Crohn's disease patients

Crohn's disease (CD) is a chronic inflammatory disorder affecting the bowel wall. Tissue-resident memory T (Trm) cells are implicated in CD, yet their characteristics remain unclear. We aimed to investigate the transcriptional profiles and functional characteristics of Trm cells in the small bowel of CD and their interactions with immune cells. Seven patients with CD and four with ulcerative colitis as controls were included. Single-cell RNA sequencing and paired T cell receptor sequencing assessed T cell subsets and transcriptional signatures in lamina propria (LP) and submucosa/muscularis propria-enriched fractions (SM/MP) from small bowel tissue samples. We detected 58,123 T cells grouped into 16 populations, including the CD4+ Trm cells with a Th17 signature and CD8+ Trm clusters. In CD, CD4+ Trm cells with a Th17 signature, termed Th17 Trm, showed significantly increased proportions within both the LP and SM/MP areas. The Th17 Trm cluster demonstrated heightened expression of tissue-residency marker genes (ITGAE, ITGA1, and CXCR6) along with elevated levels of IL17A, IL22, CCR6, and CCL20. The clonal expansion of Th17 Trm cells in CD was accompanied by enhanced transmural dynamic potential, as indicated by significantly higher migration scores. CD-prominent Th17 Trm cells displayed an increased interferon gamma (IFNγ)-related signature possibly linked with STAT1 activation, inducing chemokines (i.e., CXCL10, CXCL8, and CXCL9) in myeloid cells. Our findings underscored the elevated Th17 Trm cells throughout the small bowel in CD, contributing to disease pathogenesis through IFNγ induction and subsequent chemokine production in myeloid cells.

Identification of apolipoprotein B-reactive CDR3 motifs allows tracking of atherosclerosis-related memory CD4+T cells in multiple donors

Introduction

Atherosclerosis is a major pathological condition that underlies many cardiovascular diseases (CVDs). Its etiology involves breach of tolerance to self, leading to clonal expansion of autoreactive apolipoprotein B (APOB)-reactive CD4+T cells that correlates with clinical CVD. The T-cell receptor (TCR) sequences that mediate activation of APOB-specific CD4+T cells are unknown.

Methods

In a previous study, we had profiled the hypervariable complementarity determining region 3 (CDR3) of CD4+T cells that respond to six immunodominant APOB epitopes in most donors. Here, we comprehensively analyze this dataset of 149,065 APOB-reactive and 199,211 non-reactive control CDR3s from six human leukocyte antigen-typed donors.

Results

We identified 672 highly expanded (frequency threshold > 1.39E-03) clones that were significantly enriched in the APOB-reactive group as compared to the controls (log10 odds ratio ≥1, Fisher's test p < 0.01). Analysis of 114,755 naïve, 91,001 central memory (TCM) and 29,839 effector memory (TEM) CDR3 sequences from the same donors revealed that APOB+ clones can be traced to the complex repertoire of unenriched blood T cells. The fraction of APOB+ clones that overlapped with memory CDR3s ranged from 2.2% to 46% (average 16.4%). This was significantly higher than their overlap with the naïve pool, which ranged from 0.7% to 2% (average 1.36%). CDR3 motif analysis with the machine learning-based in-silico tool, GLIPHs (grouping of lymphocyte interactions by paratope hotspots), identified 532 APOB+ motifs. Analysis of naïve and memory CDR3 sequences with GLIPH revealed that ~40% (209 of 532) of these APOB+ motifs were enriched in the memory pool. Network analysis with Cytoscape revealed extensive sharing of the memory-affiliated APOB+ motifs across multiple donors. We identified six motifs that were present in TCM and TEM CDR3 sequences from >80% of the donors and were highly enriched in the APOB-reactive TCR repertoire.

Discussion

The identified APOB-reactive expanded CD4+T cell clones and conserved motifs can be used to annotate and track human atherosclerosis-related autoreactive CD4+T cells and measure their clonal expansion.

Mtb HLA-E-tetramer-sorted CD8+ T cells have a diverse TCR repertoire

HLA-E molecules can present self- and pathogen-derived peptides to both natural killer (NK) cells and T cells. T cells that recognize HLA-E peptides via their T cell receptor (TCR) are termed donor-unrestricted T cells due to restricted allelic variation of HLA-E. The composition and repertoire of HLA-E TCRs is not known so far. We performed TCR sequencing on CD8+ T cells from 21 individuals recognizing HLA-E tetramers (TMs) folded with two Mtb-HLA-E-restricted peptides. We sorted HLA-E Mtb TM+ and TM- CD8+ T cells directly ex vivo and performed bulk RNA-sequencing and single-cell TCR sequencing. The identified TCR repertoire was diverse and showed no conservation between and within individuals. TCRs selected from our single-cell TCR sequencing data could be activated upon HLA-E/peptide stimulation, although not robust, reflecting potentially weak interactions between HLA-E peptide complexes and TCRs. Thus, HLA-E-Mtb-specific T cells have a highly diverse TCR repertoire.

Seven-chain adaptive immune receptor repertoire analysis in rheumatoid arthritis reveals novel features associated with disease and clinically relevant phenotypes

BACKGROUND

In rheumatoid arthritis (RA), the activation of T and B cell clones specific for self-antigens leads to the chronic inflammation of the synovium. Here, we perform an in-depth quantitative analysis of the seven chains that comprise the adaptive immune receptor repertoire (AIRR) in RA.

RESULTS

In comparison to controls, we show that RA patients have multiple and strong differences in the B cell receptor repertoire including reduced diversity as well as altered isotype, chain, and segment frequencies. We demonstrate that therapeutic tumor necrosis factor inhibition partially restores this alteration but find a profound difference in the underlying biochemical reactivities between responders and non-responders. Combining the AIRR with HLA typing, we identify the specific T cell receptor repertoire associated with disease risk variants. Integrating these features, we further develop a molecular classifier that shows the utility of the AIRR as a diagnostic tool.

CONCLUSIONS

Simultaneous sequencing of the seven chains of the human AIRR reveals novel features associated with the disease and clinically relevant phenotypes, including response to therapy. These findings show the unique potential of AIRR to address precision medicine in immune-related diseases.

SARS-CoV-2-associated T-cell infiltration in the central nervous system

Objectives

Infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes coronavirus disease 2019 (COVID-19). Although an acute SARS-CoV-2 infection mainly presents with respiratory illness, neurologic symptoms and sequelae are increasingly recognised in the long-term treatment of COVID-19 patients. The pathophysiology and the neuropathogenesis behind neurologic complications of COVID-19 remain poorly understood, but mounting evidence points to endothelial dysfunction either directly caused by viral infection or indirectly by inflammatory cytokines, followed by a local immune response that may include virus-specific T cells. However, the type and role of central nervous system-infiltrating T cells in COVID-19 are complex and not fully understood.

Methods

We analysed distinct anatomical brain regions of patients who had deceased as a result of COVID-19-associated pneumonia or complications thereof and performed T cell receptor Vβ repertoire sequencing. Clonotypes were analysed for SARS-CoV-2 association using public TCR repertoire data.

Results

Our descriptive study demonstrates that SARS-CoV-2-associated T cells are found in almost all brain areas of patients with fatal COVID-19 courses. The olfactory bulb, medulla and cerebellum were brain regions showing the most SARS-CoV-2 specific sequence patterns. Neuropathological workup demonstrated primary CD8+ T-cell infiltration with a perivascular infiltration pattern.

Conclusion

Future research is needed to better define the relationship between T-cell infiltration and neurological symptoms and its long-term impact on patients' cognitive and mental health.

Single-cell and spatial architecture of primary liver cancer

Primary liver cancer (PLC) poses a leading threat to human health, and its treatment options are limited. Meanwhile, the investigation of homogeneity and heterogeneity among PLCs remains challenging. Here, using single-cell RNA sequencing, spatial transcriptomic and bulk multi-omics, we elaborated a molecular architecture of 3 PLC types, namely hepatocellular carcinoma (HCC), intrahepatic cholangiocarcinoma (ICC) and combined hepatocellular-cholangiocarcinoma (CHC). Taking a high-resolution perspective, our observations revealed that CHC cells exhibit internally discordant phenotypes, whereas ICC and HCC exhibit distinct tumor-specific features. Specifically, ICC was found to be the primary source of cancer-associated fibroblasts, while HCC exhibited disrupted metabolism and greater individual heterogeneity of T cells. We further revealed a diversity of intermediate-state cells residing in the tumor-peritumor junctional zone, including a congregation of CPE+ intermediate-state endothelial cells (ECs), which harbored the molecular characteristics of tumor-associated ECs and normal ECs. This architecture offers insights into molecular characteristics of PLC microenvironment, and hints that the tumor-peritumor junctional zone could serve as a targeted region for precise therapeutical strategies.

Sequencing the B Cell Receptor Repertoires of Antibody-Deficient Individuals With and Without Infection Susceptibility

PURPOSE

Most individuals with antibody deficiency (hypogammaglobulinemia) need immunoglobulin replacement therapy (IgG-RT) from healthy plasma donors to stay clear of infections. However, a small subset of hypogammaglobulinemic patients do not require this substitution therapy. We set out to investigate this clinical conundrum by asking whether the peripheral B cell receptor repertoires differ between antibody-deficient patients who do and do not need IgG-RT.

METHODS

We sequenced and analyzed IgG and IgM heavy chain B cell receptor repertoires from peripheral blood mononuclear cells (PBMCs) isolated from patients with low serum IgG concentrations who did or did not require IgG-RT.

RESULTS

Compared to the patients who did not need IgG-RT, those who needed IgG-RT had higher numbers of IgG antibody clones, higher IgM diversity, and less oligoclonal IgG and IgM repertoires. The patient cohorts had different heavy chain variable gene usage, and the patients who needed IgG-RT had elevated frequencies of IgG clones with higher germline identity (i.e., fewer somatic hypermutations).

CONCLUSION

Antibody-deficient patients with infection susceptibility who needed IgG-RT had more diverse peripheral antibody repertoires that were less diverged from germline and thus may not be as optimal for targeting pathogens, possibly contributing to infection susceptibility.

First-line treatment of unresectable or metastatic HER2 positive esophagogastric adenocarcinoma: liquid biomarker analysis of the phase 2 INTEGA trial

BACKGROUND

The addition of nivolumab to trastuzumab and chemotherapy in first-line unresectable or metastatic HER2 positive esophagogastric adenocarcinoma (HER2+ EGA) results in long progression-free and overall survival as shown by the INTEGA (ipilimumab or FOLFOX in combination with nivolumab and trastuzumab in HER2 positive esophagogastric adenocarcinoma) trial. This trial suggested that the chemotherapy backbone is needed in an unselected HER2+ patient population. Yet, it remains an open question if there are specific patient subsets that may benefit from an enhanced immunotherapeutic but chemotherapy-free approach.

METHODS

We analyzed blood T cell repertoire metrics determined by next-generation sequencing, circulating tumor cell (CTC) counts detected by CellSearch and their expression of HER2 and PD-L1 as potential liquid biomarkers predicting outcomes on ipilimumab versus FOLFOX (folinic acid, FOL, fluorouracil, F, oxaliplatin, OX) chemotherapy added to a backbone of trastuzumab and nivolumab in patients with HER2+ EGA in the INTEGA trial population.

RESULTS

Patients with two out of three baseline-determined liquid biomarkers-high T cell repertoire richness, absence of CTCs or HER2-expression on CTCs-made up approximately 44% of HER2+ EGA cases and did not show compromise in efficacy if treated with a chemotherapy-free regimen. Long-term responders showing a progression-free survival of >12 months were enriched in this biomarker triad, especially if treated on the chemotherapy-free arm.

CONCLUSION

Prospective validation of this liquid biomarker triad is needed to molecularly define HER2+ EGA patient subsets with different needs in the first-line systemic treatment setting.

T cell egress via lymphatic vessels is tuned by antigen encounter and limits tumor control

Antigen-specific CD8+ T cell accumulation in tumors is a prerequisite for effective immunotherapy, and yet the mechanisms of lymphocyte transit are not well defined. Here we show that tumor-associated lymphatic vessels control T cell exit from tumors via the chemokine CXCL12, and intratumoral antigen encounter tunes CXCR4 expression by effector CD8+ T cells. Only high-affinity antigen downregulates CXCR4 and upregulates the CXCL12 decoy receptor, ACKR3, thereby reducing CXCL12 sensitivity and promoting T cell retention. A diverse repertoire of functional tumor-specific CD8+ T cells, therefore, exit the tumor, which limits the pool of CD8+ T cells available to exert tumor control. CXCR4 inhibition or loss of lymphatic-specific CXCL12 boosts T cell retention and enhances tumor control. These data indicate that strategies to limit T cell egress might be an approach to boost the quantity and quality of intratumoral T cells and thereby response to immunotherapy.

A comprehensive platform for analyzing longitudinal multi-omics data

Longitudinal bulk and single-cell omics data is increasingly generated for biological and clinical research but is challenging to analyze due to its many intrinsic types of variations. We present PALMO ( https://github.com/aifimmunology/PALMO ), a platform that contains five analytical modules to examine longitudinal bulk and single-cell multi-omics data from multiple perspectives, including decomposition of sources of variations within the data, collection of stable or variable features across timepoints and participants, identification of up- or down-regulated markers across timepoints of individual participants, and investigation on samples of same participants for possible outlier events. We have tested PALMO performance on a complex longitudinal multi-omics dataset of five data modalities on the same samples and six external datasets of diverse background. Both PALMO and our longitudinal multi-omics dataset can be valuable resources to the scientific community.

Immune Activation following Irbesartan Treatment in a Colorectal Cancer Patient: A Case Study

Colorectal cancers are one of the most prevalent tumour types worldwide and, despite the emergence of targeted and biologic therapies, have among the highest mortality rates. The Personalized OncoGenomics (POG) program at BC Cancer performs whole genome and transcriptome analysis (WGTA) to identify specific alterations in an individual's cancer that may be most effectively targeted. Informed using WGTA, a patient with advanced mismatch repair-deficient colorectal cancer was treated with the antihypertensive drug irbesartan and experienced a profound and durable response. We describe the subsequent relapse of this patient and potential mechanisms of response using WGTA and multiplex immunohistochemistry (m-IHC) profiling of biopsies before and after treatment from the same metastatic site of the L3 spine. We did not observe marked differences in the genomic landscape before and after treatment. Analyses revealed an increase in immune signalling and infiltrating immune cells, particularly CD8+ T cells, in the relapsed tumour. These results indicate that the observed anti-tumour response to irbesartan may have been due to an activated immune response. Determining whether there may be other cancer contexts in which irbesartan may be similarly valuable will require additional studies.

T-Cell Receptor Repertoire Sequencing in the Era of Cancer Immunotherapy

T cells are integral components of the adaptive immune system, and their responses are mediated by unique T-cell receptors (TCR) that recognize specific antigens from a variety of biological contexts. As a result, analyzing the T-cell repertoire offers a better understanding of immune responses and of diseases like cancer. Next-generation sequencing technologies have greatly enabled the high-throughput analysis of the TCR repertoire. On the basis of our extensive experience in the field from the past decade, we provide an overview of TCR sequencing, from the initial library preparation steps to sequencing and analysis methods and finally to functional validation techniques. With regards to data analysis, we detail important TCR repertoire metrics and present several computational tools for predicting antigen specificity. Finally, we highlight important applications of TCR sequencing and repertoire analysis to understanding tumor biology and developing cancer immunotherapies.

TCR_Explore: A novel webtool for T cell receptor repertoire analysis

T cells expressing either alpha-beta or gamma-delta T cell receptors (TCR) are critical sentinels of the adaptive immune system, with receptor diversity being essential for protective immunity against a broad array of pathogens and agents. Programs available to profile TCR clonotypic signatures can be limiting for users with no coding expertise. Current analytical pipelines can be inefficient due to manual processing steps, open to data entry errors and have multiple analytical tools with unique inputs that require coding expertise. Here we present a bespoke webtool designed for users irrespective of coding expertise, coined 'TCR_Explore', enabling analysis either derived via Sanger sequencing or next generation sequencing (NGS) platforms. Further, TCR_Explore incorporates automated quality control steps for Sanger sequencing. The creation of flexible and publication ready figures are enabled for different sequencing platforms following universal conversion to the TCR_Explore file format. TCR_Explore will enhance a user's capacity to undertake in-depth TCR repertoire analysis of both new and pre-existing datasets for identification of T cell clonotypes associated with health and disease. The web application is located at https://tcr-explore.erc.monash.edu for users to interactively explore TCR repertoire datasets.

Analyzing Antibody Repertoire Using Next-Generation Sequencing and Machine Learning

Advances in high-throughput sequencing technologies have enabled comprehensive sequencing of the immune repertoire. Since repertoire analysis can help to explain the relationship between the immune system and diseases, several methods have been developed for repertoire analysis. Here, using simulated and real-world datasets, we describe how to use DeepRC, a method that applies cutting-edge machine learning techniques.

Trafficking between clonally related peripheral T-helper cells and tissue-resident T-helper cells in chronic GVHD

Chronic graft-versus-host disease (cGVHD) is an autoimmune-like syndrome. CXCR5-PD-1hi peripheral T-helper (Tph) cells have an important pathogenic role in autoimmune diseases, but the role of Tph cells in cGVHD remains unknown. We show that in patients with cGVHD, expansion of Tph cells among blood CD4+ T cells was associated with cGVHD severity. These cells augmented memory B-cell differentiation and production of immunoglobulin G via interleukin 21 (IL-21). Tph cell expansion was also observed in a murine model of cGVHD. This Tph cell expansion in the blood is associated with the expansion of pathogenic tissue-resident T-helper (Trh) cells that form lymphoid aggregates surrounded by collagen in graft-versus-host disease (GVHD) target tissues. Adoptive transfer experiments showed that Trh cells from GVHD target tissues give rise to Tph cells in the blood, and conversely, Tph cells from the blood give rise to Trh cells in GVHD target tissues. Tph cells in the blood and Trh cells in GVHD target tissues had highly overlapping T-cell receptor α and β repertoires. Deficiency of IL-21R, B-cell lymphoma 6 (BCL6), or T-bet in donor T cells markedly reduced the proportions of Tph cells in the blood and Trh cells in GVHD target tissues and reduced T-B interaction in the lymphoid aggregates. These results indicate that clonally related pathogenic Tph cells and Trh cells traffic between the blood and cGVHD target tissues, and that IL-21R-BCL6 signaling and T-bet are required for the development and expansion of Tph and Trh cells in the pathogenesis of cGVHD.

TCR repertoire profiling revealed antigen-driven CD8+ T cell clonal groups shared in synovial fluid of patients with spondyloarthritis

Spondyloarthritis (SpA) comprises a number of inflammatory rheumatic diseases with overlapping clinical manifestations. Strong association with several HLA-I alleles and T cell infiltration into an inflamed joint suggest involvement of T cells in SpA pathogenesis. In this study, we performed high-throughput T cell repertoire profiling of synovial fluid (SF) and peripheral blood (PB) samples collected from a large cohort of SpA patients. We showed that synovial fluid is enriched with expanded T cell clones that are shared between patients with similar HLA genotypes and persist during recurrent synovitis. Using an algorithm for identification of TCRs involved in immune response we discovered several antigen-driven CD8+ clonal groups associated with risk HLA-B*27 or HLA-B*38 alleles. We further show that these clonal groups were enriched in SF and had higher frequency in PB of SpA patients vs healthy donors, implying their relevance to SpA pathogenesis. Several of the groups were shared among patients with different SpAs that suggests a common immunopathological mechanism of the diseases. In summary, our results provide evidence for the role of specific CD8+ T cell clones in pathogenesis of SpA.

Combined PD-L1 and TIM3 blockade improves expansion of fit human CD8+ antigen-specific T cells for adoptive immunotherapy

Antigen-specific T cell expansion ex vivo followed by adoptive transfer enables targeting of a multitude of microbial and cancer antigens. However, clinical-scale T cell expansion from rare precursors requires repeated stimulation, which may lead to T cell dysfunction and limited therapeutic potential. We used a clinically compliant protocol to expand Epstein-Barr virus (EBV) and Wilms tumor 1 (WT1) antigen-specific CD8⁺ T cells, and leveraged T cell exhaustion-associated inhibitory receptor blockade to improve T cell expansion. Several inhibitory receptors were expressed early by ex vivo-expanded antigen-specific CD8⁺ T cells, including PD-1 and TIM3, with co-expression matching evidence of T cell dysfunction as the cultures progressed. Introduction of anti-PD-L1 and anti-TIM3 blockade in combination (but not individually) to the culture led to markedly improved antigen-specific T cell expansion without inducing T cell dysfunction. Single-cell RNA sequencing (RNA-seq) and T cell receptor (TCR) repertoire profiling revealed that double blockade does not impart specific transcriptional programs in T cells or alterations in TCR repertoires. However, combined blockade may affect gene expression in a minority of clonotypes in a donor-specific fashion. We conclude that antigen-specific CD8⁺ T cell manufacturing can be improved by using TIM3 and PD-L1/PD-1 axis blockade in combination. This approach is readily applicable to several adoptive immunotherapy strategies.

Healthy-like CD4+ Regulatory and CD4+ Conventional T-Cell Receptor Repertoires Predict Protection from GVHD Following Donor Lymphocyte Infusion

Donor lymphocyte infusion (DLI) can (re-)induce durable remission in relapsing patients after allogeneic hematopoietic stem-cell transplantation (alloHSCT). However, DLI harbors the risk of increased non-relapse mortality due to the co-occurrence of graft-versus-host disease (GVHD). GVHD onset may be caused or accompanied by changes in the clonal T-cell receptor (TCR) repertoire. To investigate this, we analyzed T cells in a cohort of 21 patients receiving DLI after alloHSCT. We performed deep T-cell receptor β (TRB) sequencing of sorted CD4+CD25+CD127low regulatory T cells (Treg cells) and CD4+ conventional T cells (Tcon cells) in order to track longitudinal changes in the TCR repertoire. GVHD following DLI was associated with less diverse but clonally expanded CD4+CD25+CD127low Treg and CD4+ Tcon TCR repertoires, while patients without GVHD exhibited healthy-like repertoire properties. Moreover, the diversification of the repertoires upon GVHD treatment was linked to steroid-sensitive GVHD, whereas decreased diversity was observed in steroid-refractory GVHD. Finally, the unbiased sample analysis revealed that the healthy-like attributes of the CD4+CD25+CD127low Treg TCR repertoire were associated with reduced GVHD incidence. In conclusion, CD4+CD25+CD127low Treg and CD4+ Tcon TRB repertoire dynamics may provide a helpful real-time tool to improve the diagnosis and monitoring of treatment in GVHD following DLI.

Allelic polymorphism controls autoreactivity and vaccine elicitation of human broadly neutralizing antibodies against influenza virus

Human broadly neutralizing antibodies (bnAbs) targeting the hemagglutinin stalk of group 1 influenza A viruses (IAVs) are biased for IGHV1-69 alleles that use phenylalanine (F54) but not leucine (L54) within their CDRH2 loops. Despite this, we demonstrated that both alleles encode for human IAV bnAbs that employ structurally convergent modes of contact to the same epitope. To resolve differences in lineage expandability, we compared F54 versus L54 as substrate within humanized mice, where antibodies develop with human-like CDRH3 diversity but are restricted to single VH genes. While both alleles encoded for bnAb precursors, only F54 IGHV1-69 supported elicitation of heterosubtypic serum bnAbs following immunization with a stalk-only nanoparticle vaccine. L54 IGHV1-69 was unproductive, co-encoding for anergic B cells and autoreactive stalk antibodies that were cleared from B cell memory. Moreover, human stalk antibodies also demonstrated L54-dependent autoreactivity. Therefore, IGHV1-69 polymorphism, which is skewed ethnically, gates tolerance and vaccine expandability of influenza bnAbs.

Follicular Lymphoma Microenvironment Characteristics Associated with Tumor Cell Mutations and MHC Class II Expression

Follicular lymphoma (FL) is a B-cell malignancy with a complex tumor microenvironment that is rich in nonmalignant immune cells. We applied single-cell RNA sequencing to characterize the diverse tumor and immune cell populations of FL and identified major phenotypic subsets of FL T cells, including a cytotoxic CD4 T-cell population. We characterized four major FL subtypes with differential representation or relative depletion of distinct T-cell subsets. By integrating exome sequencing, we observed that somatic mutations are associated with, but not definitive for, reduced MHC expression on FL cells. In turn, expression of MHCII genes by FL cells was associated with significant differences in the proportions and targetable immunophenotypic characteristics of T cells. This provides a classification framework of the FL microenvironment in association with FL genotypes and MHC expression, and informs different potential immunotherapeutic strategies based upon tumor cell MHCII expression.

SIGNIFICANCE

We have characterized the FL-infiltrating T cells, identified cytotoxic CD4 T cells as an important component that is associated with tumor cell-intrinsic characteristics, and identified sets of targetable immune checkpoints on T cells that differed from FLs with normal versus low MHC expression. See related commentary by Melnick, p. 374. This article is highlighted in the In This Issue feature, p. 369.

CD39+ tissue-resident memory CD8+ T cells with a clonal overlap across compartments mediate antitumor immunity in breast cancer

Despite being a standard treatment option in breast cancer, immune checkpoint inhibitors (ICIs) are only efficacious for a subset of patients. To gain a better understanding of the antitumor immune response in breast cancer, we examined the heterogeneity of CD8⁺ T cells in tumors, metastatic lymph nodes (mLNs), and peripheral blood from patients with early breast cancer (n = 131). Among tissue-resident memory CD8⁺ T (T_RM) cells, including virus- and tumor-specific CD8⁺ T cells, CD39 expression was observed in a tumor-specific and exhausted subpopulation in both tumors and mLNs. CD39⁺ T_RM cells from tumors and mLNs exhibited a phenotypic similarity and clonally overlapped with each other. Moreover, tumor or mLN CD39⁺ T_RM cells clonally overlapped with CD39^- T_RM and non-T_RM cells in the same compartment, implying a tissue-specific differentiation process. These inter-subpopulationally overlapping CD39⁺ T_RM clonotypes were frequently detected among effector memory CD8⁺ T cells in peripheral blood, suggesting a systemic clonal overlap. CD39⁺ T_RM cell enrichment was heterogeneous among molecular subtypes of breast cancer, which is associated with the different role of antitumor immune responses in each subtype. In vitro blockade of PD-1 and/or CTLA-4 effectively restored proliferation of CD39⁺ T_RM cells and enhanced cytokine production by CD8⁺ T cells from tumors or mLNs, particularly in the presence of CD39⁺ T_RM enrichment. This suggests that CD39⁺ T_RM cells have a capacity for functional restoration upon ICI treatment. Thus, our study indicates that CD39⁺ T_RM cells with a clonal overlap across compartments are key players in antitumor immunity in breast cancer.

GMP Manufacturing and IND-Enabling Studies of a Recombinant Hyperimmune Globulin Targeting SARS-CoV-2

Conventionally, hyperimmune globulin drugs manufactured from pooled immunoglobulins from vaccinated or convalescent donors have been used in treating infections where no treatment is available. This is especially important where multi-epitope neutralization is required to prevent the development of immune-evading viral mutants that can emerge upon treatment with monoclonal antibodies. Using microfluidics, flow sorting, and a targeted integration cell line, a first-in-class recombinant hyperimmune globulin therapeutic against SARS-CoV-2 (GIGA-2050) was generated. Using processes similar to conventional monoclonal antibody manufacturing, GIGA-2050, comprising 12,500 antibodies, was scaled-up for clinical manufacturing and multiple development/tox lots were assessed for consistency. Antibody sequence diversity, cell growth, productivity, and product quality were assessed across different manufacturing sites and production scales. GIGA-2050 was purified and tested for good laboratory procedures (GLP) toxicology, pharmacokinetics, and in vivo efficacy against natural SARS-CoV-2 infection in mice. The GIGA-2050 master cell bank was highly stable, producing material at consistent yield and product quality up to >70 generations. Good manufacturing practices (GMP) and development batches of GIGA-2050 showed consistent product quality, impurity clearance, potency, and protection in an in vivo efficacy model. Nonhuman primate toxicology and pharmacokinetics studies suggest that GIGA-2050 is safe and has a half-life similar to other recombinant human IgG1 antibodies. These results supported a successful investigational new drug application for GIGA-2050. This study demonstrates that a new class of drugs, recombinant hyperimmune globulins, can be manufactured consistently at the clinical scale and presents a new approach to treating infectious diseases that targets multiple epitopes of a virus.

Machine Learning Approaches to TCR Repertoire Analysis

Sparked by the development of genome sequencing technology, the quantity and quality of data handled in immunological research have been changing dramatically. Various data and database platforms are now driving the rapid progress of machine learning for immunological data analysis. Of various topics in immunology, T cell receptor repertoire analysis is one of the most important targets of machine learning for assessing the state and abnormalities of immune systems. In this paper, we review recent repertoire analysis methods based on machine learning and deep learning and discuss their prospects.

γ9δ2 T-Cell Expansion and Phenotypic Profile Are Reflected in the CDR3δ Repertoire of Healthy Adults

γ9δ2T cells fill a distinct niche in human immunity due to the unique physiology of the phosphoantigen-reactive γ9δ2TCR. Here, we highlight reproducible TCRδ complementarity-determining region 3 (CDR3δ) repertoire patterns associated with γ9δ2T cell proliferation and phenotype, thus providing evidence for the role of the CDR3δ in modulating in vivo T-cell responses. Features that determine γ9δ2TCR binding affinity and reactivity to the phosphoantigen-induced ligand in vitro appear to similarly underpin in vivo clonotypic expansion and differentiation. Likewise, we identify a CDR3δ bias in the γ9δ2T cell natural killer receptor (NKR) landscape. While expression of the inhibitory receptor CD94/NKG2A is skewed toward cells bearing putative high-affinity TCRs, the activating receptor NKG2D is expressed independently of the phosphoantigen-sensing determinants, suggesting a higher net NKR activating signal in T cells with TCRs of low affinity. This study establishes consistent repertoire–phenotype associations and justifies stratification for the T-cell phenotype in future research on γ9δ2TCR repertoire dynamics.

Intrahepatic CD69+Vδ1 T cells re-circulate in the blood of patients with metastatic colorectal cancer and limit tumor progression

Background

More than 50% of all patients with colorectal cancer (CRC) develop liver metastases (CLM), a clinical condition characterized by poor prognosis and lack of reliable prognostic markers. Vδ1 cells are a subset of tissue-resident gamma delta (γδ) T lymphocytes endowed with a broad array of antitumor functions and showing a natural high tropism for the liver. However, little is known about their impact in the clinical outcomes of CLM.

Methods

We isolated human γδ T cells from peripheral blood (PB) and peritumoral (PT) tissue of 93 patients undergone surgical procedures to remove CLM. The phenotype of freshly purified γδ T cells was assessed by multiparametric flow cytometry, the transcriptional profiles by single cell RNA-sequencing, the functional annotations by Gene Ontology enrichment analyses and the clonotype by γδ T cell receptor (TCR)-sequencing.

Results

The microenvironment of CLM is characterized by a heterogeneous immune infiltrate comprising different subsets of γδ tumor-infiltrating lymphocytes (TILs) able to egress the liver and re-circulate in PB. Vδ1 T cells represent the largest population of γδ TILs within the PT compartment of CLM that is greatly enriched in Vδ1 T effector (T_EF) cells expressing constitutive high levels of CD69. These Vδ1 CD69⁺ TILs express a distinct phenotype and transcriptional signature, show high antitumor potential and correlate with better patient clinical outcomes in terms of lower numbers of liver metastatic lesions and longer overall survival (OS). Moreover, intrahepatic CD69⁺ Vδ1 TILs can egress CLM tissue to re-circulate in PB, where they retain a phenotype, transcriptional signature and TCR clonal repertoires resembling their liver origin. Importantly, even the increased frequencies of the CD69⁺ terminally differentiated (T_EMRA) Vδ1 cells in PB of patients with CLM significantly correlate with longer OS. The positive prognostic score of high frequencies of CD69⁺ T_EMRA Vδ1 cells in PB is independent from the neoadjuvant chemotherapy and immunotherapy regimens administered to patients with CLM prior surgery.

Conclusions

The enrichment of tissue-resident CD69⁺ Vδ1 T_EMRA cells re-circulating at high frequencies in PB of patients with CLM limits tumor progression and represents a new important clinical tool to either predict the natural history of CLM or develop alternative therapeutic protocols of cellular therapies.

Homologous recombination repair deficient prostate cancer represents an immunologically distinct subtype

Homologous recombination repair deficiency (HRD) is observed in 10% of patients with castrate-resistant prostate cancer (PCa). Preliminary data suggest that HRD-PCa might be more responsive to immune checkpoint inhibitors (ICIs). In this study, we compare the tumor immune landscape and peripheral T cell receptor (TCR) repertoire of patients with and without HRD-PCa to gain further insight into the immunogenicity of HRD-PCa. Immunohistochemistry was performed on tumor tissue of 81 patients, including 15 patients with HRD-PCa. Peripheral TCR sequencing was performed in a partially overlapping cohort of 48 patients, including 16 patients with HRD-PCa. HRD patients more frequently had intratumoral CD3⁺, CD3⁺CD8⁻FoxP3⁻ or Foxp3⁺ TILs above median compared to patients without DNA damage repair alterations (DDRwt; CD3⁺ and Foxp3⁺: 77% vs 35%, p = .013; CD3⁺CD8⁻FoxP3⁻: 80% vs 44%, p = .031). No significant difference in CD8⁺ TILs or PD-L1 expression was observed. In peripheral blood, HRD patients displayed a more diverse TCR repertoire compared to DDRwt patients (p = .014). Additionally, HRD patients shared TCR clusters with low generation probability, suggesting patient-overlapping T cell responses. A pooled analysis of clinical data from 227 patients with molecularly characterized PCa suggested increased efficacy of ICIs in HRD-PCa. In conclusion, patients with HRD-PCa display increased TIL density and an altered peripheral TCR repertoire. Further research into the efficacy of ICIs and the presence of shared neoantigens in HRD-PCa is warranted.

Single-cell RNA sequencing coupled to TCR profiling of large granular lymphocyte leukemia T cells

T-cell large granular lymphocyte leukemia (T-LGLL) is a lymphoproliferative disease and bone marrow failure syndrome which responds to immunosuppressive therapies. We show single-cell TCR coupled with RNA sequencing of CD3+ T cells from 13 patients, sampled before and after alemtuzumab treatments. Effector memory T cells and loss of T cell receptor (TCR) repertoire diversity are prevalent in T-LGLL. Shared TCRA and TCRB clonotypes are absent. Deregulation of cell survival and apoptosis gene programs, and marked downregulation of apoptosis genes in CD8+ clones, are prominent features of T-LGLL cells. Apoptosis genes are upregulated after alemtuzumab treatment, especially in responders than non-responders; baseline expression levels of apoptosis genes are predictive of hematologic response. Alemtuzumab does not attenuate TCR clonality, and TCR diversity is further skewed after treatment. Inferences made from analysis of single cell data inform understanding of the pathophysiologic mechanisms of clonal expansion and persistence in T-LGLL.

Circulating Tregs accumulate in omental tumors and acquire adipose-resident features

Tumors that metastasize in the peritoneal cavity typically end up in the omental adipose tissue, a particularly immune-suppressive environment that includes specialized adipose-resident regulatory T cells (Tregs). Tregs rapidly accumulate in the omentum after tumor implantation and potently suppress anti-tumor immunity. However, it is unclear whether these Tregs are recruited from the circulation or derived from pre-existing adipose-resident Tregs by clonal expansion. Here we show that Tregs in tumor-bearing omenta predominantly have thymus-derived characteristics. Moreover, naïve tumor antigen-specific CD4+ T cells fail to differentiate into Tregs in tumor-bearing omenta. In fact, Tregs derived from the pre-tumor repertoire are sufficient to suppress anti-tumor immunity and promote tumor growth. However, tumor implantation in the omentum does not promote Treg clonal expansion, but instead leads to increased clonal diversity. Parabiosis experiments show that despite tissue-resident (non-circulating) characteristics of omental Tregs in naïve mice, tumor implantation promotes a rapid influx of circulating Tregs, many of which come from the spleen. Finally, we show that newly recruited Tregs rapidly acquire characteristics of adipose-resident Tregs in tumor-bearing omenta. These data demonstrate that most Tregs in omental tumors are recruited from the circulation and adapt to their environment by altering their homing, transcriptional and metabolic properties.

Targeting the tumor mutanome for personalized vaccination in a TMB low non-small cell lung cancer

BACKGROUND

Cancer is characterized by an accumulation of somatic mutations, of which a significant subset can generate cancer-specific neoepitopes that are recognized by autologous T cells. Such neoepitopes are emerging as important targets for cancer immunotherapy, including personalized cancer vaccination strategies.

METHODS

We used whole-exome and RNA sequencing analysis to identify potential neoantigens for a patient with non-small cell lung cancer. Thereafter, we assessed the autologous T-cell reactivity to the candidate neoantigens using a long peptide approach in a cultured interferon gamma ELISpot and tracked the neoantigen-specific T-cells in the tumor by T-cell receptor (TCR) sequencing. In parallel, identified gene variants were incorporated into a Modified Vaccinia Ankara-based vaccine, which was evaluated in the human leucocyte antigen A*0201 transgenic mouse model (HHD).

RESULTS

Sequencing revealed a tumor with a low mutational burden: 2219 sequence variants were identified from the primary tumor, of which 23 were expressed in the transcriptome, involving 18 gene products. We could demonstrate spontaneous T-cell responses to 5/18 (28%) mutated gene variants, and further analysis of the TCR repertoire of neoantigen-specific CD4+ and CD8+ T cells revealed TCR clonotypes that were expanded in both blood and tumor tissue. Following vaccination of HHD mice, de novo T-cell responses were generated to 4/18 (22%) mutated gene variants; T cells reactive against two variants were also evident in the autologous setting. Subsequently, we determined the major histocompatibility complex restriction of the T-cell responses and used in silico prediction tools to determine the likely neoepitopes.

CONCLUSIONS

Our study demonstrates the feasibility of efficiently identifying tumor-specific neoantigens that can be targeted by vaccination in tumors with a low mutational burden, promising successful clinical exploitation, with trials currently underway.

Haematopoietic Stem Cell Transplantation Results in Extensive Remodelling of the Clonal T Cell Repertoire in Multiple Sclerosis

Autologous haematopoietic stem cell transplantation (AHSCT) is a vital therapeutic option for patients with highly active multiple sclerosis (MS). Rates of remission suggest AHSCT is the most effective form of immunotherapy in controlling the disease. Despite an evolving understanding of the biology of immune reconstitution following AHSCT, the mechanism by which AHSCT enables sustained disease remission beyond the period of lymphopenia remains to be elucidated. Auto-reactive T cells are considered central to MS pathogenesis. Here, we analyse T cell reconstitution for 36 months following AHSCT in a cohort of highly active MS patients. Through longitudinal analysis of sorted naïve and memory T cell clones, we establish that AHSCT induces profound changes in the dominant T cell landscape of both CD4+ and CD8+ memory T cell clones. Lymphopenia induced homeostatic proliferation is followed by clonal attrition; with only 19% of dominant CD4 (p <0.025) and 13% of dominant CD8 (p <0.005) clones from the pre-transplant repertoire detected at 36 months. Recovery of a thymically-derived CD4 naïve T cell repertoire occurs at 12 months and is ongoing at 36 months, however diversity of the naïve populations is not increased from baseline suggesting the principal mechanism of durable remission from MS after AHSCT relates to depletion of putative auto-reactive clones. In a cohort of MS patients expressing the MS risk allele HLA DRB1*15:01, public clones are probed as potential biomarkers of disease. AHSCT appears to induce sustained periods of disease remission with dynamic changes in the clonal T cell repertoire out to 36 months post-transplant.

Azacitidine-induced reconstitution of the bone marrow T cell repertoire is associated with superior survival in AML patients

Hypomethylating agents (HMA) like azacitidine are licensed for the treatment of acute myeloid leukemia (AML) patients ineligible for allogeneic hematopoietic stem cell transplantation. Biomarker-driven identification of HMA-responsive patients may facilitate the choice of treatment, especially in the challenging subgroup above 60 years of age. Since HMA possesses immunomodulatory functions that constitute part of their anti-tumor effect, we set out to analyze the bone marrow (BM) immune environment by next-generation sequencing of T cell receptor beta (TRB) repertoires in 51 AML patients treated within the RAS-AZIC trial. Patients with elevated pretreatment T cell diversity (11 out of 41 patients) and those with a boost of TRB richness on day 15 after azacitidine treatment (12 out of 46 patients) had longer event-free and overall survival. Both pretreatment and dynamic BM T cell metrics proved to be better predictors of outcome than other established risk factors. The favorable broadening of the BM T cell space appeared to be driven by antigen since these patients showed significant skewing of TRBV gene usage. Our data suggest that one course of AZA can cause reconstitution to a more physiological T cell BM niche and that the T cell space plays an underestimated prognostic role in AML.

Trial registration: DRKS identifier: DRKS00004519

Heterogeneity of human anti-viral immunity shaped by virus, tissue, age, and sex

The persistence of anti-viral immunity is essential for protection and exhibits profound heterogeneity across individuals. Here, we elucidate the factors that shape maintenance and function of anti-viral T cell immunity in the body by comprehensive profiling of virus-specific T cells across blood, lymphoid organs, and mucosal tissues of organ donors. We use flow cytometry, T cell receptor sequencing, single-cell transcriptomics, and cytokine analysis to profile virus-specific CD8+ T cells recognizing the ubiquitous pathogens influenza and cytomegalovirus. Our results reveal that virus specificity determines overall magnitude, tissue distribution, differentiation, and clonal repertoire of virus-specific T cells. Age and sex influence T cell differentiation and dissemination in tissues, while T cell tissue residence and functionality are highly correlated with the site. Together, our results demonstrate how the covariates of virus, tissue, age, and sex impact the anti-viral immune response, which is important for targeting, monitoring, and predicting immune responses to existing and emerging viruses.

Clonal expansion of CD8+ T cells reflects graft-versus-leukemia activity and precedes durable remission following DLI

GVL activity in patients receiving DLI for relapse after allogeneic hematopoietic stem cell transplantation is associated with clonal expansion of CD8⁺ T cells.

Absence of clonal expansion of the CD8⁺TRB repertoire after DLI predicts relapse at a median time of 11.2 months before clinical diagnosis.

Donor lymphocyte infusion (DLI) is a standard of care for relapse of acute myeloid leukemia after allogeneic hematopoietic stem cell transplantation. Currently it is poorly understood how and when CD8⁺ αβ T cells exert graft-versus-leukemia (GVL) activity after DLI. Also, there is no reliable biomarker to monitor GVL activity of the infused CD8⁺ T cells. Therefore, we analyzed the dynamics of CD8⁺ αβ T-cell clones in patients with DLI. In this prospective clinical study of 29 patients, we performed deep T-cell receptor β (TRB ) sequencing of sorted CD8⁺ αβ T cells to track patients’ repertoire changes in response to DLI. Upon first occurrence of GVL, longitudinal analyses revealed a preferential expansion of distinct CD8⁺TRB clones (n = 14). This did not occur in samples of patients without signs of GVL (n = 11). Importantly, early repertoire changes 15 days after DLI predicted durable remission for the 36-month study follow-up. Furthermore, absence of clonal outgrowth of the CD8⁺TRB repertoire after DLI was an early biomarker that predicted relapse at a median time of 11.2 months ahead of actual diagnosis. Additionally, unbiased sample analysis regardless of the clinical outcome revealed that patients with decreasing CD8⁺TRB diversity at day 15 after DLI (n = 13) had a lower relapse incidence (P = .0040) compared with patients without clonal expansion (n = 6). In conclusion, CD8⁺TRB analysis may provide a reliable tool for predicting the efficacy of DLI and holds the potential to identify patients at risk for progression and relapse after DLI.

Response and recurrence correlates in individuals treated with neoadjuvant anti-PD-1 therapy for resectable oral cavity squamous cell carcinoma

Neoadjuvant PD-1 blockade may be efficacious in some individuals with high-risk, resectable oral cavity head and neck cancer. To explore correlates of response patterns to neoadjuvant nivolumab treatment and post-surgical recurrences, we analyzed longitudinal tumor and blood samples in a cohort of 12 individuals displaying 33% responsiveness. Pretreatment tumor-based detection of FLT4 mutations and PTEN signature enrichment favors response, and high tumor mutational burden improves recurrence-free survival. In contrast, preexisting and/or acquired mutations (in CDKN2A, YAP1, or JAK2) correlate with innate resistance and/or tumor recurrence. Immunologically, tumor response after therapy entails T cell receptor repertoire diversification in peripheral blood and intratumoral expansion of preexisting T cell clones. A high ratio of regulatory T to T helper 17 cells in pretreatment blood predicts low T cell receptor repertoire diversity in pretreatment blood, a low cytolytic T cell signature in pretreatment tumors, and innate resistance. Our study provides a molecular framework to advance neoadjuvant anti-PD-1 therapy for individuals with resectable head and neck cancer.

Single cell T cell landscape and T cell receptor repertoire profiling of AML in context of PD-1 blockade therapy

In contrast to the curative effect of allogenic stem cell transplantation in acute myeloid leukemia via T cell activity, only modest responses are achieved with checkpoint-blockade therapy, which might be explained by T cell phenotypes and T cell receptor (TCR) repertoires. Here, we show by paired single-cell RNA analysis and TCR repertoire profiling of bone marrow cells in relapsed/refractory acute myeloid leukemia patients pre/post azacytidine+nivolumab treatment that the disease-related T cell subsets are highly heterogeneous, and their abundance changes following PD-1 blockade-based treatment. TCR repertoires expand and primarily emerge from CD8⁺ cells in patients responding to treatment or having a stable disease, while TCR repertoires contract in therapy^-resistant patients. Trajectory analysis reveals a continuum of CD8⁺ T cell phenotypes, characterized by differential expression of granzyme B and a bone marrow-residing memory CD8⁺ T cell subset, in which a population with stem-like properties expressing granzyme K is enriched in responders. Chromosome 7/7q loss, on the other hand, is a cancer-intrinsic genomic marker of PD-1 blockade resistance in AML. In summary, our study reveals that adaptive T cell plasticity and genomic alterations determine responses to PD-1 blockade in acute myeloid leukemia.

Anti-PD-1/L1 lead-in before MAPK inhibitor combination maximizes antitumor immunity and efficacy

Rationally sequencing and combining PD-1/L1-and MAPK-targeted therapies may overcome innate and acquired resistance. Since increased clinical benefit of MAPK inhibitors (MAPKi) is associated with previous immune checkpoint therapy, we compare the efficacies of sequential and/or combinatorial regimens in subcutaneous murine models of melanoma driven by BrafV600, Nras, or Nf1 mutations as well as colorectal and pancreatic carcinoma driven by KrasG12C. Anti-PD-1/L1 lead-in preceding MAPKi combination optimizes response durability by promoting pro-inflammatory polarization of macrophages and clonal expansion of interferon-γhi, and CD8+ cytotoxic and proliferative (versus CD4+ regulatory) T cells that highly express activation genes. Since therapeutic resistance of melanoma brain metastasis (MBM) limits patient survival, we demonstrate that sequencing anti-PD-1/L1 therapy before MAPKi combination suppresses MBM and improves mouse survival with robust T cell clonal expansion in both intracranial and extracranial metastatic sites. We propose clinically testing brief anti-PD-1/L1 (± anti-CTLA-4) dosing before MAPKi co-treatment to suppress therapeutic resistance.

Spatial immunophenotypes predict response to anti-PD1 treatment and capture distinct paths of T cell evasion in triple negative breast cancer

Only a subgroup of triple-negative breast cancer (TNBC) responds to immune checkpoint inhibitors (ICI). To better understand lack of response to ICI, we analyze 681 TNBCs for spatial immune cell contextures in relation to clinical outcomes and pathways of T cell evasion. Excluded, ignored and inflamed phenotypes can be captured by a gene classifier that predicts prognosis of various cancers as well as anti-PD1 response of metastatic TNBC patients in a phase II trial. The excluded phenotype, which is associated with resistance to anti-PD1, demonstrates deposits of collagen-10, enhanced glycolysis, and activation of TGFβ/VEGF pathways; the ignored phenotype, also associated with resistance to anti-PD1, shows either high density of CD163+ myeloid cells or activation of WNT/PPARγ pathways; whereas the inflamed phenotype, which is associated with response to anti-PD1, revealed necrosis, high density of CLEC9A+ dendritic cells, high TCR clonality independent of neo-antigens, and enhanced expression of T cell co-inhibitory receptors.

Detection of circulating highly expanded T-cell clones in at-risk individuals for rheumatoid arthritis before the clinical onset of the disease

OBJECTIVES

To quantitatively profile the T-cell repertoire in the peripheral blood of individuals genetically at risk for RA, namely first-degree relatives of RA patients (RA-FDR) at different phases of disease development.

METHODS

Next-generation sequencing of the TCR CDR3β repertoire was performed on genomic DNA isolated from whole blood samples of RA-FDR selected at three different pre-clinical stages and of matched RA patients (n = 20/group). T-cell clones were identified by their unique sequence and their degree of expansion (frequency) within each sample was characterized. Clones with a frequency over 0.5% were considered highly expanded clones (HEC).

RESULTS

The absolute number of HEC was significantly higher in established RA patients (mean 4.65) and tended to be higher in symptomatic RA-FDR (mean 3.4) compared with asymptomatic RA-FDR (mean 1.55, P =0.003 and P =0.07, respectively). Asymptomatic individuals with high levels of ACPA did not differ from asymptomatic RA-FDR in terms of absolute number and frequency of clones. The number of HEC tended to be slightly higher at the time of RA onset (P =0.055). Neither clones shared by several patients, nor clones previously associated with RA, were preferentially present within or between the different groups. Finally, a longitudinal analysis did not allow to uncover a kinetic expansion of RA-specific clones closely correlated with disease development.

CONCLUSIONS

HEC were detected in the peripheral blood before the clinical onset of RA, in particular in the later pre-clinical phase of RA development, and their presence increased over time.

Receptor repertoires of murine follicular T helper cells reveal a high clonal overlap in separate lymph nodes in autoimmunity

Follicular T helper cells (Tfh) are a specialized subset of CD4 effector T cells that are crucial for germinal center (GC) reactions and for selecting B cells to undergo affinity maturation. Despite this central role for humoral immunity, only few data exist about their clonal distribution when multiple lymphoid organs are exposed to the same antigen (Ag) as it is the case in autoimmunity. Here, we used an autoantibody-mediated disease model of the skin and injected one auto-Ag into the two footpads of the same mouse and analyzed the T cell receptor (TCR)β sequences of Tfh located in GCs of both contralateral draining lymph nodes. We found that over 90% of the dominant GC-Tfh clonotypes were shared in both lymph nodes but only transiently. The initially dominant Tfh clonotypes especially declined after establishment of chronic disease while GC reaction and autoimmune disease continued. Our data demonstrates a dynamic behavior of Tfh clonotypes under autoimmune conditions and emphasizes the importance of the time point for distinguishing auto-Ag-specific Tfh clonotypes from potential bystander activated ones.

Entinostat induces antitumor immune responses through immune editing of tumor neoantigens

Although immune-checkpoint inhibitors (ICIs) have been a remarkable advancement in bladder cancer treatment, the response rate to single-agent ICIs remains suboptimal. There has been substantial interest in the use of epigenetic agents to enhance ICI efficacy, although precisely how these agents potentiate ICI response has not been fully elucidated. We identified entinostat, a selective HDAC1/3 inhibitor, as a potent antitumor agent in our immune-competent bladder cancer mouse models (BBN963 and BBN966). We demonstrate that entinostat selectively promoted immune editing of tumor neoantigens, effectively remodeling the tumor immune microenvironment, resulting in a robust antitumor response that was cell autonomous, dependent upon antigen presentation, and associated with increased numbers of neoantigen-specific T cells. Finally, combination treatment with anti-PD-1 and entinostat led to complete responses and conferred long-term immunologic memory. Our work defines a tumor cell-autonomous mechanism of action for entinostat and a strong preclinical rationale for the combined use of entinostat and PD-1 blockade in bladder cancer.

A comprehensive model based on temporal dynamics of peripheral T cell repertoire for predicting post-treatment distant metastasis of nasopharyngeal carcinoma

Many nasopharyngeal carcinoma (NPC) patients develop distant metastases after treatment, leading to poor outcomes. To date, there are no peripheral biomarkers suitable for all NPC patients to predict distant metastasis. Hence, we purposed to develop a noninvasive comprehensive model for predicting post-treatment distant metastasis of all NPC. Since T-cell receptor β chain (TCRB) repertoire has achieved prognostic prediction in many cancers, the clinical characteristics and parameters of TCRB repertoire of 71 cases of peripheral blood samples (pairwise pre-treatment and post-treatment samples from 40 NPC patients who without (nM, n = 21) or with (M, n = 19) post-treatment distant metastasis) were collected. The least absolute shrinkage and selection operator algorithm was used to construct a distant metastasis prediction model. In terms of TCRB repertoire parameters, the diversity of TCRB repertoire was significantly decreased in M group after treatment but not in nM group. Ascending TCRB diversity and higher similarity between pre- and post-treatment samples showed better distant metastasis-free survival (DMFS). The similarity still had robust DMFS prediction in patients with reduced TCRB diversity. More importantly, the 5-factor comprehensive model consisting of basic clinical characteristics and TCRB repertoire indices showed a higher prognostic accuracy than any one individual factor in DMFS predicting. In conclusion, treatment had different effects on the composition of TCRB repertoire in patients without and with post-treatment distant metastasis. The dynamics of TCRB diversity, the similarity of TCRB repertoires, and combinations of these factors with basic clinical characteristics could serve as noninvasive DMFS predictors for all NPC patients.

TCR repertoire characteristics predict clinical response to adoptive CTL therapy against nasopharyngeal carcinoma

The past decade has witnessed the gradual and steady progress of adoptive T cell therapy in treating various types of cancer. In combination with gemcitabine and carboplatin chemotherapy, we previously conducted a clinical trial, NCT00690872, to treat Epstein-Barr virus (EBV)-positive nasopharyngeal carcinoma (NPC) patients with autologous EBV-expanded cytotoxic T lymphocytes (CTLs). While achieving a 2-year overall survival rate of 62.9%, this trial failed to induce an anti-tumor response in a sizable fraction of patients. Thus, the identification of benchmarks capable of evaluating CTL products and predicting clinical immunotherapeutic efficacy remains an urgent need. We conducted T cell receptor (TCR) repertoire sequencing to assess EBV-expanded infusion-ready CTL products. To depict the overall repertoire landscape, we evaluated the individual repertoire diversity by Shannon entropy, and, compared the inter-patient CDR3 similarity to estimate T cells expanded by common antigens. With a recently developed bioinformatics algorithm, termed Motif Analysis, we made a machine-learning prediction of structural regions within the CDR3 of TCRβ that associate with CTL therapy prognosis. We found that long term survivors, defined as patients surviving longer than two years, had a higher CTL repertoire diversity with reduced inter-patient similarity. Furthermore, TCR Motif Analysis identified 11 structural motifs distinguishing long term survivors from short term survivors. Specifically, two motifs with a high area under the curve (AUC) values were identified as potential predictive benchmarks for efficacious CTL production. Together, these results reveal that the presence of diverse TCR sequences containing a common core motif set is associated with a favorable response to CTL immunotherapy against EBV-positive NPC.

PD-L1 targeting and subclonal immune escape mediated by PD-L1 mutations in metastatic colorectal cancer

BACKGROUND

In patients with microsatellite stable (MSS) metastatic colorectal cancer (mCRC), immune checkpoint blockade is ineffective, and combinatorial approaches enhancing immunogenicity need exploration.

METHODS

We treated 43 patients with predominantly microsatellite stable RAS/BRAF wild-type mCRC on a phase II trial combining chemotherapy with the epidermal growth factor receptor antibody cetuximab and the programmed cell death ligand 1 (PD-L1) antibody avelumab. We performed next-generation gene panel sequencing for mutational typing of tumors and liquid biopsy monitoring as well as digital droplet PCR to confirm individual mutations. Translational analyses included tissue immunohistochemistry, multispectral imaging and repertoire sequencing of tumor-infiltrating T cells. Detected PD-L1 mutations were mechanistically validated in CRISPR/Cas9-generated cell models using qRT-PCR, immunoblotting, flow cytometry, complement-dependent cytotoxicity assay, antibody-dependent cytotoxicity by natural killer cell degranulation assay and LDH release assay as well as live cell imaging of T cell mediated tumor cell killing.

RESULTS

Circulating tumor DNA showed rapid clearance in the majority of patients mirroring a high rate of early tumor shrinkage. In 3 of 13 patients expressing the high-affinity Fcγ receptor 3a (FcγR3a), tumor subclones with PD-L1 mutations were selected that led to loss of tumor PD-L1 by nonsense-mediated RNA decay in PD-L1 K162fs and protein degradation in PD-L1 L88S. As a consequence, avelumab binding and antibody-dependent cytotoxicity were impaired, while T cell killing of these variant clones was increased. Interestingly, PD-L1 mutant subclones showed slow selection dynamics reversing on avelumab withdrawal and patients with such subclones had above-average treatment benefit. This suggested that the PD-L1 mutations mediated resistance to direct antitumor effects of avelumab, while at the same time loss of PD-L1 reduced biological fitness by enhanced T cell killing limiting subclonal expansion.

CONCLUSION

The addition of avelumab to standard treatment appeared feasible and safe. PD-L1 mutations mediate subclonal immune escape to avelumab in some patients with mCRC expressing high-affinity FcγR3a, which may be a subset experiencing most selective pressure. Future trials evaluating the addition of avelumab to standard treatment in MSS mCRC are warranted especially in this patient subpopulation.

TRIAL REGISTRATION NUMBER

NCT03174405.