Next generation sequencing based assessment of the alloreactive T cell receptor repertoire in kidney transplant patients during rejection: a prospective cohort study

Clinical application of immune repertoire sequencing in solid organ transplant

Background

Measurement of T cell receptor (TCR) or B cell receptor (BCR) gene utilization may be valuable in monitoring the dynamic changes in donor-reactive clonal populations following transplantation and enabling adjustment in therapy to avoid the consequences of excess immune suppression or to prevent rejection with contingent graft damage and to indicate the development of tolerance.

Objective

We performed a review of current literature to examine research in immune repertoire sequencing in organ transplantation and to assess the feasibility of this technology for clinical application in immune monitoring.

Methods

We searched MEDLINE and PubMed Central for English-language studies published between 2010 and 2021 that examined T cell/B cell repertoire dynamics upon immune activation. Manual filtering of the search results was performed based on relevancy and predefined inclusion criteria. Data were extracted based on study and methodology characteristics.

Results

Our initial search yielded 1933 articles of which 37 met the inclusion criteria; 16 of these were kidney transplant studies (43%) and 21 were other or general transplantation studies (57%). The predominant method for repertoire characterization was sequencing the CDR3 region of the TCR β chain. Repertoires of transplant recipients were found to have decreased diversity in both rejectors and non-rejectors when compared to healthy controls. Rejectors and those with opportunistic infections were more likely to have clonal expansion in T or B cell populations. Mixed lymphocyte culture followed by TCR sequencing was used in 6 studies to define an alloreactive repertoire and in specialized transplant settings to track tolerance.

Conclusion

Methodological approaches to immune repertoire sequencing are becoming established and offer considerable potential as a novel clinical tool for pre- and post-transplant immune monitoring.

Into the multi-omics era: Progress of T cells profiling in the context of solid organ transplantation

T cells are the common type of lymphocyte to mediate allograft rejection, remaining long-term allograft survival impeditive. However, the heterogeneity of T cells, in terms of differentiation and activation status, the effector function, and highly diverse T cell receptors (TCRs) have thus precluded us from tracking these T cells and thereby comprehending their fate in recipients due to the limitations of traditional detection approaches. Recently, with the widespread development of single-cell techniques, the identification and characterization of T cells have been performed at single-cell resolution, which has contributed to a deeper comprehension of T cell heterogeneity by relevant detections in a single cell - such as gene expression, DNA methylation, chromatin accessibility, surface proteins, and TCR. Although these approaches can provide valuable insights into an individual cell independently, a comprehensive understanding can be obtained when applied joint analysis. Multi-omics techniques have been implemented in characterizing T cells in health and disease, including transplantation. This review focuses on the thesis, challenges, and advances in these technologies and highlights their application to the study of alloreactive T cells to improve the understanding of T cell heterogeneity in solid organ transplantation.

Perturbations of the T-cell immune repertoire in kidney transplant rejection

In this cross-sectional and longitudinal analysis of mapping the T-cell repertoire in kidney transplant recipients, we have investigated and validated T-cell clonality, immune repertoire chronology at rejection, and contemporaneous allograft biopsy quantitative tissue injury, to better understand the pathobiology of acute T-cell fraction, T-cell repertoire and antibody-mediated kidney transplant rejection. To follow the dynamic evolution of T-cell repertoire changes before and after engraftment and during biopsy-confirmed acute rejection, we sequenced 323 peripheral blood samples from 200 unique kidney transplant recipients, with (n=100) and without (n=100) biopsy-confirmed acute rejection. We report that patients who develop acute allograft rejection, have lower (p=0.01) T-cell fraction even before transplantation, followed by its rise after transplantation and at the time of acute rejection accompanied by high TCR repertoire turnover (p=0.004). Acute rejection episodes occurring after the first 6 months post-transplantation, and those with a component of antibody-mediated rejection, had the highest turnover; p=0.0016) of their T-cell repertoire. In conclusion, we validated that detecting repertoire changes in kidney transplantation correlates with post-transplant rejection episodes suggesting that T-cell receptor sequencing may provide recipient pre-transplant and post-transplant predictors of rejection risk.

High-Throughput Sequencing of Complementarity Determining Region 3 in the Heavy Chain of B-Cell Receptor in Renal Transplant Recipients: A Preliminary Report

Background: Graft failure resulting from rejection or any other adverse event usually originates from an aberrant and/or exaggerated immune response and is often catastrophic in renal transplantation. So, it is essential to monitor patients’ immune status for detecting a rejection/graft failure early on. Methods: We monitored the sequence change of complementary determining region 3 (CDR3) in B-cell receptor (BCR) immunoglobulin heavy-chain (IGH) immune repertoire (iR) in 14 renal transplant patients using next-generation sequencing (NGS), correlating its diversity to various clinical events occurring after transplantation. BCR-IGH-CDR3 in peripheral blood mononuclear cells was sequenced along the post-transplantation course by NGS using the iRweb server. Results: Datasets covering VDJ regions of BCR-IGH-CDR3 indicated clonal diversity (D50) variations along the post-transplant course. Furthermore, principal component analysis showed the clustering of these sequence variations. A total of 544 shared sequences were identified before transplantation. D50 remained low in three patients receiving rituximab. Among them, one’s D50 resumed after 3 m, indicating graft tolerance. The D50 rapidly increased after grafting and decreased thereafter in four patients without rejection, decreased in two patients with T-cell-mediated rejection (TCMR) and exhibited a sharp down-sliding after 3 m in two patients receiving donations after cardiac death (DCD). In another two patients with TCMR, D50 was low just before individual episodes, but either became persistently low or returned to a plateau, depending on the failure or success of the immunosuppressive treatments. Shared CDR3 clonal expansions correlated to D50 changes. Agglomerative hierarchical clustering showed a commonly shared CDR3 sequence and at least two different clusters in five patients. Conclusions: Clonal diversity in BCR-IGH-CDR3 varied depending on clinical courses of 14 renal transplant patients, including B-cell suppression therapy, TCMR, DCD, and graft tolerance. Adverse events on renal graft failure might lead to different clustering of BCR iR. However, these preliminary data need further verification in further studies for the possible applications of iR changes as genetic expression biomarkers or laboratory parameters to detect renal graft failure/rejection earlier.

Effects of Reduced-Dose Anti-Human T-Lymphocyte Globulin on Overall and Donor-Specific T-Cell Repertoire Reconstitution in Sensitized Kidney Transplant Recipients

Background

Pre-sensitized kidney transplant recipients have a higher risk for rejection following kidney transplantation and therefore receive lymphodepletional induction therapy with anti-human T-lymphocyte globulin (ATLG) whereas non-sensitized patients are induced in many centers with basiliximab. The time course of lymphocyte reconstitution with regard to the overall and donor-reactive T-cell receptor (TCR) specificity remains elusive.

Methods/Design

Five kidney transplant recipients receiving a 1.5-mg/kg ATLG induction therapy over 7 days and five patients with 2 × 20 mg basiliximab induction therapy were longitudinally monitored. Peripheral mononuclear cells were sampled pre-transplant and within 1, 3, and 12 months after transplantation, and their overall and donor-reactive TCRs were determined by next-generation sequencing of the TCR beta CDR3 region. Overall TCR repertoire diversity, turnover, and donor specificity were assessed at all timepoints.

Results

We observed an increase in the donor-reactive TCR repertoire after transplantation in patients, independent of lymphocyte counts or induction therapy. Donor-reactive CD4 T-cell frequency in the ATLG group increased from 1.14% + -0.63 to 2.03% + -1.09 and from 0.93% + -0.63 to 1.82% + -1.17 in the basiliximab group in the first month. Diversity measurements of the entire T-cell repertoire and repertoire turnover showed no statistical difference between the two induction therapies. The difference in mean clonality between groups was 0.03 and 0.07 pre-transplant in the CD4 and CD8 fractions, respectively, and was not different over time (CD4: F(1.45, 11.6) = 0.64 p = 0.496; CD8: F(3, 24) = 0.60 p = 0.620). The mean difference in R20, a metric for immune dominance, between groups was -0.006 in CD4 and 0.001 in CD8 T-cells and not statistically different between the groups and subsequent timepoints (CD4: F(3, 24) = 0.85 p = 0.479; CD8: F(1.19, 9.52) = 0.79 p = 0.418).

Conclusion

Reduced-dose ATLG induction therapy led to an initial lymphodepletion followed by an increase in the percentage of donor-reactive T-cells after transplantation similar to basiliximab induction therapy. Furthermore, reduced-dose ATLG did not change the overall TCR repertoire in terms of a narrowed or skewed TCR repertoire after immune reconstitution, comparable to non-depletional induction therapy.

ImmunoDataAnalyzer: a bioinformatics pipeline for processing barcoded and UMI tagged immunological NGS data

Background

Next-generation sequencing (NGS) is nowadays the most used high-throughput technology for DNA sequencing. Among others NGS enables the in-depth analysis of immune repertoires. Research in the field of T cell receptor (TCR) and immunoglobulin (IG) repertoires aids in understanding immunological diseases. A main objective is the analysis of the V(D)J recombination defining the structure and specificity of the immune repertoire. Accurate processing, evaluation and visualization of immune repertoire NGS data is important for better understanding immune responses and immunological behavior.

Results

ImmunoDataAnalyzer (IMDA) is a pipeline we have developed for automatizing the analysis of immunological NGS data. IMDA unites the functionality from carefully selected immune repertoire analysis software tools and covers the whole spectrum from initial quality control up to the comparison of multiple immune repertoires. It provides methods for automated pre-processing of barcoded and UMI tagged immune repertoire NGS data, facilitates the assembly of clonotypes and calculates key figures for describing the immune repertoire. These include commonly used clonality and diversity measures, as well as indicators for V(D)J gene segment usage and between sample similarity. IMDA reports all relevant information in a compact summary containing visualizations, calculations, and sample details, all of which serve for a more detailed overview. IMDA further generates an output file including key figures for all samples, designed to serve as input for machine learning frameworks to find models for differentiating between specific traits of samples.

Conclusions

IMDA constructs TCR and IG repertoire data from raw NGS reads and facilitates descriptive data analysis and comparison of immune repertoires. The IMDA workflow focus on quality control and ease of use for non-computer scientists. The provided output directly facilitates the interpretation of input data and includes information about clonality, diversity, clonotype overlap as well as similarity, and V(D)J gene segment usage. IMDA further supports the detection of sample swaps and cross-sample contamination that potentially occurred during sample preparation. In summary, IMDA reduces the effort usually required for immune repertoire data analysis by providing an automated workflow for processing raw NGS data into immune repertoires and subsequent analysis. The implementation is open-source and available on https://bioinformatics.fh-hagenberg.at/immunoanalyzer/.

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.

Prospective Tracking of Donor-Reactive T-Cell Clones in the Circulation and Rejecting Human Kidney Allografts

Background

Antigen recognition of allo-peptides and HLA molecules leads to the activation of donor-reactive T-cells following transplantation, potentially causing T-cell-mediated rejection (TCMR). Sequencing of the T-cell receptor (TCR) repertoire can be used to track the donor-reactive repertoire in blood and tissue of patients after kidney transplantation.

Methods/Design

In this prospective cohort study, 117 non-sensitized kidney transplant recipients with anti-CD25 induction were included. Peripheral mononuclear cells (PBMCs) were sampled pre-transplant and at the time of protocol or indication biopsies together with graft tissue. Next-generation sequencing (NGS) of the CDR3 region of the TCRbeta chain was performed after donor stimulation in mixed lymphocyte reactions to define the donor-reactive TCR repertoire. Blood and tissue of six patients experiencing a TCMR and six patients without rejection on protocol biopsies were interrogated for these TCRs. To elucidate common features of T-cell clonotypes, a network analysis of the TCR repertoires was performed.

Results

After transplantation, the frequency of circulating donor-reactive CD4 T-cells increased significantly from 0.86 ± 0.40% to 2.06 ± 0.40% of all CD4 cells (p < 0.001, mean dif.: -1.197, CI: -1.802, -0.593). The number of circulating donor-reactive CD4 clonotypes increased from 0.72 ± 0.33% to 1.89 ± 0.33% (p < 0.001, mean dif.: -1.168, CI: -1.724, -0.612). No difference in the percentage of donor-reactive T-cells in the circulation at transplant biopsy was found between subjects experiencing a TCMR and the control group [p = 0.64 (CD4⁺), p = 0.52 (CD8⁺)]. Graft-infiltrating T-cells showed an up to six-fold increase of donor-reactive T-cell clonotypes compared to the blood at the same time (3.7 vs. 0.6% and 2.4 vs. 1.5%), but the infiltrating TCR repertoire was not reflected by the composition of the circulating TCR repertoire despite some overlap. Network analysis showed a distinct segregation of the donor-reactive repertoire with higher modularity than the overall TCR repertoire in the blood. These findings indicate an unchoreographed process of diverse T-cell clones directed against numerous non-self antigens found in the allograft.

Conclusion

Donor-reactive T-cells are enriched in the kidney allograft during a TCMR episode, and dominant tissue clones are also found in the blood.

Trial Registration

Clinicaltrials.gov: NCT: 03422224 (https://clinicaltrials.gov/ct2/show/NCT03422224).

Integrated Analysis Toolset for Defining and Tracking Alloreactive T-cell Clones After Human Solid Organ and Hematopoietic Stem Cell Transplantation

We have developed a suite of tools for integrated analysis of T-Cell-Receptor Sequencing data to define and track alloreactive T-cells in human transplant studies. This has enabled discovery of sequences and patterns of T-cell enrichment and deletion associated with clinical outcomes such as transplant rejection and tolerance. The codebase includes user-friendly default analyses with customizable parameters which greatly accelerate computational workflows and provide robust statistics comparing post-transplant specimens to pre-transplant baseline. It also includes helper functions for robust characterization of T-cell-repertoire diversity, sample-to-sample divergence, resolution of sample-of-origin ambiguity in pooled assays, and functions to output all sequences defined as alloreactive.

Approaches for Controlling Antibody-Mediated Allograft Rejection Through Targeting B Cells

Both acute and chronic antibody-mediated allograft rejection (AMR), which are directly mediated by B cells, remain difficult to treat. Long-lived plasma cells (LLPCs) in bone marrow (BM) play a crucial role in the production of the antibodies that induce AMR. However, LLPCs survive through a T cell-independent mechanism and resist conventional immunosuppressive therapy. Desensitization therapy is therefore performed, although it is accompanied by severe side effects and the pathological condition may be at an irreversible stage when these antibodies, which induce AMR development, are detected in the serum. In other words, AMR control requires the development of a diagnostic method that predicts its onset before LLPC differentiation and enables therapeutic intervention and the establishment of humoral immune monitoring methods providing more detailed information, including individual differences in the susceptibility to immunosuppressive agents and the pathological conditions. In this study, we reviewed recent studies related to the direct or indirect involvement of immunocompetent cells in the differentiation of naïve-B cells into LLPCs, the limitations of conventional methods, and the possible development of novel control methods in the context of AMR. This information will significantly contribute to the development of clinical applications for AMR and improve the prognosis of patients who undergo organ transplantation.

Precision medicine in transplantation and hemodialysis

In kidney transplantation, precision medicine has already entered clinical practice. Donor and recipient human leucocyte antigen (HLA) regions are genotyped in two class 1 and usually three class 2 loci, and the individual degree of sensitization against alloimmune antigens is evaluated by the detection of anti-HLA donor-specific antibodies. Recently, the contribution of non-HLA mismatches to outcomes such as acute T- and B-cell-mediated rejection and even long-term graft survival was described. Tracking of specific alloimmune T- and B-cell clones by next generation sequencing and refinement of the immunogenicity of allo-epitopes specifically in the interaction with HLA and T- and B-cell receptors may further support individualized therapy. Although the choices of maintenance immunosuppression are rather limited, individualization can be accomplished by adjustment of dosing based on these risk predictors. Finally, supplementing histopathology by a transcriptomics analysis allows for a biological interpretation of the histological findings and avoids interobserver variability of results. In contrast to transplantation, the prescription of hemodialysis therapy is far from precise. Guidelines do not consider modifications by age, diet or many comorbid conditions. Patients with residual kidney function routinely receive the same treatment as those without. A major barrier hitherto is the definition of 'adequate' treatment based on urea removal. Kt/Vurea and related parameters neither reflect the severity of uremic symptoms nor predict long-term outcomes. Urea is poorly representative for numerous other compounds that accumulate in the body when the kidneys fail, yet clinicians prescribe treatment based on its measurement. Modern technology has provided the means to identify other solutes responsible for specific features of uremic illness and their measurement will be a necessary step in moving beyond the standardized prescription of hemodialysis.