Development of a diverse TCR repertoire in the thymus is needed to provide T cell-mediated protection against a wide range of pathogens. It has been suggested that, as a result of somatic rearrangements and pairing of α and β TCR chains, 10¹⁵ to 10¹⁸ unique TCRαβ sequences could theoretically form (12, 13). However, only 2 × 10¹¹ T cells exist in a human body (12, 13). With so much potential diversity and limited space for T cells in the peripheral immune system, trimming of the repertoire during positive and negative selection processes in the thymus plays an important role (14, 15). These selection processes ensure optimal self MHC-restricted T cell responses to foreign peptides and prevent the development of autoimmune diseases, respectively. A next generation sequencing study estimated that there are approximately 100 million unique TCRβ sequences in a young adult human, which declined two- to five-fold in healthy elderly people. This age-related contraction was less pronounced for memory than naïve T cells (16), reflecting the decline in thymic function with age. Another deep sequencing study has estimated that there are 40 × 10⁶ to 70 × 10⁶ unique TCRβ sequences and 60 × 10⁶ to 100 × 10⁶ TCRα sequences in human pediatric thymi (17).

Diversity of the TCR repertoire is formed at different levels, including random recombination of V, (D) and J genes for both TCRα and β chains, addition and deletion of random nucleotides at V (D) J junctional sites mediated by the enzyme terminal deoxynucleotidyl transferase (TdT) in the β chain, and pairing of TCRα and β chains (18). The processes of T cell selection in the thymus are still not completely understood, with vastly more information being available from murine studies than from humans.

Early in the T cell development process in the thymus, recombination activating enzymes RAG-1 and RAG-2 are upregulated and V, D and J genes of a β chain gene recombine, permitting β selection to occur ( Figure 1A ). In human fetal thymocytes, RAG enzymes were detected after a round of replication during the double negative (DN) and another round of replication during the double positive (DP) stage (20). Upon productive rearrangement, the formed β chain pairs with a pre-TCRα molecule and relays a signal for proliferation to the developing T cell at the late DN proliferative stage (20, 21) ( Figure 1A ). After several rounds of proliferation, RAG enzymes are expressed again and V and J genes of an α chain recombine, allowing CD4⁺ CD8⁺ DP cells to be formed. The expression of an αβ TCR allows these DP cells to undergo positive and negative selection and differentiate into single positive (SP) CD4 and CD8 T cells ( Figure 1A ). A fully formed polyclonal β chain repertoire is detected at the DP CD3^- stage, while polyclonal α chain rearrangements are only evident at the DP CD3^low stage and later (22). It is estimated that only 3% of thymocytes reach a mature state and exit the thymus (22).

For human thymocytes, the ratio of nonproductive to productive recombination events was higher prior to the proliferative phase of the DN stage, indicating that thymocytes that fail to generate a productive β chain for the first allele try another recombination with the other allele. However, the ratio of nonproductive recombinations was not increased for TCRα during the DP stage, suggesting that cells with a nonproductive TCRα may not undergo a second allele recombination, but instead die (20). However, murine studies have demonstrated that a second productive α chain rearrangement can occur (23, 24), resulting in the presence of two α chains and two alternative TCRs in one thymocyte. Observations in human T cells (25) and thymocyte selection studies (26) indicate that a similar situation pertains.

Studies of peripheral human TCR repertoires have established that CD4 cells are more diverse compared to CD8 T cells (16, 27–29), possibly due to differences between the subsets in thymic output and/or levels of clonal expansion/persistence in response to infections.

By comparing TCRβ repertoires between different individuals, the level of sharing of TCRβ sequences was found to be considerably higher than that expected based on the number of clonotypes (30). The level of TCRα sharing has been shown to be even larger than TCRβ sharing among healthy humans (31). Therefore, there may be genetic factors or selection forces that drive the sharing of certain TCR sequences. However, one limitation of most of these studies is their reliance on analysis of peripheral T cell repertoires rather than thymocytes, which are less accessible. Peripheral TCR repertoires have undergone post-thymic selection that reflects environmental factors such as antigenic exposures, potentially obscuring the role of genetic factors.

Monozygotic twin studies have been used in efforts to gain insight into the role of genetic factors in TCR repertoire generation in humans. All of these studies point to a very diverse and divergent TCR repertoire with low levels of sharing between unrelated individuals and even between identical twins. One analysis of peripheral TCR repertoires of three pairs of monozygotic twins revealed similar degrees of sharing between identical twins and unrelated individuals (32). However, CDR3 sharing between identical twins was significantly higher within the fraction of highly abundant TCRs. Analysis of out-of-frame sequences provides a window into the TCR repertoire prior to selection, as these sequences are not subjected to selection processes. In this study, V gene usage in out-of-frame sequences was more similar in identical twins than unrelated individuals, pointing to a genetically determined bias prior to selection (32). Another study involved high throughput single T cell sequencing on 15 healthy individuals, including 6 monozygotic twin pairs. Although sharing of entire TCR including α and β chains was infrequent, such sharing was greater between identical twins than between others, which could suggest that genetic components have a role in TCR repertoire selection or could point to shared environmental exposures (33).

There is controversy in the literature about the relative roles of genetic background and thymic selection in determining V (D) J gene usage (32, 34–37). Overall, a general pattern is shared throughout human populations, such that certain genes are always used more than others (38). However, genetic factors lead to subtle differences. To discover the associations between genetic variation and TCR V gene usage, Sharon et al. applied expression quantitative trait locus (eQTL) mapping in a large human cohort and discovered associations between variation in the HLA locus and TCR V gene usage (39). Other studies established increased similarity in V and J gene usage between identical twins in both peripheral (40) and thymic T cells (41), despite the overall similarity between any two individuals. A more recent study using single cell RNA sequencing elucidated further details about the process of V (D) J gene recombination and TCR α and β chain usage (20). For TCR β, an overall pattern of VDJ recombination that was established during the DN stage (after β selection) was shown to persist during subsequent developmental stages in the thymus. This pattern was correlated with the genomic positions and looping structures of different loci. Although this initial pattern formed the overall repertoire, thymic selection led to enrichment or deletion of some Vβ genes at developmental stages following β selection. On the other hand, an association between developmental timing and V-J pairing was detected for the TCRα locus. Proximal V-J pairs were formed in early stages, followed by distal pairs in later stages. Proximal pairs declined relatively as thymic selection progressed, which points to progressive recombination of the TCRα locus ( Figure 1B ). Interestingly, SP CD8 T cells were enriched for the distal pairs, suggesting a slower or less efficient commitment toward the CD8 T lineage (20).

TCRs shared between different individuals are also called “public TCRs”. Sharing of TCRβ chains in peripheral blood CD8 T cells has been shown to be 500 times greater than what is expected from TCR formation from random recombinations (30). High throughput TCRβ sequencing of 80 million TCR sequences from 666 people revealed that around 14% of TCRβs are public (42). Among DP thymocytes in syngeneic mice, about 15% of TCRs were found to be public, and most of them were detected in the preselection repertoire, suggesting recombinatorial bias (43). In another mouse study, CDR3β sequencing performed in TCRα^-/- mice, in which functional TCRαβ T cells do not form and therefore thymic selection does not occur, revealed that 11% of sequences were shared (44), again suggesting that recombinatorial bias is largely responsible for development of public TCRs. A third study involved high throughput TCRβ sequencing on a large group of mice, including syngeneic and allogeneic MHCs. Shared CDR3s were enriched among the most abundant TCRs. A group of 289 public CDR3βs constituted about 10% of the whole repertoire. Public CDR3βs were highly convergent (having high ratios of nucleotide sequences per amino acid sequence), utilized restricted V/J segments, were found in mice with different MHCs, and were shorter compared to private sequences. Computer simulations of the VDJ recombination process generated public TCRs similar to the actual dataset. These authors concluded that both recombinatorial bias and thymic selection play a role in the existence of public TCRs (45). A recent study of TCR α and β sequences in a cohort of humans showed that shared sequences are more common in productive compared to nonproductive sequences, pointing to the effect of thymic selection in enrichment of public TCRs, while gene segment usage between unrelated individuals suggested a role for recombinatorial bias as well (38).

Productive human CDR3βs have been shown to be shorter on average than out-of-frame sequences, indicating that thymic selection favors shorter sequences (27). Consistently, direct analyses of murine thymocytes have shown that mouse CDR3αs shorten as thymic selection progresses (46). The observation that public CDR3βs are shorter compared to private ones in human TCRβ repertoires suggests a possible role for thymic selection of public sequences. The observation that the same public sequences develop independently of MHC genotypes both in mice (45) and in humans (32) suggests that these sequences may be crossreactive and can interact with different peptide-MHC (p.MHC) complexes. Additionally, significant CDR3β sharing has been demonstrated between different peripheral blood human T cell subsets (47), supporting the hypothesis that shared sequences are highly crossreactive. A hypothetical increasing pattern of public TCRs during thymic and peripheral selection processes is shown in Figure 1C . Our HIS mouse studies, described in Lessons Learned From HIS Mouse Models and Figure 1D , provide experimental evidence in support of this hypothetical pattern.

Accumulating evidence indicates that public TCRs are enriched for autoreactive and cross-reactive sequences and vice versa. In a single cell TCR sequencing study, analysis of more than 3000 gluten-specific CD4 TCR clonotypes isolated from 63 celiac patients showed that gluten-specific autoreactive T cells were enriched for public TCRs (48). By comparing thymic TCRα sequences from 6 human donors to a list of sequences known for their reactivity for type 1 diabetes (T1D)-associated self-antigens and HIV antigens, Heikila et al. showed that autoreactive (T1D-associated) TCRs were generated in significantly higher numbers compared to nonself-reactive (HIV-associated) TCRs in human thymi and were not negatively selected or converted to Tregs before entering the periphery. These autoreactive TCRs were selected by different thymi, independent of their HLA type (49). In a murine study, TCRβ sequencing of T cells in thymus, periphery and central nervous system of mice with experimental autoimmune encephalomyelitis implicated public TCRs in autoimmunity. Interestingly, these autoreactive public TCRs were preferentially identified in the preselection repertoire, pointing to recombinatorial bias as the main factor in their development (50). T cells in T1D patients were shown to have shorter CDR3βs with increased deletions/reduced insertions during VDJ rearrangement compared to healthy controls, even in nonproductive sequences, implicating recombinatorial bias in the development of autoreactive T cells in T1D. CDR3s in autoreactive CD4 T cell clones in these patients were shorter than those in T cells recognizing viral antigens (51).

Flexible CDR3 loops allow TCRs to bind multiple p.MHC complexes and hence be degenerate in their antigen recognition (13, 52). This phenomenon is needed for T cells to recognize a repertoire of p.MHCs that is larger than the repertoire of possible TCRs (13). A single T1D-reactive TCR (1E6) was shown by experimental and computational analysis to be capable of binding to more than a million peptides presented on HLA-A2 (53). Available evidence indicates that shorter TCRs have a greater chance of being promiscuous and cross-reactive. Recombinatorial bias and thymic selection appear to favor these sequences, which are enriched for autoreactive TCRs. Although TCR crossreactivity allows recognition of some pathogens whose peptides are not presented on cortical thymic epithelial cells (cTECs), peptides containing TCR contact amino acid motifs not present in the human proteome were recently reported to be unrecognizable by naïve human T cells, resulting in potential holes in the TCR repertoire for recognition of potential pathogens (54).

TCR repertoires of CD4 and CD8 T cells have been found to have little CDR3β overlap and to have distinct characteristics, including the preference for certain amino acids in their CDR3βs (29). However, a recent single cell TCR sequencing study detected overlap between TCR repertoires of CD4 and CD8 T cells of the same individuals at CDR3α (7.8%), CDR3β (4.7%) and even in paired αβ (0.65%) (55). Although the overlapping repertoire at the paired αβ level is minimal, it might be enriched for cross-reactive sequences that are capable of recognizing both MHC-I and II molecules, perhaps with sufficient affinity that coreceptor expression is irrelevant. In the section on HIS mouse models below, we discuss evidence that public CDR3β sequences may be associated with such activity, as the same amino acid sequences are frequently detected among both CD4 and CD8 populations and in the context of completely disparate thymic HLA molecules (19).

Access to human thymus tissue is limited and, in contrast to mice, genetics and other variables cannot be readily manipulated in order to study the development of the human T cell repertoire. HIS mouse models involving transplantation with human fetal thymus tissue along with human hematopoietic stem cells (HSCs), however, provide an opportunity to study the factors involved in human TCR repertoire selection. Several studies have used HIS mouse models to study V(D)J recombination and CDR3 length distribution (56–59). These studies established that the human TCR repertoire in these HIS mice is diverse and has normal CDR3 length distributions.

In order to investigate the factors involved in human TCR repertoire formation, we recently performed high throughput TCRβ sequencing and single cell TCR sequencing on thymic and peripheral T cells of HIS mice generated with human fetal liver HSCs combined with either autologous or HLA-mismatched allogeneic human fetal thymus tissue (19). Measures were taken to deplete the grafted fetal thymus tissue of pre-existing thymocytes so that synchronized thymic development from progenitors could be investigated in the mice. To ensure that no human thymopoiesis occurred in the native mouse thymus of recipient NSG mice, which could confound interpretation of peripheral TCR repertoire data, we used a rapid method to thymectomize these mice (60).

Our studies revealed that the TCRβ repertoire of different human thymocyte populations is highly diverse (19) ( Figure 1D ). Formation of the human TCR repertoire was shown to be largely stochastic and was almost totally divergent between animals reconstituted simultaneously with identical HSCs, thymus, genetic backgrounds and environment. This finding may help to explain the incomplete penetrance of genetically-controlled autoimmune diseases in identical twins (61). In contrast to the original fetal thymic tissue, in which expression of the TdT enzyme was low, this enzyme was well-expressed in grafted human thymi after they matured, especially in DP thymocytes. Consequently, diversity of the TCR repertoire was significantly greater in grafted human thymi compared to the original human fetal thymus (19). Diversity was lower at the amino acid sequence level compared to the nucleotide sequence level and also in productive sequences compared to non-productive sequences for populations that had undergone selection, demonstrating that thymic selection narrows the TCR repertoire, both at the DP and the SP stages of development. Diversity of TCR repertoire decreased as thymic selection progressed and also further decreased in peripheral compared to thymic SP T cell populations (19) ( Figure 1D ).

We also examined the development of public TCRβ CDR3 sequences using the HIS mouse model. An increase in the proportion of these sequences was observed in the transition from unselected (CD69-) to positively selected (CD69+) DP thymocytes, indicating a bias for their positive selection. We also observed a trend toward an increase in the proportion of shared sequences in SP compared to DP thymocytes and in the peripheral compared to thymocyte SP repertoires (19) ( Figure 1D ). This trend suggests that both positive and negative selection processes, as well as peripheral selection, favor public sequences. The proportion of overlapping CDR3βs was not different between or among autologous (to the HSCs) and allogeneic thymi and many CDR3βs were shared between CD4 and CD8 SP thymocyte subsets, suggesting that the shared CDR3βs, like those discussed above, are cross-reactive. The shared sequences had shorter CDR3β lengths, and we observed that the length of more abundant CDR3βs decreased during thymic selection, indicating that HLA-peptide-driven selection of human thymocytes favors shorter sequences, likely due to increasing cross-reactivity of shorter CDR3s. By comparing the shared and unshared CDR3βs with a repertoire of known T1D autoreactive CDR3β sequences, we found direct evidence that the shared sequences are enriched for autoreactive sequences. Furthermore, they were greatly enriched for CDR3βs that were found in a small subset of TCRs that recognized two completely HLA-mismatched allogeneic stimulator populations, indicating that they are crossreactive to different alloantigens (19). We were able to specifically analyze the impact of selection on shared human thymocyte sequences and demonstrated that thymic selection is a significant factor driving their abundance among mature T cells.

It has been shown in a murine model that TCRβ CDR3 hydrophobicity at position 6 and position 7 (P6 and P7), the residues that interface with the p.MHC, correlated with the ability to be activated by self p.MHC (62). In our HIS mouse studies, we observed that as selection progressed, hydrophobic amino acids were enriched at position 6 and 7 of CDR3β, suggesting a role for self-peptides in human thymocyte selection. However, relatively lower hydrophobicity at P6/P7 in the shared sequences suggested that they may have weaker interactions with self-peptides than unshared sequences and that this may allow them to escape negative selection (19).

Despite the very divergent CDR3β TCR repertoires, we found that V and J gene usage in different cell populations in human thymi was overall very similar among different mice generated with different HSCs and different thymi. However, our study was not powered to detect subtle differences in V gene usage such as those that have been observed in association with HLA alleles. Our analysis of thymocyte subsets demonstrates a significant but minor role for thymic selection in determining ultimate V gene usage. Additionally, we found that the observed VJ pairing pattern correlated very well with what would be expected from stochastic VJ pairing according to their frequencies, suggesting that no major factor biases the process of TCRβ V-J recombinations (19).

High throughput methods for paired sequencing of TCR α and β chains will be needed to study selection of clonotypes in high quantities. Existing approaches for simultaneous sequencing of both TCR α and β chains are expensive and have relatively low throughput (20, 63). Development of computational methods to cluster TCR sequences according to their likely p.MHC specificity could be ground-breaking in understanding thymic TCR selection processes. GLIPH (grouping of lymphocyte interactions by paratope hotspots) was one of the first algorithms that was developed for this purpose and has the capability to group raw TCR sequences into likely p.MHC specificity groups (64). While the degeneracy of TCR recognition limits the ability to identify specific peptide antigens with this approach, the clustering may provide a fingerprint of the forces driving thymocyte selection and alloreactivity. Further progress in this area (65) is likely to greatly enhance our understanding of this process.

Alloresponses, immune responses with unique strength and diversity, are directed to non-self antigens from members of the same species due to extensive genetic polymorphisms between allogeneic donors and recipients. Of primary importance in driving the strength of the alloresponse is the MHC (66), as discussed earlier in this article. Consequences of allorecognition in transplantation include graft rejection and GVHD, two major barriers to the greater success of solid organ and hematopoietic stem cell transplantation.

Studies to track alloresponses in vivo in rodents utilized tools like endogenous superantigen-specific Vβs and alloreactive transgenic TCRs. Intrathymic deletion was thereby demonstrated as a major mechanisms of tolerance in mixed or fully allogeneic chimeric mice (67–69). Given the inability to perform such studies in vivo in humans, efforts have been made to study a small population of donor-reactive T cells via tetramers with known HLA-peptide specificities (70, 71) or to use in vitro functional assays, such as MLRs (72, 73), CML assays (74), limiting diluting assays (LDAs) (75) and enzyme-linked immunospot (ELISPOT) assays (76, 77), under particular culture conditions. However, the enormous number of putative HLA/peptide combinations potentially recognized by alloreactive TCRs precludes using the tetramer approach to study a significant fraction of the alloreactive repertoire. In vitro functional assays, on the other hand, cannot replicate the in vivo setting and have demonstrated inconsistent correlations with clinical outcomes (74, 78–80). Previous attempts to analyze repertoire diversity using conventional methods, including PCR-based CDR3 length analysis by spectratyping and flow cytometry-based V region usage, were limited by low resolution (81–83). Recently, the availability of next-generation sequencing (NGS) technologies (84–86) has opened up new possibilities for understanding this immune repertoire and to tackle fundamental questions of T cell alloimmunity in humans.

The availability of high-throughput sequencing platforms for TCRβ CDR3 sequencing allowed our group to develop a novel approach to specifically identify alloreactive T cell repertoires before transplantation (Tx). Combing CFSE-MLR and high-throughput TCRβ CDR3 DNA sequencing (TCRβ-seq) allows us to track alloreactive TCR clones in different tissues following the Tx without relying on any post-Tx functional readout (80). In normal healthy donors, we measured the clonal frequencies, specificity and repertoire diversity of alloreactive CD4 and CD8 T cells. Our data showed that the alloimmune repertoire is highly specific for a given pair of individuals, with only a few exceptional TCRβ CDR3s crossreacting on disparate HLAs (5). Most alloreactive clones are present at low frequency in the circulation, reflecting the naïve phenotype of many alloreactive T cells. The diversity of this alloreactive repertoire is quite high and increases with increased HLA mismatching. Our analysis of the contribution of naïve and memory T cells to the alloresponse (5) confirmed previous functional studies in which it was concluded that both types of T cells contribute to alloreactivity measured in vitro (87, 88). Although the diversity of memory T cell repertoires is much lower than that of naïve T cells (89), pathogen-specific memory T cells can cross-react with alloantigens (90).

We developed a statistical model to estimate the total frequency of both detected and undetected alloreactive clones within the unstimulated circulating repertoire. With this approach, we estimated that a substantial fraction of human circulating TCR repertoire (0.5%-6%) has alloreactive potential to respond to just two different allogeneic stimulators and thereby inferred that most human TCRs likely have potential alloreactivity against some allogeneic HLA alleles (5).

This CFSE-MLR + TCRβ-seq approach (91) has been utilized to investigate human alloresponses and validated in recipients of several different types of organ transplantation, including combined kidney and bone marrow transplantation (CKBMT) (80), conventional kidney transplantation (92), liver transplantation (93), and intestinal transplantation (ITx) (94–96). Through these studies, we have validated the ability of our assay to identify biologically relevant alloreactive T cell clones in pre-Tx MLR that are involved in rejection and expand in patients following transplantation.

In tolerant patients in an Immune Tolerance Network (ITN)-sponsored CKBMT study (ITN036ST), we observed significant and specific long-term declines in the number of circulating donor-reactive T cell sequences. This decline, in contrast to functional assays, distinguished the tolerant patients from a patient on the same protocol who failed to achieve tolerance (80). These studies clearly distinguished clonal deletion from anergy, a distinction which often cannot be on the basis of in vitro functional assays. Expansion of circulating donor-reactive sequences was detected in conventional kidney Tx recipients and did not appear to correlate with rejection (80, 97).

We tracked donor‐reactive TCRβ sequences from a group of tolerant (n=3) and nontolerant (n=5) patients from a randomized liver transplantation clinical trial ITN030ST and obtained evidence consistent with deletion of donor‐reactive T cells post-Tx in both groups (93). The fact that clonal reductions were greater in donor-reactive TCRβ sequences than third party-reactive sequences in all 8 subjects, regardless of successful early immunosuppression withdrawal, demonstrates an impact of the liver allograft on the donor-reactive T cell repertoire and suggests that non-deletional mechanisms may be decisive in determining outcomes of immunosuppression withdrawal following liver Tx. The difference between this result and the observation of deletion only in tolerant patients in the CKBMT study highlights the protocol- and organ-specific nature of tolerance mechanisms in humans and the need to avoid considering allograft tolerance as a single entity.

In patients receiving ITx, we found that a faster rate of donor T cell replacement by recipient T cells in the graft mucosa correlated with early rejection, which was associated with a preponderance of host-versus-graft (HvG) T cell clones within the graft, demonstrating the biological significance of the TCRs identified as donor-reactive by pre-transplant MLR (94). Given the high lymphoid cell load carried within intestinal allografts, we hypothesized that clinical outcomes in ITx are largely determined by the exchange of donor and recipient lymphoid tissue and hence the balance of graft-versus-host (GvH) and HvG-reactive T cells. Using CFSE-MLR and the high-throughput TCRβ-seq-based approach to track alloreactive T cells in the GvH and HvG directions in the graft, circulation and bone marrow (BM), we obtained data consistent with this hypothesis (96). ITx provides a unique opportunity to study two-way alloresponses in humans and we did so by applying immunogenomic tools, including single cell sequencing techniques, to enhance our understanding of immune mechanisms after ITx, as detailed in the following sections.

ITx serves as the only long-term option for patients who suffer intestinal failure and have functional disorders (98, 99). However, the success of ITx is currently limited by high rejection rates, by the risk of GVHD and by infections and posttransplant lymphoproliferative disease (PTLD) resulting from high levels of immunosuppression (100, 101). Previous reports suggest that composite allograft transplants, such as multivisceral transplantation (MVTx) and liver-intestinal transplantation (LITx), are associated with reduced rates of intestinal rejection compared to isolated intestinal transplantation (iITx) (102, 103). Our own data have confirmed this observation in association with high levels of multilineage peripheral blood chimerism in MVTx recipients (94–96, 104).

Remarkably, we detected the presence of multilineage peripheral blood chimerism in many of these recipients and, further, demonstrated that this chimerism reflects de novo lymphopoiesis from hematopoietic stem and progenitor cells (HSPCs) carried in the intestinal allograft (95). We showed that high levels of peripheral blood T cell mixed chimerism (macrochimerism: ≥4% peak donor T cells in recipient PBMCs) occurs commonly, without GVHD, in recipients of intestinal allografts, and is associated with significantly reduced graft rejection, less de novo donor-specific antibody (DSA) development and slower recipient T cell repopulation in the graft (94–96, 104). The association between blood T cell macrochimerism and slower graft T cell replacement suggests a connection between immunologic events in the graft (locally) and the blood (systemically).

We investigated the underlying mechanisms of these phenomena at both the cellular and the clonotypic levels, using our published approach of combining CFSE-MLR and high-throughput TCRβ-seq  (80) to track alloreactive T cells in the GvH and HvG directions in the graft, circulation and BM (96). These studies have provided new insights into the significance of donor T cell macrochimerism in blood, including the first direct demonstration of lymphohematopoietic graft-versus-host responses (LGVHR) in humans, in which the GvH response confined to the lymphohematopoietic system promotes engraftment without causing GVHD, a disease that requires T cell migration into epithelial tissues. We found that patients with macrochimerism have high GvH to HvG T cell clonal ratios in their allografts. These GvH clones enter the circulation, where their peak levels are associated with declines in HvG clones early post-Tx, suggesting that LGVHR may contribute to chimerism and control HvG responses systemically. Consistent with the notion that LGVHR creates hematopoietic “space” in the BM for donor hematopoietic cells, as described in murine studies (105, 106), donor-derived T cells, including GvH clones and CD34⁺ HSPCs, were simultaneously detected in the recipients’ BM >100 days post-Tx (96). Individual GvH clones were detected in either the ileal mucosa or PBMCs before detection in recipient BM, consistent with an origin in the intestinal allograft, where donor GvH-reactive T cells expand early upon entry of recipient antigen-presenting cells (APCs) into the graft. These results, combined with cytotoxic single-cell transcriptional profiles of donor T cells in recipient BM, suggest that tissue-resident GvH-reactive donor T cells develop effector function and migrate into the recipient circulation and BM, where they promote engraftment of graft-derived HSPCs and maintain mixed chimerism.

Outcomes of ITx are suboptimal in large part due to high rejection rates. We have previously demonstrated that pre-Tx MLR-defined HvG alloreactive T-cell clones are enriched in intestinal allografts during early rejection and persist long-term post-Tx despite rejection resolution (94). Recipient T cells in the mucosa eventually acquire tissue resident memory (TRM) features and seed the entire gastrointestinal tract, likely posing a constant threat of late rejection (94). CD8 intraepithelial lymphocytes (IELs) show a TRM phenotype during quiescence (CD69⁺CD103^+/-CD28^low), but regain high CD28 expression, despite maintenance of the CD103 and CD69 TRM markers, during late rejection, indicating a possible transitional phenotype of TRM and effector T (Teff) cells that correlates with rejection. Therefore, we hypothesize that HvG-reactive T cells join the TRM pool and may become Teff cells and/or be tolerized. Alternatively, HvG clones developing de novo after ITx may contribute to late rejection.

We have established a comprehensive platform to integrate T cell clonal tracking, alloreactivity and RNA profiling to allow us to directly address these hypotheses. In ongoing studies, we have carried out single cell analysis of recipient T cells in the graft mucosa (107). Single cell RNA sequencing (scRNA-seq) was performed using the 10x Genomics platform and the Seurat analysis pipeline (108) for simultaneous measurement of mRNA expression and paired V(D)J TCR α and β sequences at the single cell level. Our published protocol using Adaptive Biotechnologies TCRβ bulk sequencing to identify HvG-reactive and nonHvG-reactive TCR repertoires (80, 94) from pre-Tx recipient lymphoid tissues was applied and single cell TCRβ sequences from intestinal allograft mucosal specimens were mapped to these pre-Tx sequence sets. This methodology allows annotation of each cell, identifying CD4, CD8, HvG and nonHvG clones of each type. To further extend our ability to investigate the induction of immune tolerance, we also set up post-Tx MLR assays using recipient PBMCs post-Tx as responders against donor pre-Tx stimulators and sorted the CFSE^low CD4 and CD8 recipient T cells, allowing us to define clones by post-Tx MLR. Therefore, we can separately track the gene expression profiles of the pre-Tx and post-Txdefined HvG clones and assess their contribution to graft rejection or tolerance at the single cell level. We can further investigate the potential genes, pathways and networks involved with functional categories of HvG clones that corelate with immune tolerance and graft rejection after human ITx (107). This integrated platform could be applied to other types of SOT as well as hematopoietic cell transplantation to understand human alloresponses at a deeper level.

γδ T cells can recognize diverse antigens and exert disparate functions that include immune surveillance, modulation, and direct cytotoxicity (109). The innate- and adaptive-like features of human γδ T cells may be driven by differential γδ TCR repertoires, generally defined as innate-like Vγ9⁺δ2⁺ and “adaptive” non-Vγ9δ2, respectively (11). The repertoire can be shaped by tissue compartmentalization, age and history of antigen exposure (11, 110–112). Despite comprising a significant proportion of resident T cells in many organs, including gut and liver, γδ T cells and their possible role in transplantation outcomes are largely under‐researched (113).

Earlier functional studies were limited to in vitro systems of human peripheral blood-derived γδ T cells (114, 115), which are dominated by Vδ2 clonotypes. γδ T cells had been shown to exhibit either direct veto-type suppression of alloreactions (114), or indirect stimulation of alloreactive αβ T cell proliferation by inducing maturation of autologous dendritic cells and B cells into functional APCs (115). Tissue-resident Vδ1 T cells in liver explants have been investigated among patients who underwent liver transplantation, and were found to be polyfunctional and found to respond to both TCR and cytokine stimuli in vitro (116).

Our ongoing studies (117) suggest that γδ T cells may not only participate in host immune defense, but also likely have the potential to modulate two-way alloreactivities in the graft and hematopoietic system after human ITx. Our regimen of serial biopsies for ITx recipients allows us not only to investigate the role of intragraft γδ T cells in rejection but also to track the dynamics of their homeostatic turnover in the absence of rejection, providing a unique opportunity to study the fundamental biology of human γδ T cells. We have established a novel pipeline to integrate iRepertoire (γδ T cell primer sets) and 10x Genomics (5’cDNA library) platforms to relate functional gene profiles within individual γδ T cell clonotypes. TCR-seq of pre- and post-Tx γδ T cells in intestinal grafts, peripheral blood and BM are performed at the bulk and/or single cell levels. Clonal tracking of γδ T cells in intestinal grafts, peripheral blood and BM at both early and late time points post-Tx will provide a deeper understanding of their tissue origin, migration pattern and phenotypic maturation.

Our platform provides an opportunity to obtain a deeper understanding of the mechanisms behind γδ T cell chimerism, maturation of TRM features, and their modulatory roles on local and systemic alloresponses, facilitating the development of novel strategies to regulate γδ T cells to overcome rejection, infection and increase the utilization of ITx as a life-saving, quality of life-improving procedure. This approach can also be extended to other types of SOT and hematopoietic cell transplantation to provide immunogenomics insights into γδ T cells and their role in modulating alloresponses.