Journal Logo

Kidney Transplantation

Donor-derived Cell-free DNA Combined With Histology Improves Prediction of Estimated Glomerular Filtration Rate Over Time in Kidney Transplant Recipients Compared With Histology Alone

Huang, Edmund MD1; Gillespie, Matthew PharmD1; Ammerman, Noriko PharmD1; Vo, Ashley PharmD1; Lim, Kathlyn PharmD1; Peng, Alice MD1; Najjar, Reiad MD1; Sethi, Supreet MD1; Jordan, Stanley C. MD1; Mirocha, James MS2; Haas, Mark MD, PhD3

Author Information
doi: 10.1097/TXD.0000000000001027

Abstract

INTRODUCTION

The Banff classification of rejection in kidney transplantation relies on histologic assessment interpreted in the presence or absence of donor-specific antibodies (DSAs).1 Recognizing that rejection can be present in a functionally stable allograft, measures of allograft function, such as creatinine, glomerular filtration rate (GFR), or proteinuria, are not included in the Banff criteria for rejection. Given this, it is not surprising that allografts with similar histologic appearances can have divergent prognoses and that integration of functional and immunologic parameters with histology predicts outcomes better than histology alone.2 Prior studies have questioned the specificity of the Banff classification system and argued that histologic patterns of rejection might reflect nonimmune mediated injury in some circumstances.3 The Banff working group has considered associations with graft outcome when assessing the significance of key histologic features of rejection, and it could be argued that lesions associated with long-term stability in graft function might not be representative of immune-mediated injury.4 Therefore, the diagnostic classification of rejection might be enhanced by including determinants of allograft function and prognosis, which are best assessed using prognostic models.

Donor-derived cell-free DNA (dd-cfDNA, Allosure; CareDx, Inc.), a biomarker that can detect allograft injury, recently gained Medicare approval in 2017 for the assessment of rejection in adult kidney transplant recipients. Validation studies have indicated that higher percentages of dd-cfDNA in the blood correlate with the presence of rejection identified on biopsy,5,6 with recent literature suggesting that dd-cfDNA might have some prognostic capability.7 In this study, we investigated whether dd-cfDNA can predict the subsequent trajectory of estimated GFR (eGFR) and whether the combination of dd-cfDNA with histology can better inform on allograft prognosis than histology alone.

MATERIALS AND METHODS

This study was approved by the Cedars-Sinai Medical Center Institutional Review Board (Pro00020945). The study design, data analysis, and writing of the manuscript was performed by the authors alone and was not supported by any outside entity. Beginning in August, 2017, a total of 180 unique kidney transplant patients who were at least 1 mo posttransplant and followed at Cedars-Sinai Medical Center had assessment of dd-cfDNA within 1 mo of kidney transplant biopsy and were included in the study. All patients had dd-cfDNA assessment per manufacturer guidelines, with the following considered exclusion criteria: (1) age <18 y old, (2) pregnant status, (3) recipient of a transplant from an identical twin, (4) recipient of a nonkidney solid organ transplant, and (5) prior recipient of an allogeneic bone marrow transplant. In our center, dd-cfDNA and biopsies are assessed “for-cause” to investigate a suspicion of rejection in the setting of allograft dysfunction or presence of DSA. Biopsies were graded according to Banff 2019 criteria, which establish a minimum threshold for interstitial inflammation involving 10%–25% of unscarred cortical parenchyma (i1 lesion) in the presence of tubulitis as criteria for borderline cell-mediated rejection (CMR).8 Because there was no difference in dd-cfDNA values between borderline CMR and CMR cases, histology meeting criteria for borderline CMR was classified as CMR for this study.

All serum creatinine values were recorded beginning with the most recent assessment before dd-cfDNA measurement until graft loss (n = 5), death (n = 2), or the end of the follow-up period on November 14, 2019. eGFR was calculated using the Chronic Kidney Disease Epidemiology Collaboration creatinine equation.9 Linear mixed effects models using random slopes and intercepts and an unstructured covariance matrix were then fitted to derive a prediction model for eGFR over time. The fixed effect represents the average rate of change in eGFR by treatment group, whereas the random effect accounts for subject-specific correlation between repeated measures of eGFR within an individual. Univariate predictors and their interactions with time were tested in the model and included: recipient characteristics (age at dd-cfDNA assessment, gender, donor type [deceased versus living], presence of DSA), Banff biopsy scores (cg, ci + ct, g + ptc, inflammation in area of interstitial fibrosis and tubular atrophy [i-IFTA], i + t, and v-scores), and dd-cfDNA values (dichotomized into upper 2 quartiles versus lower 2 quartiles).

Five-fold crossvalidation was performed to identify the best-fitting model. For this procedure, the study population was randomly split into 5 equal-sized groups, and the model selection procedure was applied separately on 5 training sets incorporating a unique combination of 4 of the 5 folds. The best-fitting model identified from the training sets was then fitted on each remaining holdout (validation) set not utilized in the model training procedure. Covariates with a covariate-by-time interaction P < 0.05 on univariate analysis were entered into a multivariable linear mixed effects model. Nested models were compared using the likelihood-ratio (LR) test, Akaike Information Criterion (AIC), and Bayesian Information Criterion (BIC) to assess if inclusion of dd-cfDNA into a parsimonious model consisting of Banff biopsy scores would improve model fit. The LR test is a chi-squared test that assesses whether a model maximizes the likelihood function over another, whereas the AIC and BIC are methods that balance goodness of fit and the complexity of the model. Both the AIC and BIC favor the simplest model that yields the best fit, with a lower AIC and BIC indicating better goodness-of-fit.

A sensitivity analysis was performed on a subgroup of the study population excluding ABO-incompatible recipients and recipients with BK nephropathy, given that biopsies from both types of recipients have histologic features that overlap with but are not necessarily indicative of rejection (C4d-positivity in ABO-incompatible recipients and tubulointerstitial inflammation in BK nephropathy).

All P values were 2-tailed and a P < 0.05 was considered statistically significant. All analyses were performed using Stata version 14 (StataCorp; College Station, TX).

RESULTS

Table 1 describes the baseline characteristics of the study population. Assessment of dd-cfDNA was performed at a median of 659 d posttransplant (interquartile range [IQR], 81-2113). The majority of patients included were recipients of a deceased donor kidney. Fifty-seven patients (32% of the study population) had DSA at the time of dd-cfDNA assessment, with the vast majority consisting of class II DSA (88%). There was balanced representation of CMR, antibody-mediated rejection (AMR), and mixed CMR/AMR diagnoses.

TABLE 1.
TABLE 1.:
Baseline characteristics

The median dd-cfDNA in the study population was 0.70% (IQR, 0.28%–1.6%). Figure 1 compares the distribution of dd-cfDNA values corresponding to 4 clinical diagnostic categories: (a) no rejection, (b) isolated CMR, (c) isolated AMR, and (d) mixed CMR/AMR. There were significant differences in the distribution of dd-cfDNA between the 4 groups (P < 0.001). Pairwise comparisons indicated that the distribution of dd-cfDNA was higher among patients with isolated CMR (median, 0.80%; IQR, 0.33%–2.5%; P < 0.001), isolated AMR (median, 1.4%; IQR, 1.1%–2.7%; P < 0.001), and CMR/AMR (median, 1.5%; IQR, 0.84%–2.6%; P < 0.001) compared with those with no rejection (median, 0.35%; IQR, 0.19%–0.60%). The distribution of dd-cfDNA was higher among patients with isolated AMR (P = 0.01) and CMR/AMR (P = 0.04) compared with isolated CMR. There was no difference in the distribution of dd-cfDNA between patients with isolated AMR and those with mixed CMR/AMR (P = 0.58). With regard to cases of isolated CMR, there was no difference in dd-cfDNA between cases of borderline CMR (n = 21; median, 0.80%; IQR, 0.35%–1.3%) and cases of isolated CMR 1A or above (n = 24; median, 0.76; IQR, 0.29%–4.4%; P = 0.39).

FIGURE 1.
FIGURE 1.:
Box and whiskers plot comparing the distribution of donor-derived cell-free DNA % across histologic diagnoses. ABMR, antibody-mediated rejection; CMR, cell-mediated rejection.

Table 2 shows the distribution of Banff histologic scores corresponding to the 4 clinical diagnostic categories of no rejection, isolated CMR, isolated AMR, and mixed CMR/AMR. There was 1 patient classified as having no rejection despite having an isolated v-lesion (v1). This patient presented with a transiently increased creatinine and a dd-cfDNA value of 1.6%. There were no other histologic features of CMR or AMR identified on biopsy. On subsequent testing, the patient’s creatinine spontaneously returned to baseline and clinical suspicion of rejection was low. The patient was not treated for rejection and was thus assigned a clinical diagnosis of no rejection. Ten patients had a cg-score ≥1a, were C4d-negative, and did not have current or prior DSA, thus not meeting criteria for chronic AMR. Additionally, there were 9 recipients of an ABO-incompatible allograft found to have C4d ≥ 2 without concomitant glomerulitis, peritubular capillaritis, or DSA and were classified as not having AMR. As expected, there was a preponderance of tubulointerstitial inflammation observed among patients with CMR, whereas microvascular inflammation (glomerulitis and peritubular capillaritis) and glomerular changes (transplant glomerulopathy) were more commonly seen among patients with AMR.

TABLE 2.
TABLE 2.:
Distribution of Banff histology scores by clinical diagnosis

Figure 2 shows the eGFR trajectory over time corresponding to each individual study participant. Linear mixed effects models were fitted to model the association of recipient characteristics, Banff histologic scores, and dd-cfDNA on eGFR over time. Model discovery was performed using 5-fold cross-validation and observed that the best-fitting model consisted of a combination of the Banff biopsy scores of cg (cg = 3 versus cg < 3) and ci + ct (≥3 versus <3) with dd-cfDNA (upper 2 quartiles versus lower 2 quartiles). This model yielded a mean root mean square error of 7.6 mL/min/1.73 m2 (SE: 0.3 mL/min/1.73 m2) across the 5 validation sets. Model output applied to the entire dataset is shown in Table 3, and measures of model fit are shown in Table 4. Table 3a shows the output for univariate linear mixed effects models. In the table, the reference coefficient represents the slope in eGFR per month (Δ eGFR/mo), whereas the coefficient for the interaction represents the difference in eGFR slope between the reference and comparator groups. Δ eGFR/mo for the comparator group equals the sum of the coefficients for the reference group and the interaction. There was no association between recipient or donor characteristics and Δ eGFR/mo on univariate analysis. The covariate-by-time interactions for cg and ci + ct were significant, indicating that biopsies with higher chronicity scores were associated with a steeper decline in eGFR over time compared with biopsies with milder degrees of chronicity. There was a trend toward significance for the interaction between dd-cfDNA and time (coefficient –0.5; 95% CI, −1.0–0.05; P = 0.08), providing a signal that higher percentages of dd-cfDNA may be associated with eGFR decline. Covariate-by-time interactions for i + t and i-IFTA were also borderline significant. Table 3b shows the multivariable model output for cg and ci + ct together with their associated covariate-by-time interactions, again indicating that biopsies with higher chronicity scores were associated with a steeper decline in eGFR over time compared with biopsies with less chronicity. Table 3c shows the regression output for Banff chronicity characteristics with dd-cfDNA values together. Inferences from this model paralleled those observed from Tables 3a and b.

TABLE 3.
TABLE 3.:
(a) Univariate and (b) multivariate associations of Banff biopsy scores and (c) Banff biopsy scores and donor-derived cell-free DNA on eGFR overtime
TABLE 4.
TABLE 4.:
Comparison of model fit between nested models. A lower AIC and BIC indicate better model fit
FIGURE 2.
FIGURE 2.:
Individual eGFR trajectories over time (n = 180). eGFR, estimated glomerular filtration rate.

Table 4 shows measures of model fit for the nested multivariable models integrating Banff chronicity scores (cg and ci + ct) with dd-cfDNA and their corresponding covariate-by-time interactions versus cg and ci + ct with covariate-by-time interactions. The LR test was significant (P < 0.001), indicating that the more complex model including biopsy scores with dd-cfDNA provided a better model fit than the more parsimonious model including biopsy scores alone. This observation was supported by both penalized-likelihood criteria (AIC and BIC) favoring the more complex model of biopsy scores with dd-cfDNA. Including the presence of DSA to this model and its interaction with time did not improve fit further (LR test statistic = 3.2; df = 2; P = 0.20). Additionally, inclusion of acute inflammation (i, t, and v), microvascular inflammation (g and ptc), and i-IFTA scores to Banff chronicity scores (cg and ci + ct) did not improve model fit (data not shown).

Sensitivity Analysis

There were 9 recipients of an ABO-incompatible transplant with C4d ≥2 who did not have concomitant glomerulitis, peritubular capillaritis, or DSA and were classified as not having AMR. Additionally, there were 5 recipients who had features of BK nephropathy on biopsy, of whom 3 were recipients of an ABO-incompatible transplant. A sensitivity analysis was performed to see if exclusion of these patients resulted in any appreciable difference from what was observed in the overall population (n = 169). The univariate and multivariable regression output is shown in Table S1 SDC, http://links.lww.com/TXD/A261 and comparisons of nested models shown in Table S2, SDC, http://links.lww.com/TXD/A261. The analyses excluding ABO-incompatible recipients and recipients with BK nephropathy observed similar findings to the main analysis.

DISCUSSION

In this study, we found that both Banff biopsy characteristics of chronic transplant glomerulopathy and chronic tubulointerstitial fibrosis were predictive of eGFR over time, and the combination of dd-cfDNA with biopsy characteristics improved model prediction of eGFR over time compared with biopsy characteristics alone. These observations suggest a role for the use of biomarkers, such as dd-cfDNA, in conjunction with kidney transplant biopsies to better inform on prognosis in kidney transplant recipients.

The Banff diagnostic classification has historically relied on clinicopathologic correlations, validating key lesions as characteristic of rejection if associated with disease progression. The concept of C4d-negative AMR was introduced on the basis of studies demonstrating that expression of endothelial-associated and NK cell transcripts was associated with AMR lesions10,11 and, at least for endothelial-associated transcripts, was a better predictor of graft loss than C4d.10 The classification of C4d-negative AMR was ultimately validated and incorporated into the Banff 2013 diagnostic classification from data observing that patients with microvascular inflammation found on protocol biopsy in the absence of C4d were more likely to progress to more advanced interstitial fibrosis and tubular atrophy and lower GFR than patients without microvascular inflammation or C4d.4,12 More recently, the Banff working group has recommended revising the diagnostic classification for borderline CMR to establish a minimum threshold of interstitial inflammation involving 10%–25% of unscarred cortical parenchyma (i1 lesion) in the presence of tubulitis (at minimum t1 lesion).8 This recommendation was based on 2 reports showing that isolated tubulitis without interstitial inflammation was associated with similar graft survival compared with biopsies without rejection.13,14

Additionally, the Banff working group formally recognized molecular markers in the form of validated transcripts or gene classifiers to substitute for DSA for the diagnosis of AMR in the 2017 diagnostic classification.1 This recognition was an acknowledgment that molecular diagnostics, if thoroughly validated, could be used in conjunction with histology to classify rejection. Recent data also provide support for consideration of clinical characteristics together with histology to more precisely predict outcomes. An unsupervised archetype analysis from Aubert et al2 demonstrated that integration of functional, immunologic, and histologic characteristics could characterize distinct graft survival outcomes that could not be distinguished by histology alone. Furthermore, a risk prediction score derived from a multivariable Cox regression model incorporating demographic characteristics, DSA specificities and concentrations, and histologic scores validated in an international cohort showed good concordance with observed graft survival (C-statistic, 0.80–0.81).15

Our study provides additional support for the paradigm of a multimodal assessment of graft prognosis by demonstrating that assessment of dd-cfDNA together with histologic lesions in a predictive model can better predict eGFR over time compared with histology alone. This suggests that histologic findings might be better interpreted in the context of dd-cfDNA, and consideration of the 2 together may be more informative for making treatment decisions. This concept is supported by recent data indicating that dd-cfDNA can be used to risk stratify patients treated for CMR 1A and borderline CMR lesions on biopsy.7 In that study, patients with a dd-cfDNA <0.5% exhibited graft stability over 3–6 mo, whereas those with a dd-cfDNA ≥0.5% were more likely to develop DSA and exhibit deterioration in eGFR. Although one cannot conclude from these data that patients with low levels of dd-cfDNA need not be treated for rejection, emerging data are providing support that biomarkers, including dd-cfDNA, should be considered when interpreting biopsies for rejection.

We acknowledge that findings from our study should be validated in larger, multicenter datasets. Baseline characteristics, immunologic risk, and treatment protocols may differ between transplant centers, which may potentially have an impact on eGFR response. Additionally, this was a retrospective study that relied on clinical decisions and not prespecified protocols on dd-cfDNA and biopsy assessment. As such, it is not known whether exclusion of dd-cfDNA tests assessed without a biopsy may have biased our findings. It should also be noted that this study used cross-sectional data and does not account for changes in dd-cfDNA or biopsy findings over time. Despite these limitations, our study lends strength to the paradigm that the diagnostic classification of rejection may be enhanced when biopsies are considered alongside biomarkers, such as dd-cfDNA.

In summary, the combination of Banff cg and ci + ct scores together with dd-cfDNA provided a better model fit for eGFR over time compared with Banff scores alone. Although these findings are preliminary, our study argues that a larger-scale external validation study should be conducted. If validated, assessment of biopsies alongside dd-cfDNA could be considered as a newer paradigm for classifying rejection in kidney transplantation.

REFERENCES

1. Haas M, Loupy A, Lefaucheur C, et al. The Banff 2017 Kidney Meeting Report: revised diagnostic criteria for chronic active T cell-mediated rejection, antibody-mediated rejection, and prospects for integrative endpoints for next-generation clinical trials. Am J Transplant. 2018; 18:293–307. doi:10.1111/ajt.14625
2. Aubert O, Higgins S, Bouatou Y, et al. Archetype analysis identifies distinct profiles in renal transplant recipients with transplant glomerulopathy associated with allograft survival. J Am Soc Nephrol. 2019; 30:625–639. doi:10.1681/ASN.2018070777
3. Madill-Thomsen K, Perkowska-Ptasinska A, Bohmig GA, et al. Discrepancy analysis comparing molecular and histology diagnoses in kidney transplant biopsies. Am J Transplant. [Epub ahead of print. 2019]. doi:10.1111/ajt.15752
4. Haas M, Sis B, Racusen LC, et al. Banff meeting report writing committee. Banff 2013 Meeting Report: inclusion of c4d-negative antibody-mediated rejection and antibody-associated arterial lesions. Am J Transplant. 2014; 14:272–283. doi:10.1111/ajt.12590
5. Bloom RD, Bromberg JS, Poggio ED, et al.; Circulating Donor-Derived Cell-Free DNA in Blood for Diagnosing Active Rejection in Kidney Transplant Recipients (DART) Study Investigators Circulating Donor-Derived Cell-Free DNA in Blood for Diagnosing Active Rejection in Kidney Transplant Recipients (DART) Study Investigators. Cell-Free DNA and active rejection in kidney allografts. J Am Soc Nephrol. 2017; 28:2221–2232. doi:10.1681/ASN.2016091034
6. Huang E, Sethi S, Peng A, et al. Early clinical experience using donor-derived cell-free DNA to detect rejection in kidney transplant recipients. Am J Transplant. 2019; 19:1663–1670. doi:10.1111/ajt.15289
7. Stites E, Kumar D, Olaitan O, et al. High levels of dd-cfDNA identifies patients with TCMR 1A and borderline allograft rejection at elevated risk of graft injury. Am J Transplant. 2020. doi:10.1111/ajt.15822
8. Loupy A, Haas M, Roufosse C,. The Banff 2019 Kidney Meeting Report (I): updates on and clarification of criteria for T cell- and antibody-mediated rejection. Am J Transplant. 2020. doi:10.1111/ajt.15898
9. Levey AS, Stevens LA, Schmid CH, et al.; CKD-EPI (Chronic Kidney Disease Epidemiology Collaboration) CKD-EPI (Chronic Kidney Disease Epidemiology Collaboration). A new equation to estimate glomerular filtration rate. Ann Intern Med. 2009; 150:604–612
10. Sis B, Jhangri GS, Bunnag S, et al. Endothelial gene expression in kidney transplants with alloantibody indicates antibody-mediated damage despite lack of C4d staining. Am J Transplant. 2009; 9:2312–2323. doi:10.1111/j.1600-6143.2009.02761.x
11. Hidalgo LG, Sis B, Sellares J, et al. NK cell transcripts and NK cells in kidney biopsies from patients with donor-specific antibodies: evidence for NK cell involvement in antibody-mediated rejection. Am J Transplant. 2010; 10:1812–1822. doi:10.1111/j.1600-6143.2010.03201.x
12. Loupy A, Suberbielle-Boissel C, Hill GS, et al. Outcome of subclinical antibody-mediated rejection in kidney transplant recipients with preformed donor-specific antibodies. Am J Transplant. 2009; 9:2561–2570. doi:10.1111/j.1600-6143.2009.02813.x
13. McRae M, Bouchard-Boivin F, Béland S, et al. Impact of the current versus the previous diagnostic threshold on the outcome of patients with borderline changes suspicious for T cell-mediated rejection diagnosed on indication biopsies. Transplantation. 2018; 102:2120–2125. doi:10.1097/TP.0000000000002327
14. Nankivell BJ, P’Ng CH, Chapman JR, et al. Does tubulitis without interstitial inflammation represent borderline acute T cell mediated rejection? Am J Transplant. 2019; 19:132–144. doi:10.1111/ajt.14888
15. Loupy A, Aubert O, Orandi BJ, et al. Prediction system for risk of allograft loss in patients receiving kidney transplants: international derivation and validation study. BMJ. 2019; 366:l4923. doi:10.1136/bmj.l4923

Supplemental Digital Content

Copyright © 2020 The Author(s). Transplantation Direct. Published by Wolters Kluwer Health, Inc.