Archetypal Analysis of Injury in Kidney Transplant Biopsies Identifies Two Classes of Early AKI

Philip F Halloran, Georg A Böhmig, Jonathan Bromberg, Gunilla Einecke, Farsad A Eskandary, Gaurav Gupta, Marek Myslak, Ondrej Viklicky, Agnieszka Perkowska-Ptasinska, Katelynn S Madill-Thomsen, INTERCOMEX Investigators, Philip F Halloran, Georg A Böhmig, Jonathan Bromberg, Gunilla Einecke, Farsad A Eskandary, Gaurav Gupta, Marek Myslak, Ondrej Viklicky, Agnieszka Perkowska-Ptasinska, Katelynn S Madill-Thomsen, INTERCOMEX Investigators

Abstract

All transplanted kidneys are subjected to some degree of injury as a result of the donation-implantation process and various post-transplant stresses such as rejection. Because transplants are frequently biopsied, they present an opportunity to explore the full spectrum of kidney response-to-wounding from all causes. Defining parenchymal damage in transplanted organs is important for clinical management because it determines function and survival. In this study, we classified the scenarios associated with parenchymal injury in genome-wide microarray results from 1,526 kidney transplant indication biopsies collected during the INTERCOMEX study. We defined injury groups by using archetypal analysis (AA) of scores for gene sets and classifiers previously identified in various injury states. Six groups and their characteristics were defined in this population: No injury, minor injury, two classes of acute kidney injury ("AKI," AKI1, and AKI2), chronic kidney disease (CKD), and CKD combined with AKI. We compared the two classes of AKI, namely, AKI1 and AKI2. AKI1 had a poor function and increased parenchymal dedifferentiation but minimal response-to-injury and inflammation, instead having increased expression of PARD3, a gene previously characterized as being related to epithelial polarity and adherens junctions. In contrast, AKI2 had a poor function and increased response-to-injury, significant inflammation, and increased macrophage activity. In random forest analysis, the most important predictors of function (estimated glomerular filtration rate) and graft loss were injury-based molecular scores, not rejection scores. AKI1 and AKI2 differed in 3-year graft survival, with better survival in the AKI2 group. Thus, injury archetype analysis of injury-induced gene expression shows new heterogeneity in kidney response-to-wounding, revealing AKI1, a class of early transplants with a poor function but minimal inflammation or response to injury, a deviant response characterized as PC3, and an increased risk of failure. Given the relationship between parenchymal injury and kidney survival, further characterization of the injury phenotypes in kidney transplants will be important for an improved understanding that could have implications for understanding native kidney diseases (ClinicalTrials.gov #NCT01299168).

Keywords: archetypes; biopsy; gene expression; injury; kidney transplantation.

Conflict of interest statement

PH is a consultant to Natera, holds shares in Transcriptome Sciences Inc (TSI), a University of Alberta research company dedicated to developing molecular diagnostics, supported in part by a licensing agreement between TSI and Thermo Fisher, and by a research grant from Natera. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2022 Halloran, Böhmig, Bromberg, Einecke, Eskandary, Gupta, Myslak, Viklicky, Perkowska-Ptasinska, Madill-Thomsen and the INTERCOMEX Investigators.

Figures

Figure 1
Figure 1
Injury-based PCA, colored by injury archetype cluster. (A) PC2 vs. PC1 scores for each biopsy; (B) PC2 vs. PC3 scores; (C) moving averages of standardized injury archetype scores (window size = 400 biopsies). As there are large differences in mean scores between archetypes, all scores were standardized to a mean of 0.0 before plotting. The y-axis is in standard deviation units. Biopsies sorted by ascending time of biopsy post-transplant.
Figure 2
Figure 2
Graft survival analyses. (A) Random forests assessing variable importance in predicting low eGFR (eGFR < 30). (B) Random forests assessing variable importance in predicting 3-year death-censored graft survival. (C) Kaplan-Meier curves showing 3-year post-biopsy actuarial survival in the six AA injury archetype groups.
Figure 3
Figure 3
Kaplan-Meier curves showing actuarial 3-year death-censored graft survival in the injury AA groups after excluding all biopsies with rejection. (A) Injury AA groups were compared across time post-biopsy. We also examined the differences between, (B) short-term (≤42 days post-biopsy), and (C) long-term (>42 days post-biopsy) survival probabilities. The CKD/AKI group was removed from these plots due to the very small group size and scarcity of events.
Figure 4
Figure 4
Schematic diagram representing the relationships between sources of injury and response to injury in kidney transplant biopsies based on these analyses. Interplay between sources of injury, pre-existing limitations such as aging, and response to injury by the nephron. There are two routes to irreversible nephron shutdown, namely, the epithelial injury and through glomerulus injury. Epithelial injury should trigger the response-to-wounding, which involves epithelium, matrix, and microcirculation, and evokes innate immunity. Failure to mount a response to wounding and adopting a “PC3”-related response (e.g., PARD3) with minimal inflammation leads to failure to recover. Many sources of injury (separate from and including rejection) interact with the nephron epithelium, producing acute kidney injury (AKI). In this instance, the epithelium can be repaired and the organ can recover, or progress to nephron failure. Alternatively, aging and/or ABMR can contribute to glomerular disease and ABMR can additionally affect the microcirculation, affecting the glomerulus and again causing nephron shutdown, which eventually leads to chronic kidney disease (CKD). If this occurs, a loss of nephrons and end-stage renal disease may occur. Different sources of injury may interact to cause many forms of injury, and injury itself predicts the graft survival while the rejection status does not. Thus, defining the heterogeneity within biopsy injury is an important part of clinical management.

References

    1. Halloran PF, Famulski KS, Reeve J. Molecular assessment of disease states in kidney transplant biopsy samples. Nat Rev Nephrol. (2016) 12:534–48. 10.1038/nrneph.2016.85
    1. Reeve J, Bohmig GA, Eskandary F, Einecke G, Gupta G, Madill-Thomsen K, et al. . Generating automated kidney transplant biopsy reports combining molecular measurements with ensembles of machine learning classifiers. Am J Transplant. (2019) 19:2719–31. 10.1111/ajt.15351
    1. Reeve J, Bohmig GA, Eskandary F, Einecke G, Lefaucheur C, Loupy A, et al. . Assessing rejection-related disease in kidney transplant biopsies based on archetypal analysis of molecular phenotypes. JCI Insight. (2017) 2:e94197. 10.1172/jci.insight.94197
    1. Venner JM, Famulski KS, Reeve J, Chang J, Halloran PF. Relationships among injury, fibrosis, and time in human kidney transplants. JCI Insight. (2016) 1:e85323. 10.1172/jci.insight.85323
    1. Halloran PF, Reeve J, Akalin E, Aubert O, Bohmig GA, Brennan D, et al. . Real time central assessment of kidney transplant indication biopsies by microarrays: the INTERCOMEX study. Am J Transplant. (2017) 17:2851–62. 10.1111/ajt.14329
    1. Halloran PF, Bohmig GA, Bromberg JS, Budde K, Gupta G, Einecke G, et al. . Discovering novel injury features in kidney transplant biopsies associated with TCMR and donor aging. Am J Transplant. (2021) 21:1725–39. 10.1111/ajt.16374
    1. R Foundation for Statistical Computing . RCT. R: A Language and Environment for Statistical Computing. Vienna: R Foundation for statistical Computing; (2019) Available online at: (accessed November 2, 2021).
    1. Madill-Thomsen KS, Wiggins RC, Eskandary F, Bohmig GA, Halloran PF. The effect of cortex/medulla proportions on molecular diagnoses in kidney transplant biopsies: rejection and injury can be assessed in medulla. Am J Transplant. (2017) 17:2117–28. 10.1111/ajt.14233
    1. Loupy A, Haas M, Solez K, Racusen L, Glotz D, Seron D, et al. . The Banff 2015 kidney meeting report: current challenges in rejection classification and prospects for adopting molecular pathology. Am J Transplant. (2017) 17:28–41. 10.1111/ajt.14107
    1. Mengel M, Loupy A, Haas M, Roufosse C, Naesens M, Akalin E, et al. . Banff 2019 Meeting Report: molecular diagnostics in solid organ transplantation-Consensus for the Banff Human Organ Transplant (B-HOT) gene panel and open source multicenter validation. Am J Transplant. (2020) 20:2305–17. 10.1111/ajt.16059
    1. Halloran PF, Matas A, Kasiske BL, Madill-Thomsen KS, Mackova M, Famulski KS. Molecular phenotype of kidney transplant indication biopsies with inflammation in scarred areas. Am J Transplant. (2019) 19:1356–70. 10.1111/ajt.15178
    1. Einecke G, Reeve J, Gupta G, Bohmig GA, Eskandary F, Bromberg JS, et al. . 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. Am J Transplant. (2021) 21:1391–401. 10.1111/ajt.16161
    1. Kuhn M. caret: Classification and Regression Training. R package version 6.0-81. (2018). Available online at: (accessed November 2, 2021).
    1. Heil M, Land WG. Danger signals - damaged-self recognition across the tree of life. Front Plant Sci. (2014) 5:578. 10.3389/fpls.2014.00578
    1. Land WG, Agostinis P, Gasser S, Garg AD, Linkermann A. Transplantation and damage-associated molecular patterns (DAMPs). Am J Transplant. (2016) 16:3338–61. 10.1111/ajt.13963
    1. Famulski KS, de Freitas DG, Kreepala C, Chang J, Sellares J, Sis B, et al. . Molecular phenotypes of acute kidney injury in human kidney transplants. J Am Soc Nephrol. (2012) 23:948–58. 10.1681/ASN.2011090887
    1. Famulski KS, Broderick G, Einecke G, Hay K, Cruz J, Sis B, et al. . Transcriptome analysis reveals heterogeneity in the injury response of kidney transplants. Am J Transplant. (2007) 7:2483–95. 10.1111/j.1600-6143.2007.01980.x
    1. Einecke G, Reeve J, Mengel M, Sis B, Bunnag S, Mueller TF, et al. . Expression of B cell and immunoglobulin transcripts is a feature of inflammation in late allografts. Am J Transplant. (2008) 8:1434–43. 10.1111/j.1600-6143.2008.02232.x
    1. Mengel M, Reeve J, Bunnag S, Einecke G, Sis B, Mueller T, et al. . Molecular correlates of scarring in kidney transplants: the emergence of mast cell transcripts. Am J Transplant. (2009) 9:169–78. 10.1111/j.1600-6143.2008.02462.x
    1. Einecke G, Kayser D, Vanslambrouck JM, Sis B, Reeve J, Mengel M, et al. . Loss of solute carriers in T cell-mediated rejection in mouse and human kidneys: an active epithelial injury-repair response. Am J Transplant. (2010) 10:2241–51. 10.1111/j.1600-6143.2010.03263.x
    1. Eugster MJA, Leisch F. From spider-man to hero - archetypal analysis in R. J Stat Softw. (2009) 30:1–23. 10.18637/jss.v030.i08
    1. Thogersen JC, Morup M, Damkiaer S, Molin S, Jelsbak L. Archetypal analysis of diverse Pseudomonas aeruginosa transcriptomes reveals adaptation in cystic fibrosis airways. BMC Bioinformat. (2013) 14:279. 10.1186/1471-2105-14-279
    1. Korem Y, Szekely P, Hart Y, Sheftel H, Hausser J, Mayo A, et al. . Geometry of the gene expression space of individual cells. PLoS Comput Biol. (2015) 11:e1004224. 10.1371/journal.pcbi.1004224
    1. Ishwaran H, Kogalur UB. Fast Unified Random Forests for Survival, Regression, and Classification (RF-SRC). R package version 3.0.2. (2019). Available online at:
    1. Naik AS, Afshinnia F, Cibrik D, Hodgin JB, Wu F, Zhang M, et al. . Quantitative podocyte parameters predict human native kidney and allograft half-lives. JCI Insight. (2016) 1:e86943. 10.1172/jci.insight.86943
    1. Ronco C, Bellomo R, Kellum JA. Acute kidney injury. Lancet. (2019) 394:1949–64. 10.1016/S0140-6736(19)32563-2

Source: PubMed

3
Abonner