Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Risk factors for graft loss and death among kidney transplant recipients: A competing risk analysis

  • Jessica Pinto-Ramirez ,

    Roles Conceptualization, Data curation, Formal analysis, Methodology, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    jpinto@colombianadetrasplantes.com

    Affiliation Department of Transplant Nephrology, Colombiana de Trasplantes, Bogotá, Colombia

  • Andrea Garcia-Lopez,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Validation, Writing – original draft, Writing – review & editing

    Affiliation Department of Transplant Research, Colombiana de Trasplantes, Bogotá, Colombia

  • Sergio Salcedo-Herrera,

    Roles Conceptualization, Data curation, Methodology, Supervision

    Affiliation Department of Transplant Nephrology, Colombiana de Trasplantes, Bogotá, Colombia

  • Nasly Patino-Jaramillo,

    Roles Data curation, Formal analysis, Project administration, Writing – original draft, Writing – review & editing

    Affiliation Department of Transplant Research, Colombiana de Trasplantes, Bogotá, Colombia

  • Juan Garcia-Lopez,

    Roles Data curation, Formal analysis, Methodology, Writing – original draft

    Affiliation Departmento of Technology and Informatics, Colombiana de Trasplantes, Bogotá, Colombia

  • Jefferson Barbosa-Salinas,

    Roles Data curation, Formal analysis, Methodology

    Affiliation Departmento of Technology and Informatics, Colombiana de Trasplantes, Bogotá, Colombia

  • Sergio Riveros-Enriquez,

    Roles Data curation, Formal analysis, Methodology

    Affiliation Departmento of Technology and Informatics, Colombiana de Trasplantes, Bogotá, Colombia

  • Gilma Hernandez-Herrera,

    Roles Conceptualization, Formal analysis, Investigation, Methodology, Supervision, Validation

    Affiliation Postgraduate Program in Epidemiology, Universidad del Rosario – Universidad CES, Bogotá-Medellín, Colombia

  • Fernando Giron-Luque

    Roles Conceptualization, Investigation, Project administration, Supervision, Validation

    Affiliation Department of Transplant Surgery, Colombiana de Trasplantes, Bogotá, Colombia

Abstract

Introduction

Kidney transplantation is the best therapeutical option for CKD patients. Graft loss risk factors are usually estimated with the cox method. Competing risk analysis could be useful to determine the impact of different events affecting graft survival, the occurrence of an outcome of interest can be precluded by another. We aimed to determine the risk factors for graft loss in the presence of mortality as a competing event.

Methods

A retrospective cohort of 1454 kidney transplant recipients who were transplanted between July 1, 2008, to May 31, 2019, in Colombiana de Trasplantes, were analyzed to determine risk factors of graft loss and mortality at 5 years post-transplantation. Kidney and patient survival probabilities were estimated by the competing risk analysis. The Fine and Gray method was used to fit a multivariable model for each outcome. Three variable selection methods were compared, and the bootstrapping technique was used for internal validation as split method for resample. The performance of the final model was assessed calculating the prediction error, brier score, c-index and calibration plot.

Results

Graft loss occurred in 169 patients (11.6%) and death in 137 (9.4%). Cumulative incidence for graft loss and death was 15.8% and 13.8% respectively. In a multivariable analysis, we found that BKV nephropathy, serum creatinine and increased number of renal biopsies were significant risk factors for graft loss. On the other hand, recipient age, acute cellular rejection, CMV disease were risk factors for death, and recipients with living donor had better survival compared to deceased-donor transplant and coronary stent. The c-index were 0.6 and 0.72 for graft loss and death model respectively.

Conclusion

We developed two prediction models for graft loss and death 5 years post-transplantation by a unique transplant program in Colombia. Using a competing risk multivariable analysis, we were able to identify 3 significant risk factors for graft loss and 5 significant risk factors for death. This contributes to have a better understanding of risk factors for graft loss in a Latin-American population. The predictive performance of the models was mild.

Introduction

Kidney transplantation is the optimal renal replacement therapy for suitable patients with end stage renal failure [1]. Identifying risk factors for graft failure in kidney transplant recipients is useful for recognizing those patients at high risk and anticipating potential therapeutic interventions to improve graft survival [2, 3]. Risk factors in the field of organ transplantation are typically assessed using time-to-event outcomes, for instance, when recording time-to-death or time-to-graft loss. Survival analyses are key statistical tools in transplantation research [4]. This analyses are the most used methods to estimate the incidence of an outcome of interest, often censoring for a competing event [4]. For example, death is competing event to graft loss because a patient may die before losing the graft, as such no opportunity for graft loss in that case. Thus, competing events are present when another event precludes the event of interest. In this condition, the Kaplan-Meier (KM) approach is not suitable because this method assumes that censored patients are at the same risk as patients who remain in the study. In general, this leads to an overestimation of the cumulative incidence of the event of interest [57]. To solve this limitation, Kalbfleisch and Prentice introduced the Cumulative Incidence Function (CIF) [8]. The CIF calculates all events’ probability as the sum of the event of interest’s probabilities and those of the competing risks. The competing risk analysis (CRA) allows using a modified chi-squared test to compare CIF curves between groups and the Fine and Gray model with subdistribution hazards (sdHR) [9, 10]. Thus, patients are followed until observing the first of multiple event types in the CRA. Adjusting this fact, the inferred incidence of the event of interest is lower, and the sum of calculated incidences across all event types sums up to 100% [4, 10].

The CRA may provide further insights into the effect of interventions on the separate endpoints, comparing CIF curves to explore the association between covariates and the absolute risk. Indeed, CRA may be essential for clinical decision-making and prognostic research questions [11]. Despite this advantage, competing risk models have not been used frequently by researchers [2, 10]. In particular, the advantage of using CRA method was highlighted in a study evaluating race, age, and survival among patients undergoing dialysis, where accounting for transplant as competing risk brought to light a greater disparity in death on dialysis among younger black patients (related to disparity in access to transplantation) [12]. Other studies used the CRA method to evaluate risk of mortality and subsequent graft survival in older recipients after sustaining fracture [13], as well as the risk of graft loss and mortality in older recipients (age ≥60 y) receiving older kidneys (age ≥80 y) versus remaining on dialysis [14]. So far, there are no Latin-American studies related to competing risks analysis in transplantation. In this context and given the substantial variability of the identified risk factors for graft loss across different transplant populations, our transplant program (Colombiana de Trasplantes—CT) aims to use a well characterized Latin-American cohort of kidney transplant recipients with long term follow-up to determine the risk factors for graft loss in the presence of mortality as a competing risk event.

Materials and methods

Study design and population

We conducted a retrospective cohort study at Colombiana de Trasplantes. To give a context it must be said that Colombiana de Trasplantes is a transplant network in Colombia with 4 centers with around 21% of the annual national kidney transplant activity. We included first time, kidney transplant recipients aged ≥ 18 years who were transplanted between July 1, 2008, to May 31, 2019. Patients with primary renal graft thrombosis (arterial or venous) were excluded. Recipients were followed up to graft failure, death, or end of follow up at 5 years post transplantation, whichever was earliest.

Kidney donors

Informed consent was obtained from organ/tissue deceased donors (DD) by a family interview where both family and the donation team go through information related to answer inquiries about encephalic death, emotional support, and the possibility of organ donation. The main causes of death in our DD were cerebrovascular/stroke, followed by head trauma, anoxia, CNS tumor, or others. Organ donation and tissue consent form is provided in S1 Table in S1 File.

In the case of the live kidney donors, our transplant team provides kidney donation and nephrectomy informed consent and the affidavit from the live kidney donors. Overall, less than 1% of our kidney donors were previously registered as organ donors. According to Colombian law, the total of donor medical costs is covered including organ donor maintenance and procurement with an average of 5000 USD and without any economic contribution to the family donor.

Immunosuppression and follow-up protocol

All patients received standard induction therapy with alemtuzumab, basiliximab or antihuman thymocyte immunoglobulin according to immunologic risk or transplant clinical guidelines. All patients received a fixed-dose of methylprednisolone perioperatively for 3 days with a transition to fixed-dose oral prednisone from day 4 to day 7 in the postoperative period. One-week post transplantation steroids were withdrawn. Chronic immunosuppression consists of Calcineurin inhibitors-based therapy and mycophenolate mofetil. Patients are monitored closely in the first 4 weeks post transplantation, and they return for follow-up monthly thereafter.

The acute rejection was classified under parameters described by Banff (2015) [15]. Biopsy was performed on those patients with increase of serum creatinine by >20% from baseline. Our center does not perform biopsies per protocol.

Treatment for acute cellular rejection was started once the histological diagnosis was confirmed as follows: Methylprednisolone: 500 mg. IV / day in infusion for three days. Oral prednisolone from the fourth day at a dose of 0.5 mg / Kg / day divided into two doses and for two weeks. After completing the two weeks, a weekly decrease of 10 mg / day was made until reaching the previous dose that was received before the rejection episode. Serum creatinine was done 5 days after finishing the boluses. A response to corticosteroid treatment was defined with a decrease in serum creatinine greater than and equal to 20% of the patient’s baseline creatinine.

Histocompatibility tests performed in our center correspond to Human Leukocyte Antigens (HLA), Panel Reactive Antibodies (PRA) I and II, flow cytometric crossmatch and anti-HLA antibodies (the latter only in living donors when the crossmatch is positive). HLA matching is when the recipient and the donor shared the same HLA antigens (HLA-A,-B,-DR antigens) [16].

We do not perform routine preimplant biopsies. The decision to take or not the organs from expanded criteria donors is made by macroscopic evaluation of the graft and, if required, it is sent for histological evaluation. Indices like KDPI / KDRI are not considered for taking or allocating organs [17]. Organ allocation is made according to the allocation criteria for kidney transplantation in Colombia [18].

Outcomes

Our primary outcome of interest was graft loss, not including death with function. Graft loss was defined by center reported as permanent return to dialysis or retransplantation. Death was defined as mortality from any cause and was ascertained by review of the Colombiana de Trasplantes database which records patient’s death and supplemented with the National register Master File. Patients were censored at 5 years of follow-up since the last follow-up date if they were transferred to another transplant center or lost to follow-up. Thus, survival analyses were performed using a competing risk approach, where graft loss and mortality were treated as competing events.

Statistical analysis

Descriptive analysis.

Descriptive statistics were used to report the population characteristics. Frequencies and percentages were used for categorical variables. The numerical variables were reported according to its distribution using mean and standard deviation for normally distributed variables, and median and inter-quartile range (IQR) for non-normally distributed variables. Multiple imputation was not considered as there were few missing values (5.9% of the total number of cases), and those values were at random. According to this, we performed a complete case analysis in the univariable and multivariable models.

Predictors.

Prespecified variables based on published literature and those available in our data, were collected as potential risk factors for graft loss. Data collected included demographics, medical history and clinic characteristics of kidney transplant recipient and donor. Definition of predictor measurement is provided in Supplementary material (S1 Table in S1 File).

Incidence estimates.

The overall incidence of graft loss and/or death at 5 years post transplantation was calculated by Competing risk analysis method (CRA) using cumulative incidence function (CIF) where mortality was treated as a competing risk with graft loss. Log Rank test for graft loss and death were compared in the entire population and in specific patient population including living and deceased donor.

Comparisons between the two groups (graft loss yes/no and death yes/no) were performed using modified χ2 test. The subdistribution Hazard Ratio (sdHR) also known as Fine and Gray model was calculated for each independent variable and the two outcomes.

Variable selection and prediction.

Variables with p value <0.25 in an univariable analysis and those with clinical importance were selected to perform further analysis. Variable selection to build the final model for graft loss was performed comparing three methods:

  1. 1. Full model: contains all the predictors selected in previous analysis and no variable selection was done.
  2. 2. Backward selection based on the Akaike information criterion (AIC).
  3. 3. Backward selection based on the Bayesian Information Criterion (BIC).

The model was selected on the model´s better performance.

Bootstrapping technique was used for internal validation as split method for resample a single dataset to create many simulated samples. The prediction models were trained on B bootstrap samples with replacement. The models were assessed in the observations that were not included in the bootstrap sample. This allowed us to calculate the prediction error, brier score and c-index. Calibration plot was used to compare the predicted probability with the observed probability.

The Fine and Gray model directly models the covariate effect on CIF, and it reports the sdHR. To model the impact of covariates on graft loss, we used the Fine and Gray method [9] for performing competing risk regression. The association between the primary outcome and the independent variables were assessed by the sdHR.

The model development and report was based on The Transparent Reporting of a multivariable prediction model for Individual Prognosis Or Diagnosis (TRIPOD) [19]. More details of modeling process can be found in Supplementary material (Modelling process in S1 File).

Analysis was performed using the Software R version 3.6.3. Library used to perform competing risk analysis was cmprsk [20].

Ethics considerations

This study was approved by the ethics research committee of the institution, acting in concordance with local and national regulations, as well as with the Helsinki declaration. Confidentiality of all patients was secured all the time during the execution of the research. None of the transplant donors was from a vulnerable population and all donors or next of kin provided written informed consent that was freely given.

Results

Patient characteristics

A total of 1454 out of 1621 recipients met selection criteria. Exclusions took place in 167 (113 pediatric transplants and 54 kidney transplants with graft thrombosis). In gender distribution most of patients were male, the overall mean age was 43.58 ± 13 years. Table 1 summarizes the clinical and demographic characteristics.

thumbnail
Table 1. Clinical and demographic characteristics of kidney transplant recipients.

https://doi.org/10.1371/journal.pone.0269990.t001

Overall cumulative incidence

During the follow-up period graft loss occurred in 169 patients (11.6%) and death occurred in 137 (9.4%). Cumulative incidence for graft loss and death was 15.8% and 13.8% respectively. Fig 1 displays the combined cumulative incidence for the entire cohort. Significant differences in estimates of both outcomes were found when analyzing live and deceased transplant separately, where deceased transplant (17.1% and 16.3% for death and graft loss respectively) had greater cumulative incidences (deceased transplant 17.1% and 16.3% for death and graft loss respectively vs live 5.4% and 15% for death and graft loss respectively). Fig 2 shows the difference between type of transplant in the cumulative incidence of graft loss and death.

thumbnail
Fig 1. Cumulative incidence of graft loss and death estimated by the method for competing risk.

https://doi.org/10.1371/journal.pone.0269990.g001

thumbnail
Fig 2. Log Rank of cumulative incidence of risk of death by type of transplant, and Log Rank of cumulative incidence of risk graft loss by type of transplant.

https://doi.org/10.1371/journal.pone.0269990.g002

Risk factors for cumulative incidence of graft loss and death with functioning graft

We fit the Fine and Gray competing risk survival regression model for identifying the potential determinants of graft loss using covariates with significant association and those with clinical importance. The covariates that had a significant impact on the graft loss were stroke, cold ischemia time, qualitative PRA II, BKV nephropathy, number of allograft biopsies, acute cellular rejection, serum creatinine at 12 months and number of hospital readmissions. The covariates that had a significant impact on death were recipient age, diabetes mellitus, type of donor, expanded criteria donor, CMV disease, cold ischemia time, coronary stent, and number of hospital readmissions.

Table 2 provides the sdHR of risk factors estimated by using the final multivariate Fine and Gray model. The risk of graft failure was noted to increase in the presence of nephropathy due BK virus, higher rates of serum creatinine at 12 months post transplantation, and greater number of kidney biopsies. Significant risk factors associated with cumulative incidence of death were recipient age, deceased donor, CMV disease, coronary stent, and acute cellular rejection.

thumbnail
Table 2. Factors associated with graft loss using death as a competing risk in kidney transplant recipients in a final Fine and Gray model.

https://doi.org/10.1371/journal.pone.0269990.t002

Variable selection to build the final model for graft loss was performed comparing three methods:

  1. 4. Full model: contains all the predictors selected in previous analysis and no variable selection was done. 2. Backward selection based on the AIC. 3. Backward selection based on the BIC.

The performance and prediction error of the three models were evaluated using Bootstrap cross-validation, showing similar results for the AIC and BIC models. The C-index for the full model was 0.57, for the AIC model was 0.6 and, for the BIC model was 0.6.

The same process for variable selection and performance assessment was performed for death model. Similar results were obtained for the AIC and BIC models. The C-index for the full model was 0.78, for the AIC model was 0.72 and, for the BIC model was 0.72. Fig 3 provides the prediction errors and calibration plot of the final Fine and Gray model for graft loss. Fig 4 provides the prediction errors and calibration plot of the final Fine and Gray model for death.

thumbnail
Fig 3. Prediction errors and calibration plot of the final Fine and Gray model for graft loss.

https://doi.org/10.1371/journal.pone.0269990.g003

thumbnail
Fig 4. Prediction errors and calibration plot of the final Fine and Gray model for death.

https://doi.org/10.1371/journal.pone.0269990.g004

Discussion

Kidney transplantation is the best therapy available for most patients with end- stage kidney disease [21]. We developed two predictive models of risk of graft loss and risk of death in kidney transplant patients. Risk prediction models are useful for identifying kidney recipients at high risk of graft failure, and optimize clinical care, decision-making and resource allocation; is a challenging issue in kidney transplantation [2]. Our objective was characterized Latin-American cohort of kidney transplant recipients with long term follow-up and to predict the risk factors for graft loss in the presence of mortality as a competing risk event. We were able to identify 3 significant risk factors for graft loss and 5 significant risk factors for death. This contributes to have a better understanding of risk factors for graft loss in a Latin-American population.

Graft survival is one of the most critical concerns in kidney transplant recipients, and our ability to accurately monitor the cumulative incidence of graft loss its importance. Risk prediction models are useful for identifying kidney recipients at high risk of graft failure, thus optimizing clinical care. Therefore, using competing risks methods that provide more accurate estimates, we sought to identify risk factors leading to graft loss, considering death as a competing risk in kidney transplant recipients [4].

Particularly, one study highlights the advantage of using CRA method assessing both the probabilities of death and graft loss.

Risk factors for graft loss

Late failure of kidney transplants remains an important clinical problem [2, 22]. Renal allograft loss is multifactorial [23]. In the United States, 5469 kidney transplants developed end-stage kidney failure in 2008 (data provided by Jon Snyder from USRDS), making kidney transplant failure the fourth leading cause of end-stage renal disease. The reasons for failure are not well understood. Some have postulated that late deterioration reflects dysregulated fibrosis, drug toxicity or progressive “chronic allograft nephropathy” [2426]. In our study we found that BKV nephropathy, serum creatinine at 12 months and increased number of renal allograft biopsies were significant risk factors for graft loss. Sellarés et al., attributed causes of graft failure in the biopsy-for-cause population to antibody-mediated rejection (ABMR), probable ABMR or mixed rejection, with nonadherence recorded in nearly half. There was evidence of ABMR in 18 of 19 nonadherent patients who failed. There were also three groups of nonrejection causes of failures: glomerulonephritis, BKV nephropathy and failure in the context of an intercurrent illness. The results emphasize the burden of ABMR and mixed rejection and its interaction with nonadherence in observed failures, making these key targets for further progress. The results also illustrate the range of clinical courses leading to failure and the sometimes-complex relationships to the indication biopsy findings [22].

On the other hand, renal allograft biopsy (RAB) is still the best approach to diagnose renal transplant complications [27]. We found that kidney recipients with more significant requirements to perform RAB had greater risk of graft loss. According to our guidelines, biopsy was performed on those patients with increase of serum creatinine by >20% from baseline, generally when acute or chronic renal allograft rejection is suspected, antibody-mediated rejection, polyoma virus nephropathy, glomerular diseases, atrophy u other. Thus, greater number of RAB may be associated with renal allograft disfunction or detection of other lesions that may influence graft loss [2729]. However, as a limitation, we do not have an electronic database of all the histological findings of dysregulated fibrosis and C4d of the renal biopsies of our patients. Of those factors related to graft loss, there is a high magnitude of association with BKV nephropathy [30, 31]. Previous studies have associated the BKV nephropathy with premature loss of kidney function [3235], graft loss and alteration of renal histology [2123, 36]. The reactivation of the virus may occur with the use of immunosuppression. Polyomavirus BK virus reactivation in kidney transplant recipients can lead to BK polyoma virus-associated nephropathy (BKPyVAN), which is associated with graft dysfunction in >90% and graft loss in over 50% of the affected individuals [37].

Our results also showed an association with higher serum creatinine level at 12 months. This factor has been widely described as predictor of graft loss [3, 3845].

The identification of risk factors for graft loss in the long term has been provided by several studies, however, there is substantial variability in data collection, the methods used for model development and included predictors [2]. Among others, the most described predictors are: chronic dysfunction [38, 42, 46], episodes of acute rejection [3, 3841], death with functional graft [38, 46], glomerulonephritis [38], donor age [47], hypertension [47, 48], diabetes [41, 47], type of immunosuppression [47], delayed graft function [47], recipient age [3], race [3], albumin [3], proteinuria [3, 42, 47], low-density lipoprotein (LDL) cholesterol levels [48] and higher BMI [49]. However, some of them included in the analysis but that finally were not significant.

Risk factors for death

Identification and quantification of the relevant factors for death can improve patients’ individual risk assessment and decision-making [50]. In this study we confirm risk factors for death like recipient age, deceased donor, CMV disease, CMV disease, coronary stent, acute cellular rejection. Our findings show that older recipients are more likely to die, which is consistent with several published studies that report youngest age groups demonstrated a clear trend toward lower mortality compared with those ≥60–65 years [5053]. However, it must be said that long term patient survival in the elderly has been shown to be significantly better in transplant patients compared with remaining on the waitlist [5457]. Similarly to what happens with large series (Collaborative Transplant Study [58] and UNOS Register [59]), it is observed that living-donor kidney transplantation provides better outcomes than deceased-donor transplantation. Besides, it is associated with shorter transplant waiting list period and better early outcome [60]. We have found that CMV disease represents a risk of death in our population. This is one of the most important infectious complications in transplant recipient leading to significant morbidity and mortality [61]. Various direct and indirect detrimental effects occur because of CMV infection on patients and grafts. Indirect effects may include rejection, immunosuppression resulting in infections by other microorganisms, graft dysfunction, and poor survival of the kidney graft [62].

Cardiovascular disease (CVD) is frequent after kidney transplantation, is a major cause of morbidity and of death with functioning graft in recipients [63, 64]. We found as a risk factor in the model that coronary disease, specifically coronary stent placement, as a risk factor for death. OPTN/SRTR 2017 Annual Data Report: Kidney, Death with a functioning graft is the leading cause of graft loss in kidney transplant recipients, and a major cause of death is cardiovascular disease, accounting for about one third of known causes [65]. Another of the factors related in the model with the death of kidney transplant patients that we found was the presence of Acute Rejection (AR). Clayton et al., proposed that AR and its treatment may directly or indirectly affect longer-term outcomes for kidney transplant recipients, they found AR was also associated with death with a functioning graft (HR, 1.22; 95% CI, 1.08 to 1.36), and with death due to cardiovascular disease (HR, 1.30; 95% CI, 1.11 to 1.53) and concluded AR is associated with increased risks of longer-term graft failure and death, particularly death from cardiovascular disease and cancer. The results suggest AR remains an important short-term outcome to monitor in kidney transplantation [66].

Previous studies have attempted to identify and integrate risk factors for death into predictive models, including the pre-transplant variables gender, race, body mass index (BMI), time on dialysis, cause of end-stage renal disease, panel reactive antibodies, HLA mismatches, comorbidities such as diabetes, and heart failure, and donor age. In some models, the post-transplant factors Delayed graft function (DGF), and graft function were included [50], however, in our study population these were not significant. We think that in the case of diabetes the sample size was not sufficient.

Strengths and limitations

Unlike most previous studies, the main strength of this study is that our analysis includes a competing risk model. Many papers have pointed out the important issue of competing events in kidney transplantation [2, 4, 5, 7]. This method allowed us to determine graft loss risk factors differentiating those who increase recipient mortality. We believe that this integral view is best suited to a rational and patient-centered risk assessment.

Further, our cohort is the largest reporting risk factors for graft loss and death by a unique transplant program in Colombia and contributes to have a better understanding of Latin-American population as most of previous studies have been reported by transplant programs that treat mainly Caucasian patients. Other strengths include consistent data collection with a high degree of completeness and several variables.

Potential limitations attendant with the nature of data collection. The retrospective nature of our study prohibited adjusting for unmeasured confounding factors that may explain the association between independent factors and adverse graft outcomes. Besides, donor age was no considered in our analysis due to no available information.

On the other hand, variable selection with backward regression is not ideal. A fundamental problem with stepwise regression is that some real explanatory variables that have causal effects on the dependent variable may happen to not be statistically significant, while nuisance variables may be coincidentally significant. As a result, the model may fit the data well in-sample but do poorly out-of-sample. Unfortunately, penalized methods for Fine-Gray models have some limitations and the output from the crrp () function does not include convenient parameters such as a p-value. In addition, this package has not been maintained since its first commit in 2015.

We did not perform external validation, and this could be useful to assess the generalizability to other similar populations.

Conclusion

In summary, we found that stroke, BKV nephropathy, serum creatinine at 12 months and increased number of renal allograft biopsies were significant risk factors for graft loss. On the other hand, recipient age, acute cellular rejection, CMV disease were risk factors for death, and recipients with living donor had better survival compared to deceased-donor transplant and coronary stent. This contributes to have a better understanding of Latin-American population. However, the predictive performance of the models was mild.

Acknowledgments

We are grateful with Colombiana de Trasplantes to make this study possible.

References

  1. 1. Scheffner I., Gietzelt M., Abeling T., Marschollek M., and Gwinner W., “Patient Survival after Kidney Transplantation: Important Role of Graft-sustaining Factors as Determined by Predictive Modeling Using Random Survival Forest Analysis,” Transplantation, pp. 1095–1107, 2020. pmid:31403555
  2. 2. Kaboré R., Haller M. C., Harambat J., Heinze G., and LeffondrÉ K., “Risk prediction models for graft failure in kidney transplantation: A systematic review,” Nephrol. Dial. Transplant., vol. 32, no. February, pp. ii68–ii76, 2017. pmid:28206633
  3. 3. Shabir S. et al., “Predicting 5-year risk of kidney transplant failure: A prediction instrument using data available at 1 year posttransplantation,” Am. J. Kidney Dis., vol. 63, no. 4, pp. 643–651, 2014. pmid:24387794
  4. 4. El Ters M., Smith B. H., Cosio F. G., and Kremers W. K., “Competing Risk Analysis in Renal Allograft Survival: A New Perspective to an Old Problem,” Transplantation, vol. Publish Ah, 2020.
  5. 5. Lau B., Cole S. R., and Gange S. J., “Competing risk regression models for epidemiologic data,” Am. J. Epidemiol., vol. 170, no. 2, pp. 244–256, 2009. pmid:19494242
  6. 6. Schuster N. A., Hoogendijk E. O., Kok A. A. L., Twisk J. W. R., and Heymans M. W., “Ignoring competing events in the analysis of survival data may lead to biased results: a nonmathematical illustration of competing risk analysis,” J. Clin. Epidemiol., vol. 122, pp. 42–48, 2020. pmid:32165133
  7. 7. Gjertson D. W., Dabrowska D. M., Cui X., and Cecka J. M., “Four causes of cadaveric kidney transplant failure: A competing risk analysis,” Am. J. Transplant., vol. 2, no. 1, pp. 84–93, 2002. pmid:12095062
  8. 8. Prentice R. L., Kalbfleisch J. D., V Peterson A., Flournoy N., Farewell V. T., and Breslow N. E., “The Analysis of Failure Times in the Presence of Competing Risks,” Biometrics, vol. 34, no. 4, pp. 541–554, Feb. 1978. pmid:373811
  9. 9. Fine J. P. and Gray R. J., “A Proportional Hazards Model for the Subdistribution of a Competing Risk,” J. Am. Stat. Assoc., vol. 94, no. 446, pp. 496–509, Jun. 1999.
  10. 10. Sapir-Pichhadze R. et al., “Survival Analysis in the Presence of Competing Risks: The Example of Waitlisted Kidney Transplant Candidates,” Am. J. Transplant., vol. 16, no. 7, pp. 1958–1966, 2016. pmid:26751409
  11. 11. Wolbers M. et al., “Competing risks analyses: Objectives and approaches,” Eur. Heart J., vol. 35, no. 42, pp. 2936–2941, 2014. pmid:24711436
  12. 12. Kucirka L. M. et al., “Association of race and age with survival among patients undergoing dialysis,” JAMA—J. Am. Med. Assoc., vol. 306, no. 6, pp. 620–626, 2011. pmid:21828325
  13. 13. Salter M. L. et al., “Fractures and Subsequent Graft Loss and Mortality among Older Kidney Transplant Recipients,” J. Am. Geriatr. Soc., vol. 67, no. 8, pp. 1680–1688, 2019. pmid:31059126
  14. 14. Arcos E., José Pérez-Sáez M., Comas J., Lloveras J., Tort J., and Pascual J., “Assessing the Limits in Kidney Transplantation: Use of Extremely Elderly Donors and Outcomes in Elderly Recipients,” Transplantation, vol. 104, no. 1, pp. 176–183, 2020. pmid:30985579
  15. 15. Loupy A. et al., “The Banff 2015 Kidney Meeting Report: Current Challenges in Rejection Classification and Prospects for Adopting Molecular Pathology,” Am. J. Transplant., vol. 17, no. 1, pp. 28–41, 2017. pmid:27862883
  16. 16. Montgomery R. A., Tatapudi V. S., Leffell M. S., and Zachary A. A., “HLA in transplantation,” Nat. Rev. Nephrol. 2018 149, vol. 14, no. 9, pp. 558–570, Jul. 2018. pmid:29985463
  17. 17. Nino-Torres L., Garcia-Lopez A., Giron-Luque F., and Nino-Murcia A., “Retrospective Analysis of the Kidney Donor Profile Index to Predict Patient and Graft Survival at 5 Years Posttransplantation in a Colombian Cohort,” Transplant. Proc., vol. 7, pp. 1–7, 2021.
  18. 18. Y. Arias, M. A. Salinas Nova, and J. Montaño, Criterios de Asignación para Trasplante Renal en Colombia. 2018.
  19. 19. Collins G. S., Reitsma J. B., Altman D. G., and Moons K. G. M., “Transparent reporting of a multivariable prediction model for individual prognosis or diagnosis (TRIPOD): The TRIPOD statement,” Ann. Intern. Med., vol. 162, no. 1, pp. 55–63, 2015. pmid:25560714
  20. 20. Gray B., “Package ‘cmprsk.’” p. 13, 2020.
  21. 21. El-Zoghby Z. M. et al., “Identifying Specific Causes of Kidney Allograft Loss,” Am. J. Transplant., vol. 9, no. 3, pp. 527–535, Mar. 2009. pmid:19191769
  22. 22. Sellarés J. et al., “Understanding the causes of kidney transplant failure: The dominant role of antibody-mediated rejection and nonadherence,” Am. J. Transplant., vol. 12, no. 2, pp. 388–399, Feb. 2012. pmid:22081892
  23. 23. Naesens M. et al., “The histology of kidney transplant failure: A long-term follow-up study,” Transplantation, vol. 98, no. 4, pp. 427–435, 2014. pmid:25243513
  24. 24. Goldfarb D. A., “The natural history of chronic allograft nephropathy.,” J. Urol., vol. 173, no. 6, p. 2106, 2005. pmid:15879858
  25. 25. Solez K. et al., “International standardization of criteria for the histologic diagnosis of renal allograft rejection: The Banff working classification of kidney transplant pathology,” Kidney Int., vol. 44, no. 2, pp. 411–422, 1993. pmid:8377384
  26. 26. Nankivell B. J. and Chapman J. R., “Chronic allograft nephropathy: Current concepts and future directions,” Transplantation, vol. 81, no. 5, pp. 643–654, 2006. pmid:16534463
  27. 27. Metter C. and Torrealba J. R., “Pathology of the kidney allograft,” Semin. Diagn. Pathol., vol. 37, no. 3, pp. 148–153, 2020. pmid:32249077
  28. 28. Salcedo-Herrera S., Pinto Ramirez J. L., García-Lopez A., Amaya-Nieto J., and Girón-Luque F., “Acute Rejection in Kidney Transplantation and Early Beginning of Tacrolimus,” Transplant. Proc., vol. 51, no. 6, pp. 1758–1762, 2019. pmid:31399163
  29. 29. Sakai K., Oguchi H., Muramatsu M., and Shishido S., “Protocol graft biopsy in kidney transplantation.,” Nephrology (Carlton)., vol. 23 Suppl 2, pp. 38–44, Jul. 2018. pmid:29968403
  30. 30. Findlay M. D. et al., “Risk factors and outcome of stroke in renal transplant recipients,” Clin. Transplant., vol. 30, no. 8, pp. 918–924, 2016. pmid:27218240
  31. 31. Huang S.-T. et al., “The Risk of Stroke in Kidney Transplant Recipients with End-Stage Kidney Disease,” Int. J. Environ. Res. Public Health, vol. 16, no. 3, p. 326, Jan. 2019. pmid:30682846
  32. 32. Ramos E. et al., “Clinical Course of Polyoma Virus Nephropathy in 67 Renal Transplant Patients,” J. Am. Soc. Nephrol., vol. 13, no. 8, pp. 2145–2151, 2002. pmid:12138148
  33. 33. Leeaphorn N., Thongprayoon C., Chon W. J., Cummings L. S., Mao M. A., and Cheungpasitporn W., “Outcomes of kidney retransplantation after graft loss as a result of BK virus nephropathy in the era of newer immunosuppressant agents,” Am. J. Transplant., vol. 20, no. 5, pp. 1334–1340, May 2020. pmid:31765056
  34. 34. Vasudev B., Hariharan S., Hussain S. A., Zhu Y.-R., Bresnahan B. A., and Cohen E. P., “BK virus nephritis: Risk factors, timing, and outcome in renal transplant recipients,” Kidney Int., vol. 68, no. 4, pp. 1834–1839, 2005. pmid:16164661
  35. 35. Favi E. et al., “Incidence, risk factors, and outcome of BK polyomavirus infection after kidney transplantation,” World J. Clin. cases, vol. 7, no. 3, pp. 270–290, Feb. 2019. pmid:30746369
  36. 36. Gago M., Cornell L. D., Kremers W. K., Stegall M. D., and Cosio F. G., “Kidney Allograft Inflammation and Fibrosis, Causes and Consequences,” Am. J. Transplant., vol. 12, no. 5, pp. 1199–1207, May 2012. pmid:22221836
  37. 37. Myint T. M., Chong C. H. Y., Wyld M., Nankivell B., Kable K., and Wong G., “Polyoma BK Virus in Kidney Transplant Recipients: Screening, Monitoring, and Management,” Transplantation, vol. 106, no. 1, pp. E76–E89, 2022. pmid:33908382
  38. 38. Tanabe K., Takahashi K., and Toma H., “Causes of long-term graft failure in renal transplantation,” World J. Urol., vol. 14, no. 4, pp. 230–235, 1996. pmid:8873436
  39. 39. Meier-Kriesche H. U. et al., “Increased impact of acute rejection on chronic allograft failure in recent era.,” Transplantation, vol. 70, no. 7, pp. 1098–1100, Oct. 2000. pmid:11045649
  40. 40. Pallardó Mateu L. M., Sancho Calabuig A., Capdevila Plaza L., and Franco Esteve A., “Acute rejection and late renal transplant failure: risk factors and prognosis.,” Nephrol. Dial. Transplant. Off. Publ. Eur. Dial. Transpl. Assoc.—Eur. Ren. Assoc., vol. 19 Suppl 3, pp. iii38–42, Jun. 2004. pmid:15192134
  41. 41. Cole E. H., Johnston O., Rose C. L., and Gill J. S., “Impact of acute rejection and new-onset diabetes on long-term transplant graft and patient survival.,” Clin. J. Am. Soc. Nephrol., vol. 3, no. 3, pp. 814–821, May 2008. pmid:18322046
  42. 42. Fellström B. et al., “Risk factors for reaching renal endpoints in the Assessment of Lescol in Renal Transplantation (ALERT) trial,” Transplantation, vol. 79, no. 2, pp. 205–212, 2005. pmid:15665769
  43. 43. Hariharan S., McBride M. A., Cherikh W. S., Tolleris C. B., Bresnahan B. A., and Johnson C. P., “Post-transplant renal function in the first year predicts long-term kidney transplant survival.,” Kidney Int., vol. 62, no. 1, pp. 311–318, Jul. 2002. pmid:12081593
  44. 44. Salvadori M. et al., “Estimated one-year glomerular filtration rate is the best predictor of long-term graft function following renal transplant.,” Transplantation, vol. 81, no. 2, pp. 202–206, Jan. 2006. pmid:16436963
  45. 45. Rumsfeld J. S. et al., “Recurrent glomerulonephritis after kidney transplantation: risk factors and allograft outcomes,” Transplantation, vol. 13, no. 2, pp. 1–5, Oct. 2019.
  46. 46. Morales J. M. et al., “Risk factors for graft loss and mortality after renal transplantation according to recipient age: A prospective multicentre study,” Nephrol. Dial. Transplant., vol. 27, no. SUPPL.4, 2012. pmid:23258810
  47. 47. Khalkhali H. R., Ghafari A., Hajizadeh E., and Kazemnejad A., “Risk factors of long-term graft loss in renal transplant recipients with chronic allograft dysfunction.,” Exp. Clin. Transplant. Off. J. Middle East Soc. Organ Transplant., vol. 8, no. 4, pp. 277–282, Dec. 2010. pmid:21143092
  48. 48. Soltanian A. R., Mahjub H., Taghizadeeh-Afshari A., Gholami G., and Sayyadi H., “Identify survival predictors of the first kidney transplantation: A retrospective cohort study,” Iran. J. Public Health, vol. 44, no. 5, pp. 683–689, 2015. pmid:26284210
  49. 49. Papalia T., Greco R., Lofaro D., Maestripieri S., Mancuso D., and Bonofiglio R., “Impact of body mass index on graft loss in normal and overweight patients: Retrospective analysis of 206 renal transplants,” Clin. Transplant., vol. 24, no. 6, pp. 25–30, 2010. pmid:20482558
  50. 50. Abeling T. et al., “Risk factors for death in kidney transplant patients: analysis from a large protocol biopsy registry,” Nephrol. Dial. Transplant., vol. 34, no. 7, pp. 1171–1181, Jul. 2019. pmid:29860340
  51. 51. Fabrizii V. et al., “Patient and Graft Survival in Older Kidney Transplant Recipients: Does Age Matter?,” J. Am. Soc. Nephrol., vol. 15, no. 4, pp. 1052–1060, 2004. pmid:15034109
  52. 52. Mohamed Ali A. A. et al., “Renal transplantation in the elderly: South Indian experience.,” Int. Urol. Nephrol., vol. 43, no. 1, pp. 265–271, Mar. 2011. pmid:21203840
  53. 53. Gerbase-DeLima M., de Marco R., Monteiro F., Tedesco-Silva H., Medina-Pestana J. O., and Mine K. L., “Impact of Combinations of Donor and Recipient Ages and Other Factors on Kidney Graft Outcomes,” Front. Immunol., vol. 11, no. May, pp. 1–9, 2020. pmid:32528472
  54. 54. Rao P. S., Merion R. M., Ashby V. B., Port F. K., a Wolfe R., and Kayler L. K., “Renal transplantation in elderly patients older than 70 years of age: results from the Scientific Registry of Transplant Recipients.,” Transplantation, vol. 83, no. 8, pp. 1069–1074, 2007. pmid:17452897
  55. 55. Schold J., Srinivas T. R., Sehgal A. R., and Meier-Kriesche H.-U., “Half of kidney transplant candidates who are older than 60 years now placed on the waiting list will die before receiving a deceased-donor transplant,” Clin. J. Am. Soc. Nephrol., vol. 4, no. 7, pp. 1239–1245, Jul. 2009. pmid:19541814
  56. 56. Schold J. D., Srinivas T. R., Kayler L. K., and Meier-Kriesche H. U., “The overlapping risk profile between dialysis patients listed and not listed for renal transplantation.,” Am. J. Transplant. Off. J. Am. Soc. Transplant. Am. Soc. Transpl. Surg., vol. 8, no. 1, pp. 58–68, Jan. 2008. pmid:17979999
  57. 57. Hernández D. et al., “Mortality in Elderly Waiting-List Patients Versus Age-Matched Kidney Transplant Recipients: Where is the Risk?,” Kidney Blood Press. Res., vol. 43, no. 1, pp. 256–275, 2018. pmid:29490298
  58. 58. Opelz G., Wujciak T., Döhler B., Scherer S., and Mytilineos J., “HLA compatibility and organ transplant survival. Collaborative Transplant Study.,” Rev. Immunogenet., vol. 1, no. 3, pp. 334–342, 1999. pmid:11256424
  59. 59. Cecka J. M., “The UNOS renal transplant registry.,” Clin. Transpl., pp. 1–18, 2001. pmid:12211771
  60. 60. Basiri A. et al., “Living or deceased-donor kidney transplant: the role of psycho-socioeconomic factors and outcomes associated with each type of transplant,” Int. J. Equity Health, vol. 19, no. 1, p. 79, 2020. pmid:32487079
  61. 61. Minz R. W. et al., “Cytomegalovirus Infection in Postrenal Transplant Recipients: 18 Years’ Experience From a Tertiary Referral Center,” Transplant. Proc., vol. 52, no. 10, pp. 3173–3178, 2020. pmid:32624232
  62. 62. Freeman R. B., “The ‘Indirect’ effects of cytomegalovirus infection: Minireview,” Am. J. Transplant., vol. 9, no. 11, pp. 2453–2458, 2009. pmid:19843027
  63. 63. Ojo A. O., “Cardiovascular complications after renal transplantation and their prevention,” Transplantation, vol. 82, no. 5, pp. 603–611, 2006. pmid:16969281
  64. 64. Birdwell K. A., “Post-Transplant Cardiovascular Disease,” vol. 16, pp. 1878–1889, 2021. pmid:34556500
  65. 65. Hart A. et al., “OPTN/SRTR 2017 Annual Data Report: Kidney,” Am. J. Transplant, vol. 19, pp. 19–123, 2019. pmid:30811893
  66. 66. Clayton P. A., McDonald S. P., Russ G. R., and Chadban S. J., “Long-term outcomes after acute rejection in kidney transplant recipients: An ANZDATA analysis,” J. Am. Soc. Nephrol., vol. 30, no. 9, pp. 1697–1707, 2019. pmid:31308074