The estimated glomerular filtration rate (eGFR) quantifies kidney graft function and is measured repeatedly after transplantation. Kidney graft rejection is diagnosed by performing biopsies on a regular basis (protocol biopsies at time of stable eGFR) or by performing biopsies due to clinical cause (indication biopsies at time of declining eGFR). The diagnostic value of the eGFR evolution as biomarker for rejection is not well established. To this end, we built a joint model which combines characteristics of transition models and shared parameter models to carry over information from one biopsy to the next, taking into account the longitudinal information of eGFR collected in between. From our model, applied to data of University Hospitals Leuven (870 transplantations, 2 635 biopsies), we conclude that a negative deviation from the mean eGFR slope increases the probability of rejection in indication biopsies, but that, on top of the biopsy history, there is little benefit in using the eGFR profile for diagnosing rejection. Methodologically, our model fills a gap in the biomarker literature by relating a frequently (repeatedly) measured continuous outcome with a less frequently (repeatedly) measured binary indicator. The developed joint transition model is flexible and applicable to multiple other research settings.
BillingsleyP (1986) Probability and measure. Chichester: Wiley.
2.
ChenH, HuangY and ZhangN (2015) Joint modeling of a linear mixed effects model for selfesteem from mean ages 13 to 22 and a generalized linear model for anxiety disorder at mean age 33. Journal of Medical Statistics and Informatics, 3, 1.
3.
CoxC (1998) Delta method. In Encyclopedia of Biostatistics, edited by P Armitage and T Colton, pages 1125–1126. New York, NY: John Wiley and Sons.
4.
de RooijM (2018) Transitional modeling of experimental longitudinal data with missing values. Advances in Data Analysis and Classification, 12, 107–30.
5.
DigglePJ, HeagertyP, LiangKY and ZegerS (2002) Analysis of longitudinal data. Oxford: Oxford University Press.
6.
DjamaliA, KaufmanDB, EllisTM, ZhongW, MatasA and SamaniegoM (2014) Diagnosis and management of antibody-mediated rejection: Current status and novel approaches. American Journal of Transplantation, 14, 255–71.
7.
EineckeG, ReeveJ, GuptaG, BoehmigGA, EskandaryF, BrombergJS, Budde, et al (2020) Factors associated with kidney graft survival in pure antibody-mediated rejection at the time of indication biopsy: Importance of parenchymal injury but not disease activity. American Journal of Transplantation, 21, 1391–1401.
8.
ElashoffRM, LiG and LiN (2008) A joint model for longitudinal measurements and survival data in the presence of multiple failure types. Biometrics, 64, 762–71.
9.
FassettRG, VenuthurupalliSK, GobeGC, CoombesJS, CooperMA and HoyWE (2011) Biomarkers in chronic kidney disease: A review. Kidney International, 80, 806–21.
10.
FerrerL, RondeauV, DignamJ, PicklesT, Jacqmin-GaddaH and Proust-LimaC (2016) Joint modelling of longitudinal and multistate processes: Application to clinical progressions in prostate cancer. Statistics in Medicine, 35, 3933–48.
11.
FitzmauriceG and MolenberghsG (2009) Advances in longitudinal data analysis: An historical perspective. In Longitudinal Data Analysis, edited by G Fitzmaurice, MG Davidian, G Verbeke and G Molenberghs, pages3–30. New York, NY: Chapman and Hall/CRC.
12.
GueorguievaR (2001) A multivariate generalized linear mixed model for joint modelling of clustered outcomes in the exponential family. Statistical Modelling, 1, 177–93.
13.
HalloranPF, MatasA, KasiskeBL, Madill-ThomsenKS, MackovaM and FamulskiKS (2019) Molecular phenotype of kidney transplant indication biopsies with inflammation in scarred areas. American Journal of Transplantation, 19, 1356–70.
14.
HariharanS, McbrideMA, CherikhWS, TollerisCB, BresnahanBA and JohnsonCP (2002) Post-transplant renal function in the first year predicts long-term kidney transplant survival. Kidney international, 62, 311–18.
15.
HollisE, ShehataM, KhalifaF, Abou El-GharM, El-DiastyT and El-BazA (2017) Towards non-invasive diagnostic techniques for early detection of acute renal transplant rejection: A review. The Egyptian Journal of Radiology and Nuclear Medicine, 48, 257–69.
16.
HorrocksJ and van Den HeuvelMJ (2009) Prediction of pregnancy: A joint model for longitudinal and binary data. Bayesian Analysis, 4, 523–38.
17.
IvanovaA, MolenberghsG and VerbekeG (2016) Mixed models approaches for joint modeling of different types of responses. Journal of Biopharmaceutical Statistics, 26, 601–18.
18.
JaffaMA, GebregziabherM and JaffaAA (2015) Analysis of multivariate longitudinal kidney function outcomes using generalized linear mixed models. Journal of Translational Medicine, 13, 1–12.
19.
JoostenSA, SijpkensYWJ, Van Kooten C and PaulLC (2005) Chronic renal allograft rejection: Pathophysiologic considerations. Kidney International, 68, 1–13.
20.
KasiskeBL, AndanyMA and DanielsonB (2002) A thirty percent chronic decline in inverse serum creatinine is an excellent predictor of late renal allograft failure. American Journal of Kidney Diseases, 39, 762–68.
21.
KimS, ZengD, ChamblessL and LiY (2012) Joint models of longitudinal data and recurrent events with informative terminal event. Statistics in Biosciences, 4, 262–81.
22.
LesaffreE and VerbekeG (1998) Local influence in linear mixed models. Biometrics, 54, 570–82.
23.
LiL, LuoS, HuB and GreeneT (2017) Dynamic prediction of renal failure using longitudinal biomarkers in a cohort study of chronic kidney disease. Statistics in Biosciences, 9, 357–78.
24.
LiangKY and ZegerSL (1989) A class of logistic regression models for multivariate binary time series. Journal of the American Statistical Association, 84, 447–51.
25.
LiuL and HuangX (2009) Joint analysis of correlated repeated measures and recurrent events processes in the presence of death, with application to a study on acquired immune deficiency syndrome. Journal of the Royal Statistical Society: Series C (Applied Statistics), 58, 65–81.
26.
LoDJ, KaplanB and KirkAD (2014) Biomarkers for kidney transplant rejection. Nature Reviews Nephrology, 10, 215.
27.
Lopez-GiacomanS and MaderoM (2015) Biomarkers in chronic kidney disease, from kidney function to kidney damage. World Journal of Nephrology, 4, 57–73.
28.
LoupyA, AubertO, OrandiBJ, NaesensM, BouatouY, RaynaudM, DivardG, et al (2019) Prediction system for risk of allograft loss in patients receiving kidney transplants: International derivation and validation study. BMJ, 366, 4923.
29.
LoupyA, HaasM, RoufosseC, NaesensM, AdamB, AfrouzianM, AkalinE, et al (2020) The Banff 2019 Kidney Meeting Report (I): Updates on and clarification of criteria for T cell -and antibody-mediated rejection. American Journal of Transplantation, 20, 2318–31.
30.
LoupyA, VernereyD, TinelC, AubertO, van HuyenJPD, RabantM, VerineJ, et al (2015) Subclinical rejection phenotypes at 1 year post-transplant and outcome of kidney allografts. Journal of the American Society of Nephrology, 26, 1721–31.
31.
MolenberghsG and VerbekeG (2006) Models for discrete longitudinal data. Springer Science and Business Media.
32.
PetersB, AnderssonY, StegmayrB, MölneJ, JensenG, DahlbergP, Holm-GunnarssonI, et al (2014) A study of clinical complications and risk factors in 1001 native and transplant kidney biopsies in Sweden. Acta Radiologica, 55, 890–96.
33.
RacusenLC (2006) Protocol transplant biopsies in kidney allografts: Why and when are they indicated?Clinical Journal of the American Society of Nephrology, 1, 144–47.
34.
RakhmawatiTW, MolenberghsG, VerbekeG and FaesC (2016) Local influence diagnostics for incomplete overdispersed longitudinal counts. Journal of Applied Statistics, 43, 1722–37.
35.
RakhmawatiTW, MolenberghsG, VerbekeG and FaesC (2016) Local influence diagnostics for hierarchical count data models with overdispersion and excess zeros. Biometrical Journal, 58, 1390–1408.
36.
RakhmawatiTW, MolenberghsG, VerbekeG and FaesC (2017) Local influence diagnostics for generalized linear mixed models with overdispersion. Journal of Applied Statistics, 44, 620–41.
37.
RizopoulosD (2011) Dynamic predictions and prospective accuracy in joint models for longitudinal and time-to-event data. Biometrics, 67, 819–29.
38.
RizopoulosD (2012) Joint models for longitudinal and time-to-event data: With applications in R. CRC Press.
39.
SalvadoriM and TsalouchosA (2017) Biomarkers in renal transplantation: An updated review. World Journal of Transplantation, 7, 161–78.
40.
SerroyenJ, MolenberghsG, VerbekeG and DavidianM (2009) Nonlinear models for longitudinal data. The American Statistician, 63, 378–88.
41.
Van LoonE, BernardsJ, Van CraenenbroeckAH and NaesensM (2020) The causes of kidney allograft failure: More than alloimmunity. A viewpoint article. Transplantation, 104, e46–e56.
42.
Van LoonE, GazutS, YazdaniS, LerutE, de LoorH, CoemansM, NoëlLH, et al (2019) Development and validation of a peripheral blood mRNA assay for the assessment of antibody-mediated kidney allograft rejection: A multicentre, prospective study. EBioMedicine, 46, 463–72.
43.
VerbekeG, FieuwsS, MolenberghsG and DavidianM (2014) The analysis of multivariate longitudinal data: A review. Statistical Methods in Medical Research, 23, 42–59.
44.
VigliettiD, LoupyA, AubertO, BestardO, Van HuyenJPD, TaupinJL, GlotzD, et al (2018) Dynamic prognostic score to predict kidney allograft survival in patients with antibody-mediated rejection. Journal of the American Society of Nephrology, 29, 606–19.
45.
WilkinsonA (2006) Protocol transplant biopsies: Are they really needed?Clinical Journal of the American Society of Nephrology, 1, 130–37.
