Distinctive properties of identical twins' TCR repertoires revealed by high-throughput sequencing

Age-Related Dynamics and Spectral Characteristics of the TCRβ Repertoire in Healthy Children: Implications for Immune Aging

T-cell receptor (TCR) diversity is crucial for adaptive immunity, yet baseline characterizations in pediatric populations remain sparse. We sequenced the TCRβ chain of 325 healthy Chinese children aged 0-18, categorized into six age groups. We also analyzed cellular composition and TCRβ associations using flow cytometry in 81 of these samples. Our results indicate a decrease in TCRβ diversity with age, characterized by an increase in high-frequency clonotypes and notable changes in CDR3 length and V(D)J gene usage. These changes are influenced by early life vaccinations and antigen exposures. Additionally, we found a significant association between reduced TCRβ diversity and a decrease in CD4+ T naïve cells. We also developed a predictive model that identifies specific TCRβ features as potential biomarkers for biological age, validated by their significant correlation with changes in the immune repertoire. These findings enhance our understanding of age-related variations in the TCRβ repertoire among children, providing resourceful information for research on pediatric TCR in health and disease.

HTS and scRNA-seq revealed that the location and RSS quality of the mammalian TRBV and TRBJ genes impact biased rearrangement

The quality of Recombination signal sequences (RSSs), location, and genetics of mammalian V, D, and J genes synergistically affect the recombination frequency of genes; however, the specific regulatory mechanism and efficiency have not been elucidated. By taking advantage of single-cell RNA-sequencing (scRNA-seq) and high-throughput sequencing (HTS) to investigate V(D)J rearrangement characteristics in the CDR3 repertoire, we found that the distal and proximal V genes (or J genes) "to D" gene were involved in rearrangement significantly more frequently than the middle V genes (or J genes) in the TRB locus among various species, including Primates (human and rhesus monkey), Rodentia (BALB/c, C57BL/6, and Kunming mice), Artiodactyla (buffalo), and Chiroptera (Rhinolophus affinis). The RSS quality of the V and J genes affected their frequency in rearrangement to varying degrees, especially when the V-RSSs with recombination signal information content (RIC) score < -45 significantly reduced the recombination frequency of the V gene. The V and J genes that were "away from D" had the dual advantages of recombinant structural accessibility and relatively high-quality RSSs, which promoted their preferential utilization in rearrangement. The quality of J-RSSs formed during mammalian evolution was apparently greater than that of V-RSSs, and the D-J distance was obviously shorter than that of V-D, which may be one of the reasons for guaranteeing that the "D-to-J preceding V-to-DJ rule" occurred when rearranged. This study provides a novel perspective on the mechanism and efficiency of V-D-J rearrangement in the mammalian TRB locus, as well as the biased utilization characteristics and application of V and J genes in the initial CDR3 repertoire.

Human TCR repertoire in cancer

BACKGROUND

T cells, the "superstar" of the immune system, play a crucial role in antitumor immunity. T-cell receptors (TCR) are crucial molecules that enable T cells to identify antigens and start immunological responses. The body has evolved a unique method for rearrangement, resulting in a vast diversity of TCR repertoires. A healthy TCR repertoire is essential for the particular identification of antigens by T cells.

METHODS

In this article, we systematically summarized the TCR creation mechanisms and analysis methodologies, particularly focusing on the application of next-generation sequencing (NGS) technology. We explore the TCR repertoire in health and cancer, and discuss the implications of TCR repertoire analysis in understanding carcinogenesis, cancer progression, and treatment.

RESULTS

The TCR repertoire analysis has enormous potential for monitoring the emergence and progression of malignancies, as well as assessing therapy response and prognosis. The application of NGS has dramatically accelerated our comprehension of TCR diversity and its role in cancer immunity.

CONCLUSIONS

To substantiate the significance of TCR repertoires as biomarkers, more thorough and exhaustive research should be conducted. The TCR repertoire analysis, enabled by advanced sequencing technologies, is poised to become a crucial tool in the future of cancer diagnosis, monitoring, and therapy evaluation.

Establishment of a humanized mouse model using steady-state peripheral blood-derived hematopoietic stem and progenitor cells facilitates screening of cancer-targeted T-cell repertoires

Background

Cancer-targeted T-cell receptor T (TCR-T) cells hold promise in treating cancers such as hematological malignancies and breast cancers. However, approaches to obtain cancer-reactive TCR-T cells have been unsuccessful.

Methods

Here, we developed a novel strategy to screen for cancer-targeted TCR-T cells using a special humanized mouse model with person-specific immune fingerprints. Rare steady-state circulating hematopoietic stem and progenitor cells were expanded via three-dimensional culture of steady-state peripheral blood mononuclear cells, and then the expanded cells were applied to establish humanized mice. The human immune system was evaluated according to the kinetics of dendritic cells, monocytes, T-cell subsets, and cytokines. To fully stimulate the immune response and to obtain B-cell precursor NAML-6- and triple-negative breast cancer MDA-MB-231-targeted TCR-T cells, we used the inactivated cells above to treat humanized mice twice a day every 7 days. Then, human T cells were processed for TCR β-chain (TRB) sequencing analysis. After the repertoires had been constructed, features such as the fraction, diversity, and immune signature were investigated.

Results

The results demonstrated an increase in diversity and clonality of T cells after treatment. The preferential usage and features of TRBV, TRBJ, and the V-J combination were also changed. The stress also induced highly clonal expansion. Tumor burden and survival analysis demonstrated that stress induction could significantly inhibit the growth of subsequently transfused live tumor cells and prolong the survival of the humanized mice.

Conclusions

We constructed a personalized humanized mouse model to screen cancer-targeted TCR-T pools. Our platform provides an effective source of cancer-targeted TCR-T cells and allows for the design of patient-specific engineered T cells. It therefore has the potential to greatly benefit cancer treatment.

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.

Distinguishing between monozygotic twins' blood samples through immune repertoire sequencing

Monozygotic (MZ) twins with highly similar genomic DNA sequences can not be distinguished by conventional forensic DNA testing. The immune repertoire (IR) reflects an individual's immune history, which is unique between individuals, has been applied to individualized treatment in precision medicine. However, the application of IR in forensic genetics has not been reported to date. In this study, the diversity in the complementary determining region 3 (CDR3) of both the T-cell receptor β chain (TCRβ) and B-cell receptor heavy chain (also known as immunoglobulin heavy chain, IGH) in four pairs of MZ twins were analyzed. The results showed that the amino acid sequences length distribution frequency of TCRβ CDR3 had 4-10 differences, and the nucleic acid sequences length distribution frequency of TCRβ CDR3 had 2-7 differences between MZ twins. The shared difference of four pairs of MZ twins focused on the length distribution frequency of 34 bp nucleotide sequences in TCRβ. By analyzing the usage frequency of V and J genes in TCRβ and IGH CDR3 DNA sequence rearrangements, we also found that there were biases between each pair of MZ twins, and the usage frequency of TRBJ2-3 showed common differences between each pair of MZ twins. Furthermore, each pair of MZ twins had its own unique V-J genes combination mode in TCRβ and IGH CDR3 DNA sequences. This study, for the first time, suggested that IR can be used as a potential biological marker to distinguish MZ twins.

Targeted long-read sequencing facilitates phased diploid assembly and genotyping of the human T cell receptor alpha, delta, and beta loci

T cell receptors (TCRs) recognize peptide fragments presented by the major histocompatibility complex (MHC) and are critical to T cell-mediated immunity. Recent data have indicated that genetic diversity within TCR-encoding gene regions is underexplored, limiting understanding of the impact of TCR loci polymorphisms on TCR function in disease, even though TCR repertoire signatures (1) are heritable and (2) associate with disease phenotypes. To address this, we developed a targeted long-read sequencing approach to generate highly accurate haplotype resolved assemblies of the TCR beta (TRB) and alpha/delta (TRA/D) loci, facilitating the genotyping of all variant types, including structural variants. We validate our approach using two mother-father-child trios and 5 unrelated donors representing multiple populations. This resulted in improved genotyping accuracy and the discovery of 84 undocumented V, D, J, and C alleles, demonstrating the utility of this framework for improving our understanding of TCR diversity and function in disease.

Evolution and modulation of antigen-specific T cell responses in melanoma patients

Analyzing antigen-specific T cell responses at scale has been challenging. Here, we analyze three types of T cell receptor (TCR) repertoire data (antigen-specific TCRs, TCR-repertoire, and single-cell RNA + TCRαβ-sequencing data) from 515 patients with primary or metastatic melanoma and compare it to 783 healthy controls. Although melanoma-associated antigen (MAA) -specific TCRs are restricted to individuals, they share sequence similarities that allow us to build classifiers for predicting anti-MAA T cells. The frequency of anti-MAA T cells distinguishes melanoma patients from healthy and predicts metastatic recurrence from primary melanoma. Anti-MAA T cells have stem-like properties and frequent interactions with regulatory T cells and tumor cells via Galectin9-TIM3 and PVR-TIGIT -axes, respectively. In the responding patients, the number of expanded anti-MAA clones are higher after the anti-PD1(+anti-CTLA4) therapy and the exhaustion phenotype is rescued. Our systems immunology approach paves the way for understanding antigen-specific responses in human disorders.

Previous studies have characterized the diversity and dynamics of the T cell receptor (TCR) repertoire in patients with solid cancer. Here, by analyzing TCR repertoire data from multiple datasets, the authors report that melanoma-associated antigen-specific TCRs can be used to separate metastatic melanoma patients from healthy controls and to follow anti-tumor responses in patients treated with immunotherapy.

Common Trajectories of Highly Effective CD19-Specific CAR T Cells Identified by Endogenous T-cell Receptor Lineages

Current chimeric antigen receptor-modified (CAR) T-cell products are evaluated in bulk, without assessing functional heterogeneity. We therefore generated a comprehensive single-cell gene expression and T-cell receptor (TCR) sequencing data set using pre- and postinfusion CD19-CAR T cells from blood and bone marrow samples of pediatric patients with B-cell acute lymphoblastic leukemia. We identified cytotoxic postinfusion cells with identical TCRs to a subset of preinfusion CAR T cells. These effector precursor cells exhibited a unique transcriptional profile compared with other preinfusion cells, corresponding to an unexpected surface phenotype (TIGIT+, CD62Llo, CD27-). Upon stimulation, these cells showed functional superiority and decreased expression of the exhaustion-associated transcription factor TOX. Collectively, these results demonstrate diverse effector potentials within preinfusion CAR T-cell products, which can be exploited for therapeutic applications. Furthermore, we provide an integrative experimental and analytic framework for elucidating the mechanisms underlying effector development in CAR T-cell products.

SIGNIFICANCE

Utilizing clonal trajectories to define transcriptional potential, we find a unique signature of CAR T-cell effector precursors present in preinfusion cell products. Functional assessment of cells with this signature indicated early effector potential and resistance to exhaustion, consistent with postinfusion cellular patterns observed in patients. This article is highlighted in the In This Issue feature, p. 2007.

T-Cell Receptor Repertoire Sequencing and Its Applications: Focus on Infectious Diseases and Cancer

The immune system is a dynamic feature of each individual and a footprint of our unique internal and external exposures. Indeed, the type and level of exposure to physical and biological agents shape the development and behavior of this complex and diffuse system. Many pathological conditions depend on how our immune system responds or does not respond to a pathogen or a disease or on how the regulation of immunity is altered by the disease itself. T-cells are important players in adaptive immunity and, together with B-cells, define specificity and monitor the internal and external signals that our organism perceives through its specific receptors, TCRs and BCRs, respectively. Today, high-throughput sequencing (HTS) applied to the TCR repertoire has opened a window of opportunity to disclose T-cell repertoire development and behavior down to the clonal level. Although TCR repertoire sequencing is easily accessible today, it is important to deeply understand the available technologies for choosing the best fit for the specific experimental needs and questions. Here, we provide an updated overview of TCR repertoire sequencing strategies, providers and applications to infectious diseases and cancer to guide researchers’ choice through the multitude of available options. The possibility of extending the TCR repertoire to HLA characterization will be of pivotal importance in the near future to understand how specific HLA genes shape T-cell responses in different pathological contexts and will add a level of comprehension that was unthinkable just a few years ago.

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.

SARS-CoV-2 (COVID-19)-specific T cell and B cell responses in convalescent rheumatoid arthritis: Monozygotic twins pair case observation

Rheumatoid arthritis (RA) patients present higher risk of SARS-CoV-2 infection (COVID-19), and proper management of the disease in this population requires a better understanding of how the immune system controls the virus. We analyzed the T cell and B cell phenotypes, and their repertoire in a pair of monozygotic twins with RA mismatched for COVID-19 infection. Twin- was not infected, while Twin+ was infected and effectively controlled the infection. We found no significant changes on the αβ T cell composition, while γδ T cells and B cells presented considerable expansion of memory population in Twin+ and robust T/B cell responses to several SARS-CoV-2 peptides. T cell receptor β/γ-chain and immunoglobulin heavy chain next-generation sequencing depicted a remarkable higher diversity in Twin+ compared with Twin-, despite no significant changes being found in variable/joining family usage. Repertoire overlap analyses showed that, although being identical twins, very few clones were shared between them, indicating that COVID-19 may lead to deep changes on the immune cell repertoire in RA patients. Altogether, our results indicate that RA patients may develop robust and persistent COVID-19-specific T/B cell responses; γδ T cells and B cells may play a key role in the management of COVID-19 in RA, and the infection may lead to a profound reshaping of immune cell receptor specificities.

Combining genotypes and T cell receptor distributions to infer genetic loci determining V(D)J recombination probabilities

Every T cell receptor (TCR) repertoire is shaped by a complex probabilistic tangle of genetically determined biases and immune exposures. T cells combine a random V(D)J recombination process with a selection process to generate highly diverse and functional TCRs. The extent to which an individual’s genetic background is associated with their resulting TCR repertoire diversity has yet to be fully explored. Using a previously published repertoire sequencing dataset paired with high-resolution genome-wide genotyping from a large human cohort, we infer specific genetic loci associated with V(D)J recombination probabilities using genome-wide association inference. We show that V(D)J gene usage profiles are associated with variation in the TCRB locus and, specifically for the functional TCR repertoire, variation in the major histocompatibility complex locus. Further, we identify specific variations in the genes encoding the Artemis protein and the TdT protein to be associated with biasing junctional nucleotide deletion and N-insertion, respectively. These results refine our understanding of genetically-determined TCR repertoire biases by confirming and extending previous studies on the genetic determinants of V(D)J gene usage and providing the first examples of trans genetic variants which are associated with modifying junctional diversity. Together, these insights lay the groundwork for further explorations into how immune responses vary between individuals.

T Cell Repertoire Abnormality in Immunodeficiency Patients with DNA Repair and Methylation Defects

Both DNA damage response and methylation play a crucial role in antigen receptor recombination by creating a diverse repertoire in developing lymphocytes, but how their defects relate to T cell repertoire and phenotypic heterogeneity of immunodeficiency remains obscure. We studied the TCR repertoire in patients with the mutation in different genes (ATM, DNMT3B, ZBTB24, RAG1, DCLRE1C, and JAK3) and uncovered distinct characteristics of repertoire diversity. We propose that early aberrancies in thymus T cell development predispose to the heterogeneous phenotypes of the immunodeficiency spectrum. Shorter CDR3 lengths in ATM-deficient patients, resulting from a decreased number of nucleotide insertions during VDJ recombination in the pre-selected TCR repertoire, as well as the increment of CDR3 tyrosine residues, lead to the enrichment of pathology-associated TCRs, which may contribute to the phenotypes of ATM deficiency. Furthermore, patients with DNMT3B and ZBTB24 mutations who exhibit discrepant phenotypes present longer CDR3 lengths and reduced number of known pathology-associated TCRs.

Analysis of T-Cell Receptor Repertoire in Transplantation: Fingerprint of T Cell-mediated Alloresponse

T cells play a key role in determining allograft function by mediating allogeneic immune responses to cause rejection, and recent work pointed their role in mediating tolerance in transplantation. The unique T-cell receptor (TCR) expressed on the surface of each T cell determines the antigen specificity of the cell and can be the specific fingerprint for identifying and monitoring. Next-generation sequencing (NGS) techniques provide powerful tools for deep and high-throughput TCR profiling, and facilitate to depict the entire T cell repertoire profile and trace antigen-specific T cells in circulation and local tissues. Tailing T cell transcriptomes and TCR sequences at the single cell level provides a full landscape of alloreactive T-cell clones development and biofunction in alloresponse. Here, we review the recent advances in TCR sequencing techniques and computational tools, as well as the recent discovery in overall TCR profile and antigen-specific T cells tracking in transplantation. We further discuss the challenges and potential of using TCR sequencing-based assays to profile alloreactive TCR repertoire as the fingerprint for immune monitoring and prediction of rejection and tolerance.

High Throughput Human T Cell Receptor Sequencing: A New Window Into Repertoire Establishment and Alloreactivity

Recent advances in high throughput sequencing (HTS) of T cell receptors (TCRs) and in transcriptomic analysis, particularly at the single cell level, have opened the door to a new level of understanding of human immunology and immune-related diseases. In this article, we discuss the use of HTS of TCRs to discern the factors controlling human T cell repertoire development and how this approach can be used in combination with human immune system (HIS) mouse models to understand human repertoire selection in an unprecedented manner. An exceptionally high proportion of human T cells has alloreactive potential, which can best be understood as a consequence of the processes governing thymic selection. High throughput TCR sequencing has allowed assessment of the development, magnitude and nature of the human alloresponse at a new level and has provided a tool for tracking the fate of pre-transplant-defined donor- and host-reactive TCRs following transplantation. New insights into human allograft rejection and tolerance obtained with this method in combination with single cell transcriptional analyses are reviewed here.

Generation of self-reactive, shared T-cell receptor α chains in the human thymus

The T-cell receptor (TCR) repertoire is generated in a semistochastic process of gene recombination and pairing of TCRα to TCRβ chains with the estimated total TCR diversity of >10⁸. Despite this high diversity, similar or identical TCR chains are found to recur in immune responses. Here, we analyzed the thymic generation of TCR sequences previously associated with recognition of self- and nonself-antigens, represented by sequences associated with autoimmune diabetes and HIV, respectively. Unexpectedly, in the CD4⁺ compartment TCRα chains associated with the recognition of self-antigens were generated in significantly higher numbers than TCRα chains associated with the recognition of nonself-antigens. The analysis of the circulating repertoire further showed that these chains are not lost in negative selection nor predominantly converted to the regulatory T-cell lineage. The high abundance of self-reactive TCRα chains in multiple individuals suggests that the human thymus has a predilection to generate self-reactive TCRα chains independently of the HLA-type and that the individual risk of autoimmunity may be modulated by the TCRβ repertoire associated with these chains.

VJ Segment Usage of TCR-Beta Repertoire in Monozygotic Cystic Fibrosis Twins

Sixteen monozygotic cystic fibrosis (CF) twin pairs of whom 14 pairs were homozygous for the most common p.Phe508del CFTR mutation were selected from the European Cystic Fibrosis Twin and Sibling Study Cohort. The monozygotic twins were examined in their T cell receptor (TCR) repertoire in peripheral blood by amplicon sequencing of the CDR3 variable region of the ß-chain. The recruitment of TCR J and V genes for recombination and selection in the thymus showed a strong genetic influence in the CF twin cohort as indicated by the shortest Jensen-Shannon distance to the twin individual. Exceptions were the clinically most discordant and/or most severely affected twin pairs where clonal expansion probably caused by recurrent pulmonary infections overshadowed the impact of the identical genomic blueprint. In general the Simpson clonality was low indicating that the population of TCRß clonotypes of the CF twins was dominated by the naïve T-cell repertoire. Intrapair sharing of clonotypes was significantly more frequent among monozygotic CF twins than among pairs of unrelated CF patients. Complete nucleotide sequence identity was observed in about 0.11% of CDR3 sequences which partially should represent persisting fetal clones derived from the same progenitor T cells. Complete amino acid sequence identity was noted in 0.59% of clonotypes. Of the nearly 40,000 frequent amino acid clonotypes shared by at least two twin siblings 99.8% were already known within the immuneACCESS database and only 73 had yet not been detected indicating that the CDR3ß repertoire of CF children and adolescents does not carry a disease-specific signature but rather shares public clones with that of the non-CF community. Clonotypes shared within twin pairs and between unrelated CF siblings were highly abundant among healthy non-CF people, less represented in individuals with infectious disease and uncommon in patients with cancer. This subset of shared CF clonotypes defines CDR3 amino acid sequences that are more common in health than in disease.

Distinct organization of adaptive immunity in the long-lived rodent Spalax galili

A balanced immune response is a cornerstone of healthy aging. Here, we uncover distinctive features of the long-lived blind mole-rat (Spalax spp.) adaptive immune system, relative to humans and mice. The T-cell repertoire remains diverse throughout the Spalax lifespan, suggesting a paucity of large long-lived clones of effector-memory T cells. Expression of master transcription factors of T-cell differentiation, as well as checkpoint and cytotoxicity genes, remains low as Spalax ages. The thymus shrinks as in mice and humans, while interleukin-7 and interleukin-7 receptor expression remains high, potentially reflecting the sustained homeostasis of naive T cells. With aging, immunoglobulin hypermutation level does not increase and the immunoglobulin-M repertoire remains diverse, suggesting shorter B-cell memory and sustained homeostasis of innate-like B cells. The Spalax adaptive immune system thus appears biased towards sustained functional and receptor diversity over specialized, long-lived effector-memory clones-a unique organizational strategy that potentially underlies this animal's extraordinary longevity and healthy aging.

CDR3 and V genes show distinct reconstitution patterns in T cell repertoire post-allogeneic bone marrow transplantation

Restoration of T cell repertoire diversity after allogeneic bone marrow transplantation (allo-BMT) is crucial for immune recovery. T cell diversity is produced by rearrangements of germline gene segments (V (D) and J) of the T cell receptor (TCR) α and β chains, and selection induced by binding of TCRs to MHC-peptide complexes. Multiple measures were proposed for this diversity. We here focus on the V-gene usage and the CDR3 sequences of the beta chain. We compared multiple T cell repertoires to follow T cell repertoire changes post-allo-BMT in HLA-matched related donor and recipient pairs. Our analyses of the differences between donor and recipient complementarity determining region 3 (CDR3) beta composition and V-gene profile show that the CDR3 sequence composition does not change during restoration, implying its dependence on the HLA typing. In contrast, V-gene usage followed a time-dependent pattern, initially following the donor profile and then shifting back to the recipients' profile. The final long-term repertoire was more similar to that of the recipient's original one than the donor's; some recipients converged within months, while others took multiple years. Based on the results of our analyses, we propose that donor-recipient V-gene distribution differences may serve as clinical biomarkers for monitoring immune recovery.

Characterization of human T cell receptor repertoire data in eight thymus samples and four related blood samples

T cell receptor (TCR) is a heterodimer consisting of TCRα and TCRβ chains that are generated by somatic recombination of multiple gene segments. Nascent TCR repertoire undergoes thymic selections where non-functional and potentially autoreactive receptors are removed. During the last years, the development of high-throughput sequencing technology has allowed a large scale assessment of TCR repertoire and multiple analysis tools are now also available.

In our recent manuscript, Human thymic T cell repertoire is imprinted with strong convergence to shared sequences[1], we show highly overlapping thymic TCR repertoires in unrelated individuals. In the current Data in Brief article, we provide a more detailed characterization of the basic features of these thymic and related peripheral blood TCR repertoires. The thymus samples were collected from eight infants undergoing corrective cardiac surgery, two of whom were monozygous twins [2]. In parallel with the surgery, a small aliquot of peripheral blood was drawn from four of the donors. Genomic DNA was extracted from mechanically released thymocytes and circulating leukocytes. The sequencing of TCRα and TCRβ repertoires was performed at ImmunoSEQ platform (Adaptive Biotechnologies). The obtained repertoire data were analysed applying relevant features from immunoSEQ® 3.0 Analyzer (Adaptive Biotechnologies) and a freely available VDJTools software package for programming language R [3].

The current data analysis displays the basic features of the sequenced repertoires including observed TCR diversity, various descriptive TCR diversity measures, and V and J gene usage. In addition, multiple methods to calculate repertoire overlap between two individuals are applied. The raw sequence data provide a large database of reference TCRs in healthy individuals at an early developmental stage. The data can be exploited to improve existing computational models on TCR repertoire behaviour as well as in the generation of new models.

CD8+ T-Cell Repertoire in Human Leukocyte Antigen Class I-Mismatched Alloreactive Immune Response

In transplantation, direct allorecognition is a complex interplay between T-cell receptors (TCR) and HLA molecules and their bound peptides expressed on antigen-presenting cells. In analogy to HLA mismatched hematopoietic stem cell transplantation (HSCT), the TCR CDR3β repertoires of alloreactive cytotoxic CD8⁺ responder T cells, defined by the cell surface expression of CD137 and triggered in vitro by HLA mismatched stimulating cells, were analyzed in different HLA class I mismatched combinations. The same HLA mismatched stimulatory cells induced very different repertoires in distinct but HLA identical responders. Likewise, stimulator cells derived from HLA identical donors activated CD8⁺ cells expressing very different repertoires in the same mismatched responder. To mimic in vivo inflammation, expression of HLA class l antigens was upregulated in vitro on stimulating cells by the inflammatory cytokines TNFα and IFNβ. The repertoires differed whether the same responder cells were stimulated with cells treated or not with both cytokines. In conclusion, the selection and expansion of alloreactive cytotoxic T-cell clonotypes expressing a very diverse repertoire is observed repeatedly despite controlling for HLA disparities and is significantly influenced by the inflammatory status. This makes prediction of alloreactive T-cell repertoires a major challenge in HLA mismatched HSCT.

Immune fingerprinting through repertoire similarity

Immune repertoires provide a unique fingerprint reflecting the immune history of individuals, with potential applications in precision medicine. However, the question of how personal that information is and how it can be used to identify individuals has not been explored. Here, we show that individuals can be uniquely identified from repertoires of just a few thousands lymphocytes. We present "Immprint," a classifier using an information-theoretic measure of repertoire similarity to distinguish pairs of repertoire samples coming from the same versus different individuals. Using published T-cell receptor repertoires and statistical modeling, we tested its ability to identify individuals with great accuracy, including identical twins, by computing false positive and false negative rates < 10-6 from samples composed of 10,000 T-cells. We verified through longitudinal datasets that the method is robust to acute infections and that the immune fingerprint is stable for at least three years. These results emphasize the private and personal nature of repertoire data.

Predicting Cross-Reactivity and Antigen Specificity of T Cell Receptors

Adaptive immune recognition is mediated by specific interactions between heterodimeric T cell receptors (TCRs) and their cognate peptide-MHC (pMHC) ligands, and the methods to accurately predict TCR:pMHC interaction would have profound clinical, therapeutic and pharmaceutical applications. Herein, we review recent developments in predicting cross-reactivity and antigen specificity of TCR recognition. We discuss current experimental and computational approaches to investigate cross-reactivity and antigen-specificity of TCRs and highlight how integrating kinetic, biophysical and structural features may offer valuable insights in modeling immunogenicity. We further underscore the close inter-relationship of these two interconnected notions and the need to investigate each in the light of the other for a better understanding of T cell responsiveness for the effective clinical applications.

Preselection TCR repertoire predicts CD4+ and CD8+ T-cell differentiation state

T cells must display diversity regarding both the cell state and T-cell receptor (TCR) repertoire to provide effective immunity against pathogens; however, the generation and evolution of cellular T-cell heterogeneity in the adaptive immune system remains unclear. In the present study, a combination of multiplex PCR and immune repertoire sequencing (IR-seq) was used for a standardized analysis of the TCR β-chain repertoire of CD4⁺ naive, CD4⁺ memory, CD8⁺ naive and CD8⁺ memory T cells. We showed that the T-cell subsets could be distinguished from each another with regard to the TCR β-chain (TCR-β) diversity, CDR3 length distribution and TRBV usage, which could be observed both in the preselection and in the post-selection repertoire. Moreover, the Dβ-Jβ and Vβ-Dβ combination patterns at the initial recombination step, template-independent insertion of nucleotides and inter-subset overlap were consistent between the pre- and post-selection repertoires, with a remarkably positive correlation. Taken together, these results support differentiation of the CD4⁺ and CD8⁺ T-cell subsets prior to thymic selection, and these differences survived both positive and negative selection. In conclusion, these findings provide deeper insight into the generation and evolution of TCR repertoire generation.

Human thymic T cell repertoire is imprinted with strong convergence to shared sequences

A highly diverse repertoire of T cell antigen receptors (TCR) is created in the thymus by recombination of gene segments and the insertion or deletion of nucleotides at the junctions. Using next-generation TCR sequencing we define here the features of recombination and selection in the human TCRα and TCRβ locus, and show that a strikingly high proportion of the repertoire is shared by unrelated individuals. The thymic TCRα nucleotide repertoire was more diverse than TCRβ, with 4.1 × 10⁶ vs. 0.81 × 10⁶ unique clonotypes, and contained nonproductive clonotypes at a higher frequency (69.2% vs. 21.2%). The convergence of distinct nucleotide clonotypes to the same amino acid sequences was higher in TCRα than in TCRβ repertoire (1.45 vs. 1.06 nucleotide sequences per amino acid sequence in thymus). The gene segment usage was biased, and generally all individuals favored the same genes in both TCRα and TCRβ loci. Despite the high diversity, a large fraction of the repertoire was found in more than one donor. The shared fraction was bigger in TCRα than TCRβ repertoire, and more common in in-frame sequences than in nonproductive sequences. Thus, both biases in rearrangement and thymic selection are likely to contribute to the generation of shared repertoire in humans.

The human T-cell receptor repertoire in health and disease and potential for omics integration

The adaptive immune system arose 600 million years ago in a cold-blooded fish. Over countless generations, our antecedents tuned the function of the T-cell receptor (TCR). The TCR system is arguably the most complex known to science. The TCR evolved hypervariability to fight the hypervariability of pathogens and cancers that look to consume our resources. This review describes the genetics and architecture of the human TCR and highlights surprising new discoveries over the past years that have disproved very old dogmas. The standardization of TCR sequencing data is discussed in preparation for big data bioinformatics and predictive analysis. We next catalogue new signatures and phenomenon discovered by TCR next generation sequencing (NGS) in health and disease and work that remain to be done in this space. Finally, we discuss how TCR NGS can add to immunodiagnostics and integrate with other omics platforms for both a deeper understanding of TCR biology and its use in the clinical setting.

T cell receptor sequence clustering and antigen specificity

There has been increasing interest in the role of T cells and their involvement in cancer, autoimmune and infectious diseases. However, the nature of T cell receptor (TCR) epitope recognition at a repertoire level is not yet fully understood. Due to technological advances a plethora of TCR sequences from a variety of disease and treatment settings has become readily available. Current efforts in TCR specificity analysis focus on identifying characteristics in immune repertoires which can explain or predict disease outcome or progression, or can be used to monitor the efficacy of disease therapy. In this context, clustering of TCRs by sequence to reflect biological similarity, and especially to reflect antigen specificity have become of paramount importance. We review the main TCR sequence clustering methods and the different similarity measures they use, and discuss their performance and possible improvement. We aim to provide guidance for non-specialists who wish to use TCR repertoire sequencing for disease tracking, patient stratification or therapy prediction, and to provide a starting point for those aiming to develop novel techniques for TCR annotation through clustering.

Identification of Disease-associated Traits and Clonotypes in the T Cell Receptor Repertoire of Monozygotic Twins Affected by Inflammatory Bowel Diseases

BACKGROUND AND AIMS

Intestinal inflammation in inflammatory bowel diseases [IBD] is thought to be T cell mediated and therefore dependent on the interaction between the T cell receptor [TCR] and human leukocyte antigen [HLA] proteins expressed on antigen presenting cells. The collection of all TCRs in one individual, known as the TCR repertoire, is characterised by enormous diversity and inter-individual variability. It was shown that healthy monozygotic [MZ] twins are more similar in their TCR repertoire than unrelated individuals. Therefore MZ twins, concordant or discordant for IBD, may be useful to identify disease-related and non-genetic factors in the TCR repertoire which could potentially be used as disease biomarkers.

METHODS

Employing unique molecular barcoding that can distinguish between polymerase chain reaction [PCR] artefacts and true sequence variation, we performed deep TCRα and TCRβ repertoire profiling of the peripheral blood of 28 MZ twin pairs from Denmark and Germany, 24 of whom were discordant and four concordant for IBD.

RESULTS

We observed disease- and smoking-associated traits such as sharing, diversity and abundance of specific clonotypes in the TCR repertoire of IBD patients, and particularly in patients with active disease, compared with their healthy twins.

CONCLUSIONS

Our findings identified TCR repertoire features specific for smokers and IBD patients, particularly when signs of disease activity were present. These findings are a first step towards the application of TCR repertoire analyses as a valuable tool to characterise inflammatory bowel diseases and to identify potential biomarkers and true disease causes.

Immunotyping of clinically divergent p.Phe508del homozygous monozygous cystic fibrosis twins

Blood of the three clinically most concordant and most discordant p.Phe508del homozygous monozygous twin pairs of the European Cystic Fibrosis Twin and Sibling Study was examined in two postzygotic attributes that generate diversity between monozygous twins, i.e. the repertoire of the CDR3 region of the T-cell receptor ß chains and the DNA methylation at 450,000 genomic CpG sites. Methylation patterns in peripheral blood of twins changed at selected cell-type-independent positions and the immune cells of the twins showed individual profiles of the T cell receptor repertoire reflecting the plasticity of the immune system of genetically identical humans with cystic fibrosis to cope with the environment.

Lineage Tracking the Generation of T Regulatory Cells From Microbial Activated T Effector Cells in Naïve Mice

Regulatory T cells (Tregs) are essential for the maintenance of gut homeostasis by suppressing conventional CD4⁺ helper T cells (Tconvs) that are activated by microbial antigens. Although thymus is the major source of the peripheral Tregs, peripheral conversion from Tconvs to Tregs have also been shown to occur under various experimental conditions. It remains less clear about the frequency of lineage conversion from Tconvs to Tregs in naïve animals. Here we used a newly established reporter system to track a group of post expansion Tregs (eTregs), which exhibited a stronger suppressive ability than the non-lineage marked Tregs. Notably, microbial antigens are the primary driver for the formation of eTregs. TCR repertoire analysis of Peyer's patch T cells revealed that eTregs are clonally related to Tconvs, but not to the non-lineage tracked Tregs. Adoptive transfer of Tconvs into lymphopenic hosts demonstrated a conversion from Tconvs to eTregs. Thus, our lineage tracking method was able to capture the lineage conversion from microbial activated effector T cells to Tregs in naïve animals. This study suggests that a fraction of clonally activated T cells from the natural T cell repertoire exhibits lineage conversion to Tregs in response to commensal microbes under homeostatic conditions.

Determinants governing T cell receptor α/β-chain pairing in repertoire formation of identical twins

The T cell repertoire in each individual includes T cell receptors (TCRs) of enormous sequence diversity through the pairing of diverse TCR α- and β-chains, each generated by somatic recombination of paralogous gene segments. Whether the TCR repertoire contributes to susceptibility to infectious or autoimmune diseases in concert with disease-associated major histocompatibility complex (MHC) polymorphisms is unknown. Due to a lack in high-throughput technologies to sequence TCR α-β pairs, current studies on whether the TCR repertoire is shaped by host genetics have so far relied only on single-chain analysis. Using a high-throughput single T cell sequencing technology, we obtained the largest paired TCRαβ dataset so far, comprising 965,523 clonotypes from 15 healthy individuals including 6 monozygotic twin pairs. Public TCR α- and, to a lesser extent, TCR β-chain sequences were common in all individuals. In contrast, sharing of entirely identical TCRαβ amino acid sequences was very infrequent in unrelated individuals, but highly increased in twins, in particular in CD4 memory T cells. Based on nucleotide sequence identity, a subset of these shared clonotypes appeared to be the progeny of T cells that had been generated during fetal development and had persisted for more than 50 y. Additional shared TCRαβ in twins were encoded by different nucleotide sequences, implying that genetic determinants impose structural constraints on thymic selection that favor the selection of TCR α-β pairs with entire sequence identities.

An overview of immunoinformatics approaches and databases linking T cell receptor repertoires to their antigen specificity

Recent advances in molecular and bioinformatic methods have greatly improved our ability to study the formation of an adaptive immune response towards foreign pathogens, self-antigens, and cancer neoantigens. T cell receptors (TCR) are the key players in this process that recognize peptides presented by major histocompatibility complex (MHC). Owing to the huge diversity of both TCR sequence variants and peptides they recognize, accumulation and complex analysis of large amounts of TCR-antigen specificity data is required for understanding the structure and features of adaptive immune responses towards pathogens, vaccines, cancer, as well as autoimmune responses. In the present review, we summarize recent efforts on gathering and interpreting TCR-antigen specificity data and outline the critical role of tighter integration with other immunoinformatics data sources that include epitope MHC restriction, TCR repertoire structure models, and TCR/peptide/MHC structural data. We suggest that such integration can lead to the ability to accurately annotate individual TCR repertoires, efficiently estimate epitope and neoantigen immunogenicity, and ultimately, in silico identify TCRs specific to yet unstudied antigens and predict self-peptides related to autoimmunity.

No difference in TCRβ repertoire of CD4+ naive T cell between patients with primary biliary cholangitis and healthy control subjects

Primary biliary cholangitis (PBC) is considered as a model of organ-specific autoimmune disease based on the serological findings of anti-mitochondrial antibodies (AMA), infiltrates of T cells, and selective destruction of epithelial cells in the liver. T-cell-mediated autoimmune mechanisms are considered to be involved in the pathogenesis of primary biliary cholangitis (PBC). In this context, we used a combination of multiplex-PCR, Illumina sequencing and IMGT/HighV-QUEST for a standardized analysis of the T cell receptor β-chain (TCRβ) repertoire of CD4^＋naive T cells in PBC patients compared with healthy volunteers. Nonfunctional TCRs were used to study the pre-selection TCR repertoire, as they are not subject to functional selection (positive and negative selection). Functional TCRs were used to study the post-selection TCR repertoire. The results showed that there was not significant difference between PBC patients and healthy volunteers in TCRβ diversity, CDR3 length distributions, degree of sequence sharing, and usage frequency of TRBV and TRBJ segments, no matter in Pre-selection or Post-selection repertoires. In conclusion, early events in thymic T cell development and repertoire generation are not abnormality in PBC patients. The breakdown of self-tolerance to autoantigen may be derived from other immunological dysregulation or environmental agents.

Balancing sensitivity and specificity in distinguishing TCR groups by CDR sequence similarity

Background

Repertoire sequencing is enabling deep explorations into the cellular immune response, including the characterization of commonalities and differences among T cell receptor (TCR) repertoires from different individuals, pathologies, and antigen specificities. In seeking to understand the generality of patterns observed in different groups of TCRs, it is necessary to balance how well each pattern represents the diversity among TCRs from one group (sensitivity) vs. how many TCRs from other groups it also represents (specificity). The variable complementarity determining regions (CDRs), particularly the third CDRs (CDR3s) interact with major histocompatibility complex (MHC)-presented epitopes from putative antigens, and thus encode the determinants of recognition.

Results

We here systematically characterize the predictive power that can be obtained from CDR3 sequences, using representative, readily interpretable methods for evaluating CDR sequence similarity and then clustering and classifying sequences based on similarity. An initial analysis of CDR3s of known structure, clustered by structural similarity, helps calibrate the limits of sequence diversity among CDRs that might have a common mode of interaction with presented epitopes. Subsequent analyses demonstrate that this same range of sequence similarity strikes a favorable specificity/sensitivity balance in distinguishing twins from non-twins based on overall CDR3 repertoires, classifying CDR3 repertoires by antigen specificity, and distinguishing general pathologies.

Conclusion

We conclude that within a fairly broad range of sequence similarity, matching CDR3 sequences are likely to share specificities.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2864-8) contains supplementary material, which is available to authorized users.

Human T cell repertoire: what happens in thymus does not stay in thymus

The T cell receptor (TCR) repertoire is diverse, thus allowing recognition of a wide range of pathogens by T cells. In humans, the study of the formation of TCR repertoires is problematic because of the difficulty in performing investigations in vivo. In this issue of the JCI, Khosravi-Maharlooei and colleagues describe a new humanized mouse model that allows direct investigations on this topic. Using high-throughput and single-cell TCR-complementarity-determining region 3 β (TCR-CDR3β) sequencing, the authors were able to demonstrate that human thymic selection is a major driver of TCR sequence sharing, also implicating a preferential selection of shared cross-reactive CDR3βs during repertoire formation.

Broad cross-reactivity of the T-cell repertoire achieves specific and sufficiently rapid target searching

The molecular recognition of T-cell receptors is the hallmark of the adaptive immunity. Given the finiteness of the T-cell repertoire, individual T-cell receptors are necessary to be cross-reactive to multiple antigenic peptides. In this study, we quantify the variability of the cross-reactivity by using a string model that estimates the binding affinity between two sequences of amino acids. We examine sequences of 10,000 human T-cell receptors and 10,000 antigenic peptides, and obtain a full spectrum of cross-reactivity of the receptor-peptide binding. Then, we find that the cross-reactivity spectrum is broad. Some T-cells are reactive to 1000 peptides, but some T-cells are reactive to only one or two peptides. Since the degree of cross-reactivity has a correlation with the (un)binding affinity of receptors, we further investigate how the broad cross-reactivity affects the target searching of T-cells. High cross-reactive T-cells may not require many trials for searching correct targets, but they may spend long time to unbind from incorrect targets. In contrast, low cross-reactive T-cells may not spend long time to ignore incorrect targets, but they require many trials for screening correct targets. We evaluate this hypothesis, and show that the broad cross-reactivity of the natural T-cell repertoire can balance the trade-off between the rapid screening and unbinding penalty.

Crossreactive public TCR sequences undergo positive selection in the human thymic repertoire

We investigated human T-cell repertoire formation using high throughput TCRβ CDR3 sequencing in immunodeficient mice receiving human hematopoietic stem cells (HSCs) and human thymus grafts. Replicate humanized mice generated diverse and highly divergent repertoires. Repertoire narrowing and increased CDR3β sharing was observed during thymocyte selection. While hydrophobicity analysis implicated self-peptides in positive selection of the overall repertoire, positive selection favored shorter shared sequences that had reduced hydrophobicity at positions 6 and 7 of CDR3βs, suggesting weaker interactions with self-peptides than unshared sequences, possibly allowing escape from negative selection. Sharing was similar between autologous and allogeneic thymi and occurred between different cell subsets. Shared sequences were enriched for allo-crossreactive CDR3βs and for Type 1 diabetes-associated autoreactive CDR3βs. Single-cell TCR-sequencing showed increased sharing of CDR3αs compared to CDR3βs between mice. Our data collectively implicate preferential positive selection for shared human CDR3βs that are highly cross-reactive. While previous studies suggested a role for recombination bias in producing "public" sequences in mice, our study is the first to demonstrate a role for thymic selection. Our results implicate positive selection for promiscuous TCRβ sequences that likely evade negative selection, due to their low affinity for self-ligands, in the abundance of "public" human TCRβ sequences.

Localization-associated immune phenotypes of clonally expanded tumor-infiltrating T cells and distribution of their target antigens in rectal cancer

The degree and type of T cell infiltration influence rectal cancer prognosis regardless of classical tumor staging. We asked whether clonal expansion and tumor infiltration are restricted to selected-phenotype T cells; which clones are accessible in peripheral blood; and what the spatial distribution of their target antigens is. From five rectal cancer patients, we isolated paired tumor-infiltrating T cells (TILs) and T cells from unaffected rectum mucosa (T_UM) using 13-parameter FACS single cell index sorting. TCRαβ sequences, cytokine, and transcription factor expression were determined with single cell sequencing. TILs and T_UM occupied distinct phenotype compartments and clonal expansion predominantly occurred within CD8⁺ T cells. Expanded TIL clones identified by paired TCRαβ sequencing and exclusively detectable in the tumor showed characteristic PD-1 and TIM-3 expression. TCRβ repertoire sequencing identified 49 out of 149 expanded TIL clones circulating in peripheral blood and 41 (84%) of these were PD-1^- TIM-3^-. To determine whether clonal expansion of predominantly tumor-infiltrating T cell clones was driven by antigens uniquely presented in tumor tissue, selected TCRs were reconstructed and incubated with cells isolated from corresponding tumor or unaffected mucosa. The majority of clones exclusively detected in the tumor recognized antigen at both sites. In summary, rectal cancer is infiltrated with expanded distinct-phenotype T cell clones that either i) predominantly infiltrate the tumor, ii) predominantly infiltrate the unaffected mucosa, or iii) overlap between tumor, unaffected mucosa, and peripheral blood. However, the target antigens of predominantly tumor-infiltrating TIL clones do not appear to be restricted to tumor tissue.

Shorter TCR β-Chains Are Highly Enriched During Thymic Selection and Antigen-Driven Selection

The adaptive immune system uses several strategies to generate a repertoire of T cell receptors (TCR) with sufficient diversity to recognize the universe of potential pathogens. However, it remains unclear how differences in the T cell receptor (TCR) contribute to heterogeneity in T cell state. In this study, we used polychromatic flow cytometry to isolate highly pure CD4+/CD8+ naive and memory T cells, and applied deep sequencing to characterize corresponding TCR β-chain (TCRβ) complementary-determining region 3 (CDR3) repertoires. We find that shorter TCRβ CDR3s with fewer insertions were highly enriched during thymic selection. Antigen-experienced T cells (memory T cells) harbor shorter CDR3s vs. naive T cells. Moreover, the public TCRβ CDR3 clonotypes within cell subsets or interindividual tend to have shorter CDR3 length and a significantly larger size compared with "private" clonotypes. Taken together, shorter CDR3s highly enriched during thymic selection and antigen-driven selection, and further enriched in public T-cell responses. These results indicated that it may be evolutionary pressures drive short CDR3s to recognize most of antigen in nature.

Using T Cell Receptor Repertoires to Understand the Principles of Adaptive Immune Recognition

Adaptive immune recognition is mediated by antigen receptors on B and T cells generated by somatic recombination during lineage development. The high level of diversity resulting from this process posed technical limitations that previously limited the comprehensive analysis of adaptive immune recognition. Advances over the last ten years have produced data and approaches allowing insights into how T cells develop, evolutionary signatures of recombination and selection, and the features of T cell receptors that mediate epitope-specific binding and T cell activation. The size and complexity of these data have necessitated the generation of novel computational and analytical approaches, which are transforming how T cell immunology is conducted. Here we review the development and application of novel biological, theoretical, and computational methods for understanding T cell recognition and discuss the potential for improved models of receptor:antigen interactions.

Precise tracking of vaccine-responding T cell clones reveals convergent and personalized response in identical twins

T cell receptor (TCR) repertoire data contain information about infections that could be used in disease diagnostics and vaccine development, but extracting that information remains a major challenge. Here we developed a statistical framework to detect TCR clone proliferation and contraction from longitudinal repertoire data. We applied this framework to data from three pairs of identical twins immunized with the yellow fever vaccine. We identified 600 to 1,700 responding TCRs in each donor and validated them using three independent assays. While the responding TCRs were mostly private, albeit with higher overlap between twins, they could be well-predicted using a classifier based on sequence similarity. Our method can also be applied to samples obtained postinfection, making it suitable for systematic discovery of new infection-specific TCRs in the clinic.

TCR repertoire profiling of tumors, adjacent normal tissues, and peripheral blood predicts survival in nasopharyngeal carcinoma

The T-cell immune responses in nasopharyngeal carcinoma patients have been extensively investigated recently for designing adoptive immunotherapy or immune checkpoint blockade therapy. However, the distribution characteristics of T cells associated with NPC pathogenesis are largely unknown. We performed deep sequencing for TCR repertoire profiling on matched tumor/adjacent normal tissue from 15 NPC patients and peripheral blood from 39 NPC patients, 39 patients with other nasopharyngeal diseases, and 33 healthy controls. We found that a lower diversity of TCR repertoire in tumors than paired tissues or a low similarity between the paired tissues was associated with a poor prognosis in NPC. A more diverse TCR repertoire was identified in the peripheral blood of NPC patients relative to the controls; this was related to a significant decrease in the proportion of high-frequency TCR clones in NPC. Higher diversity in peripheral blood of NPC patients was associated with a worse prognosis. Due to the peculiarity of the Vβ gene usage patterns in the peripheral blood of NPC patients, 15 Vβ genes were selected to distinguish NPC patients from controls by the least absolute shrinkage and selection operator analysis. We identified 11 clonotypes shared by tumors and peripheral blood samples from different NPC patients, defined as "NPC-associated" that might have value in adoptive immunotherapy. In conclusion, we here report the systematic and overall characteristics of the TCR repertoire in tumors, adjacent normal tissues, and peripheral blood of NPC patients. The data obtained may be relevant to future clinical studies in the setting of immunotherapy for NPC patients.

The landscape and diagnostic potential of T and B cell repertoire in Immunoglobulin A Nephropathy

Immunoglobulin A Nephropathy (IgAN) is the most common glomerulonephritis worldwide. The pathologic hallmark of IgAN is immune complex deposited in glomerular mesangium, which induces inflammation and affects the kidney's normal functions. The exact pathogenesis of IgAN, however, remains obscure. Further, in current clinical practice, the diagnosis relies on needle biopsy of renal tissue. Therefore, a non-invasive method for diagnosis and prognosis surveillance of the disease is highly desirable. To this end, we investigated the T cell receptor beta chain (TCRB) and immunoglobulin heavy chain (IGH) repertoire in circulating lymphocytes and compared them with kidney infiltrating lymphocytes using immune repertoire high throughput sequencing. We found that some features of TCRB and IGH in renal tissues were remarkably different from that in the blood, including decreased repertoire diversity, increased IgA and IgG frequency, and more antigen-experienced B cells. The complementarity-determining region 3 (CDR3) length of circulating TCRB and IGH in IgAN patients was significantly shorter than that in healthy controls, which is the result of both VDJ rearrangement and clonal selection. The IgA1 frequency in the blood of IgAN patients is significantly higher than that in other Nephropathy (NIgAN) patients and healthy control. Importantly we identified a set of TCRB and IGH clones, which can be used to distinguish IgAN from NIgAN and healthy controls with high accuracy. These results indicated that the TCRB and IGH repertoire can potentially serve as non-invasive biomarkers for the diagnosis of IgAN. The characteristics of the kidney infiltrating and circulating lymphocytes repertoires shed light on IgAN detection, treatment and surveillance.

High-throughput T cell receptor sequencing identifies clonally expanded CD8+ T cell populations in alopecia areata

Alopecia areata (AA) is an autoimmune disease in which cytotoxic T cells specifically target growing hair follicles. We used high-throughput TCR sequencing in the C3H/HeJ mouse model of AA and in human AA patients to gain insight into pathogenic T cell populations and their dynamics, which revealed clonal CD8+ T cell expansions in lesional skin. In the C3H/HeJ model, we observed interindividual sharing of TCRβ chain protein sequences, which strongly supports a model of antigenic drive in AA. The overlap between the lesional TCR repertoire and a population of CD8+NKG2D+ T cells in skin-draining lymph nodes identified this subset as pathogenic effectors. In AA patients, treatment with the oral JAK inhibitor tofacitinib resulted in a decrease in clonally expanded CD8+ T cells in the scalp but also revealed that many expanded lesional T cell clones do not completely disappear from either skin or blood during treatment with tofacitinib, which may explain in part the relapse of disease after stopping treatment.

The expanding role of systems immunology in decoding the T cell receptor repertoire

T cells play a crucial role in the immune system's defense against many infectious diseases, including persistent infections for which no effective vaccines currently exist. The T cell component of the adaptive immune system is highly complex involving a constantly evolving landscape of various inter-related T cell populations. These T cell populations are characterized by their phenotypic and functional properties as well as the collection, or repertoire, of T cell receptors (TCR) that mediate T cell recognition of antigenic peptides derived from pathogens. Understanding the various processes and factors that impact the development and evolution of the broader T cell repertoire available to recognize and respond to pathogens and the characteristics of antigen-experienced T cell repertoires associated with effective immune control of pathogens is critical to the rational design of T cell-based vaccines and therapies. In this article we discuss, using examples of recent research, the promise that systems immunology approaches, involving quantitative analysis and mathematical and computational modeling of immunological data, hold for decoding the complex TCR repertoire system in the current era of advancing technologies.

Memory B Cells Activate Brain-Homing, Autoreactive CD4+ T Cells in Multiple Sclerosis

Multiple sclerosis is an autoimmune disease that is caused by the interplay of genetic, particularly the HLA-DR15 haplotype, and environmental risk factors. How these etiologic factors contribute to generating an autoreactive CD4⁺ T cell repertoire is not clear. Here, we demonstrate that self-reactivity, defined as “autoproliferation” of peripheral Th1 cells, is elevated in patients carrying the HLA-DR15 haplotype. Autoproliferation is mediated by memory B cells in a HLA-DR-dependent manner. Depletion of B cells in vitro and therapeutically in vivo by anti-CD20 effectively reduces T cell autoproliferation. T cell receptor deep sequencing showed that in vitro autoproliferating T cells are enriched for brain-homing T cells. Using an unbiased epitope discovery approach, we identified RASGRP2 as target autoantigen that is expressed in the brain and B cells. These findings will be instrumental to address important questions regarding pathogenic B-T cell interactions in multiple sclerosis and possibly also to develop novel therapies.

•

Autoproliferation of CD4⁺ T cells and B cells is involved in multiple sclerosis

•

The main genetic factor of MS, HLA-DR15, plays a central role in autoproliferation

•

Memory B cells drive autoproliferation of Th1 brain-homing CD4⁺ T cells

•

Autoproliferating T cells recognize antigens expressed in B cells and brain lesions

Predicting the spectrum of TCR repertoire sharing with a data-driven model of recombination

Despite the extreme diversity of T-cell repertoires, many identical T-cell receptor (TCR) sequences are found in a large number of individual mice and humans. These widely shared sequences, often referred to as "public," have been suggested to be over-represented due to their potential immune functionality or their ease of generation by V(D)J recombination. Here, we show that even for large cohorts, the observed degree of sharing of TCR sequences between individuals is well predicted by a model accounting for the known quantitative statistical biases in the generation process, together with a simple model of thymic selection. Whether a sequence is shared by many individuals is predicted to depend on the number of queried individuals and the sampling depth, as well as on the sequence itself, in agreement with the data. We introduce the degree of publicness conditional on the queried cohort size and the size of the sampled repertoires. Based on these observations, we propose a public/private sequence classifier, "PUBLIC" (Public Universal Binary Likelihood Inference Classifier), based on the generation probability, which performs very well even for small cohort sizes.