However, PCA shown closer proximity of the combined rejection instances to AMR. we apply a novel gene expression platform in sequential formalin-fixed paraffin-embedded (FFPE) samples from nonhuman primate (NHP) renal transplants. We analyzed 34 previously-described gene transcripts related to AMR in humans in 197 archival NHP samples, including 102 from recipients that developed chronic AMR, 80 from recipients without AMR, and 15 normal native nephrectomies. Three endothelial genes (were manufactured for the mRNA sequences of 38 genes (Integrated DNA Systems, Coralville, IA). These included a previously-described AMR 34-gene-set comprised of 18 endothelial, 6 NK cell, and 10 inflammation-related genes, as well as 4 housekeeping genes (12). Probe sequences are provided in Table S1. Gene manifestation was then quantified with the NanoString? nCounter? Gene Manifestation assay (NanoString Systems, Seattle, WA) as per manufacturer instructions. To assess reproducibility, eight samples were randomly selected for duplicate analysis in individual runs. Quality control assessment and normalization of natural NanoString? gene expression results were performed with nSolver? Analysis Software Version 3.0 (NanoString Technologies, Seattle, WA) using the manufacturer-recommended default parameters. Retrospective Analysis of Human Microarray Data Set A publicly-available human cDNA microarray data set (GSE36059) was retrieved from your NCBI Gene Expression Omnibus (www.ncbi.nlm.nih.gov/geo) to compare gene set overall performance in human samples. The data set included gene expression results (Affymetrix? Human Genome U133 Plus 2.0 Array) for 403 renal allograft indication biopsies with diagnostic labels of AMR (n=65), TCMR (n=35), mixed rejection (n=22), and non-rejection (n=281), as per Banff 2009 criteria (22). The natural data files were imported and normalized with BRB-ArrayTools Version 4.5.0 (23). Mean values were utilized for genes with repeat data points. Data Analysis Post-normalization statistical analysis and visualization were performed with R version 3.3.2 (R Foundation for Statistical Computing, Vienna, Austria). Normalized transcript counts (NanoString? data) and log intensity values (microarray data) were converted to z-scores (quantity of standard deviations away from the population mean for each gene) for individual gene analysis. Mean z-scores were utilized for aggregate gene set analysis. Gene expression and correlation warmth map analyses (heatmap.2 function in gplots package) were performed with unsupervised hierarchical clustering by Euclidean distance. Spearmans rank correlation coefficients (cor function in stats package) and unsupervised principal component analysis (PCA; prcomp function in stats package) were used to characterize inter-variable associations. Mann-Whitney U-tests (wilcox.test function in stats package) were utilized for class comparison analyses. Receiver operating characteristic (ROC) curve analysis (roc function in pROC package) was utilized for assessment of diagnostic overall performance. Youdens J-statistic BMS-986165 (point on ROC curve farthest from diagonal index collection) was utilized for defining diagnostic thresholds (24). Individual gene rating, gene set construction, and gene set ranking were achieved with repeated 10-fold cross-validation analysis (train function in caret package) using three repeats and Naive Bayes model. P-values less than 0.05 were considered statistically significant. RESULTS RNA and Quality Control The imply RNA yield obtained from three 20-m sections per FFPE block was 6138 ng (range: 212C66508 ng) with a imply concentration of 153.5 ng/L (5.3C1662.7 ng/L) and a mean A260/A280 RNA purity ratio of 1 1.85 (1.54C2.20). No quality control or normalization flags were encountered BMS-986165 during nSolver? analysis for any of the 197 NHP samples included. The eight samples analyzed in duplicate exhibited excellent reproducibility with a mean correlation coefficient of 0.990 (range: 0.953C0.999). Individual Gene Expression vs. Diagnosis Following histological evaluation, the samples were assigned one of the following diagnostic labels according to Banff 2015 criteria: AMR (n=38), suspicious for AMR (n=15), mixed rejection (n=27), TCMR (n=41), borderline (n=21), no rejection (n=32), other (n=8), and normal native nephrectomy (n=15). Warmth map analysis revealed general grouping of diagnostic groups based on individual gene expression patterns (Physique 1). Sixty-seven (74%) of the 91 samples diagnosed as normal, no rejection, other, borderline, or suspicious for AMR clustered within the larger No Rejection group indicated in Physique 1. Thirty-five (54%) of the 65 samples diagnosed as real AMR or mixed rejection clustered within the AMR group. Five (33%) of the 15 samples diagnosed as suspicious for AMR also ACVRL1 clustered within the AMR group. Sixteen (39%) of the 41 samples diagnosed as TCMR BMS-986165 clustered within the TCMR group. Endothelium-associated transcripts generally exhibited higher expression in the AMR group. Inflammation-related transcripts showed a tendency for higher expression in the.