Abstract
Background and objectives T cell mediated rejection (TCMR) is the most frequent type of acute rejection and is associated with kidney allograft failure. Almost 40% of TCMR episodes fail to respond to anti-rejection therapy. FOXP3 is a specification factor for regulatory T cells and our single center study of 83 kidney allograft recipients showed that urinary cell FOXP3 mRNA level is diagnostic of TCMR and predicts TCMR reversibility and allograft survival. The objective of the current study is to determine whether our original findings could be replicated in an independent cohort of 480 kidney allograft recipients enrolled in the multicenter Clinical Trials of Organ Transplantation (CTOT)-04.
Design, setting, participants, and measurements We measured levels of FOXP3 mRNA and levels of mRNA for CD25, CD3E, and perforin, and 18S rRNA in 3559 urines from 480 kidney allograft recipients prospectively enrolled in CTOT-04. RNA was isolated from the urinary cells and preamplification-enhanced real-time quantitative PCR assays were used to measure mRNAs.
Results 18S rRNA normalized levels of mRNA for FOXP3 (P=0.01, Kruskal-Wallis test), CD25 (P=0.01), CD3E (P<0.0001), and perforin (P<0.0001) distinguished patients with TCMR biopsies from those with No Rejection biopsies and those with stable graft function. FOXP3 mRNA level, but not the levels of mRNA for CD25, CD3E, or perforin, predicted TCMR reversal (AUC=0.764; 95% confidence interval, 0.611 to 0.917; P=0.008). Multivariable logistic regression analysis showed that FOXP3 mRNA level remains predictive after adjustment for potential cofounders. Kaplan-Meier survival curve analysis showed that FOXP3 mRNA level (P=0.0306) but not the levels of mRNA for CD25, CD3E, or perforin, is associated with kidney allograft survival.
Conclusion In the independent CTOT-04 cohort, we demonstrate that urinary cell level of FOXP3 mRNA is diagnostic of TCMR, predicts its reversibility, and is prognostic of kidney allograft survival following an episode of TCMR.
Introduction
T cell mediated rejection (TCMR) is the most frequent type of acute rejection (1–4). Anti-rejection therapy has evolved over the years, but Bouatou et al found in a comprehensive study of 256 kidney allograft recipients with kidney allograft biopsies showing TCMR that 40% of TCMR episodes fail to respond to anti-rejection treatment (5). Moreover, the histologic hallmarks of TCMR – Banff tubulitis score and the interstitial inflammation score at the time of TCMR biopsy diagnosis– were not determinants of kidney allograft loss whereas the inflammation score and peritubular capillaritis score observed in the biopsy performed 3 months after anti-rejection treatment were independent predictors of kidney allograft failure (5). The need for post-treatment parameters to improve prognostic accuracy is problematic from the perspective of clinical decision making at the time of TCMR diagnosis. Also, the need for multiple biopsies – one to diagnose TCMR and a second one to better prognosticate TCMR outcome – is challenging in view of complications related to the invasive biopsy procedure; the well-documented inter-observer variability in the grading of biopsies is yet another problem (6, 7). Because molecular or cellular mechanisms for recalcitrant TCMR have not been deciphered, a mechanistic approach to treatment of TCMR is currently lacking, and often the “one-size-fits all” strategy is followed in the clinical setting.
To fill the existing gaps and improve our understanding of mechanisms for non-responsiveness to anti-rejection treatment, we investigated whether urinary cell FOXP3 mRNA profiles predict TCMR reversal. Our approach was informed by the following considerations: (i) our earlier single center study of 83 kidney allograft recipients demonstrating that urinary cell level of FOXP3 mRNA predicts TCMR reversibility and identifies patients at risk for kidney allograft failure (8); (ii) naturally occurring CD4+FOXP3+ regulatory T cells (Tregs) play a non-redundant role in immune homeostasis and in the maintenance of self-tolerance (9, 10); (iii) preclinical data showing that Tregs prevent rejection and promote transplant tolerance (11–13); and (iv) consideration of adoptive Treg therapy in the clinic to control autoimmunity or promote tolerance to allografts (14, 15). We also considered that it is important to investigate the reproducibility of the associations we previously reported for urinary cell FOXP3 mRNA levels (8) in view of the reproducibility crisis in the replication of scientific findings (16, 17).
In the current investigation, we measured levels of mRNAs in 3505 urine specimens from 480 kidney transplant recipients enrolled in the multicenter Clinical Trials in Organ Transplantation 04 (CTOT-04) study in order to determine whether urinary cell FOXP3 gene expression profile is diagnostic of TCMR, predicts TCMR reversibility, and is prognostic of kidney allograft survival following an episode of TCMR.
Materials and Methods
Study subjects
In the parent CTOT-04 study, 485 kidney allograft recipients were prospectively enrolled at 5 transplant sites. Urine was prospectively collected on days 3, 7, 15, and 30 and in months 2, 3, 4, 5, 6, 9 and 12 as well as the time of each kidney allograft biopsy and 2 weeks thereafter. Urine cell pellets were prepared at each clinical site and shipped to the Gene Expression Monitoring (GEM) core at Weill Cornell Medicine. The GEM core isolated RNA from the urinary cell pellets, reverse transcribed it to cDNA, and checked for transcript quality thresholds – 100 copies of TGFB1 mRNA and 5×107 copies of 18S rRNA per one microgram of total RNA – prior to downstream data analysis. Absolute levels of mRNA for CD3E, perforin, granzyme B, proteinase inhibitor-9 (PI-9), CD103, interferon inducible protein-10 (IP-10), CXCR3, TGFB1, and 18S rRNA (CTOT-04 Prespecified mRNA Panel) were measured. A total of 3559 urine specimens from 485 kidney allograft recipients passed the transcript quality thresholds. The primary objectives of the parent CTOT-04 study were to determine whether urinary cell levels of a prespecified mRNA panel, measured at time of biopsy, are diagnostic of TCMR and whether the levels in sequential samples predict future development of TCMR. The parent CTOT-04 study did not investigate whether the measured transcripts predict TCMR reversal or kidney allograft survival. Urinary cell levels of FOXP3 mRNA and mRNA for CD25 were not measured in the parent CTOT-04 study.
In the current study, cDNA prepared from the same 3559 urine cell pellets were retrieved for the measurement of FOXP3 mRNA and CD25 mRNA. Prior to measurement of mRNAs, the cDNA was again assessed for transcript quality thresholds and 3505 of the 3559 cDNAs (98.5%) prepared from the urine specimens from 480 of 485 kidney allograft recipients from the parent CTOT-04 study cohort met the quality threshold, and were used to measure urinary cell mRNA levels using preamplification-enhanced real-time quantitative polymerase chain reaction (customized RT-QPCR) assays developed in our laboratory (4). Table S1 lists the oligonucleotide primers and TaqMan probes we designed for the absolute quantification of mRNAs. Additional details of the customized RT-QPCR assay developed in our laboratory and used to measure transcript levels have been reported (4).
In the current investigation, we applied the same criterion we used in our earlier study (8) to classify an episode of TCMR as reversible or nonreversible. Specifically, a TCMR episode was classified as reversible if the serum creatinine level returned to within 15 percent of the pre-rejection level within 4 weeks after initiation of anti-rejection therapy.
Statistical methods
Copy numbers for each mRNA were analyzed after normalization with 18S rRNA copies (x10−6) and log10-transformation. Kruskal-Wallis and Mann-Whitney statistical tests were used to compare 18S rRNA normalized mRNA levels across diagnoses. Receiver-operating-characteristic (ROC) curves were used to determine the predictive accuracy of urinary cell mRNA levels and the optimal sensitivity and specificity combination was determined based on Youden’s index (18). Multivariable logistic regression was used to evaluate urinary cell levels of mRNA after controlling for potential confounders.
Kaplan-Meier survival curves were used to evaluate graft survival rates and to determine the ability of urinary cell mRNA level, measured at the time of biopsy, to anticipate graft failure. Time to event was calculated from time of TCMR biopsy until graft failure or censoring. Subjects were censored if they died prior to experiencing graft failure or were lost to follow-up. Log-rank tests were used to compare survival curves across strata. Long-term follow-up information on graft outcomes, beyond the 3-year duration of the CTOT-04 study, was obtained from the Organ Procurement and Transplantation Network (OPTN). Multivariable Cox proportional hazards regression was used to evaluate urinary cell mRNA as a predictor of graft failure after controlling for potential confounders.
All analyses were performed with SAS software, version 9.4 (SAS Institute) except for the creation of ROC curves using the pROC R package (19). The datasets and code generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.
Study approval
The institutional review board at Weill Cornell Medicine approved the study, “Multicenter Study: Use of PCR to Evaluate Immune Regulatory Molecules”, Protocol Number 9608002317.
Results
Study groups
Figure 1 shows the distribution of 3505 urine samples collected from the 480 patients and with urinary cell mRNA data. Among the CTOT-04 cohort of 480 kidney allograft recipients, 218 underwent kidney allograft biopsies, which were classified by the on-site pathologist using the Banff classification schema (20). Among 262 patients who did not undergo a biopsy, 199 patients contributed 1524 urine samples and met the criteria for stable graft function based on average study recorded serum creatinine values ≤ 2.0 mg/dl, no treatment for acute rejection, and no evidence of cytomegalovirus or polyomavirus type BK infection during the first 12 months post-transplantation (Stable).
Urinary cell mRNAs diagnostic of TCMR
We compared urinary cell levels of mRNA in 43 urine samples matched to TCMR biopsies from 34 patients, 162 samples matched to No Rejection biopsies from 126 patients, and 1524 samples collected longitudinally from 199 Stable patients. Table 1 is a summary of baseline characteristics of the kidney allograft recipients included in this analysis.
Violin plots with in-laid box-and-whisker plots in Figure 2 illustrate the distribution of log10-transformed ratios of mRNA to 18S rRNA copies (x10−6) per microgram of total RNA. We replicated our earlier findings that FOXP3 mRNA levels and the levels of mRNA for CD25, CD3E, and perforin are all diagnostic of TCMR; the levels of all 4 mRNAs were significantly higher in urine samples matched to TCMR biopsies than in urine samples matched to No Rejection biopsies. The levels of these mRNAs in urine matched to TCMR biopsies were also higher than in urine samples collected longitudinally from Stable patients. Table S2 shows the median and lower and upper quartiles of the log10-transformed 18S rRNA normalized ratios of mRNAs for the diagnostic groups of TCMR, No Rejection, Stable and other biopsy diagnoses.
Urinary cell mRNA levels stratified by reversibility of TCMR
Among 43 biopsy-confirmed TCMR episodes from 34 patients with biopsy-matched urine samples, 39 episodes from 33 patients were classifiable as reversible (n=24) or nonreversible (n=15); the remaining 4 were not classified due to missing creatinine values (n=2), close proximity to BKV nephropathy diagnosis (n=1), or close proximity to an earlier biopsy (n=1). Tables S2 and S3 list recipient and donor characteristics and biopsy-associated variables, respectively, by TCMR reversible status.
We compared urinary cell levels of mRNAs in urine matched to reversible TCMR biopsies to those matched to nonreversible TCMR biopsies (Table 2). Urinary cell levels of mRNA for FOXP3 (P=0.0096, Mann-Whitney Test), but not the levels of mRNA for CD25 (P=0.1531), CD3E (P=0.1887), or perforin (P=0.4322), were significantly higher in urine matched to reversible TCMR biopsies than in samples matched to nonreversible TCMR biopsies.
We also compared levels of mRNAs in urine matched to reversible or nonreversible TCMR biopsies to levels in urine matched to No Rejection biopsies (Table 2). These comparisons identified additional differential gene expression patterns. Urinary cell levels of mRNA for Treg specification factor FOXP3 (P<0.0001, Mann-Whitney Test) and CD25 (P=0.0066) were significantly higher in urine matched to reversible TCMR biopsies than in urine matched to No Rejection biopsies. In contrast, levels of FOXP3 mRNA (P=0.7212) and CD25 mRNA (P=0.7623) were not different between urine matched to nonreversible TCMR biopsies and urine matched to No Rejection biopsies. The levels of mRNA encoding cytotoxic perforin and pan T cell specific CD3E were significantly higher both in urine matched to reversible and nonreversible TCMR biopsies compared to urine matched to No Rejection biopsies (Table 2).
ROC curve analysis of reversal of TCMR
ROC curve analysis yielded an AUC of 0.764 (95% confidence interval [CI], 0.611 to 0.917; P=0.008) for the 18S rRNA normalized values of FOXP3 mRNA (Figure 3A). The cut-point that maximized Youden’s index (18) was 0.06, yielding a sensitivity of 92% (95% CI, 73% to 99%) and specificity of 87% (95% CI, 60% to 98%). Multivariable logistic regression analyses showed that FOXP3 mRNA level remains predictive of TCMR reversal after individual adjustment for serum creatinine measured at time of biopsy, time from transplant to biopsy, and type of organ donor (Table 3).
As previously observed in our single center study (8), levels of mRNA for CD25 (AUC=0.642; 95% CI, 0.463 to 0.820; P=0.100) (Figure 3B), CD3E (AUC=0.631; 95% CI, 0.440 to 0.821; P=0.173) (Figure 3C), and perforin (AUC=0.578; 95% CI, 0.383 to 0.772; P=0.399) (Figure 3D) did not predict TCMR reversal in the multicenter CTOT-04 cohort.
The urinary cell three-gene signature of log10-transformed 18S rRNA normalized values of CD3E and IP10 mRNAs and log10-transformed values of 18S rRNA ([x10−6]) was developed and validated in the parent CTOT-04 study to be diagnostic of TCMR (4). The three-gene TCMR diagnostic signature did not discriminate reversible from nonreversible TCMR episodes (AUC=0.589; 95% CI, 0.396 to 0.781; P=0.253) (Figure 3E). Serum creatinine measured at the time of biopsy did not predict TCMR reversal (AUC=0.611; 95% CI, 0.428 to 0.795; P=0.308) (Figure 3F). The use of thymoglobulin as anti-rejection treatment for TCMR (AUC=0.479; 95% CI, 0.321 to 0.637; P=0.792) (data not shown) did not discriminate between reversible and nonreversible TCMR.
Survival analyses of kidney allografts
We examined whether the functional criterion used in this study to classify an episode of TCMR as reversible or nonreversible is associated with kidney allograft survival. Kaplan-Meier survival analysis showed that survival rates of recipients with reversible TCMR are similar to those with No Rejection biopsies, whereas almost all graft losses were in those with nonreversible TCMR biopsies (Figure 4A, P<0.0001, Log-rank test).
Kidney allograft recipients with FOXP3 mRNA levels above the median FOXP3 mRNA level of -1.251 had significantly better graft survival compared to those with levels below the median (Figure 4B, P=0.0306). Multivariable Cox proportionate hazards regression analysis showed that FOXP3 mRNA level remains significantly predictive of kidney allograft outcomes after adjustment for serum creatinine measured at time of biopsy, but not after adjustment for time from transplant to biopsy or type of transplant (Table 4).
Urinary cell levels of mRNA for CD25 (P=0.4371), CD3E (P=0.8062), and perforin (P=0.9839), measured at the time of TCMR biopsy and stratified by their median levels, were not associated with kidney allograft survival. Neither the use of thymoglobulin for anti-rejection therapy (Figure 4C, P=0.8678) nor the level of serum creatinine measured at time of biopsy was associated with kidney allograft survival (Figure 4D, P=0.7384).
Discussion
FOXP3+ Tregs play a pivotal role in preventing autoimmunity and maintaining immune homeostasis (9, 10). Preclinical studies suggest that Tregs prevent or delay the onset of allograft rejection and may induce tolerance (11–13). In our earlier single center clinical study, we found that urinary cell FOXP3 mRNA level is diagnostic of TCMR and predicts TCMR reversal; we also found that urinary cell FOXP3 mRNA level is associated with kidney allograft failure (8). In the current investigation, we have replicated these findings using urine samples collected from an external cohort of kidney allograft recipients enrolled in the multicenter CTOT-04 study. Our replication of earlier findings is significant from a biological perspective regarding the potential role of Tregs in regulating kidney allograft rejection, and is reassuring in the context of the existing reproducibility crisis in replicating scientific findings (16, 17).
Our profiling identified molecular features consistent with the hypothesis that the balance between T effectors and Tregs is tipped towards Treg deficiency at the onset of nonreversible TCMR. Whereas the levels of mRNA for cytopathic molecule perforin, pan T cell marker CD3E, and T cell activation marker CD25 were not significantly different between urine matched to reversible TCMR biopsies and urine matched to nonreversible TCMR biopsies, level of mRNA for Treg specification factor FOXP3 was significantly lower in urine matched to nonreversible TCMR biopsies compared to urine matched to reversible TCMR biopsies.
Recombinant IL-2 or Tregs have been used to prevent graft versus host disease (GVHD) in recipient hematopoietic stem cell recipients and to prevent acute rejection in solid organ graft rejection (14, 15, 21). The findings from the current study suggest that infusion of IL-2 or Tregs may be a consideration for those with TCMR associated with a low urinary cell FOXP3 mRNA. A switch from calcineurin inhibitor (CNI)-based immunosuppression to rapamycin-based immunosuppression could also be considered in this setting given that CNIs are potent inhibitors of IL-2 gene transcription, whereas rapamycin can reprogram Treg metabolism and improve its function (22–24).
The overall survival of kidney allografts in the CTOT-04 study cohort was similar to the US kidney graft survival rates (25) suggesting that our study participants are representative of the US kidney transplant population. Kidney graft survival was significantly inferior in patients with biopsy-confirmed TCMR, with most graft failures occurring in those with nonreversible TCMR biopsies. Importantly, urinary cell FOXP3 mRNA level, measured at time of TCMR, foretold allograft failure. Laboratory parameters such as serum creatinine, measured at the time of allograft biopsy, or the use T cell depleting agent antithymocyte globulin as anti-rejection therapy did not predict kidney allograft failure in this investigation.
A number of features of our study are worthy of emphasis. First, urinary cell FOXP3 level was assessed in what is likely to be the largest number of urine samples from the largest number of kidney allograft recipients from 5 academic transplant sites. Second, urine samples were collected under the direction of an independent Statistical and Clinical Coordinating Center ensuring protocol-based sample collection, and data entry and integrity (4). Third, the oligonucleotide primers we designed for the quantification of mRNAs in the customized RT-QPCR assays obviated DNA amplification and ensured that only mRNAs are measured in the assays. Fourth, our incorporation of preamplification prior to mRNA measurement compensated for the low RNA yield from urinary cells and facilitated successful measurement of multiple mRNAs from the same cDNA sample. Fifth, the use of a customized BAK amplicon in the RT-QPCR assays allowed for absolute quantification of transcript numbers and avoided pitfalls associated with relative quantification (26). Finally, our study included a clinically relevant endpoint of graft failure following TCMR.
Our study has a number of shortcomings. We characterized TCMR reversal based on functional recovery rather than by histological confirmation. This may not be a significant limitation since graft outcome in this study was strongly associated with the functional criterion used to classify an episode of TCMR. Another limitation is that we did not assess the functional activity of Tregs and we inferred Treg deficiency based on FOXP3 mRNA abundance. We note that mRNA expression patterns in themselves have helped inform therapeutic decisions (27). Our profiling of urinary cells rather than the allograft itself might be considered a limitation as well. This is less of a concern in the light of RNA sequencing demonstrating that kidney allograft biopsy transcriptional signatures are enriched in urinary cells (28).
In sum, we have replicated our single center findings that urinary cell FOXP3 mRNA level is diagnostic of TCMR, predicts TCMR reversal and is associated with graft survival following an episode of TCMR. It has been reported that 40% of TCMR episodes fail to respond to state-of-the-art anti-rejection therapy (5). Our finding that recalcitrant TCMR is associated with low FOXP3 abundance suggest that infusion of IL-2, a trophic factor for Tregs, infusion of Tregs or rapamycin that may reprogram Treg metabolism may represent novel treatment options for TCMR exemplified by a urinary cell FOXP3 mRNA level consistent with Treg deficiency.
Data Availability
The data that support the findings of this study are available on request from the corresponding author. The data are not publicly available due to privacy or ethical restrictions.
Funding
This investigation was supported by awards from the National Institute of Allergy and Infectious Diseases, National Institutes of Health (RO1 AI072790 and R37 AI051652 to MS). The Clinical Trials in Organ Transplantation Study-04 (ClinicalTrials.gov NCT00337220) was supported by an award (UO1AI63589) from the National Institute of Allergy and Infectious Diseases, National Institutes of Health to Abraham Shaked, University of Pennsylvania School of Medicine, Philadelphia, PA.
Author contributions
MS conceived and coordinated the study. DL, PA, JES, and MS wrote the manuscript. DD made critical contributions to the preparation of the manuscript including data collection and kidney allograft survival analysis. RD designed the oligonucleotide primers and TaqMan probes, and developed and supervised the performance of preamplification-enhanced real-time quantitative PCR assays. TM, ML, JL, and FM made significant contributions to data analysis. VKS coordinated data acquisition and management. DL performed the statistical analyses. JES provided critical statistical guidance to DL.
Competing interests
The authors of this manuscript have no conflicts of interest to disclose.
Data and materials availability
The data that support the findings of this study are available on request from the corresponding author. The data are not publicly available due to privacy or ethical restrictions.
Supplementary Materials
Materials and Methods
Acknowledgments
The authors gratefully acknowledge Dr. Nancy D. Bridges, Transplantation Branch, National Institute of Allergy and Infectious Diseases, Bethesda, MD, USA for her careful review of the manuscript and key insights. The authors thank Ms. Christina Chang and Ms. Christine Hoang (Weill Cornell Medicine, New York, NY) for their superb technical assistance in performing the RT-QPCR assays. The authors thank the United Network for Organ Sharing for providing the clinical data for the calculation of survival curves of kidney allograft recipients. The data reported here have been supplied by UNOS as the contractor for the Organ Procurement and Transplantation Network (OPTN). The OPTN data system includes data on all donors, wait-listed candidates, and transplant recipients in the US, submitted by all members of OPTN. The Health Resources and Services Administration (HRSA), U.S. Department of Health and Human Services provides oversight to the activities of the OPTN contractor. The interpretation and reporting of these data are the responsibility of the authors and in no way should be seen as an official policy of or interpretation by the OPTN or the U.S. Government.