Abstract
Background and Aims
TOMOGRAM, multicenter study founded by DESCARTES ERA/EDTA WG, aims to identify transcriptomic and genomic signatures of operational tolerance (OT) in recently identified cohort of OT kidney transplant recipients.
Method
RNA sequencing of peripheral blood was evaluated in 15 OT patients recently identified by TOMOGRAM consortium in 8 European countries, 23 stable patients (≥ 15 years on immunosuppression, STA), 14 CABMR patients (≥ 1 year, CR), 14 non-transplant CNI-treated patients and 14 healthy controls (HC). Differential expression was performed using DESEq2 and gene annotation analysis using Enrichr. Besides immunosuppression unadjusted model, robust negative-binomial regression model was created to adjust for immunosuppression intake. The models was trained on homogeneous group of STA patients.
Results
Using model unadjusted for immunosuppression, no differences in transcriptomic profiles between OT, STA and HC groups were identified. Nine transcripts were upregulated and 2 downregulated in OT compared CR group. The number of deregulated transcripts substantially increased when the model was adjusted for immunosuppression. Gene annotation analysis of top ranked deregulated 1109 transcripts (FC>2, adjusted p value <0.0001) showed deregulation of biological processes related to interferon-γ-mediated signaling pathway (p=1.4*10-5), response to cytokine (p=1.5*10-5), type I interferon signaling pathway (p=0.00036), regulation of I-kappaB kinase/NF-kappaB signaling (p=0.0021), cytokine-mediated signaling pathway (p=0.019) and neutrophil mediated immunity (p=0.033). While interferon-γ-mediated and type I interferon signaling were related to transcripts increased in CR, neutrophils associated transcripts were increased in OT. Analysis of cell types transcripts showed enrichment of CD19 B cells (p=1.6*10-9) in CR, while CD56NK cells (p=2.5*10-11) and CD8 T cells (p=1.6*10-11) transcripts predominated in OT. To reveal probability of operational tolerance inside STA group, 13 transcripts able to discriminate OT and CR cohorts with high AUC (>0.89) were used in PCA analysis (ADGRG3, ATG2A, GDPD5, IL16, MX2, SLA2, PRKD2, SLIRP, GNLY, SRCAP, ARGHAP9, IGHM, CD5). The high probability of OT signature was found in a single STA patient.
Conclusion
Contrary to previous reports which pointed out towards naïve B cell signatures, unique OT patients exhibit other specific immunosuppression-independent transcriptomic profiles.
Positive natural selection in primate genes of the type I interferon response