Abstract
Chronic allograft dysfunction (CAD), characterized histologically by interstitial fibrosis and tubular atrophy, is the major cause of kidney allograft loss. Here, using single nuclei RNA sequencing and transcriptome analysis, we identified the origin, functional heterogeneity, and regulation of fibrosis-forming cells in kidney allografts with CAD. A robust technique was used to isolate individual nuclei from kidney allograft biopsies and successfully profiled 23,980 nuclei from five kidney transplant recipients with CAD and 17,913 nuclei from three patients with normal allograft function. Our analysis revealed two distinct states of fibrosis in CAD; low and high extracellular matrix (ECM) with distinct kidney cell subclusters, immune cell types, and transcriptional profiles. Imaging mass cytometry analysis confirmed increased ECM deposition at the protein level. Proximal tubular cells transitioned to an injured mixed tubular (MT1) phenotype comprised of activated fibroblasts and myofibroblast markers, generated provisional ECM which recruited inflammatory cells, and served as the main driver of fibrosis. MT1 cells in the high ECM state achieved replicative repair evidenced by dedifferentiation and nephrogenic transcriptional signatures. MT1 in the low ECM state showed decreased apoptosis, decreased cycling tubular cells, and severe metabolic dysfunction, limiting the potential for repair. Activated B, T and plasma cells were increased in the high ECM state, while macrophage subtypes were increased in the low ECM state. Intercellular communication between kidney parenchymal cells and donor-derived macrophages, detected several years post-transplantation, played a key role in injury propagation. Thus, our study identified novel molecular targets for interventions aimed to ameliorate or prevent allograft fibrogenesis in kidney transplant recipients.