# Circulating Donor-Specific Anti-HLA Antibodies Associate With Immune Activation Independent of Kidney Transplant Histopathological Findings

Elisabet van Loon, Baptiste Lamarthée, Thomas Barba, Sandra Claes, Maarten Coemans, Henriette de Loor, Marie-Paule Emonds, Priyanka Koshy, Dirk Kuypers, Paul Proost, et al.

► To cite this version:

Elisabet van Loon, Baptiste Lamarthée, Thomas Barba, Sandra Claes, Maarten Coemans, et al. Circulating Donor-Specific Anti-HLA Antibodies Associate With Immune Activation Independent of Kidney Transplant Histopathological Findings. Frontiers in Immunology, 2022, 13, 10.3389/fimmu.2022.818569. [hal-04963260](https://hal.science/hal-04963260v1)

HAL Id: hal-04963260

https://hal.science/hal-04963260v1

Submitted on 24 Feb 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Circulating Donor-Specific Anti-HLA Antibodies Associate With Immune Activation Independent of Kidney Transplant Histopathological Findings

Elisabet Van Loon<sup>1,2</sup>, Baptiste Lamarthée<sup>1</sup>, Thomas Barba<sup>3</sup>, Sandra Claes<sup>4</sup>, Maarten Coemans<sup>1,5</sup>, Henriette de Loor<sup>1</sup>, Marie-Paule Emonds<sup>1,6</sup>, Priyanka Koshy<sup>7</sup>, Dirk Kuypers<sup>1,2</sup>, Paul Proost<sup>8</sup>, Aleksandar Senev<sup>1,6</sup>, Ben Sprangers<sup>1,2,8</sup>, Claire Tinel<sup>1</sup>, Olivier Thaunat<sup>3</sup>, Amaryllis H. Van Craenenbroeck<sup>1,2</sup>, Dominique Schols<sup>4</sup> and Maarten Naesens<sup>1,2*</sup>

OPEN ACCESS

**Edited by:**
Mahzad Akbarpour,
University of Chicago Medicine,
United States

**Reviewed by:**
Anna Valujskikh,
Cleveland Clinic, United States
Nuala Mooney,
Institut National de la Santé et de la
Recherche Médicale (INSERM),
France

*Correspondence:
Maarten Naesens
maarten.naesens@uzleuven.be

**Specialty section:**
This article was submitted to
Alloimmunity and Transplantation,
a section of the journal
Frontiers in Immunology

**Received:** 19 November 2021
**Accepted:** 10 January 2022
**Published:** 23 February 2022

**Citation:**
Van Loon E, Lamarthée B, Barba T, Claes S, Coemans M, de Loor H, Emonds M-P, Koshy P, Kuypers D, Proost P, Senev A, Sprangers B, Tinel C, Thaunat O, Van Craenenbroeck AH, Schols D and Naesens M (2022) Circulating Donor-Specific Anti-HLA Antibodies Associate With Immune Activation Independent of Kidney Transplant Histopathological Findings. *Front. Immunol.* 13:818569.
doi: 10.3389/fimmu.2022.818569

<sup>1</sup> Department of Microbiology, Immunology and Transplantation, Nephrology and Kidney Transplantation Research Group, Katholieke Universiteit (KU) Leuven, Leuven, Belgium, <sup>2</sup> Department of Nephrology and Kidney Transplantation, University Hospitals Leuven, Leuven, Belgium, <sup>3</sup> Department of Transplantation, Nephrology and Clinical Immunology, Edouard Herriot Hospital Lyon, Hospices Civils de Lyon, Lyon, France, <sup>4</sup> Department of Microbiology, Immunology and Transplantation, Laboratory of Virology and Chemotherapy, Rega Institute, Katholieke Universiteit (KU) Leuven, Leuven, Belgium, <sup>5</sup> Leuven Biostatistics and Statistical Bioinformatics Centre, Department of Public Health and Primary Care, Katholieke Universiteit (KU) Leuven, Leuven, Belgium, <sup>6</sup> Histocompatibility and Immunogenetics Laboratory, Red Cross-Flanders, Mechelen, Belgium, <sup>7</sup> Department of Imaging and Pathology, Katholieke Universiteit (KU) Leuven, Leuven, Belgium, <sup>8</sup> Laboratory of Molecular Immunology, Department of Microbiology, Immunology and Transplantation, Rega Institute for Medical Research, Katholieke Universiteit (KU) Leuven, Leuven, Belgium

Despite the critical role of cytokines in allograft rejection, the relation of peripheral blood cytokine profiles to clinical kidney transplant rejection has not been fully elucidated. We assessed 28 cytokines through multiplex assay in 293 blood samples from kidney transplant recipients at time of graft dysfunction. Unsupervised hierarchical clustering identified a subset of patients with increased pro-inflammatory cytokine levels. This patient subset was hallmarked by a high prevalence (75%) of donor-specific anti-human leukocyte antigen antibodies (HLA-DSA) and histological rejection (70%) and had worse graft survival compared to the group with low cytokine levels (HLA-DSA in 1.7% and rejection in 33.7%). Thirty percent of patients with high pro-inflammatory cytokine levels and HLA-DSA did not have histological rejection. Exploring the cellular origin of these cytokines, we found a corresponding expression in endothelial cells, monocytes, and natural killer cells in single-cell RNASeq data from kidney transplant biopsies. Finally, we confirmed secretion of these cytokines in HLA-DSA-mediated cross talk between endothelial cells, NK cells, and monocytes. In conclusion, blood pro-inflammatory cytokines are increased in kidney transplant patients with HLA-DSA, even in the absence of histology of rejection. These observations challenge the concept that histology is the gold standard for identification of ongoing allo-immune activation after transplantation.

**Keywords:** cytokines, chemokines, donor-specific anti-HLA antibodies, allograft rejection, kidney transplantation
# INTRODUCTION

Despite marked improvement in short-term graft outcome through reduced incidence of acute T-cell-mediated rejection (TCMR) with the current immunosuppressive armamentarium, long-term kidney allograft survival remains suboptimal (1, 2). An important cause for this long-term graft failure is acute rejection, especially antibody-mediated rejection (ABMR) (3, 4). ABMR is initiated by donor-specific antibodies (DSA), either to human leukocyte antigens (HLA) or less commonly, but also less easily detectable, to other donor-recipient mismatched antigens (5, 6). Part of the explanation for the impaired graft survival associated with ABMR (7–10) is the lack of proven effective treatments to prevent or treat ABMR (4, 11, 12).

In the search for biomarkers and therapeutic targets for acute rejection, several research groups have suggested to study cytokines, as they play a crucial role in the pathophysiology of rejection (13–31). However, studies focusing on only one or few molecules miss the complexity of the interactions and the interplay between the full landscape of cytokines. Also, details on the specificity of some cytokines for allo-immune processes, different rejection subtypes and associated histological lesions is lacking.

Therefore, we evaluated the complex landscape of 28 cytokines, chemokines, and growth factors in the blood of patients with a broad range of acute rejection types and histological lesions and analyzed the cellular origin of the most relevant cytokines, through publicly available single-cell RNASeq data and through *in vitro* models of HLA-DSA-mediated activation of NK cells, monocytes, and endothelial cells.

# MATERIALS AND METHODS

## Patient Population

We included all adequate for-cause (“indication”) biopsies performed within the first year after transplantation between August 07, 2012, and July 13, 2016, from consenting patients who received a single kidney transplant at the University Hospitals Leuven, Belgium. All patients gave written informed consent for collection and analysis in the Kidney Transplant Biobank, approved by the local ethical committee (S53364 and S61971). Details on the data collection are described in the **Supplementary Appendix**.

## Histopathology

Histological lesions were semiquantitatively scored according to the Banff consensus (32), as reported previously (33). Diagnosis of the phenotypes of ABMR, TCMR, and borderline changes was

established using the Banff 2019 criteria (34). More details are described in the **Supplementary Appendix**.

## Detection of Circulating Anti-HLA Antibodies

The follow-up of anti-HLA antibodies was systematically monitored in one histocompatibility laboratory (HILA—Belgian Red Cross Flanders); details on this assessment were previously published (35), and further details are provided in the **Supplementary Appendix**.

## Cytokine Quantification

Peripheral blood serum samples were analyzed for 28 cytokines, chemokines, and growth factors with Bio-Plex immunoassay using Luminex magnetic beads following the manufacturer's instructions (Bio-Rad, M50-0KCAFOY; Temse, Belgium) using a 27-multiplex panel and an additional single-plex one for CXCL9. Details are provided in **Supplementary Table S1** and **Supplementary Appendix**.

## Single-Cell RNA Sequencing Data Analysis

Previously published human single-cell data from two ABMR biopsies and five healthy references corresponding to transplant surveillance biopsies were used. The associated raw counts or matrices were downloaded from the Gene Expression Omnibus (GEO, GSE145927, https://www.ncbi.nlm.nih.gov/geo) (36) and Kidney Precision Medicine Project (https://atlas.kpmp.org/repository). Details on the gene expression analyses are described in the **Supplementary Appendix**. The CellChat R package (37) was used to analyze the cell-cell communication between these cells.

## NK Cell and Monocyte Sorting and Coculture With Endothelial Cells

Peripheral blood mononuclear cells were isolated from the blood of healthy volunteers from the Etablissement Français du Sang (Lyon, France) or from the Belgian Red Cross Flanders by Ficoll gradient centrifugation (Eurobio, Courtaboeuf, France). Cell sorting methods are detailed in **Supplementary Appendix**. Purified NK cells and non-classical monocytes were then cocultured with glomerular endothelial cells. In a first experiment, NK cells were cocultured with the human conditionally immortalized glomerular endothelial cell line (38) (ciGENC; HLA-A2), in combination with either anti-HLA-A2 DSA-containing serum or control human serum. After 4 h of coculture, supernatants were collected. In a second experiment, NK cells and non-classical monocytes were cocultured with primary glomerular endothelial cells (GENC, Cell Systems, USA), after incubation with either anti-HLA-A, -B, -C purified antibody (BD Biosciences, San Jose, CA, USA; cat #560187), or control isotype (BD Biosciences, cat #553447). After 24 h of coculture, supernatants were collected. Details on the two experiments are specified in the **Supplementary Appendix**.

Levels of 27 cytokines, chemokines, and growth factors in the supernatants were assessed using the same 27-multiplex analysis as used for the clinical samples, following the manufacturer's

<footnote>**Abbreviations:** ABMR, antibody-mediated rejection; APC, antigen-presenting cells; DSA, donor-specific antibodies; ECa, activated endothelial cells; FcγR, Fc-gamma receptor; FDR, false discovery rate; HLA, human leukocyte antigen; NK, natural killer; PCA, principal component analysis; scRNASeq, single-cell RNA sequencing; TCMR, T cell-mediated rejection.</footnote>
instructions (Bio-Rad, M50-0KCAF0Y; Bio-Rad, Nazareth, Belgium).

## Statistical Analysis

We report descriptive statistics using mean and standard deviation (or median and interquartile range for skewed distributions) for continuous variables or numbers, and percentages for discrete variables. Protein levels were log10-transformed for analysis. We restricted the Pearson correlation analyses to the first sample per patient, to avoid dependent measurements. Principal component analysis (PCA) and heatmap analysis were done using Euclidean distance to obtain the distance matrix and complete agglomeration method for clustering in the heatmap analysis. In logistic regression analysis, we corrected for repeat sampling by using mixed models with a random intercept. We controlled for multiple testing using the false discovery rate (FDR). Cox regression analyses and Kaplan-Meier survival curves with log-rank testing were used to assess associations with death-censored graft failure. Patients were censored at the time of death with a functioning graft or at the time of the last follow-up date. Time to failure was counted starting at the day of sample. We assessed the differences in protein levels in the *in vitro* experiments using paired T-tests on the log10-transformed protein levels. We used SAS (version 9.4; SAS Institute, Cary, NC), R (version 3.6.3, R core team, 2014, packages corrplot, FactoMineR, heatmap3), and GraphPad Prism (version 8; GraphPad Software, San Diego, CA) for statistical analysis and data presentation. Figures were created with BioRender.com.

## RESULTS

### Blood Cytokine Profiles Associate With HLA-DSA

In the first samples of each patient (N = 192; **Figure 1A, Table 1** and **Supplementary Table S2**), strong positive correlations were demonstrated among cytokines CXCL10/IP10, PDGF-BB, IL-6, GM-CSF, IL-9, CCL5/RANTES, basic-FGF, CCL4/MIP-1β, TNF-α, IFN-γ, CCL2/MCP-1, and IL-17α, followed by IL-15, IL-1β, and IL-2, and negative correlations with IL-1Ra, IL-7, and G-CSF (**Figure 1B**). PCA based on the 28 blood proteins in first samples revealed two distinct clusters (**Figure 1C**). The pro-inflammatory cytokines IFN-γ, TNF-α, IL-6, IL-9, IL-17α, CXCL10/IP-10, CCL4/MIP1β, CCL5/RANTES, basic-FGF, and GM-CSF (**Figures 1D, E**) were the strongest contributors to differentiate a cluster of 20/192 (10.4%) samples, which was hallmarked by positivity for HLA-DSA (15/20; 75%), whereas the majority of the other cluster were HLA-DSA negative (169/172; 98.3%) (**Figure 1C**). Unsupervised hierarchical clustering confirmed this differentiation in the levels of pro-inflammatory cytokines, corresponding to the HLA-DSA status, and the histological lesions scores (cluster I) (**Figure 1F**). In 14/15 HLA-DSA+ patients in cluster I, HLA-DSA were present at time of the biopsy (12 persistent pretransplant, 1 *de novo*, and 1 both pretransplant and *de novo* HLA-DSA), whereas in 1 patient

the HLA-DSA were present pretransplant but had resolved at the first biopsy. In the other five patients (25%) of cluster I, no HLA-DSA was detected at the time of transplantation or at any time post-transplantation. No differences in specificity (class) or strength (MFI) of the HLA-DSA were noted between the HLA-DSA+ samples in cluster I (N = 15 HLA-DSA) vs. cluster II (N = 3 HLA-DSA) (**Supplementary Table S3**). When comparing patients in cluster I vs. II, differences were noted in HLA-DSA, repeat transplantation, induction therapy, and recipient sex, all reflecting the background risk and management of HLA-DSA, but no differences in C-reactive protein or polyomavirus or cytomegalovirus viremia (**Table 2**). When repeating the analyses in all samples (N = 293), including the repeat sample-biopsy pairs, a similar cluster of HLA-DSA+ samples resulted from the PCA and heatmap analyses (**Supplementary Figure S1**). Follow-up samples from patients whose first biopsy was clustered in cluster I clustered together in the pro-inflammatory cluster.

As the clustering approach indicated a strong association between the global protein levels and HLA-DSA, we next evaluated which proteins individually associated with HLA-DSA using logistic mixed models. All proteins except IL-4, IL-5, IL-10, IL-12p70, IL-13, CCL11/Eotaxin, and VEGF significantly associated with HLA-DSA (**Figure 2A**). A high number of samples had levels below the detection limit in both clusters for IL-5, IL-10, IL-12p70, IL-13, and VEGF (**Supplementary Table S1** and **Supplementary Figure S2**).

### Cytokine Profiles in Relation to Histological Inflammation

Next to the association with HLA-DSA, also rejection associated with the cytokine clusters (**Table 2**). Rejection was present in the concomitant biopsy in 14/20 (70%) patients of cluster I vs. 58/172 (33.7%) of cluster II (p < 0.0001). ABMR (including mixed rejection) was present in 9/20 (45%) cases in cluster I vs. 6/172 (3.5%) cases in cluster II (p < 0.0001). Glomerulitis, peritubular capillaritis, thrombi, C4d deposition in the peritubular capillaries, and interstitial inflammation scores differed between the clusters (**Table 2**). In non-rejecting samples in cluster I (N = 6/20), the biopsies showed no inflammatory lesions, except for isolated tubulitis in three of these biopsies.

When considering all biopsies, significant associations were found between cytokines and diagnosis of ABMR (N = 26) vs. no ABMR (N = 267) (**Figure 2C** and **Supplementary Figure S3**). However, after adjustment for HLA-DSA, none of the cytokines remained significantly associated with ABMR (**Figure 2D**). Conversely, the associations between HLA-DSA and the individual proteins remained largely unaltered after adjustment for ABMR (**Figure 2B**). This illustrates that HLA-DSA status is a stronger determinant of the cytokine profiles than ABMR histology. Compared to ABMR, fewer cytokines associated with TCMR (N = 71) vs. no TCMR (N = 222), and none associated with borderline changes vs. no rejection (**Supplementary Figure S3**). When considering pure TCMR (N = 54) vs. no TCMR (N = 223), none of the cytokines remained associated after FDR correction (**Supplementary**
<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Visualization of the structural relationships between blood protein levels and clinical and histological characteristics. (A) Study design. (B) Correlation matrix of the 28 proteins using Pearson correlations of the log10-transformed protein levels. Colors indicate correlation coefficient (ρ); dots indicate p-values (only p-values <0.05 are represented with a circle). Proteins are ordered as defined by hierarchical clustering. (C) Principal component analysis (PCA) of first biopsies per patient (N = 192) shows two clusters distinct in their 28 blood protein levels. Colors indicate the presence of HLA-DSA. (D) Contributions of the 28 cytokines to the principal component analysis. (E) Top 10 contributing cytokines to axis 1 of the principal component analysis (PCA1). (F) Heatmap analysis of histological lesions and blood proteins of first biopsies per patient (N = 192). Reordering of dendrograms based on hierarchical clustering. Two distinct clusters, cluster I (N = 20) and II (N = 172) can be distinguished. (G) Kaplan-Meier survival curve illustrating survival probability for cluster I and II, counted from the day of the first biopsy (N = 192). ABMR, antibody-mediated rejection; TCMR, T cell-mediated rejection; HLA-DSA, anti-HLA donor-specific antibodies; HLA-Abs, anti-HLA antibodies. Mixed rejection was defined as ABMR concomitant with TCMR or borderline changes. All protein levels are log10 transformed. Humoral score = sum of glomerulitis; peritubular capillaritis, intimal arteritis and C4d deposition in the peritubular capillaries. Cellular score = sum of tubulitis, interstitial inflammation, and intimal arteritis.</figcaption></figure>

**Figure S3**). These findings were corroborated by the associations between protein levels and inflammatory histological lesions. Basic-FGF, IL-1β, IL-9, IL-2, IL-17α, CCL4/MIP1β, CCL5/RANTES, CXCL9/MIG, IFN-γ, and TNF-α positively and IL-7 negatively associated with microvascular inflammation, the hallmark of ABMR (**Supplementary Figure S4**). Proteins with a positive association with tubulo-interstitial inflammation, the hallmark of TCMR, were CXCL9/MIG and CXCL10/IP10 (**Supplementary Figure S4**). When considering intimal arteritis, which can occur both in ABMR and in TCMR,

upregulated proteins were IL-6, CXCL8/IL-8, CXCL9/MIG, CXCL10/IP10, IFN-γ, and TNF-α; downregulated proteins were IL-7, IL-13, and G-CSF (**Supplementary Figure S4**).

## Association of Blood Cytokines With Future Rejection and Graft Survival

Of the 14 patients in cluster I with rejection in the first biopsy, 12 (85.7%) had recurring rejections in the follow-up biopsies. Of the 6 patients in cluster I without rejection in the first biopsy, 2 patients developed transplant glomerulopathy in follow-up
<table><caption>**TABLE 1** | Sample characteristics of first and all biopsies.</caption><thead><tr><th></th><th>First biopsies (N = 192)</th><th>Mean ± std</th><th>All biopsies (N = 293)</th></tr><tr><th></th><th></th><th>Or median (IQR)</th><th>Or no (%)</th></tr></thead><tbody><tr><td><strong>Clinical parameters</strong></td><td></td><td></td><td></td></tr><tr><td>Serum C reactive protein, mg/L</td><td>9.0 (1.8-24.7)</td><td></td><td>5.3 (1.4-18.3)</td></tr><tr><td>eGFR, MDRD, mL/min</td><td>19.6 ± 13.1</td><td></td><td>19.8 ± 12.2</td></tr><tr><td>Proteinuria, g/g creatinine</td><td>0.3 (0.2-1.1)</td><td></td><td>0.3 (0.2-0.8)</td></tr><tr><td>Time after transplantation (days)</td><td>11 (7-24)</td><td></td><td>18 (9-51)</td></tr><tr><td>HLA-DSA at time of biopsy</td><td>17 (8.9%)</td><td></td><td>29 (9.9%)</td></tr><tr><td>HLA-DSA at time of transplant</td><td>17 (8.9%)</td><td></td><td>34 (11.6%)</td></tr><tr><td><em>De novo</em> HLA-DSA</td><td>2 (1.0%)</td><td></td><td>6 (2.1%)</td></tr><tr><td><strong>Histological diagnosis</strong></td><td></td><td></td><td></td></tr><tr><td>Any rejection</td><td>72 (37.5%)</td><td></td><td>106 (36.2%)</td></tr><tr><td>- ABMR (pure)</td><td>5 (2.6%)</td><td></td><td>7 (2.4%)</td></tr><tr><td>- TCMR (pure)</td><td>35 (18.2%)</td><td></td><td>54 (18.4%)</td></tr><tr><td>- Borderline changes (pure)</td><td>22 (11.5%)</td><td></td><td>26 (8.9%)</td></tr><tr><td>- Mixed rejection</td><td>10 (5.2%)</td><td></td><td>19 (6.5%)</td></tr><tr><td><strong>Histological lesion scores</strong></td><td></td><td></td><td></td></tr><tr><td>Glomerulitis &gt; 0</td><td>33 (17.2%)</td><td></td><td>53 (18.1%)</td></tr><tr><td>Peritubular capillaritis &gt; 0<sup>a</sup></td><td>19 (9.9%)</td><td></td><td>37 (12.6%)</td></tr><tr><td>Tubulitis &gt;0</td><td>128 (66.7%)</td><td></td><td>200 (68.3%)</td></tr><tr><td>Interstitial inflammation &gt;0</td><td>51 (26.6%)</td><td></td><td>73 (24.9%)</td></tr><tr><td>Intimal arteritis &gt;0</td><td>38 (19.8%)</td><td></td><td>62 (21.2%)</td></tr><tr><td>C4d ptc &gt;1<sup>a</sup></td><td>17 (8.9%)</td><td></td><td>24 (8.2%)</td></tr><tr><td>Trombi (yes)</td><td>10 (5.2%)</td><td></td><td>14 (4.8%)</td></tr><tr><td>ABMRh (yes)</td><td>23 (12.0%)</td><td></td><td>41 (14.0%)</td></tr><tr><td>Arteriolar hyalinosis &gt;1</td><td>25 (13.0%)</td><td></td><td>44 (15.0%)</td></tr><tr><td>Vascular intimal thickening &gt;1</td><td>40 (20.8%)</td><td></td><td>70 (23.9%)</td></tr><tr><td>Interstitial fibrosis &gt;1</td><td>7 (3.7%)</td><td></td><td>25 (8.5%)</td></tr><tr><td>Tubular atrophy &gt;1</td><td>4 (2.1%)</td><td></td><td>11 (3.8%)</td></tr><tr><td>Transplant glomerulopathy &gt;0</td><td>0</td><td></td><td>3 (1.0%)</td></tr><tr><td>Mesangial matrix expansion &gt;0</td><td>10 (5.2%)</td><td></td><td>19 (6.5%)</td></tr></tbody></table>

biopsies without previous acute rejections detected on pathology; 2 patients developed rejection (one TCMR, one borderline); and 2 never developed rejection. In total, 25/192 patients experienced death-censored graft failure during a mean follow-up time of 5.1 (± 1.9) years. Graft failure rate differed between cluster I and cluster II (univariable hazard ratio [HR] 2.81, 95% CI 1.12–7.06, p = 0.03) (Figure 1G). The association of cluster I with graft failure was independent of recipient sex, repeat transplantation, induction therapy, and ABMRh (adjusted HR 3.31, 95% CI 1.09–10.03, p = 0.03).

## Cellular Origin of the Cytokine Profiles: Single-Cell RNA Sequencing of Kidney Transplant Biopsies

Based on the observed association between cytokines and HLA-DSA, we hypothesized that these cytokines originated from the transplanted kidneys as the allograft is the site where the DSA-mediated injury occurs. To determine the cellular origin of the cytokine profiles in the allograft, we used publicly available single-cell data from kidney biopsies. Briefly, scRNASeq was performed on two renal allograft biopsies with DSA-mediated injury, namely, ABMR (one with class I and II anti-HLA-DSA (B8, DQA1) and one with class II HLA-DSA [DR1; DR53 and DQ2 (36)]), and five reference biopsies. After quality control and filtering

(Supplementary Figure S5), 33,216 cells were detected and clustering revealed 14 clusters corresponding to the main renal cell subtypes but also infiltrating immune cells (Figures 3A, B). After subclustering using well-established markers (39), we evaluated the expression of genes corresponding to the cytokines across the kidney structural cells and infiltrated immune cells (Figures 3A, B). The cytokines of interest were expressed by activated endothelial cells (ECa), antigen-presenting cells (APCs), and lymphocytes. More specifically, ECa expressed CSF3, FGF2, IL15, CXCL9, and CCL2. APCs expressed IL1RN, IL8, CCL3, IL1B, CXCL10, IL10, TNF, and CCL3 whereas lymphocytes expressed CCL4, CSF2, IL2, IFNG, IL4, IL13, and CCL5 (Figure 3C). Of note, IL9 and IL17A were not detected in this dataset. Given the heterogeneity of APC and lymphocyte populations, we subclustered these two populations. In APCs, we identified CD19+MS4A1+ B cells, the three main monocytic populations—classical CD14+CD16-, intermediate CD14+CD16+, and non-classical CD14-CD16+ monocytes—and two dendritic cell (DC) populations (Figure 3D). Non-classical monocytes expressed the most relevant cytokines with IL1RN, CCL3, IL1B, CXCL10, TNF, IL15, and CCL4. In lymphocytes, we detected CD4 T cells, CD8 T cells, NKT cells, NK cells, and plasmacytoid DC, whereas CD4 and CD8 T cells expressed IL4, IL2, and IL13, NK and NKT cells expressed IFNG, CCL3, CCL4, and CCL5 (Figure 3E).

<footnote><sup>a</sup>Missing data for peritubular capillaritis (N = 5 in first biopsies, and N = 6 in all biopsies), C4d deposition in ptc (N = 6 in first biopsies, N = 9 in all biopsies).</footnote>
<table><caption>**TABLE 2** | Comparison of clinical characteristics between cluster I (N = 20) and cluster II (N = 172).</caption><thead><tr><th></th><th>Cluster I</th><th>Cluster II</th><th>p-value</th></tr></thead><tbody><tr><td>Recipient age (years)</td><td>52.0 ± 11.7</td><td>55.9 ± 11.6</td><td>0.16</td></tr><tr><td>Donor age (years)</td><td>47.2 ± 14.8</td><td>54.3 ± 15.2</td><td>0.05</td></tr><tr><td>Recipient sex (male)</td><td>7 (35.0%)</td><td>126 (73.3%)</td><td><strong>0.0004</strong></td></tr><tr><td>Donor sex (male)<sup>a</sup></td><td>12 (60.0%)</td><td>98 (58.0%)</td><td>0.86</td></tr><tr><td>Cold ischemia time (hours)<sup>a</sup></td><td>11.2 ± 5.8</td><td>11.9 ± 5.5</td><td>0.59</td></tr><tr><td>Type of donor</td><td></td><td></td><td>0.55</td></tr><tr><td>- Circulatory death</td><td>3 (15.0%)</td><td>39 (22.7%)</td><td></td></tr><tr><td>- Brain death</td><td>15 (75.0%)</td><td>124 (72.1%)</td><td></td></tr><tr><td>- Living</td><td>2 (10.0%)</td><td>9 (5.2%)</td><td></td></tr><tr><td>Repeat transplantation (yes)</td><td>6 (30.0%)</td><td>17 (9.9%)</td><td><strong>0.009</strong></td></tr><tr><td>Immunosuppressive regimen (TAC-MMF-CS)</td><td>20 (100%)</td><td>155 (90.1%)</td><td>0.14</td></tr><tr><td>Induction therapy</td><td>12 (60.0%)</td><td>54 (31.4%)</td><td><strong>0.01</strong></td></tr><tr><td>Polyoma viremia (positive)<sup>a</sup></td><td>0 (0%)</td><td>15 (8.7%)</td><td>0.20</td></tr><tr><td>CMV viremia (positive)<sup>a</sup></td><td>0 (0%)</td><td>9 (5.2%)</td><td>0.34</td></tr><tr><td>Serum C reactive protein, mg/L</td><td>8.7 (1.8–25.3)</td><td>10.3 (2.9–21.1)</td><td>0.66</td></tr><tr><td>eGFR, MDRD, mL/min</td><td>17.7 ± 11.9</td><td>19.8 ± 13.2</td><td>0.49</td></tr><tr><td>Proteinuria, g/g creatinine</td><td>0.3 (0.2–1.1)</td><td>0.5 (0.2–1.2)</td><td>0.30</td></tr><tr><td>Time after transplantation (days)</td><td>11 (7–24)</td><td>10 (7–23)</td><td>0.53</td></tr><tr><td>HLA-DSA at time of biopsy</td><td>14 (70.0%)</td><td>3 (1.7%)</td><td><strong>&lt;0.0001</strong></td></tr><tr><td>HLA-DSA at time of transplant</td><td>14 (70.0%)</td><td>3 (1.7%)</td><td><strong>&lt;0.0001</strong></td></tr><tr><td>De novo HLA-DSA</td><td>2 (10.0%)</td><td>0 (0%)</td><td><strong>&lt;0.0001</strong></td></tr><tr><td>Any rejection</td><td>14 (70.0%)</td><td>58 (33.7%)</td><td><strong>&lt;0.0001</strong></td></tr><tr><td>- Pure ABMR</td><td>1 (5.0%)</td><td>4 (2.3%)</td><td></td></tr><tr><td>- Pure TCMR</td><td>2 (10.0%)</td><td>33 (19.2%)</td><td></td></tr><tr><td>- Pure borderline changes</td><td>3 (15.0%)</td><td>19 (11.1%)</td><td></td></tr><tr><td>- Mixed rejection</td><td>8 (40.0%)</td><td>2 (1.2%)</td><td></td></tr><tr><td>Acute histological lesions</td><td></td><td></td><td></td></tr><tr><td>- Glomerulitis &gt; 0</td><td>10 (50.0%)</td><td>23 (13.4%)</td><td><strong>&lt;0.0001</strong></td></tr><tr><td>- Peritubular capillaritis &gt; 0<sup>a</sup></td><td>6 (30.0%)</td><td>13 (7.8%)</td><td><strong>0.002</strong></td></tr><tr><td>- Tubulitis &gt;0</td><td>16 (80.0%)</td><td>112 (65.1%)</td><td>0.18</td></tr><tr><td>- Interstitial inflammation &gt;0</td><td>11 (55.0%)</td><td>40 (23.3%)</td><td><strong>0.002</strong></td></tr><tr><td>- Intimal arteritis &gt;0</td><td>7 (35.0%)</td><td>31 (18.0%)</td><td>0.07</td></tr><tr><td>- C4d ptc &gt;1<sup>a</sup></td><td>8 (40.0%)</td><td>9 (5.4%)</td><td><strong>&lt;0.0001</strong></td></tr><tr><td>- Thrombi (yes)</td><td>4 (20.0%)</td><td>6 (3.5%)</td><td><strong>0.002</strong></td></tr><tr><td>- ABMRh (yes)</td><td>9 (45.0%)</td><td>14 (8.1%)</td><td><strong>&lt;0.0001</strong></td></tr></tbody></table>

TAC, tacrolimus; MMF, mycophenolate mofetil; CS, corticosteroids.

Significant differences are indicated in bold.

## Cell-Cell Communication

As cytokines and chemokines are the mediators of cell-cell communication, we used the CellChat package to infer the intercellular communication from the scRNASeq data. Twenty signaling pathways showed significant communication between the allograft cells. Among these were the CXCL, CCL, CX3C, PDGF, VEGF, and CSF signaling pathways (Figure 4A). ECa were important senders but also receivers of CXCL and CCL signaling (Figure 4B). Lymphocytes were important senders of CCL signaling, whereas APCs were receiving CX3C and CSF signals from ECa in these biopsies. When looking at all signals received by ECa, sent by ECa, APCs, and lymphocytes, the atypical chemokine receptor (ACKR1), also known as DARC, was the identified receptor for all ligand-receptor pairs (Figure 4C).

## Cytokine Secretion in Antibody-Mediated NK Cell and Monocyte Activation

Given that the most clinically relevant cytokines corresponded to expressed genes on activated glomerular endothelial cells, non-classical monocytes, and NK cells, we next assessed the cytokine

production by these cells. Indeed, endothelial injury is the hallmark of DSA-mediated injury in kidney allografts and both NK cells and monocytes are key players in ABMR (40–45). Both NK cells and non-classical monocytes express the Fcγ receptor, through which antibody-dependent cellular cytotoxicity can be induced against a target cell, i.e., the endothelial cell in the setting of ABMR. We tested the specificity of the studied proteins for these cell types in an *in vitro* model of antibody-mediated NK cell and/or monocyte activation in interaction with glomerular endothelial cells. After 4 h of coculture with conditionally immortalized glomerular endothelial cells (ciGENC), IFN-γ, TNF-α, IL-6, IL-8, IL-9, IL-15, CCL2/MCP-1, CCL3/MIP1α, CCL4/MIP1β, CCL5/RANTES, CCL11/Eotaxin, and basic-FGF were significantly higher in the supernatants of the antibody-mediated NK cell activation condition, compared to the control condition (Figure 5A and Supplementary Figure S6). After 24 h of coculture with primary GENC, IFN-γ, IL-4, IL-6, CXCL8/IL-8, CXCL10/IP10, CCL2/MCP-1, CCL3/MIP1α, CCL4/MIP1β, and basic-FGF were significantly higher in the supernatants of the antibody-mediated monocyte activation condition, compared to

<footnote><sup>a</sup>Missing data: donor sex N = 3; cold ischemia time N = 3; peritubular capillaritis N = 5; C4d deposition ptc N = 6; polyoma viremia N = 6; CMV viremia N = 40.</footnote>
<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /><figcaption>**FIGURE 2** | Associations of the cytokines with HLA-DSA and antibody-mediated rejection in all samples (N = 293). (A) Forest plots of odds ratios for association of the cytokines with HLA-DSA (N = 38 vs. 255), (B) with HLA-DSA when adjusted for presence of ABMR, (C) with ABMR vs. no ABMR (N = 26 vs. N =267), and (D) with ABMR when adjusted for presence of HLA-DSA. Blue color indicates associations that were significant after false discovery rate correction (FDR p-value <0.05). ABMR, antibody-mediated rejection incl. mixed rejection; HLA-DSA, anti-HLA donor-specific antibodies; FDR, false discovery rate.</figcaption></figure>

control (Figure 5B and Supplementary Figure S7). After 24 h of coculture of NK cells, non-classical monocytes and GENC, IFN-γ, IL-5, CXCL8/IL-8, IL-9, IL-15, CXCL10/IP10, CCL2/MCP-1, CCL4/MIP1β, and GM-CSF were significantly higher in the NK cell and monocyte activation condition, compared to control, suggesting a synergistic effect between these two cell types (Figure 5B and Supplementary Figure S7). When comparing *in vitro* models of GENC alone, with an HLA class I-specific antibody vs. non-specific isotype antibody, no significant cytokine production was noted (Supplementary Figure S8). The cell–cell contact of immune cells with GENC was required for cytokine production, illustrated by the higher cytokine levels in cocultures of immune cells with GENC compared to immune cells alone (Supplementary Figure S9).

DISCUSSION

In this representative cohort of consecutive samples from kidney
transplant recipients at time of allograft dysfunction, we
demonstrate that pro-inflammatory cytokines in the peripheral
blood associate with presence of HLA-DSA and impaired kidney
transplant survival. These pro-inflammatory blood signals associate
with HLA-DSA, rather than with rejection phenotypes as identified
in kidney transplant biopsies. This illustrates that HLA-DSA can
lead to immune activation that is not always reflected in histology

of the concomitant kidney transplant biopsies. Single-cell data
analysis demonstrated that the most contributing cytokines to the
clustering of the clinical samples corresponded to genes expressed
on activated glomerular endothelial cells, non-classical monocytes,
and NK cells. Subsequent *in vitro* models confirmed that the
cytokines observed in patients’ serum can be induced by HLA-
DSA-mediated cross talk between glomerular endothelial cells, NK
cells, and monocytes.

To our knowledge, this is the first study demonstrating the
strong association of pro-inflammatory blood cytokine profiles
with presence of HLA-DSA. One previous study found a panel of
6 pro-inflammatory cytokines (TNF-α, IFN-γ, IL-1β, IL-6, IL-8/
CXCL8, and IL-17) to be associated with angiotensin II type 1
receptor antibodies, a non-HLA antibody, but not with presence
of HLA-DSA (13). Most of the cytokines that associated with
HLA-DSA in our study are known as potent mediators of the
alloimmune responses and associated injury. The function of
cytokines and chemokines is characterized by pleiotropism,
redundancy, synergism, and antagonism (46), clearly illustrated
by the strong correlations we observed between the proteins. The
signals that stand out primarily reflect IFN-γ-related
inflammation including CXCR3 and CCR5 signaling. Levels of
CXCL9/MIG and CXCL10/IP10, ligands to CXCR3, and levels of
the MIP-1 family proteins (especially CCL4/MIP-1β and CCL5/
RANTES), ligands to CCR5, are strongly associated with
diagnosis of presence of HLA-DSA. The importance of these
<figure><img src="image_8.png" /></figure>

receptors in rejection has been demonstrated in literature, as well as their potential as therapeutic targets (16, 17, 27, 47–52). The presence of both CXCR3 and CCR5 appears to identify subsets of NK cells and T cells in blood with a strong predilection for homing to inflammatory sites (53, 54), and CXCL9/MIG and CXCL10/IP10 have been well characterized in kidney allograft rejection (20, 55).

Next, the scRNASeq analysis on kidney transplant biopsies suggested that, although the different cytokines are highly collinear, they do not originate from the same cells. CXCL9/MIG is mainly expressed by activated endothelial cells, whereas CXCL10/IP10 is expressed by both activated endothelial cells and CD16+ monocytes. These results are in line with reports showing that CD16+ non-classical
<figure><img src="image_9.png" /><figcaption>**FIGURE 4** | Cell–cell communications derived from scRNASeq data of kidney allograft biopsies. (A) Cytokine and chemokine signaling pathways that were significant communication pathways in the kidney allograft cells. (B) Most important senders, receivers, mediators, and influencers of the CXCL, CCL, and CX3C signaling network. (C) Circos plot demonstrating all signals received by activated endothelial cells (ECa), sent by activated endothelial cells, antigen-presenting cells (APCs), and lymphocytes, illustrating the ligand–receptor pairs.</figcaption></figure>

monocytes present an increased capacity for antibody-
dependent cellular cytotoxicity (56). Other relevant identified
cytokines are the pro-inflammatory cytokines IL-6, IL-8, IL-9,
and IL-17α.

As the scRNASeq analysis suggested expression of the
relevant pro-inflammatory cytokine profiles in kidney
endothelial cells, monocytes, and NK cells, we next performed
*in vitro* studies to evaluate whether these cytokines were indeed
<figure><img src="image_10.png" /><figcaption>**FIGURE 5** | Cytokine profiles in supernatants of NK cells and/or monocytes cocultured with HLA class I antibody-activated glomerular endothelial cells or isotype antibodies. (A) NK cell experiment using NK cells from 7 human donors, cocultured with human antibody serum (anti-HLA-A2-containing serum or control pooled serum that was negative for anti-HLA A2) and ciGENC (immortalized glomerular endothelial cells), harvested after 4 h. (B) NK cell and monocyte experiment using NK cells and CD16+ monocytes from 5 human donors, cocultured with murine antibody (anti-HLA-A, -B, -C, or control isotype) and GENC (primary glomerular endothelial cells, Cell Systems, USA); harvesting of supernatants was done after 24 h of coculture. Violin plots of the cytokine concentrations (pg/mL) of the significantly differently expressed cytokines in both conditions, as assessed by a paired parametric T-test of the log10 transformed concentrations.</figcaption></figure>

produced by these cells upon the presence of HLA-DSA. The strongest increases in cytokine levels were observed when combining both NK cells and monocytes with the glomerular endothelial cells, suggesting a synergistic effect in their cytokine secretion. These observations suggest that cytokine production in the presence of HLA-DSA is FcyR-mediated; however, they do not provide evidence of causality. Future experiments using FcyR blocking or using only the F(Ab')2 fragments are needed to prove causality. In a previous study by Wei et al. (57) studying

monocytes and DSA-mediated injury of endothelial cells, using F(ab')2 alone induced only upregulation of the antigens on the endothelial cells, leading to increased tethering and adhesion of monocytes. Only in the presence of the Fc region was the extra signal for activation of the monocytes provided, inducing cytokine secretion, suggesting that FcyR-mediated signaling is essential for the cytokine production in DSA-mediated injury.

Since the endothelial cells present in the *in vitro* conditions could contribute to the cytokine profiles observed in these
supernatants (58, 59), supernatants of glomerular endothelial cells alone in the presence of a specific anti-HLA antibody vs. isotype antibody were also analyzed. No significant differences in cytokine secretion were noted between the two conditions. Together with the single-cell data, this suggests that endothelial cells express cytokines, but production is only observed when interacting with immune cells. For this interaction with immune cells through chemokines, the ACKR1 receptor, or DARC, was identified as a common receptor on activated endothelial cells for multiple incoming chemokine signals of the CXCL and CCL families. This receptor is known to regulate the dynamics of inflammatory chemokine presentation on vascular endothelial cells and chemokine concentrations in plasma (60). The importance of this receptor in ABMR has been reported previously (28, 61, 62).

Although older studies did demonstrate that HLA class I antibodies can directly trigger cytokine secretion from endothelial cells (15, 63), our data are in accordance with the study of Wei et al., where HLA class I antibodies did not induce cytokine production directly from endothelial cells derived from aortic rings, but cell-cell contact of monocytes with antibody-activated endothelial cells was required for cytokine production (57). The observation that this is also true for glomerular endothelial cells used in our experiments, which are more relevant for kidney microvascular injury, is novel. The discrepancy with previous studies (15, 63) could be explained by differences in vascular origin of the endothelial cells or in the type/concentration of the HLA antibodies.

Our present work challenges the vision that kidney transplant histology is the gold standard for identification of ongoing alloimmune processes. The immune activation detectable in the peripheral blood of patients with HLA-DSA could be an explanation for the bad outcome observed in patients with HLA-DSA, even in the absence of histological lesions in the biopsy (35). This could be related to the observer dependency of kidney transplant biopsy reading and the risk for sampling error, in addition to the potential that classic light microscopy is not sensitive enough for detecting kidney transplant inflammation. As illustrated previously (64), increased intrarenal rejection-associated gene transcripts can be observed in biopsy samples of DSA-positive/ABMR-negative patients, potentially proceeding to later overt histological injury (65). Given that the pathological evaluation of biopsies is observer-dependent and subject to sampling error, inflammatory lesions could be missed or remain below the thresholds of the Banff classification, yet be reflected by inflammatory changes detectable in the graft/circulation. Markers for inflammatory changes in the graft below the threshold of the Banff classification have been suggested before, by studying intragraft transcriptional changes ("molecular microscope") (66) or donor-derived cell-free DNA (67) as markers for severity of graft injury to complement the histological assessment using Banff classification. Our study illustrates that also blood inflammatory cytokine profiles could help in objectifying deleterious inflammation in patients at high risk of rejection but with negative histology. This is further supported by our finding that patients with HLA-DSA and enhanced pro-inflammatory signals had worse outcome than

HLA-DSA-positive patients without pro-inflammatory signals in the blood. All of this seems to be in contrast to the study of Parajuli et al., who did not observe different outcomes in patients with and without HLA-DSA if the index biopsy was negative for histological inflammation (68). This previous study was however focused on first biopsies in the absence of rejection, the study design is not comparable to our strategy with a relatively short follow-up time, and the interpretation of these data could be hampered by inherent selection bias.

There are some limitations to this study. The single-center nature of our study population, consisting largely of white Europeans with low to intermediate immunologic risk, treated mainly with a tacrolimus-based regimen, limits the generalizability of these findings to other and higher-risk populations and warrants external validation. Second, this study consisted of biopsies performed at time of graft dysfunction with possible confounders and inherently excluding cases with subclinical rejection. All samples were taken in the first year after transplantation, and the majority in the first 3 months after transplantation, reflecting the real-life timing of the majority of indication biopsies but obviating generalization of the results to biopsies later after transplantation and chronic rejection phenotypes. Induction therapy was different between the two clusters, with more induction therapy given to sensitized patients, and could have confounded cytokine levels early after transplantation. However, we would expect the induction therapy to cause a lesser immune activation, whereas here the opposite was observed. T cell sensitization was not measured, and its effect could therefore not be explored in this study. Also, we used a preselected multiplex assay of cytokines, chemokines, and growth factors, but other potentially important proteins were not studied. Downstream analyses were confined to the main clinical association, namely, DSA-mediated injury. This implies that other rejection types that could also contribute to cytokine expression and secretion were not studied in the downstream analyses. Nevertheless, we did not find indications for their contribution in the clinical associations. With regard to the *in vitro* analysis, we chose to focus on NK cells and non-classical monocytes, while other cell types capable of antibody recognition like CD16+ CD8 T cells (69), NKT cells, γδT cells, neutrophils, and macrophages could also play a role in the chemokine profiles associated with DSA. Moreover, we did not study the effect of complement, which can cause endothelial injury from antibodies independent of immune cells (70). We also did not study the effect of missing self, which could enhance NK cell activation in addition to DSA (43, 71). Finally, the *in vitro* experiments in this study are not adequate to define the exact cellular source of the secreted cytokines.

This study did not intend to explore the biomarker potential of these cytokines. Although biomarker potential is not excluded, different study design and analyses would be needed to assess this. Furthermore, one could question the value of peripheral blood cytokines as clinical biomarkers as we can easily measure the gold standard, which are the circulating HLA-DSA.

In conclusion, this study demonstrates that blood pro-inflammatory cytokines are increased in kidney transplant
patients with HLA-DSA and associate with worse graft survival, even in the absence of histology of rejection. These observations illustrate that HLA-DSA can lead to immune activation that is not always reflected in the histology of the concomitant kidney biopsy, challenging the concept that histology is the gold standard for identification of ongoing allo-immune activation after transplantation.

## DATA AVAILABILITY STATEMENT

All of the individual de-identified participant data that underlie the results reported in this article (text, tables, figures, and appendices) can be made available on a collaborative basis following institutional review board approval. Proposals should be directed at the corresponding author.

## ETHICS STATEMENT

The studies involving human participants were reviewed and approved by the local ethical committee (S53364 and S61971). The patients/participants provided their written informed consent to participate in this study.

## AUTHOR CONTRIBUTIONS

EVL and MN conceived and designed the study. EVL, MC, HdL, M-PE, PK, DK, AS, BS, AVC, and MN collected the clinical data

## REFERENCES

1. Hart A, Smith JM, Skeans MA, Gustafson SK, Stewart DE, Cherikh WS, et al. Opttn/Srtr 2015 Annual Data Report: Kidney. *Am J Transplant* (2017) 17 Suppl 1:21–116. doi: 10.1111/ajt.14124

2. Coemans M, Süsal C, Döhler B, Anglicheau D, Giral M, Bestard O, et al. Analyses of the Short- and Long-Term Graft Survival After Kidney Transplantation in Europe Between 1986 and 2015. *Kidney Int* (2018) 94 (5):964–73. doi: 10.1016/j.kint.2018.05.018

3. Loupy A, Hill GS, Jordan SC. The Impact of Donor-Specific Anti-HLA Antibodies on Late Kidney Allograft Failure. *Nat Rev Nephrol* (2012) 8 (6):348–57. doi: 10.1038/nrneph.2012.81

4. Loupy A, Lefaucheur C. Antibody-Mediated Rejection of Solid-Organ Allografts. *N Engl J Med* (2018) 379(12):1150–60. doi: 10.1056/NEJMra1802677

5. Reindl-Schwaighofer R, Heinzel A, Kainz A, van Setten J, Jelencsics K, Hu K, et al. Contribution of Non-HLA Incompatibility Between Donor and Recipient to Kidney Allograft Survival: Genome-Wide Analysis in a Prospective Cohort. *Lancet* (2019) 393(10174):910–7. doi: 10.1016/S0140-6736(18)32473-5

6. Steers NJ, Li Y, Drace Z, D'Addario JA, Fischman C, Liu L, et al. Genomic Mismatch at. *N Engl J Med* (2019) 380(20):1918–28. doi: 10.1056/NEJMoa1803731

7. El-Zoghby ZM, Stegall MD, Lager DJ, Kremers WK, Amer H, Gloor JM, et al. Identifying Specific Causes of Kidney Allograft Loss. *Am J Transplant* (2009) 9(3):527–35. doi: 10.1111/j.1600-6143.2008.02519.x

8. Sellarés J, de Freitas DG, Mengel M, Reeve J, Einecke G, Sis B, et al. Understanding the Causes of Kidney Transplant Failure: The Dominant Role of Antibody-Mediated Rejection and Nonadherence. *Am J Transplant* (2012) 12(2):388–99. doi: 10.1111/j.1600-6143.2011.03840.x

and samples. SC performed the Luminex analyses. EVL and BL performed the scRNASeq analyses. TB, OT, and BL performed the *in vitro* experiments. EVL, BL, and MN did the statistical analyses and interpreted the data. EVL and MN wrote the draft of the report. All authors revised the report. All authors contributed to the article and approved the submitted version.

## FUNDING

EVL holds a fellowship grant (1143919N) from The Research Foundation Flanders (F.W.O.). MN is a senior clinical investigator of F.W.O. (1844019N) and is supported by an F.W.O TBM grant (grant no. T004417N), an F.W.O. junior project grant (grant no. G087620N), an ERACoSysMed H2020 grant (grant no. JTC2_29), and a C3 internal grant from the KU Leuven (grant no. C32/17/049).

## ACKNOWLEDGMENTS

We thank the clinicians and surgeons, nursing staff, and the patients of University Hospitals Leuven.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2022.818569/full#supplementary-material

9. Naesens M, Kuypers DR, De Vusser K, Evenepoel P, Claes K, Bammens B, et al. The Histology of Kidney Transplant Failure: A Long-Term Follow-Up Study. *Transplantation* (2014) 98(4):427–35. doi: 10.1097/TP.0000000000000183

10. Loupy A, Vernerey D, Tinel C, Aubert O, Duong van Huyen JP, Rabant M, et al. Subclinical Rejection Phenotypes at 1 Year Post-Transplant and Outcome of Kidney Allografts. *J Am Soc Nephrol* (2015) 26(7):1721–31. doi: 10.1681/ASN.2014040399

11. Roberts DM, Jiang SH, Chadban SJ. The Treatment of Acute Antibody-Mediated Rejection in Kidney Transplant Recipients-a Systematic Review. *Transplantation* (2012) 94(8):775–83. doi: 10.1097/TP.0b013e31825d1587

12. Wan SS, Ying TD, Wyburn K, Roberts DM, Wyld M, Chadban SJ. The Treatment of Antibody-Mediated Rejection in Kidney Transplantation: An Updated Systematic Review and Meta-Analysis. *Transplantation* (2018) 102 (4):557–68. doi: 10.1097/TP.00000000000002049

13. Pearl MH, Grotts J, Rossetti M, Zhang Q, Gjertson DW, Weng P, et al. Cytokine Profiles Associated With Angiotensin Ii Type 1 Receptor Antibodies. *Kidney Int Rep* (2019) 4(4):541–50. doi: 10.1016/j.ekir.2018.12.011

14. Muczynski KA, Leca N, Anderson AE, Kieran N, Anderson SK. Multicolor Flow Cytometry and Cytokine Analysis Provides Enhanced Information on Kidney Transplant Biopsies. *Kidney Int Rep* (2018) 3(4):956–69. doi: 10.1016/j.ekir.2018.02.012

15. Reyes-Vargas E, Pavlov IY, Martins TB, Schwartz JJ, Hill HR, Delgado JC. Binding of Anti-HLA Class I Antibody to Endothelial Cells Produce an Inflammatory Cytokine Secretory Pattern. *J Clin Lab Anal* (2009) 23(3):157–60. doi: 10.1002/jcla.20310

16. Schnickel GT, Bastani S, Hsieh GR, Shefizadeh A, Bhatia R, Fishbein MC, et al. Combined CXCR3/CCR5 Blockade Attenuates Acute and Chronic Rejection. *J Immunol* (2008) 180(7):4714–21. doi: 10.4049/jimmunol.180.7.4714
17. Akashi S, Sho M, Kashizuka H, Hamada K, Ikeda N, Kuzumoto Y, et al. A Novel Small-Molecule Compound Targeting CXCR5 and CXCR3 Prevents Acute and Chronic Allograft Rejection. *Transplantation* (2005) 80(3):378-84. doi: 10.1097/01.tp.0000166338.99933.e1

18. Salom RN, Maguire JA, Hancock WW. Endothelial Activation and Cytokine Expression in Human Acute Cardiac Allograft Rejection. *Pathology* (1998) 30 (1):24-9. doi: 10.1080/00313029800169625

19. Blydt-Hansen TD, Sharma A, Gibson IW, Wiebe C, Sharma AP, Langlois V, et al. Validity and Utility of Urinary CXCL10/Cr Immune Monitoring in Pediatric Kidney Transplant Recipients. *Am J Transplant* (2021) 21(4):1545-55. doi: 10.1111/ajt.16336

20. Tinel C, Devresse A, Vermorel A, Sauvaget V, Marx D, Avettand-Fenoel V, et al. Development and Validation of an Optimized Integrative Model Using Urinary Chemokines for Noninvasive Diagnosis of Acute Allograft Rejection. *Am J Transplant* (2020) 20(12):3462-76. doi: 10.1111/ajt.15959

21. Hricik DE, Nickerson P, Formica RN, Poggio ED, Rush D, Newell KA, et al. Multicenter Validation of Urinary CXCL9 as a Risk-Stratifying Biomarker for Kidney Transplant Injury. *Am J Transplant* (2013) 13(10):2634-44. doi: 10.1111/ajt.12426

22. Jackson JA, Kim EJ, Begley B, Cheeseman J, Harden T, Perez SD, et al. Urinary Chemokines CXCL9 and CXCL10 Are Noninvasive Markers of Renal Allograft Rejection and BK Viral Infection. *Am J Transplant* (2011) 11 (10):2228-34. doi: 10.1111/j.1600-6143.2011.03680.x

23. Rosenblum JM, Shimoda N, Schenk AD, Zhang H, Kish DD, Keslar K, et al. CXC Chemokine Ligand (CXCL) 9 and CXCL10 Are Antagonistic Costimulation Molecules During the Priming of Alloreactive T Cell Effectors. *J Immunol* (2010) 184(7):3450-60. doi: 10.4049/jimmunol.0903831

24. Schaub S, Nickerson P, Rush D, Mayr M, Hess C, Golian M, et al. Urinary CXCL9 and CXCL10 Levels Correlate With the Extent of Subclinical Tubulitis. *Am J Transplant* (2009) 9(6):1347-53. doi: 10.1111/j.1600-6143.2009.02645.x

25. Krüger B, Böger CA, Obed A, Farkas S, Hoffmann U, Banas B, et al. RANTES/ CCL5 Polymorphisms as a Risk Factor for Recurrent Acute Rejection. *Clin Transplant* (2007) 21(3):385-90. doi: 10.1111/j.1399-0012.2006.00654.x

26. Panzer U, Steinmetz OM, Reinking RR, Meyer TN, Fehr S, Schneider A, et al. Compartment-Specific Expression and Function of the Chemokine IP-10/ CXCL10 in a Model of Renal Endothelial Microvascular Injury. *J Am Soc Nephrol* (2006) 17(2):454-64. doi: 10.1681/ASN.2005040364

27. Hancock WW, Wang L, Ye Q, Han R, Lee I. Chemokines and Their Receptors as Markers of Allograft Rejection and Targets for Immunosuppression. *Curr Opin Immunol* (2003) 15(5):479-86. doi: 10.1016/S0952-7915(03)00103-1

28. Segerer S, Regele H, MacK M, Kain R, Cartron JP, Colin Y, et al. The Duffy Antigen Receptor for Chemokines Is Up-Regulated During Acute Renal Transplant Rejection and Crescentic Glomerulonephritis. *Kidney Int* (2000) 58(4):1546-56. doi: 10.1046/j.1523-1755.2000.00316.x

29. Nelson PJ, Krensky AM. Chemokines, Chemokine Receptors, and Allograft Rejection. *Immunity* (2001) 14(4):377-86. doi: 10.1016/S1074-7613(01) 00118-2

30. Haskell CA, Ribeiro S, Horuk R. Chemokines in Transplant Rejection. *Curr Opin Investig Drugs* (2002) 3(3):399-405.

31. Mühlbacher J, Doberer K, Kozakowski N, Regele H, Camovic S, Haindl S, et al. Non-Invasive Chemokine Detection: Improved Prediction of Antibody-Mediated Rejection in Donor-Specific Antibody-Positive Renal Allograft Recipients. *Front Med (Lausanne)* (2020) 7:114. doi: 10.3389/fmed.2020.000114

32. Roufosse C, Simmonds N, Clahsen-van Groningen M, Haas M, Henriksen KJ, Horsfield C, et al. A 2018 Reference Guide to the Banff Classification of Renal Allograft Pathology. *Transplantation* (2018) 102(11):1795-814. doi: 10.1097/ TP.00000000000002366

33. Senev A, Coemans M, Lerut E, Van Sandt V, Daniëls L, Kuypers D, et al. Histological Picture of Antibody-Mediated Rejection Without Donor-Specific Anti-HLA Antibodies: Clinical Presentation and Implications for Outcome. *Am J Transplant* (2019) 19(3):763-80. doi: 10.1111/ajt.15074

34. Loupy A, Haas M, Roufosse C, Naesens M, Adam B, Afrouzian M, et al. The Banff 2019 Kidney Meeting Report (I): Updates on and Clarification of Criteria for T Cell- and Antibody-Mediated Rejection. *Am J Transplant* (2020) 20(9):2318-31. doi: 10.1111/ajt.15898

35. Senev A, Lerut E, Van Sandt V, Coemans M, Callemeyn J, Sprangers B, et al. Specificity, Strength, and Evolution of Pretransplant Donor-Specific HLA

Antibodies Determine Outcome After Kidney Transplantation. *Am J Transplant* (2019) 19(11):3100-13. doi: 10.1016/j.humimm.2019.07.028

36. Malone AF, Wu H, Fronick C, Fulton R, Gaut JP, Humphreys BD. Harnessing Expressed Single Nucleotide Variation and Single Cell RNA Sequencing to Define Immune Cell Chimerism in the Rejecting Kidney Transplant. *J Am Soc Nephrol* (2020) 31(9):1977-86. doi: 10.1681/ASN.2020030326

37. Jin S, Guerrero-Juarez CF, Zhang L, Chang I, Ramos R, Kuan CH, et al. Inference and Analysis of Cell-Cell Communication Using Cellchat. *Nat Commun* (2021) 12(1):1088. doi: 10.1038/s41467-021-21246-9

38. Satchell SC, Tasman CH, Singh A, Ni L, Geelen J, von Ruhland CJ, et al. Conditionally Immortalized Human Glomerular Endothelial Cells Expressing Fenestrations in Response to VEGF. *Kidney Int* (2006) 69(9):1633-40. doi: 10.1038/s41467-019-10861-2

39. Lake BB, Chen S, Hoshi M, Plongthongkum N, Salamon D, Knoten A, et al. A Single-Nucleus RNA-Sequencing Pipeline to Decipher the Molecular Anatomy and Pathophysiology of Human Kidneys. *Nat Commun* (2019) 10 (1):2832. doi: 10.1038/s41467-019-10861-2

40. Charreau B. Cellular and Molecular Crosstalk of Graft Endothelial Cells During AMR: Effector Functions and Mechanisms. *Transplantation* (2021) 105(11):e156-67. doi: 10.1097/TP.00000000000003741

41. Resch T, Fabritius C, Ebner S, Ritschl P, Kotsch K. The Role of Natural Killer Cells in Cellulor. *Transplantation* (2015) 99(7):1335-40. doi: 10.1097/TP.0000000000000757

42. Yazdani S, Callemeyn J, Gazut S, Lerut E, de Loor H, Wevers M, et al. Natural Killer Cell Infiltration Is Discriminative for Antibody-Mediated Rejection and Predicts Outcome After Kidney Transplantation. *Kidney Int* (2019) 95 (1):188-98. doi: 10.1016/j.kint.2018.08.027

43. Koenig A, Chen CC, Marçais A, Barba T, Mathias V, Sicard A, et al. Missing Self Triggers NK Cell-Mediated Chronic Vascular Rejection of Solid Organ Transplants. *Nat Commun* (2019) 10(1):5350. doi: 10.1038/s41467-019-13113-5

44. Hirohashi T, Chase CM, Della Pelle P, Sebastian D, Alessandrini A, Madsen JC, et al. A Novel Pathway of Chronic Allograft Rejection Mediated by NK Cells and Alloantibody. *Am J Transplant* (2012) 12(2):313-21. doi: 10.1111/ j.1600-6143.2011.03836.x

45. Calvani J, Terada M, Lesaffre C, Eloudzeri M, Lamarthée B, Burger C, et al. *In Situ* Multiplex Immunofluorescence Analysis of the Inflammatory Burden in Kidney Allograft Rejection: A New Tool to Characterize the Alloimmune Response. *Am J Transplant* (2020) 20(4):942-53. doi: 10.1111/ajt.15699

46. Punt J, Stanford SA, Jones PP, Owen JA. *Kuby Immunology. 8th ed.* New York: Macmillan Education (2019).

47. Cascieri MA, Springer MS. The Chemokine/Chemokine-Receptor Family: Potential and Progress for Therapeutic Intervention. *Curr Opin Chem Biol* (2000) 4(4):420-7. doi: 10.1016/S1367-5931(00)00113-7

48. Melter M, Exeni A, Reinders ME, Fang JC, McMahon G, Ganz P, et al. Expression of the Chemokine Receptor CXCR3 and Its Ligand IP-10 During Human Cardiac Allograft Rejection. *Circulation* (2001) 104(21):2558-64. doi: 10.1161/hc4601.098010

49. Schnickel GT, Hsieh GR, Garcia C, Shefizadeh A, Fishbein MC, Ardehali A. Role of CXCR3 and CCR5 in Allograft Rejection. *Transplant Proc* (2006) 38 (10):3221-4. doi: 10.1016/j.transproceed.2006.10.164

50. Bickerstaff A, Nozaki T, Wang JJ, Pelletier R, Hadley G, Nadasdy G, et al. Acute Humoral Rejection of Renal Allografts in CCR5(-/-) Recipients. *Am J Transplant* (2008) 8(3):557-66. doi: 10.1111/j.1600-6143.2007.02125.x

51. Fischereder M, Luckow B, Hocher B, Wüthrich RP, Rothenpieler U, Schneeberger H, et al. CC Chemokine Receptor 5 and Renal-Transplant Survival. *Lancet* (2001) 357(9270):1758-61. doi: 10.1016/S0140-6736(00) 04898-4

52. Sullivan JS. CC Chemokine Receptor 5 and Renal-Transplant Survival. *Lancet* (2001) 358(9289):1269-70. doi: 10.1016/S0140-6736(01)06360-7

53. Qin S, Rottman JB, Myers P, Kassam N, Weinblatt M, Loetscher M, et al. The Chemokine Receptors CXCR3 and CCR5 Mark Subsets of T Cells Associated With Certain Inflammatory Reactions. *J Clin Invest* (1998) 101(4):746-54. doi: 10.1172/JCI1422

54. Carlin LE, Hemann EA, Zacharias ZR, Heusel JW, Legge KL. Natural Killer Cell Recruitment to the Lung During Influenza a Virus Infection Is Dependent on CXCR3, CCR5, and Virus Exposure Dose. *Front Immunol* (2018) 9:781. doi: 10.3389/fimmu.2018.00781
55. Rabant M, Amrouche L, Lebreton X, Aulagnon F, Benon A, Sauvaget V, et al. Urinary C-X-C Motif Chemokine 10 Independently Improves the Noninvasive Diagnosis of Antibody-Mediated Kidney Allograft Rejection. *J Am Soc Nephrol* (2015) 26(11):2840-51. doi: 10.1681/ASN.2014080797

56. Yeaپ WH, Wong KL, Shimasaki N, Teo EC, Quek JK, Yong HX, et al. CD16 Is Indispensable for Antibody-Dependent Cellular Cytotoxicity by Human Monocytes. *Sci Rep* (2016) 6:34310. doi: 10.1038/srep34310

57. Wei X, Valenzuela NM, Rossetti M, Sosa RA, Nevarez-Mejia J, Fishbein GA, et al. Antibody-Induced Vascular Inflammation Skews Infiltrating Macrophages to a Novel Remodeling Phenotype in a Model of Transplant Rejection. *Am J Transplant* (2020) 20(10):2686-702. doi: 10.1111/ajt.15934

58. Krishnaswamy G, Kelley J, Yerra L, Smith JK, Chi DS. Human Endothelium as a Source of Multifunctional Cytokines: Molecular Regulation and Possible Role in Human Disease. *J Interferon Cytokine Res* (1999) 19(2):91-104. doi: 10.1089/107999099314234

59. Krausgruber T, Fortelny N, Fife-Gernedl V, Senekowitsch M, Schuster LC, Lercher A, et al. Structural Cells Are Key Regulators of Organ-Specific Immune Responses. *Nature* (2020) 583(7815):296-302. doi: 10.1038/s41586-020-2424-4

60. Nibbs RJ, Graham GJ. Immune Regulation by Atypical Chemokine Receptors. *Nat Rev Immunol* (2013) 13(11):815-29. doi: 10.1038/nri3544

61. Halloran PF, Venner JM, Famulski KS. Comprehensive Analysis of Transcript Changes Associated With Allograft Rejection: Combining Universal and Selective Features. *Am J Transplant* (2017) 17(7):1754-69. doi: 10.1111/ajt.14200

62. Kläger J, Eskandary F, Böhmig GA, Kozakowski N, Kainz A, Colin Aronovicz Y, et al. Renal Allograft Darcness in Subclinical Acute and Chronic Active ABMR. *Transpl Int* (2021) 34(8):1494-505. doi: 10.1111/tri.13904

63. Naemi FM, Carter V, Kirby JA, Ali S. Anti-Donor HLA Class I Antibodies: Pathways to Endothelial Cell Activation and Cell-Mediated Allograft Rejection. *Transplantation* (2013) 96(3):258-66. doi: 10.1097/TP.0b013e3182985504

64. Hayde N, Broin P, Bao Y, de Boccardo G, Lubetzky M, Ajaimy M, et al. Increased Intragraft Rejection-Associated Gene Transcripts in Patients With Donor-Specific Antibodies and Normal Biopsies. *Kidney Int* (2014) 86 (3):600-9. doi: 10.1038/ki.2014.75

65. Naesens M, Khatri P, Li L, Sigdel TK, Vitalone MJ, Chen R, et al. Progressive Histological Damage in Renal Allografts Is Associated With Expression of Innate and Adaptive Immunity Genes. *Kidney Int* (2011) 80(12):1364-76. doi: 10.1038/ki.2011.245

66. Loupy A, Lefaucheur C, Vernerey D, Chang J, Hidalgo LG, Beuscart T, et al. Molecular Microscope Strategy to Improve Risk Stratification in Early

Antibody-Mediated Kidney Allograft Rejection. *J Am Soc Nephrol* (2014) 25 (10):2267-77. doi: 10.1681/ASN.201311149

67. Stites E, Kumar D, Olaitan O, John Swanson S, Leca N, Weir M, et al. High Levels of Dd-Cfdna Identify Patients With TCMR 1A and Borderline Allograft Rejection at Elevated Risk of Graft Injury. *Am J Transplant* (2020) 20(9):2491-8. doi: 10.1111/ajt.15822

68. Parajuli S, Joachim E, Alagusundaramoorthy S, Aziz F, Blazel J, Garg N, et al. Donor-Specific Antibodies in the Absence of Rejection Are Not a Risk Factor for Allograft Failure. *Kidney Int Rep* (2019) 4(8):1057-65. doi: 10.1016/j.ekir.2019.04.011

69. Jacquemont L, Tilly G, Yap M, Doan-Ngoc TM, Danger R, Guérif P, et al. Terminally Differentiated Effector Memory CD8. *J Am Soc Nephrol* (2020) 31 (4):876-91. doi: 10.1681/ASN.2019080847

70. Loupy A, Lefaucheur C, Vernerey D, Prugger C, Duong van Huyen JP, Mooney N, et al. Complement-Binding Anti-HLA Antibodies and Kidney-Allograft Survival. *N Engl J Med* (2013) 369(13):1215-26. doi: 10.1056/NEJMoa1302506

71. Callemeyn J, Senev A, Coemans M, Lerut E, Sprangers B, Kuypers D, et al. Missing Self-Induced Microvascular Rejection of Kidney Allografts: A Population-Based Study. *J Am Soc Nephrol* (2021) 32(8):2070-82. doi: 10.1681/ASN.2020111558

**Conflict of Interest:** The 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.

**Publisher's Note:** All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Copyright © 2022 Van Loon, Lamarthée, Barba, Claes, Coemans, de Loor, Emonds, Koshy, Kuypers, Proost, Senev, Sprangers, Tinel, Thaunat, Van Craenenbroeck, Schols and Naesens. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.