A Computational Gene Expression Score for Predicting Immune Injury in Renal Allografts

BACKGROUND: Whole genome microarray meta-analyses of 1030 kidney, heart, lung and liver allograft biopsies identified a common immune response module (CRM) of 11 genes that define acute rejection (AR) across different engrafted tissues. We evaluated if the CRM genes can provide a molecular microscope to quantify graft injury in acute rejection (AR) and predict risk of progressive interstitial fibrosis and tubular atrophy (IFTA) in histologically normal kidney biopsies. METHODS: Computational modeling was done on tissue qPCR based gene expression measurements for the 11 CRM genes in 146 independent renal allografts from 122 unique patients with AR (n = 54) and no-AR (n = 92). 24 demographically matched patients with no-AR had 6 and 24 month paired protocol biopsies; all had histologically normal 6 month biopsies, and 12 had evidence of progressive IFTA (pIFTA) on their 24 month biopsies. Results were correlated with demographic, clinical and pathology variables. RESULTS: The 11 gene qPCR based tissue CRM score (tCRM) was significantly increased in AR (5.68 ± 0.91) when compared to STA (1.29 ± 0.28; p < 0.001) and pIFTA (7.94 ± 2.278 versus 2.28 ± 0.66; p = 0.04), with greatest significance for CXCL9 and CXCL10 in AR (p <0.001) and CD6 (p<0.01), CXCL9 (p<0.05), and LCK (p<0.01) in pIFTA. tCRM was a significant independent correlate of biopsy confirmed AR (p < 0.001; AUC of 0.900; 95% CI = 0.705–903). Gene expression modeling of 6 month biopsies across 7/11 genes (CD6, INPP5D, ISG20, NKG7, PSMB9, RUNX3, and TAP1) significantly (p = 0.037) predicted the development of pIFTA at 24 months. CONCLUSIONS: Genome-wide tissue gene expression data mining has supported the development of a tCRM-qPCR based assay for evaluating graft immune inflammation. The tCRM score quantifies injury in AR and stratifies patients at increased risk of future pIFTA prior to any perturbation of graft function or histology.