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

scRepertoire: An R-based toolkit for single-cell immune receptor analysis

Single-cell sequencing is an emerging technology in the field of immunology and oncology that allows researchers to couple RNA quantification and other modalities, like immune cell receptor profiling at the level of an individual cell. A number of workflows and software packages have been created to process and analyze single-cell transcriptomic data. These packages allow users to take the vast dimensionality of the data generated in single-cell-based experiments and distill the data into novel insights. Unlike the transcriptomic field, there is a lack of options for software that allow for single-cell immune receptor profiling. Enabling users to easily combine mRNA and immune profiling, scRepertoire was built to process data derived from 10x Genomics Chromium Immune Profiling for both T-cell receptor (TCR) and immunoglobulin (Ig) enrichment workflows and subsequently interacts with a number of popular R packages for single-cell expression, such as Seurat. The scRepertoire R package and processed data are open source and available on GitHub and provides in-depth tutorials on the capability of the package.

scRepertoire: An R-based toolkit for single-cell immune receptor analysis

Single-cell sequencing is an emerging technology in the field of immunology and oncology that allows researchers to couple RNA quantification and other modalities, like immune cell receptor profiling at the level of an individual cell. A number of workflows and software packages have been created to process and analyze single-cell transcriptomic data. These packages allow users to take the vast dimensionality of the data generated in single-cell-based experiments and distill the data into novel insights. Unlike the transcriptomic field, there is a lack of options for software that allow for single-cell immune receptor profiling. Enabling users to easily combine mRNA and immune profiling, scRepertoire was built to process data derived from 10x Genomics Chromium Immune Profiling for both T-cell receptor (TCR) and immunoglobulin (Ig) enrichment workflows and subsequently interacts with the popular Seurat R package. The scRepertoire R package and processed data are open source and available on GitHub and provides in-depth tutorials on the capability of the package.

Staphylococcus aureus Isolated from Skin from Atopic-Dermatitis Patients Produces Staphylococcal Enterotoxin Y, Which Predominantly Induces T-Cell Receptor Vα-Specific Expansion of T Cells

While investigating the virulence traits of Staphylococcus aureus adhering to the skin of atopic-dermatitis (AD) patients, we identified a novel open reading frame (ORF) with structural similarity to a superantigen from genome sequence data of an isolate from AD skin. Concurrently, the same ORF was identified in a bovine isolate of S. aureus and designated SElY (H. K. Ono, Y. Sato'o, K. Narita, I. Naito, et al., Appl Environ Microbiol 81:7034-7040, 2015, https://doi.org/10.1128/AEM.01873-15). Recombinant SElY_bov had superantigen activity in human peripheral blood mononuclear cells. It further demonstrated emetic activity in a primate animal model, and it was proposed that SElY be renamed SEY (H. K. Ono, S. Hirose, K. Narita, M. Sugiyama, et al., PLoS Pathog 15:e1007803, 2019, https://doi.org/10.1371/journal.ppat.1007803). Here, we investigated the prevalence of the sey gene in 270 human clinical isolates of various origins in Japan. Forty-two strains were positive for the sey gene, and the positive isolates were from patients with the skin diseases atopic dermatitis and impetigo/staphylococcal scalded skin syndrome (SSSS), with a detection rate of ∼17 to 22%. There were three variants of SEY (SEY₁, SEY₂, and SEY₃), and isolates producing SEY variants formed three distinct clusters corresponding to clonal complexes (CCs) 121, 59, and 20, respectively. Most sey ⁺ isolates produced SEY in broth culture. Unlike SEY_bov, the three recombinant SEY variants exhibited stability against heat treatment. SEY predominantly activated human T cells with a particular T-cell receptor (TCR) Vα profile, a unique observation since most staphylococcal enterotoxins exert their superantigenic activities through activating T cells with specific TCR Vβ profiles. SEY may act to induce localized inflammation via skin-resident T-cell activation, facilitating the pathogenesis of S. aureus infection in disrupted epithelial barriers.

T Cell Receptor Repertoire Sequencing

The status of T cell receptors (TCRs) repertoire is associated with the occurrence and progress of various diseases and can be used in monitoring the immune responses, predicting the prognosis of disease and other medical fields. High-throughput sequencing promotes the studying in TCR repertoire. The chapter focuses on the whole process of TCR profiling, including DNA extraction, library construction, high-throughput sequencing, and how to analyze data.

NKp46-expressing human gut-resident intraepithelial Vδ1 T cell subpopulation exhibits high antitumor activity against colorectal cancer

γδ T cells account for a large fraction of human intestinal intraepithelial lymphocytes (IELs) endowed with potent antitumor activities. However, little is known about their origin, phenotype, and clinical relevance in colorectal cancer (CRC). To determine γδ IEL gut specificity, homing, and functions, γδ T cells were purified from human healthy blood, lymph nodes, liver, skin, and intestine, either disease-free, affected by CRC, or generated from thymic precursors. The constitutive expression of NKp46 specifically identifies a subset of cytotoxic Vδ1 T cells representing the largest fraction of gut-resident IELs. The ontogeny and gut-tropism of NKp46+/Vδ1 IELs depends both on distinctive features of Vδ1 thymic precursors and gut-environmental factors. Either the constitutive presence of NKp46 on tissue-resident Vδ1 intestinal IELs or its induced expression on IL-2/IL-15-activated Vδ1 thymocytes are associated with antitumor functions. Higher frequencies of NKp46+/Vδ1 IELs in tumor-free specimens from CRC patients correlate with a lower risk of developing metastatic III/IV disease stages. Additionally, our in vitro settings reproducing CRC tumor microenvironment inhibited the expansion of NKp46+/Vδ1 cells from activated thymic precursors. These results parallel the very low frequencies of NKp46+/Vδ1 IELs able to infiltrate CRC, thus providing insights to either follow-up cancer progression or to develop adoptive cellular therapies.

Histone deacetylase inhibition promotes intratumoral CD8+ T-cell responses, sensitizing murine breast tumors to anti-PD1

Histone deacetylase (HDAC) inhibitors impair tumor cell proliferation and alter gene expression. However, the impact of these changes on anti-tumor immunity is poorly understood. Here, we showed that the class I HDAC inhibitor, entinostat (ENT), promoted the expression of immune-modulatory molecules, including MHCII, costimulatory ligands, and chemokines on murine breast tumor cells in vitro and in vivo. ENT also impaired tumor growth in vivo-an effect that was dependent on both CD8⁺ T cells and IFNγ. Moreover, ENT promoted intratumoral T-cell clonal expansion and enhanced their functional activity. Importantly, ENT sensitized normally unresponsive tumors to the effects of PD1 blockade, predominantly through increases in T-cell proliferation. Our findings suggest that class I HDAC inhibitors impair tumor growth by enhancing the proliferative and functional capacity of CD8⁺ T cells and by sensitizing tumor cells to T-cell recognition.

A High-Avidity T-cell Receptor Redirects Natural Killer T-cell Specificity and Outcompetes the Endogenous Invariant T-cell Receptor

T-cell receptor (TCR) gene transfer redirects T cells to target intracellular antigens. However, the potential autoreactivity generated by TCR mispairing and occurrence of graft-versus-host disease in the allogenic setting due to the retention of native TCRs remain major concerns. Natural killer T cells (NKT) have shown promise as a platform for adoptive T-cell therapy in cancer patients. Here, we showed their utility for TCR gene transfer. We successfully engineered and expanded NKTs expressing a functional TCR (TCR NKTs), showing HLA-restricted antitumor activity in xenogeneic mouse models in the absence of graft-versus-mouse reactions. We found that TCR NKTs downregulated the invariant TCR (iTCR), leading to iTCR⁺TCR⁺ and iTCR^-TCR⁺ populations. In-depth analyses of these subsets revealed that in iTCR^-TCR⁺ NKTs, the iTCR, although expressed at the mRNA and protein levels, was retained in the cytoplasm. This effect resulted from a competition for binding to CD3 molecules for cell-surface expression by the transgenic TCR. Overall, our results highlight the feasibility and advantages of using NKTs for TCR expression for adoptive cell immunotherapies. NKT-low intrinsic alloreactivity that associated with the observed iTCR displacement by the engineered TCR represents ideal characteristics for "off-the-shelf" products without further TCR gene editing.

sumrep: A Summary Statistic Framework for Immune Receptor Repertoire Comparison and Model Validation

The adaptive immune system generates an incredible diversity of antigen receptors for B and T cells to keep dangerous pathogens at bay. The DNA sequences coding for these receptors arise by a complex recombination process followed by a series of productivity-based filters, as well as affinity maturation for B cells, giving considerable diversity to the circulating pool of receptor sequences. Although these datasets hold considerable promise for medical and public health applications, the complex structure of the resulting adaptive immune receptor repertoire sequencing (AIRR-seq) datasets makes analysis difficult. In this paper we introduce sumrep, an R package that efficiently performs a wide variety of repertoire summaries and comparisons, and show how sumrep can be used to perform model validation. We find that summaries vary in their ability to differentiate between datasets, although many are able to distinguish between covariates such as donor, timepoint, and cell type for BCR and TCR repertoires. We show that deletion and insertion lengths resulting from V(D)J recombination tend to be more discriminative characterizations of a repertoire than summaries that describe the amino acid composition of the CDR3 region. We also find that state-of-the-art generative models excel at recapitulating gene usage and recombination statistics in a given experimental repertoire, but struggle to capture many physiochemical properties of real repertoires.

Germline-Encoded Affinity for Cognate Antigen Enables Vaccine Amplification of a Human Broadly Neutralizing Response against Influenza Virus

Antibody paratopes are formed by hypervariable complementarity-determining regions (CDRH3s) and variable gene-encoded CDRs. The latter show biased usage in human broadly neutralizing antibodies (bnAbs) against both HIV and influenza virus, suggesting the existence of gene-endowed targeting solutions that may be amenable to pathway amplification. To test this, we generated transgenic mice with human CDRH3 diversity but simultaneously constrained to individual user-defined human immunoglobulin variable heavy-chain (VH) genes, including IGHV1-69, which shows biased usage in human bnAbs targeting the hemagglutinin stalk of group 1 influenza A viruses. Sequential immunization with a stalk-only hemagglutinin nanoparticle elicited group 1 bnAbs, but only in IGHV1-69 mice. This VH-endowed response required minimal affinity maturation, was elicited alongside pre-existing influenza immunity, and when IGHV1-69 B cells were diluted to match the frequency measured in humans. These results indicate that the human repertoire could, in principle, support germline-encoded bnAb elicitation using a single recombinant hemagglutinin immunogen.

Multiplexed activation of endogenous genes by CRISPRa elicits potent antitumor immunity

Immunotherapy has transformed cancer treatment. However, current immunotherapy modalities face various limitations. Here, we developed MAEGI, a new form of immunotherapy that elicits anti-tumor immunity through multiplexed activation of endogenous genes. We leveraged CRISPR activation (CRISPRa) to directly augment the in situ expression of endogenous genes, thereby the presentation of tumor antigens, leading to dramatic anti-tumor immune responses. Deploying this as a cell-based vaccination strategy showed efficacy in both prophylactic and therapeutic settings. Intratumoral adeno-associated virus delivery of CRISPRa libraries elicited strong anti-tumor immunity across multiple cancer types. Precision targeting of mutated gene sets eradicated a large fraction of established tumors at both local and distant sites. This treatment modality led to alterations of the tumor microenvironment, marked by enhanced T cell infiltration and anti-tumor immune signatures. Multiplexed endogenous gene activation is a versatile and highly scalable strategy to elicit potent immune responses against cancer, distinct from existing cancer therapies.

Analysis of TCR β CDR3 sequencing data for tracking anti-tumor immunity

Anti-tumor T cells are the soldiers in the body's war against cancer. Effector T cells can detect and eliminate cells expressing their cognate antigen via activation through engagement of the T cell receptor (TCR) with its cognate peptide:MHC complex. Owing to the recent success of immunotherapy in the treatment of many different cancer types, research efforts have shifted toward identifying and tracking anti-tumor T cell responses upon treatment in cancer patients. While traditional methods, such as ELISpot and flow cytometric intracellular staining have had limited success, likely owing to the inability to get viable biospecimens or the lower magnitude of tumor-specific T cell responses relative to virus-specific responses, new techniques that utilize next generation sequencing enable T cell response tracking independent of cytokine production or cell viability. The TCR, which confers T cell antigen-specificity, can be used as a molecular barcode to track T cell clonotypic dynamics across biological compartments and over time in cancer patients undergoing treatment. Because this method does not require viable cells, these T cell clonotypes can also be tracked in archival tumor tissue and flash frozen cell pellets. While exciting, quantitative TCR sequencing (TCRseq) technologies have been met with the conundrum of how to properly analyze and interpret the data. Here we provide a comprehensive guide on how to acquire, analyze, and interpret TCRseq data, as well as special considerations that should be taken prior to experimental setup.

Hepatitis E Virus (HEV)-Specific T Cell Receptor Cross-Recognition: Implications for Immunotherapy

T cell immunotherapy is a concept developed for the treatment of cancer and infectious diseases, based on cytotoxic T lymphocytes to target tumor- or pathogen-specific antigens. Antigen-specificity of the T cell receptors (TCRs) is an important selection criterion in the developmental design of immunotherapy. However, off-target specificity is a possible autoimmunity concern if the engineered antigen-specific T cells are cross-reacting to self-peptides in-vivo. In our recent work, we identified several hepatitis E virus (HEV)-specific TCRs as potential candidates to be developed into T cell therapy to treat chronic hepatitis E. One of the identified TCRs, targeting a HLA-A2-restricted epitope at the RNA-dependent RNA polymerase (HEV-1527: LLWNTVWNM), possessed a unique multiple glycine motif in the TCR-β CDR3, which might be a factor inducing cross-reactivity. The aim of our study was to explore if this TCR could cross-recognize self-peptides to underlay autoimmunity. Indeed, we found that this HEV-1527-specific TCR could also cross-recognize an apoptosis-related epitope, Nonmuscle Myosin Heavy Chain 9 (MYH9-478: QLFNHTMFI). While this TCR had dual specificities to both viral epitope and a self-antigen by double Dextramer binding, it was selectively functional against HEV-1527 but not activated against MYH9-478. The consecutive glycine motif in β chain may be the reason promoting TCR binding promiscuity to recognize a secondary target, thereby facilitating cross-recognition. In conclusion, candidate TCRs for immunotherapy development should be screened for autoimmune potential, especially when the TCRs exhibit unique sequence pattern.

TCR repertoire and CDR3 motif analyses depict the role of αβ T cells in Ankylosing spondylitis

Background

Ankylosing spondylitis (AS) is a chronic inflammatory disease with worldwide high prevalence. Although AS is strongly associated with HLA-B27 MHC-I antigen presentation, the role played by αβ T cells in AS remains elusive.

Methods

Utilizing TCRβ repertoire sequencing and bioinformatics tools developed in house, we analyzed overall TCR repertoire structures and antigen-recognizing CDR3 motifs in AS patients with different disease activities.

Findings

We found that disease progression is associated with both CD4+ and CD8+ T cell oligo-clonal expansion, which suggests that αβ T cell activation may mediate AS disease progression. By developing a bioinformatics platform to dissect antigen-specific responses, we discovered a cell population consisting of both CD4+ and CD8+ T cells expressing identical TCRs, herein termed CD4/8 T cells. CD4/8 clonotypes were highly enriched in the spondyloarthritic joint fluid of patients, and their expansion correlated with the activity of disease.

Interpretation

These results provide evidence on the T cell clone side to reveal the potential role of CD4/8 T cells in the etiology of AS development.

High-Throughput Immunogenetics Reveals a Lack of Physiological T Cell Clusters in Patients With Autoimmune Cytopenias

Autoimmune cytopenias (AIC) such as immune thrombocytopenia or autoimmune hemolytic anemia are claimed to be essentially driven by a dysregulated immune system. Using next-generation immunosequencing we profiled 59 T and B cell repertoires (TRB and IGH) of 25 newly diagnosed patients with primary or secondary (lymphoma-associated) AIC to test the hypothesis if these patients present a disease-specific immunological signature that could reveal pathophysiological clues and eventually be exploited as blood-based biomarker. Global TRB and IGH repertoire metrics as well as VJ gene usage distribution showed uniform characteristics for all lymphoma patients (high clonality and preferential usage of specific TRBV- and TRBJ genes), but no AIC-specific signature. Since T cell immune reactions toward antigens are unique and polyclonal, we clustered TCRβ clones in-silico based on target recognition using the GLIPH (grouping of lymphocyte interactions by paratope hotspots) algorithm. This analysis revealed a considerable lack of physiological T cell clusters in patients with primary AIC. Interestingly, this signature did not discriminate between the different subentities of AIC and was also found in an independent cohort of 23 patients with active autoimmune hepatitis. Taken together, our data suggests that the identified T cell cluster signature could represent a blood biomarker of autoimmune conditions in general and should be functionally validated in future studies.

Analysis of the TCR Repertoire in HIV-Exposed but Uninfected Infants

Maternal human immunodeficiency virus (HIV) infection has been shown to leave profound and lasting impacts on the HIV-exposed uninfected (HEU) infant, including increased mortality and morbidity, immunological changes, and developmental delays compared to their HIV-unexposed (HU) counterparts. Exposure to HIV or antiretroviral therapy may influence immune development, which could increase morbidity and mortality. However, a direct link between the increased mortality and morbidity and the infant's immune system has not been identified. To provide a global picture of the neonatal T cell repertoire in HEU versus HU infants, the diversity of the T cell receptor beta chain (TRB) expressed in cord blood samples from HEU infants was determined using next-generation sequencing and compared to healthy (HU) infants collected from the same community. While the TRB repertoire of HU infants was broadly diverse, in line with the expected idea of a naïve T cell repertoire, samples of HEU infants showed a significantly reduced TRB diversity. This study is the first to demonstrate differences in TRB diversity between HEU and HU cord blood samples and provides evidence that maternal HIV, in the absence of transmission, influences the adaptive immune system of the unborn child.

Polyfunctional tumor-reactive T cells are effectively expanded from non-small cell lung cancers, and correlate with an immune-engaged T cell profile

Non-small cell lung cancer (NSCLC) is the second most prevalent type of cancer. With the current treatment regimens, the mortality rate remains high. Therefore, better therapeutic approaches are necessary. NSCLCs generally possess many genetic mutations and are well infiltrated by T cells (TIL), making TIL therapy an attractive option. Here we show that T cells from treatment naive, stage I-IVa NSCLC tumors can effectively be isolated and expanded, with similar efficiency as from normal lung tissue. Importantly, 76% (13/17) of tested TIL products isolated from NSCLC lesions exhibited clear reactivity against primary tumor digests, with 0.5%-30% of T cells producing the inflammatory cytokine Interferon (IFN)-γ. Both CD4⁺ and CD8⁺ T cells displayed tumor reactivity. The cytokine production correlated well with CD137 and CD40L expression. Furthermore, almost half (7/17) of the TIL products contained polyfunctional T cells that produced Tumor Necrosis Factor (TNF)-α and/or IL-2 in addition to IFN-γ, a hallmark of effective immune responses. Tumor-reactivity in the TIL products correlated with high percentages of CD103⁺CD69⁺CD8⁺ T cell infiltrates in the tumor lesions, with PD-1^hiCD4⁺ T cells, and with FoxP3⁺CD25⁺CD4⁺ regulatory T cell infiltrates, suggesting that the composition of T cell infiltrates may predict the level of tumor reactivity. In conclusion, the effective generation of tumor-reactive and polyfunctional TIL products implies that TIL therapy will be a successful treatment regimen for NSCLC patients.

Single-Cell Transcriptome Analysis of T Cells

Single-cell RNA-seq (scRNA-seq) has provided novel routes to investigate the heterogeneous populations of T cells and is rapidly becoming a common tool for molecular profiling and identification of novel subsets and functions. This chapter offers an experimental and computational workflow for scRNA-seq analysis of T cells. We focus on the analyses of scRNA-seq data derived from plate-based sorted T cells using flow cytometry and full-length transcriptome protocols such as Smart-Seq2. However, the proposed pipeline can be applied to other high-throughput approaches such as UMI-based methods. We describe a detailed bioinformatics pipeline that can be easily reproduced and discuss future directions and current limitations of these methods in the context of T cell biology.

BRepertoire: a user-friendly web server for analysing antibody repertoire data

Antibody repertoire analysis by high throughput sequencing is now widely used, but a persisting challenge is enabling immunologists to explore their data to discover discriminating repertoire features for their own particular investigations. Computational methods are necessary for large-scale evaluation of antibody properties. We have developed BRepertoire, a suite of user-friendly web-based software tools for large-scale statistical analyses of repertoire data. The software is able to use data preprocessed by IMGT, and performs statistical and comparative analyses with versatile plotting options. BRepertoire has been designed to operate in various modes, for example analysing sequence-specific V(D)J gene usage, discerning physico-chemical properties of the CDR regions and clustering of clonotypes. Those analyses are performed on the fly by a number of R packages and are deployed by a shiny web platform. The user can download the analysed data in different table formats and save the generated plots as image files ready for publication. We believe BRepertoire to be a versatile analytical tool that complements experimental studies of immune repertoires. To illustrate the server’s functionality, we show use cases including differential gene usage in a vaccination dataset and analysis of CDR3H properties in old and young individuals. The server is accessible under http://mabra.biomed.kcl.ac.uk/BRepertoire.

IL-7-dependent compositional changes within the γδ T cell pool in lymph nodes during ageing lead to an unbalanced anti-tumour response

How the age-associated decline of immune function leads to increased cancer incidence is poorly understood. Here, we have characterised the cellular composition of the γδ T-cell pool in peripheral lymph nodes (pLNs) upon ageing. We find that ageing has minimal cell-intrinsic effects on function and global gene expression of γδ T cells, and γδTCR diversity remains stable. However, ageing alters TCRδ chain usage and clonal structure of γδ T-cell subsets. Importantly, IL-17-producing γδ17 T cells dominate the γδ T-cell pool of aged mice-mainly due to the selective expansion of Vγ6+ γδ17 T cells and augmented γδ17 polarisation of Vγ4+ T cells. Expansion of the γδ17 T-cell compartment is mediated by increased IL-7 expression in the T-cell zone of old mice. In a Lewis lung cancer model, pro-tumourigenic Vγ6+ γδ17 T cells are exclusively activated in the tumour-draining LN and their infiltration into the tumour correlates with increased tumour size in aged mice. Thus, upon ageing, substantial compositional changes in γδ T-cell pool in the pLN lead to an unbalanced γδ T-cell response in the tumour that is associated with accelerated tumour growth.

T cell receptor next-generation sequencing reveals cancer-associated repertoire metrics and reconstitution after chemotherapy in patients with hematological and solid tumors

The dynamics of immunoaging and the onset of immunoparesis in healthy individuals and cancer patients has been controversially discussed. Moreover, the role of chemotherapy on T cell regeneration needs further elucidation in light of novel immunotherapies that have become standard of care for many elderly cancer patients. We used next-generation immunosequencing to study T cell receptor (TCR) repertoire metrics on 346 blood samples from healthy individuals and cancer patients producing a dataset with around 8.8 million TCR reads. This analysis showed that decline of T cell diversity and increase in T cell clonality is a continuous process beginning in healthy individuals over 40 years of age. Untreated patients with both hematological and solid tumors showed blood TCR repertoires with significantly lower diversity and higher clonality as compared to healthy individuals across all decades. Loss in T cell diversity was essentially driven by a loss in richness in aging healthy individuals, while in cancer patients a loss in repertoire evenness was an additional contributing factor. Interestingly, chemotherapy did not impair the regeneration of blood TCR repertoire diversity to pre-treatment age-specific levels. Surprisingly, even patients over the age of 70 years receiving highly T cell toxic therapies reestablished their pre-treatment T cell diversity suggesting rebound thymic activity rather than recovery of T cell counts by peripheral expansion only. Taken together, these data suggest that human TCR repertoire metrics gradually deteriorate in the aging individual, but age-specific TCR metrics are restored after T cell depleting therapy even in elderly cancer patients.

Resident memory CD8 T cells persist for years in human small intestine

Bartolomé-Casado et al. demonstrate that human gut contains large numbers of resident memory CD8 T cells that survive for years. Intestinal CD8 Trm cells have a clonally expanded immune repertoire that is stable over time and exhibit enhanced protective capabilities.

Resident memory CD8 T (Trm) cells have been shown to provide effective protective responses in the small intestine (SI) in mice. A better understanding of the generation and persistence of SI CD8 Trm cells in humans may have implications for intestinal immune-mediated diseases and vaccine development. Analyzing normal and transplanted human SI, we demonstrated that the majority of SI CD8 T cells were bona fide CD8 Trm cells that survived for >1 yr in the graft. Intraepithelial and lamina propria CD8 Trm cells showed a high clonal overlap and a repertoire dominated by expanded clones, conserved both spatially in the intestine and over time. Functionally, lamina propria CD8 Trm cells were potent cytokine producers, exhibiting a polyfunctional (IFN-γ⁺ IL-2⁺ TNF-α⁺) profile, and efficiently expressed cytotoxic mediators after stimulation. These results suggest that SI CD8 Trm cells could be relevant targets for future oral vaccines and therapeutic strategies for gut disorders.

Bioinformatic methods for cancer neoantigen prediction

Tumor cells accumulate aberrations not present in normal cells, leading to presentation of neoantigens on MHC molecules on their surface. These non-self neoantigens distinguish tumor cells from normal cells to the immune system and are thus targets for cancer immunotherapy. The rapid development of molecular profiling platforms, such as next-generation sequencing, has enabled the generation of large datasets characterizing tumor cells. The simultaneous development of algorithms has enabled rapid and accurate processing of these data. Bioinformatic software tools encoding the algorithms can be strung together in a workflow to identify neoantigens. Here, with a focus on high-throughput sequencing, we review state-of-the art bioinformatic tools along with the steps and challenges involved in neoantigen identification and recognition.

HLA class II-Restricted CD8+ T cells in HIV-1 Virus Controllers

A paradigm shifting study demonstrated that induction of MHC class E and II-restricted CD8+ T cells was associated with the clearance of SIV infection in rhesus macaques. Another recent study highlighted the presence of HIV-1-specific class II-restricted CD8+ T cells in HIV-1 patients who naturally control infection (virus controllers; VCs). However, questions regarding class II-restricted CD8+ T cells ontogeny, distribution across different HIV-1 disease states and their role in viral control remain unclear. In this study, we investigated the distribution and anti-viral properties of HLA-DRB1*0701 and DQB1*0501 class II-restricted CD8+ T cells in different HIV-1 patient cohorts; and whether class II-restricted CD8+ T cells represent a unique T cell subset. We show that memory class II-restricted CD8+ T cell responses were more often detectable in VCs than in chronically infected patients, but not in healthy seronegative donors. We also demonstrate that VC CD8+ T cells inhibit virus replication in both a class I- and class II-dependent manner, and that in two VC patients the class II-restricted CD8+ T cells with an anti-viral gene signature expressed both CD4+ and CD8+ T cell lineage-specific genes. These data demonstrated that anti-viral memory class II-restricted CD8+ T cells with hybrid CD4+ and CD8+ features are present during natural HIV-1 infection.

High-throughput analysis of the human thymic Vδ1+ T cell receptor repertoire

γδ T cells are a relatively rare subset of lymphocytes in the human peripheral blood, but they play important roles at the interface between the innate and the adaptive immune systems. The γδ T cell lineage is characterized by a signature γδ T cell receptor (γδTCR) that displays extensive sequence variability originated by DNA rearrangement of the corresponding V(D)J loci. Human γδ T cells comprise Vγ9Vδ2 T cells, the major subset in the peripheral blood; and Vδ1⁺ T cells, the predominant subpopulation in the post-natal thymus and in peripheral tissues. While less studied, Vδ1⁺ T cells recently gathered significant attention due to their anti-cancer and anti-viral activities. In this study we applied next-generation sequencing (NGS) to analyse the γδTCR repertoire of highly (FACS-)purified Vδ1⁺ T cells from human thymic biopsies. Our analysis reveals unsuspected aspects of thymically rearranged and expressed (at the mRNA level) TRG and TRD genes, thus constituting a data resource that qualifies previous conclusions on the TCR repertoire of γδ T cells developing in the human thymus.

Design Type(s)	transcription profiling design • cellular data analysis objective
Measurement Type(s)	gamma delta T-lymphocyte
Technology Type(s)	RNA sequencing
Factor Type(s)	sex • gene
Sample Characteristic(s)	Homo sapiens • thymus

Machine-accessible metadata file describing the reported data (ISA-Tab format)

Application of immune repertoire sequencing in cancer immunotherapy

With the prominent breakthrough in the field of tumor immunology, diverse cancer immunotherapies have attracted great attention in the last decade. The immune checkpoint inhibitors, adoptive cell therapies, and therapeutic cancer vaccines have already achieved impressive clinical success. However, the fact that only a small subset of patients with specific tumor types can benefit from these treatments limits the application of cancer immunotherapy. To seek out the molecular mechanisms behind this challenge and to select cancer precision medicine for different individuals, researchers apply the immune repertoire sequencing (IRS) to evaluate genetic responses of each patient to current immunotherapies. This review summarizes the technical advances and recent applications of IRS in cancer immunotherapy, indicates the limitations of this technique, and predicts future perspectives both in basic studies and clinical trials.

RNA sequencing identifies clonal structure of T-cell repertoires in patients with adult T-cell leukemia/lymphoma

The diversity of T-cell receptor (TCR) repertoires, as generated by somatic DNA rearrangements, is central to immune system function. High-throughput sequencing technologies now allow examination of antigen receptor repertoires at single-nucleotide and, more recently, single-cell resolution. The TCR repertoire can be altered in the context of infections, malignancies or immunological disorders. Here we examined the diversity of TCR clonality and its association with pathogenesis and prognosis in adult T-cell leukemia/lymphoma (ATL), a malignancy caused by infection with human T-cell leukemia virus type-1 (HTLV-1). We analyzed 62 sets of high-throughput RNA sequencing data from 59 samples of HTLV-1−infected individuals—asymptomatic carriers (ACs), smoldering, chronic, acute and lymphoma ATL subtypes—and three uninfected controls to evaluate TCR distribution. Based on these TCR profiles, CD4-positive cells and ACs showed polyclonal patterns, whereas ATL patients showed oligo- or monoclonal patterns (with 446 average clonotypes across samples). Expression of TCRα and TCRβ genes in the dominant clone differed among the samples. ACs, CD4-positive samples and smoldering patients showed significantly higher TCR diversity compared with chronic, acute and lymphoma subtypes. CDR3 sequence length distribution, amino acid conservation and gene usage variability for ATL patients resembled those of peripheral blood cells from ACs and healthy donors. Thus, determining monoclonal architecture and clonal diversity by RNA sequencing might be useful for prognostic purposes and for personalizing ATL diagnosis and assessment of treatments.

LKB1 expressed in dendritic cells governs the development and expansion of thymus-derived regulatory T cells

Liver Kinase B1 (LKB1) plays a key role in cellular metabolism by controlling AMPK activation. However, its function in dendritic cell (DC) biology has not been addressed. Here, we find that LKB1 functions as a critical brake on DC immunogenicity, and when lost, leads to reduced mitochondrial fitness and increased maturation, migration, and T cell priming of peripheral DCs. Concurrently, loss of LKB1 in DCs enhances their capacity to promote output of regulatory T cells (Tregs) from the thymus, which dominates the outcome of peripheral immune responses, as suggested by increased resistance to asthma and higher susceptibility to cancer in CD11c^ΔLKB1 mice. Mechanistically, we find that loss of LKB1 specifically primes thymic CD11b⁺ DCs to facilitate thymic Treg development and expansion, which is independent from AMPK signalling, but dependent on mTOR and enhanced phospholipase C β1-driven CD86 expression. Together, our results identify LKB1 as a critical regulator of DC-driven effector T cell and Treg responses both in the periphery and the thymus.

Preferential Infiltration of Unique Vγ9Jγ2-Vδ2 T Cells Into Glioblastoma Multiforme

Glioblastoma multiforme (GBM) is clinically highly aggressive as a result of evolutionary dynamics induced by cross-talk between cancer cells and a heterogeneous group of immune cells in tumor microenvironment. The brain harbors limited numbers of immune cells with few lymphocytes and macrophages; thus, innate-like lymphocytes, such as γδ T cells, have important roles in antitumor immunity. Here, we characterized GBM-infiltrating γδ T cells, which may have roles in regulating the GBM tumor microenvironment and cancer cell gene expression. V(D)J repertoires of tumor-infiltrating and blood-circulating γδ T cells from four patients were analyzed by next-generation sequencing-based T-cell receptor (TCR) sequencing in addition to mutation and immune profiles in four GBM cases. In all tumor tissues, abundant innate and effector/memory lymphocytes were detected, accompanied by large numbers of tumor-associated macrophages and closely located tumor-infiltrating γδ T cells, which appear to have anti-tumor activity. The immune-related gene expression analysis using the TCGA database showed that the signature gene expression extent of γδ T cells were more associated with those of cytotoxic T and Th1 cells and M1 macrophages than those of Th2 cells and M2 macrophages. Although the most abundant γδ T cells were Vγ9Vδ2 T cells in both tumor tissues and blood, the repertoire of intratumoral Vγ9Vδ2 T cells was distinct from that of peripheral blood Vγ9Vδ2 T cells and was dominated by Vγ9Jγ2 sequences, not by canonical Vγ9JγP sequences that are mostly commonly found in blood γδ T cells. Collectively, unique GBM-specific TCR clonotypes were identified by comparing TCR repertoires of peripheral blood and intra-tumoral γδ T cells. These findings will be helpful for the elucidation of tumor-specific antigens and development of anticancer immunotherapies using tumor-infiltrating γδ T cells.

Dengue-specific CD8+ T cell subsets display specialized transcriptomic and TCR profiles

Accumulating evidence demonstrates that CD8+ T cells contribute to protection from severe dengue virus (DENV) disease and vaccine efficacy. Nevertheless, molecular programs associated with DENV-specific CD8+ T cell subsets have not been defined. Here, we studied the transcriptomic profiles of human DENV-specific CD8+ T cells isolated after stimulation with DENV epitopes from donors who had been infected with DENV multiple times and would therefore be expected to have significant levels of adaptive immunity. We found that DENV-specific CD8+ T cells mainly consisted of effector memory subsets, namely CD45RA-CCR7- effector memory (Tem) and CD45RA+CCR7- effector memory re-expressing CD45RA (Temra) cells, which enacted specific gene expression profiles upon stimulation with cognate antigens. DENV-specific CD8+ T cell subsets in general, and Temra cells in particular, were fully activated and polyfunctional, yet associated with relatively narrow transcriptional responses. Furthermore, we found that DENV-specific CD8+ Tem and Temra cells showed some unique T cell receptor features in terms of overlap and variable (V) gene usage. This study provides a transcriptomic definition of DENV-specific activated human CD8+ T cell subsets and defines a benchmark profile that vaccine-specific responses could aim to reproduce.

Comprehensive Analysis of TCR-β Repertoire in Patients with Neurological Immune-mediated Disorders

In this study we characterized the TCR repertoire profiles in patients with chronic progressive inflammatory neurological disorders including HAM/TSP, associated with human T-cell lymphotropic virus type I (HTLV-I) infection, and multiple sclerosis (MS), an inflammatory, demyelinating disease of the CNS of unknown etiology. We hypothesized that a T-cell receptor (TCR) clonal repertoire ‘signature’ could distinguish HAM/TSP patients from healthy controls, as well as from patients with a more heterogeneous CNS-reactive inflammatory disease such as MS. In this study, we applied an unbiased molecular technique – unique molecular identifier (UMI) library-based strategy to investigate with high accuracy the TCR clonal repertoire by high throughput sequencing (HTS) technology. cDNA-TCR β-chain libraries were sequenced from 2 million peripheral mononuclear cells (PBMCs) in 14 HAM/TSP patients, 34 MS patients and 20 healthy controls (HC). While HAM/TSP patients showed a higher clonal T-cell expansion compared to MS and HC, increase of the TCR clonal expansion was inversely correlated with the diversity of TCR repertoire in all subjects. In addition, longitudinal analysis of TCR repertoires from HAM/TSP patients demonstrated a correlation of the TCR clonal expansion with HTLV-I proviral load. Surprisingly, MS patients showed a higher diversity of TCR repertoires than other groups. Despite higher TCR clonal expansions in HAM/TSP patients, no disease-specific TCRs were shared among patients. Only non-shared or “private” TCR repertoires was observed. While no clones that shared the same CDR3 amino acid sequences were seen in either HC or MS patients, there was a cluster of related CDR3 amino acid sequences observed for 18 out of 34 MS patients when evaluated by phylogenetic tree analysis. This suggests that a TCR-repertoire signature may be identified in a subset of patients with MS.

Vedolizumab is associated with changes in innate rather than adaptive immunity in patients with inflammatory bowel disease

OBJECTIVE

Vedolizumab, a monoclonal antibody directed against the integrin heterodimer α4β7, is approved for the treatment of Crohn's disease and ulcerative colitis. The efficacy of vedolizumab has been suggested to result from inhibition of intestinal T cell trafficking although human data to support this conclusion are scarce. We therefore performed a comprehensive analysis of vedolizumab-induced alterations in mucosal and systemic immunity in patients with inflammatory bowel disease (IBD), using anti-inflammatory therapy with the TNFα antibody infliximab as control.

DESIGN

Immunophenotyping, immunohistochemistry, T cell receptor profiling and RNA sequencing were performed using blood and colonic biopsies from patients with IBD before and during treatment with vedolizumab (n=18) or, as control, the anti-TNFα antibody infliximab (n=20). Leucocyte trafficking in vivo was assessed using single photon emission computed tomography and endomicroscopy.

RESULTS

Vedolizumab was not associated with alterations in the abundance or phenotype of lamina propria T cells and did not affect the mucosal T cell repertoire or leucocyte trafficking in vivo. Surprisingly, however, α4β7 antibody treatment was associated with substantial effects on innate immunity including changes in macrophage populations and pronounced alterations in the expression of molecules involved in microbial sensing, chemoattraction and regulation of the innate effector response. These effects were specific to vedolizumab, not observed in response to the TNFα antibody infliximab, and associated with inhibition of intestinal inflammation.

CONCLUSION

Our findings suggest that modulation of innate immunity contributes to the therapeutic efficacy of vedolizumab in IBD.

TRIAL REGISTRATION NUMBER

NCT02694588.

CapTCR-seq: hybrid capture for T-cell receptor repertoire profiling

Mature T-cell lymphomas consisting of an expanded clonal population of T cells that possess common rearrangements of the T-cell receptor (TCR) encoding genes can be identified and monitored using molecular methods of T-cell repertoire analysis. We have developed a hybrid-capture method that enriches DNA sequencing libraries for fragments encoding rearranged TCR genes from all 4 loci in a single reaction. We use this method to describe the TCR repertoires of 63 putative lymphoma clinical isolates, 7 peripheral blood mononuclear cell (PBMC) populations, and a collection of tumor infiltrating lymphocytes. Dominant Variable (V) and Joining (J) gene pair rearrangements in cancer cells were confirmed by polymerase chain reaction (PCR) amplification and Sanger sequencing; clonality assessment of clinical isolates using BIOMED-2 methods showed agreement for 73% and 77% of samples at the β and γ loci, respectively, whereas β locus V and J allele prevalence in PBMCs were well correlated with results from commercial PCR-based DNA sequencing assays (r 2 = 0.94 with Adaptive ImmunoSEQ, 0.77-0.83 with Invivoscribe LymphoTrack TRB Assay). CapTCR-seq allows for rapid, high-throughput and flexible characterization of dominant clones within TCR repertoire that will facilitate quantitative analysis of patient samples and enhance sensitivity of tumor surveillance over time.

Adipose Tissue is Enriched for Activated and Late-Differentiated CD8+ T Cells and Shows Distinct CD8+ Receptor Usage, Compared With Blood in HIV-Infected Persons

BACKGROUND

Adverse viral and medication effects on adipose tissue contribute to the development of metabolic disease in HIV-infected persons, but T cells also have a central role modulating local inflammation and adipocyte function. We sought to characterize potentially proinflammatory T-cell populations in adipose tissue among persons on long-term antiretroviral therapy and assess whether adipose tissue CD8 T cells represent an expanded, oligoclonal population.

METHODS

We recruited 10 HIV-infected, non-diabetic, overweight or obese adults on efavirenz, tenofovir, and emtricitabine for >4 years with consistent viral suppression. We collected fasting blood and subcutaneous abdominal adipose tissue to measure the percentage of CD4 and CD8 T cells expressing activation, exhaustion, late differentiation/senescence, and memory surface markers. We performed T-cell receptor (TCR) sequencing on sorted CD8 cells. We compared the proportion of each T-cell subset and the TCR repertoire diversity, in blood versus adipose tissue.

RESULTS

Adipose tissue had a higher percentage of CD3CD8 T cells compared with blood (61.0% vs. 51.7%, P < 0.01) and was enriched for both activated CD8HLA-DR T cells (5.5% vs. 0.9%, P < 0.01) and late-differentiated CD8CD57 T cells (37.4% vs. 22.7%, P < 0.01). Adipose tissue CD8 T cells displayed distinct TCRβ V and J gene usage, and the Shannon Entropy index, a measure of overall TCRβ repertoire diversity, was lower compared with blood (4.39 vs. 4.46; P = 0.05).

CONCLUSIONS

Adipose tissue is enriched for activated and late-differentiated CD8 T cells with distinct TCR usage. These cells may contribute to tissue inflammation and impaired adipocyte fitness in HIV-infected persons.

High CD8 T-Cell Receptor Clonality and Altered CDR3 Properties Are Associated With Elevated Isolevuglandins in Adipose Tissue During Diet-Induced Obesity

Adipose tissue (AT) CD4⁺ and CD8⁺ T cells contribute to obesity-associated insulin resistance. Prior studies identified conserved T-cell receptor (TCR) chain families in obese AT, but the presence and clonal expansion of specific TCR sequences in obesity has not been assessed. We characterized AT and liver CD8⁺ and CD4⁺ TCR repertoires of mice fed a low-fat diet (LFD) and high-fat diet (HFD) using deep sequencing of the TCRβ chain to quantify clonal expansion, gene usage, and CDR3 sequence. In AT CD8⁺ T cells, HFD reduced TCR diversity, increased the prevalence of public TCR clonotypes, and selected for TCR CDR3 regions enriched in positively charged and less polarized amino acids. Although TCR repertoire alone could distinguish between LFD- and HFD-fed mice, these properties of the CDR3 region of AT CD8⁺ T cells from HFD-fed mice led us to examine the role of negatively charged and nonpolar isolevuglandin (isoLG) adduct-containing antigen-presenting cells within AT. IsoLG-adducted protein species were significantly higher in AT macrophages of HFD-fed mice; isoLGs were elevated in M2-polarized macrophages, promoting CD8⁺ T-cell activation. Our findings demonstrate that clonal TCR expansion that favors positively charged CDR3s accompanies HFD-induced obesity, which may be an antigen-driven response to isoLG accumulation in macrophages.

Complement C5a Fosters Squamous Carcinogenesis and Limits T Cell Response to Chemotherapy

Complement is a critical component of humoral immunity implicated in cancer development; however, its biological contributions to tumorigenesis remain poorly understood. Using the K14-HPV16 transgenic mouse model of squamous carcinogenesis, we report that urokinase (uPA)+ macrophages regulate C3-independent release of C5a during premalignant progression, which in turn regulates protumorigenic properties of C5aR1+ mast cells and macrophages, including suppression of CD8+ T cell cytotoxicity. Therapeutic inhibition of C5aR1 via the peptide antagonist PMX-53 improved efficacy of paclitaxel chemotherapy associated with increased presence and cytotoxic properties of CXCR3+ effector memory CD8+ T cells in carcinomas, dependent on both macrophage transcriptional programming and IFNγ. Together, these data identify C5aR1-dependent signaling as an important immunomodulatory program in neoplastic tissue tractable for combinatorial cancer immunotherapy.

T-cell receptor-α repertoire of CD8+ T cells following allogeneic stem cell transplantation using next-generation sequencing

Alloreactivity or opportunistic infections following allogeneic stem cell transplantation are difficult to predict and contribute to post-transplantation mortality. How these immune reactions result in changes to the T-cell receptor repertoire remains largely unknown. Using next-generation sequencing, the T-cell receptor alpha (TRα) repertoire of naïve and memory CD8⁺ T cells from 25 patients who had received different forms of allogeneic transplantation was analyzed. In parallel, reconstitution of the CD8⁺/CD4⁺ T-cell subsets was mapped using flow cytometry. When comparing the influence of anti-T-cell therapy, a delay in the reconstitution of the naïve CD8⁺ T-cell repertoire was observed in patients who received in vivo T-cell depletion using antithymocyte globulin or post-transplantation cyclophosphamide in case of haploidentical transplantation. Sequencing of the TRα identified a repertoire consisting of more dominant clonotypes (>1% of reads) in these patients at 6 and 18 months post transplantation. When comparing donor and recipient, approximately 50% and approximately 80% of the donors’ memory repertoire were later retrieved in the naïve and memory CD8⁺ T-cell receptor repertoire of the recipients, respectively. Although there was a remarkable expansion of single clones observed in the recipients’ memory CD8⁺ TRα repertoire, no clear association between graft-versus-host disease or cytomegalovirus infection and T-cell receptor diversity was identified. A lower TRα diversity was observed in recipients of a cytomegalovirus-seropositive donor (P=0.014). These findings suggest that CD8⁺ T-cell reconstitution in transplanted patients is influenced by the use of T-cell depletion or immunosuppression and the donor repertoire.

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.

Meta-analysis of RNA sequencing datasets reveals an association between TRAJ23, psoriasis, and IL-17A

Numerous studies of relatively few patients have linked T cell receptor (TCR) genes to psoriasis but have yielded dramatically conflicting results. To resolve these discrepancies, we have chosen to mine RNA-Seq datasets for patterns of TCR gene segment usage in psoriasis. A meta-analysis of 3 existing and 1 unpublished datasets revealed a statistically significant link between the relative expression of TRAJ23 and psoriasis and the psoriasis-associated cytokine IL-17A. TRGV5, a TCR-γ segment, was also associated with psoriasis but correlated instead with IL-36A, other IL-36 family members, and IL-17C (not IL-17A). In contrast, TRAJ39 was strongly associated with healthy skin. T cell diversity measurements and analysis of CDR3 sequences were also conducted, revealing no psoriasis-associated public CDR3 sequences. Finally, in comparison with the expression of TCR-αβ genes, the expression of TCR-γδ genes was relatively low but mildly elevated in psoriatic skin. These results have implications for the development of targeted therapies for psoriasis and other autoimmune diseases. Also, the techniques employed in this study have applications in other fields, such as cancer immunology and infectious disease.

A transcriptionally and functionally distinct PD-1+ CD8+ T cell pool with predictive potential in non-small-cell lung cancer treated with PD-1 blockade

Evidence from mouse chronic viral infection models suggests that CD8⁺ T cell subsets characterized by distinct expression levels of the receptor PD-1 diverge in their state of exhaustion and potential for reinvigoration by PD-1 blockade. However, it remains unknown whether T cells in human cancer adopt a similar spectrum of exhausted states based on PD-1 expression levels. We compared transcriptional, metabolic and functional signatures of intratumoral CD8⁺ T lymphocyte populations with high (PD-1^T), intermediate (PD-1^N) and no PD-1 expression (PD-1^-) from non-small-cell lung cancer patients. PD-1^T T cells showed a markedly different transcriptional and metabolic profile from PD-1^N and PD-1^- lymphocytes, as well as an intrinsically high capacity for tumor recognition. Furthermore, while PD-1^T lymphocytes were impaired in classical effector cytokine production, they produced CXCL13, which mediates immune cell recruitment to tertiary lymphoid structures. Strikingly, the presence of PD-1^T cells was strongly predictive for both response and survival in a small cohort of non-small-cell lung cancer patients treated with PD-1 blockade. The characterization of a distinct state of tumor-reactive, PD-1-bright lymphocytes in human cancer, which only partially resembles that seen in chronic infection, provides potential avenues for therapeutic intervention.

Recasting Human Vδ1 Lymphocytes in an Adaptive Role

γδ T cells are unconventional lymphocytes commonly described as 'innate-like' in function, which can respond in both a T cell receptor (TCR)-independent and also major histocompatibility complex (MHC)-unrestricted TCR-dependent manner. While the relative importance of TCR recognition had remained unclear, recent studies revealed that human Vδ1 T cells display unexpected parallels with adaptive αβ T cells. Vδ1 T cells undergo profound and highly focussed clonal expansion from an initially diverse and private TCR repertoire, most likely in response to specific immune challenges. Concomitantly, they differentiate from a Vδ1 T cell naïve (Tnaïve) to a Vδ1 T cell effector (Teffector) phenotype, marked by the downregulation of lymphoid homing receptors and upregulation of peripheral homing receptors and effector markers. This suggests that an adaptive paradigm applies to Vδ1 T cells, likely involving TCR-dependent but MHC-unrestricted responses to microbial and non-microbial challenges.

CD154 Costimulation Shifts the Local T-Cell Receptor Repertoire Not Only During Thymic Selection but Also During Peripheral T-Dependent Humoral Immune Responses

CD154 is a transmembrane cytokine expressed transiently on activated CD4 T cells upon T-cell receptor (TCR) stimulation that interacts with CD40 on antigen-presenting cells. The signaling via CD154:CD40 is essential for B-cell maturation and germinal center formation and also for the final differentiation of CD4 T cells during T-dependent humoral immune responses. Recent data demonstrate that CD154 is critically involved in the selection of T-cell clones during the negative selection process in the thymus. Whether CD154 signaling influences the TCR repertoire during peripheral T-dependent humoral immune responses has not yet been elucidated. To find out, we used CD154-deficient mice and assessed the global TCRβ repertoire in T-cell zones (TCZ) of spleens by high-throughput sequencing after induction of a Th2 response to the multiepitopic antigen sheep red blood cells. Qualitative and quantitative comparison of the splenic TCZ-specific TCRβ repertoires revealed that CD154 deficiency shifts the distribution of Vβ-Jβ genes after antigen exposure. This data led to the conclusion that costimulation via CD154:CD40 during the interaction of T cells with CD40-matured B cells contributes to the recruitment of T-cell clones into the immune response and thereby shapes the peripheral TCR repertoire.

RepSeq Data Representativeness and Robustness Assessment by Shannon Entropy

High-throughput sequencing (HTS) has the potential to decipher the diversity of T cell repertoires and their dynamics during immune responses. Applied to T cell subsets such as T effector and T regulatory cells, it should help identify novel biomarkers of diseases. However, given the extreme diversity of TCR repertoires, understanding how the sequencing conditions, including cell numbers, biological and technical sampling and sequencing depth, impact the experimental outcome is critical to proper use of these data. Here, we assessed the representativeness and robustness of TCR repertoire diversity assessment according to experimental conditions. By comparative analyses of experimental datasets and computer simulations, we found that (i) for small samples, the number of clonotypes recovered is often higher than the number of cells per sample, even after removing the singletons; (ii) high-sequencing depth for small samples alters the clonotype distributions, which can be corrected by filtering the datasets using Shannon entropy as a threshold; and (iii) a single sequencing run at high depth does not ensure a good coverage of the clonotype richness in highly polyclonal populations, which can be better covered using multiple sequencing. Altogether, our results warrant better understanding and awareness of the limitation of TCR diversity analyses by HTS and justify the development of novel computational tools for improved modeling of the highly complex nature of TCR repertoires.

IRProfiler - a software toolbox for high throughput immune receptor profiling

Background

The study of the huge diversity of immune receptors, often referred to as immune repertoire profiling, is a prerequisite for diagnosis, prognostication and monitoring of hematological disorders. In the era of high-throughput sequencing (HTS), the abundance of immunogenetic data has revealed unprecedented opportunities for the thorough profiling of T-cell receptors (TR) and B-cell receptors (BcR). However, the volume of the data to be analyzed mandates for efficient and ease-to-use immune repertoire profiling software applications.

Results

This work introduces Immune Repertoire Profiler (IRProfiler), a novel software pipeline that delivers a number of core receptor repertoire quantification and comparison functionalities on high-throughput TR and BcR sequencing data. Adopting 5 alternative clonotype definitions, IRProfiler implements a series of algorithms for 1) data filtering, 2) calculation of clonotype diversity and expression, 3) calculation of gene usage for the V and J subgroups, 4) detection of shared and exclusive clonotypes among multiple repertoires, and 5) comparison of gene usage for V and J subgroups among multiple repertoires. IRProfiler has been implemented as a toolbox of the Galaxy bioinformatics platform, comprising 6 tools. Theoretical and experimental evaluation has shown that the tools of IRProfiler are able to scale well with respect to the size of input dataset(s). IRProfiler has been utilized by a number of recently published studies concerning hematological disorders.

Conclusion

IRProfiler is made freely available via 3 distribution channels, including the Galaxy Tool Shed. Despite being a new entry in a crowded ecosystem of immune repertoire profiling software, IRProfiler founds its added value on its support for alternative clonotype definitions in conjunction with a combination of properties stemming from its user-centric design, namely ease-of-use, ease-of-access, exploitability of the output data, and analysis flexibility.

Human γδ T Cell Receptor Repertoires in Peripheral Blood Remain Stable Despite Clearance of Persistent Hepatitis C Virus Infection by Direct-Acting Antiviral Drug Therapy

Human γδ T cells can contribute to clearance of hepatitis C virus (HCV) infection but also mediate liver inflammation. This study aimed to understand the clonal distribution of γδ T cells in peripheral blood of chronic HCV patients and following HCV clearance by interferon-free direct-acting antiviral drug therapies. To this end, γδ T cell receptor (TCR) repertoires were monitored by mRNA-based next-generation sequencing. While the percentage of Vγ9⁺ T cells was higher in patients with elevated liver enzymes and a few expanded Vδ3 clones could be identified in peripheral blood of 23 HCV-infected non-cirrhotic patients, overall clonality and complexity of γδ TCR repertoires were largely comparable to those of matched healthy donors. Monitoring eight chronic HCV patients before, during and up to 1 year after therapy revealed that direct-acting antiviral (DAA) drug therapies induced only minor alterations of TRG and TRD repertoires of Vγ9⁺ and Vγ9^- cells. Together, we show that peripheral γδ TCR repertoires display a high stability (1) by chronic HCV infection in the absence of liver cirrhosis and (2) by HCV clearance in the course of DAA drug therapy.

Analyzing Immunoglobulin Repertoires

Somatic assembly of T cell receptor and B cell receptor (BCR) genes produces a vast diversity of lymphocyte antigen recognition capacity. The advent of efficient high-throughput sequencing of lymphocyte antigen receptor genes has recently generated unprecedented opportunities for exploration of adaptive immune responses. With these opportunities have come significant challenges in understanding the analysis techniques that most accurately reflect underlying biological phenomena. In this regard, sample preparation and sequence analysis techniques, which have largely been borrowed and adapted from other fields, continue to evolve. Here, we review current methods and challenges of library preparation, sequencing and statistical analysis of lymphocyte receptor repertoire studies. We discuss the general steps in the process of immune repertoire generation including sample preparation, platforms available for sequencing, processing of sequencing data, measurable features of the immune repertoire, and the statistical tools that can be used for analysis and interpretation of the data. Because BCR analysis harbors additional complexities, such as immunoglobulin (Ig) (i.e., antibody) gene somatic hypermutation and class switch recombination, the emphasis of this review is on Ig/BCR sequence analysis.

The transcription factors Runx3 and ThPOK cross-regulate acquisition of cytotoxic function by human Th1 lymphocytes

Cytotoxic CD4 (CD4_CTX) T cells are emerging as an important component of antiviral and antitumor immunity, but the molecular basis of their development remains poorly understood. In the context of human cytomegalovirus infection, a significant proportion of CD4 T cells displays cytotoxic functions. We observed that the transcriptional program of these cells was enriched in CD8 T cell lineage genes despite the absence of ThPOK downregulation. We further show that establishment of CD4_CTX-specific transcriptional and epigenetic programs occurred in a stepwise fashion along the Th1-differentiation pathway. In vitro, prolonged activation of naive CD4 T cells in presence of Th1 polarizing cytokines led to the acquisition of perforin-dependent cytotoxic activity. This process was dependent on the Th1 transcription factor Runx3 and was limited by the sustained expression of ThPOK. This work elucidates the molecular program of human CD4_CTX T cells and identifies potential targets for immunotherapy against viral infections and cancer.