Activation of counter-regulatory mechanisms in a rat renal acute rejection model

Background Microarray analysis provides a powerful approach to identify gene expression alterations following transplantation. In patients the heterogeneity of graft specimens, co-morbidity, co-medications and the challenges in sample collection and preparation complicate conclusions regarding the underlying mechanisms of graft injury, rejection and immune regulation. Results We used a rat kidney transplantation model with strict transplant and sample preparation procedures to analyze genome wide changes in gene expression four days after syngeneic and allogeneic transplantation. Both interventions were associated with substantial changes in gene expression. After allogeneic transplantation, genes and pathways related to transport and metabolism were predominantly down-regulated consistent with rejection-mediated graft injury and dysfunction. Up-regulated genes were primarily related to the acute immune response including antigen presentation, T-cell receptor signaling, apoptosis, interferon signaling and complement cascades. We observed a cytokine and chemokine expression profile consistent with activation of a Th1-cell response. A novel finding was up-regulation of several regulatory and protective genes after allogeneic transplantation, specifically IL10, Bcl2a1, C4bpa, Ctla4, HO-1 and the SOCS family. Conclusion Our data indicate that in parallel with the predicted activation of immune response and tissue injury pathways, there is simultaneous activation of pathways for counter regulatory and protective mechanisms that would balance and limit the ongoing inflammatory/immune responses. The pathophysiological mechanisms behind and the clinical consequences of alterations in expression of these gene classes in acute rejection, injury and dysfunction vs. protection and immunoregulation, prompt further analyses and open new aspects for therapeutic approaches.


Background
Microarray analyses have been used to link changes at the level of gene expression to different kidney diseases to obtain markers for diagnosis and prognosis. Microarray analyses of human patients have been used to analyze post-transplant events [1][2][3][4][5][6]. A recent review of microarray-based studies and the search for biomarkers in organ transplantation is given by Kurian et al. [7]. Microarray analysis of kidney transplant biopsies with acute rejection identified gene expression patterns that distinguish three distinct subtypes of acute rejection that, although indistinguishable by histology, were marked by differences in mechanisms of immune activation and cell proliferation [8]. Interestingly, another study of acute rejection biopsies showed no evidence for up-regulation of cytotoxic T-cell effector molecules despite the fact that these have been considered markers for acute renal rejection [9]. In parallel, gene expression signatures of peripheral blood lymphocytes (PBLs) from transplant patients have also been shown to be capable of classifying patients with acute rejection [6]. Surprisingly, there was essentially no overlap in rejection-diagnostic genes up-regulated in PBL and genes up-regulated in the biopsies from the same patients, indicating that the blood must be considered a very different immune compartment in this setting. As a consequence of these initial studies, the hope of developing validated diagnostic biomarkers for transplantation is tempered by results that challenge an existing dogma on immune response mechanisms and demand additional studies to deconvolute.
There are a number of limitations inherent in clinical studies of gene expression in transplantation. Patients typically represent major differences in gender, age, co-morbidities, clinical histories, immunosuppressive regimes, race/ethnicities, and genetics and it is always difficult to collect samples at specific times relative to the procedures. These differences significantly influence gene expression and that operant variability compromises the results obtained by microarray analysis as well as the confidence to identify the principal underlying immune mechanisms. The use of experimental transplantation models in rats or mice is an opportunity to mitigate these variables and work in a controlled system. Several studies have made use of experimental transplantation to study changes in gene expression related to acute rejection after murine heart, kidney and lung transplantation [10][11][12][13]. While several studies have been performed using heart transplantation in rats [14,15], to our knowledge, a study of acute kidney transplant rejection in rats has not been reported.
In previous studies [16,17], we showed in well characterized rat renal transplantation models, no evidence of tissue necrosis at day 4 while histological changes consistent with acute rejection including activated lymphocyte infil-tration were found consistently. In those studies we showed differences in gene expression for selected sodium and water transporters that were relevant to tissue injury and renal compensation. In the present study, we made use of these transplantation models and continued our studies using day 4 as the measuring point to represent a time at which there is full representation of acute rejection mechanisms but before the widespread renal tissue injury that can confuse the gene signatures. Thus, changes in gene expression related solely to acute rejection in allogeneic transplantation or linked to the surgical procedure in syngeneic transplantation can be evaluated. Native, untreated kidneys were used as the control group. No immunosuppression was used to specifically analyze the influence of acute rejection on gene expression independent of the effects of immunosuppressive regimens. Thereby, maximum homogeneity in procedures, sample preparation, timing, and handling was achieved. In contrast to the design of renal transplantation studies performed in mice [12,13], we used bilaterally nephrectomized recipient rats, such that the animal's transplant was required to maintain the full metabolic and synthetic function of the organism as is the case in the clinical transplant situation.

Correlation and clustering of the samples
The global experiment structure was inspected utilizing different methods intended to validate whether the groups were homogeneous and whether the within-group variance was smaller than the between-group variance. Furthermore, group distances could be roughly estimated. The correlations within the "acute group" (allogeneic kidney transplant with acute rejection), "syngeneic group" (syngeneic transplant as control for surgical changes) and "control group" (untreated, native kidneys) showed a homogenous correlation of 0.97-0.99 (Fig. 1A). The correlation between the control group and the syngeneic group was 0.95-0.97, between control and the acute group was 0.84-0.87, and between the acute group and the syngeneic group was 0.88-0.92. The representative scatter plots of all genes between different samples showed that the observed correlations are not dominated by a small number of highly expressed genes (Fig. 1B). The syngeneic group was more similar to the control group than the acute group was to the control group. Thus, the processes of acute rejection were the most relevant factors in driving the observed changes in gene expression.

Identification of differentially expressed genes
Significantly regulated and differentially expressed genes were identified using class comparisons as implemented in BRB ArrayTools for all the possible group comparisons (acute to control, acute to syngeneic and syngeneic to con-trol). Genes with a parametric p value < 0.001 were classified as significantly differentially expressed ( Fig. 2A). Thus, 3871 genes were up-regulated and 3483 were downregulated for the acute group compared to the control. Comparison of the acute to the syngeneic group revealed 2668 genes were up-and 3236 genes were down-regulated. In contrast, comparing the control to syngeneic groups for the impact of the transplant surgery revealed 564 genes were down-regulated and 1291 were up-regulated. The Venn diagram in Figure 2B shows the overlap between the lists of significantly altered genes as a function of the three experimental groups. For further analysis we used the gene lists from syngeneic transplants compared to controls and from allogeneic transplants compared to controls. Specifically, we decided not to report on our comparison of the allogeneic transplants to the syngeneic transplants because the results of the functional and pathway analysis were nearly identical with the results obtained using the comparison of allogeneic transplants with controls. There were marginal differences in the number of genes related to a distinct function or pathway. However, the complete output files with all regulated genes are provided as additional file 1.
Genes with protective or regulatory properties up-regulated after allogeneic transplantation like HO-1 (22 fold), Bcl2a1 (59 fold), C4bpa (219 fold), IL10 (120 fold) or members of a gene family that were highly induced like TIMP1 (18 fold), MMP7 (15 fold), TAP1/2 (24/14 fold) were analyzed by real time PCR ( Fig. 2C and 2D). FoxP3 and CD25 are expressed by regulatory T-cells [18]. The results for expression of genes that were not significantly regulated after syngeneic transplantation according to the array were set to 1. The PCR results confirmed the up-regulated genes after allogeneic transplantation identified by the expression arrays. The changes in expression varied between 10-fold for Tap2 up to ~500-fold for C4bp. It is important to note that there is no corresponding probe set on the Rat 230 2.0 microarray for FoxP3. Therefore the FoxP3 was also chosen for Real-Time PCR analysis. Interestingly the expression of FoxP3 was highly induced after allogeneic transplantation (~14 fold).
The genes analyzed by real time PCR were not identified to be differently expressed after syngeneic transplantation compared to control. However, the real time PCR results showed up-regulated expression of MMP7 (~21-fold) and C4bp (~13-fold) after syngeneic transplantation.

Genes involved in immune response
A selection of transcripts determined to be up-regulated as a function of acute rejection is shown in Table 1. The recruitment of cells related to the immune response is indicated by up-regulation of CD8 (229-fold), the marker for cytotoxic T-cells and by the up-regulation of granzyme-A (67-fold), granzyme-B (103-fold) and perforin (25-fold), cytotoxic effector molecules released by CD8+ T-cells. Evidence for the involvement and activation of macrophages is the elevated expression of CD163, the marker for macrophages [19] and ficolin B, which is expressed by activated macrophages [20].
T-cells and macrophages release factors like cytokines for differentiation and proliferation and chemokines for recruitment of other immunologically active cells. Thus, Correlation and clustering of the samples Figure 1 Correlation and clustering of the samples. A) Correlation analyses of the samples showed the expected results with higher correlation between control group (ctr) and syngeneic group (sTX) samples compared to acute rejection group (aTX). B) The representative scatter plots of all genes between different samples shows, that the observed correlation is not dominated by a small number of high-expressing genes.
Differently expressed genes between the control, syngeneic and allogeneic group Figure 2 Differently expressed genes between the control, syngeneic and allogeneic group. Class comparisons were performed using BRB array tools to identify genes with significant changes in expression levels between the syngeneic group (sTX_con) or the allogeneic group (aTX_con) compared to the control group and between allogeneic and syngeneic group (aTX_sTX; Fig. 2A). Genes with a p-value < 0.001 were classified as significant regulated. The complete lists with the significant regulated genes are provided with this manuscript (see additional file 1). The Venn diagram in Fig. 2B shows the number of overlapping genes. The changes in expression compared to control for selected genes were analyzed by real time PCR using specific primer pairs or Taq-Man gene expression assays. Relative changes were evaluated using the 2 -ΔΔ Ct method. The changes in gene expression after allogeneic transplantation are shown in Fig. 2C and the results for the syngeneic transplantation group in Fig. 2D. The black columns show the corresponding expression levels on the arrays (ratio of mean signal intensities). MMP7. IL10, C4bpa, Tap1 and Tap2 were not differently expressed after syngeneic transplantation on the array. There was no corresponding probe set for FoxP3 on the array. Immune cells related genes the expression pattern of cytokines and chemokines like Cxcl9, Xcl1, Cxcl10 or Ccl3 shown in Table 1 is consistent with their well established importance in the signaling mechanisms activated after allogeneic transplantation (for review see [21]).

Functional annotation using gene ontology terms
To identify activated or suppressed functions, the lists of significantly up-or down-regulated genes after syngeneic transplantation and allogeneic transplantation were analyzed to identify over-represented GO-terms. The third GO-level gives the best compromise between specificity and list coverage [22]. After syngeneic transplantation, 18 GO-terms were enriched in the up-regulated and 1 in the down-regulated group of genes on GO-level 3. We also analyzed over-represented GO-terms without filtering on the third level. Table 2 shows the ten most significant over-represented genes related to biological process (BP), cellular component (CC) and molecular function (MF). Due to the limitation of space the complete lists are also provided with this manuscript (see additional file 2).
Decreased expression of metabolism-related genes after syngeneic transplantation suggests the suppression of BPrelated genes. Increased expression of genes related to "cell proliferation" and "cell cycle" within the BP category, "extra cellular matrix" and "collagen" within the CC category and "collagen binding" within the MF category reveal the activation of tissue injury recovery mechanisms related to the transplant surgery itself.
After allogeneic transplantation, 58 GO-terms were significantly enriched in the up-regulated and 29 in the downregulated genes on GO-level 3. Without the application of any filter, 250 GO-terms were enriched in the up-and 119 in the down-regulated set of genes. Table 3 shows the ten most significant over-represented genes related to BP, CC and MF. As expected, the increased expression of genes related to the BP categories "immune response", "defense response" and the CC category "T-cell receptor complex" and the MF category called "cytokine binding" indicates the activation of the immune response. These GO-terms are consistent with infiltration of the graft by T-cells and describe signaling mechanisms involved in activation, differentiation and recruitment of immune cells. The BP 'cell death' group describes the downstream processes following T-cell activation. The decreased expression of genes related to the BP categories "metabolism", "transport", "excretion" and the MF category "primary active transporter activity" indicates the challenges posed to the graft to maintain renal function. A major part of renal function is mediated by several transporters and transport systems all primarily or secondarily coupled to energy consumption. The down-regulation of genes that are involved in metabolism or transport leads to a decreased tubular function of the kidney.

Pathways affected after syngeneic and allogeneic transplantation
The functional annotation with DAVID identified pathways related to the Kyoto Encyclopedia of Genes and Genomes (KEGG) [23] that were affected after syngeneic or allogeneic transplantation. The pathways affected after syngeneic transplantation are listed in Table 4. The majority of down-regulated genes are related to metabolism pathways. Increased expression of genes involved in "cell cycle" and "extra cellular matrix (ECM)-receptor interaction" are suggestive of active recovery processes following the surgery.
Pathways with enrichment of down-or up-regulated genes after allogeneic transplantation are listed in Table 4. The majority of the pathways with down-regulated genes are also related to metabolism. In all cases the numbers of genes related to metabolism pathways were higher after allogeneic than syngeneic transplantation indicating a stronger suppression of these pathways in the context of acute rejection. But some pathways, however, were affected only after allogeneic transplantation, like 'oxidative phosphorylation', 'citrate cycle' or 'metabolism of xenobiotics by cytochrome p450' indicating an acute rejection mediated down-regulation of genes related to these pathways. The majority of pathways with enrich-  To illustrate the pathways activated by acute rejection related to the immune response, we used IPA to map net-works of up-and down-regulated genes after allogeneic transplantation to functional pathways for antigen-presentation (Fig. 3), TCR-signaling (Fig. 4), apoptosis (Fig.  5), IFNγ (Fig. 6) and complement cascades (Fig. 7). Within the antigen processing and presentation pathway  Table 2 The numbers indicate the quantity of genes corresponding with the GO-terms.   nearly all genes were affected (Fig. 3). The immunoproteasome alters the activity of the proteasome by replacement of the β1, β2 and β5 subunits of the 20S core proteasome by LMP2, LMP7 and Psmb10 [24]. The expression of these subunits was strongly induced by acute rejection. The antigens processed by the proteasome are transported by molecules associated with the antigen processing (TAP) complex [25]. Thus, we saw strong induction of TAP1 (24fold) and TAP2 (15-fold; Fig. 2C).

05, Fisher's exact test, with Benjamini and Hochberg correction) related to biological process (BP), cellular component (CC) and molecular function (MF) in the set of genes down-or up-regulated after allogeneic transplantation (aTX) compared to control are shown in
The recognition of allogeneic MHC by the T-cell's cognate antigen receptor (TCR) will induce activation signals and transcriptional activity (Fig. 4). The gene expression changes indicate an induction of NFATc1 and NF-κB related pathways. Within these pathways, genes with counter regulatory properties like CTLA-4 or cbl were also up-regulated. Fig. 6 illustrates the signaling pathway for IFNγ. IFNγ acts as a pro-inflammatory cytokine with altered expression after allogeneic transplantation (Table 1). Nearly all the genes within this pathway showed altered expression in the setting of acute rejection. The expression of signal transducer and activator of transcription (STAT) 1 and STAT 2, the downstream effectors of IFNγ signaling, was highly induced. But the increased expression of the suppressor of cytokine signaling (SOCS) 1 and of the protein tyrosine phosphatase non-receptor type 2 (TC-PtP), both negative regulators of IFNγ signaling, indicate the simultaneous activation of regulatory mechanisms to limit the IFNγ signaling.

Potential novel markers for acute rejection
Several genes in the list of up-regulated genes after allogeneic transplantation might serve as novel markers for  acute rejection (Table 5). One candidate is the ETS transcription factor, Spic (SPI-C; 113-fold). SPI-C is expressed during B-lymphocyte development and in activated macrophages [26]. SPI-C also interacts with STAT6 promoting IL4 dependent IgE expression [27]. No elevated expression level of IL4 and only a marginal up-regulation of STAT6 after allogeneic transplantation were found so that the function of SPI-C remains unclear. Phospholipase-A2, group IID (98-fold) and phospholipase-A2, group IIA (75-fold) have been shown to be up-regulated in human broncho-epithelial and nasal epithelial cells after treatment with IFNγ [28]. The authors postulated a possible involvement of both in cytokine-mediated inflammation.

Discussion
The MHC antigens are the main barrier for an acceptance of the graft by the host organism and antigen presentation is essential for the activation of T-cells. T-cells are activated when the TCR recognizes allo-MHC class I molecules in the case of CD8+ and MHC class II molecules in the case of CD4+ T-cells [29]. Host T-cells can recognize MHC antigens by either indirect or direct presentation [30]. During indirect presentation, peptides derived from donor MHC class I or MHC class II molecules are processed and presented to CD4+ T-cells by host antigen presenting cells expressing the recipient's MHC class II molecules. The direct antigen presentation pathway is unique to transplantation. In this pathway the T-cell receptors directly recognize an intact MHC molecule expressed on donor cells. In this context, it is important that for our studies with acute rejection, nearly all the genes associated with the antigen processing and presentation pathways were affected (Fig. 3). Our data also indicates that T-cell cytotoxicity is the main immunological pathway activated in response to acute rejection in allogeneic transplantation. The cytotoxic molecules are perforin/granzyme for CD8+ T-cells and the Fas/Fas-ligand system for CD4+ T-cells [31,32]. In this context, the absolute changes in expression of genes involved in antigen processing and presentation by MHC class I molecules indicates that this part of the antigen presentation is intensely up-regulated. That is also consistent with the up-regulation of the immunoproteasome and TAP1/2 transcripts (Fig. 3) likely to be driven by the upregulation of interferon γ [33,34], which we have also documented. IFNγ can also induce the expression of the chemokines, CXCL9, CXCL10 and CXCL11 [35]. These data are consistent with other studies showing the changes in gene expression related to IFNγ or cytotoxic T lymphocytes in murine transplantation models [12,13].

Antigen presentation pathway
The induction of tolerance to the graft would be the treatment of choice after transplantation. Increased expression of protective genes like metallothionein-1 or α2βcrystallin were reported with a mouse cardiac transplant model [36,37]. In this same context, we observed increased expression of a number of genes known to have protective or regulatory properties (Figs. 4, 5, 6, 7) creating a much more comprehensive list than reported in any previous studies. Our study also indicates the endogenous activation of counter regulatory mechanisms within several pathways that have been linked to acute rejection. For example, genes with anti-apoptotic functions include Bcell leukemia/lymphoma 2 related protein A1 (Bfl-1) (Fig.  5) and the Inhibitor of Apoptosis Proteins (IAP) [38,39]. The expression of Heme Oxygenase 1 (HO1; 23-fold) was induced after allogeneic transplantation. The expression of HO1 can suppress graft rejection and leads to long term graft survival in some models [40]. The Cytotoxic T-Lymphocyte-Associated protein 4 (CTLA-4) is a negative regulator of T-cells (Fig. 4) [41]. The function of IL10 is to limit TCR signaling pathway Figure 4 TCR signaling pathway. Ingenuity Pathways Knowledge Base was used to display pathways containing multiple differently expressed genes. Significant up-regulation of genes in the canonical TCR signaling pathway was observed after allogeneic transplantation. Up-regulated genes are highlighted in red the down-regulated genes are highlighted in green. Figure 5 Apoptosis signaling. In the same way as described in figure  3 significant regulated genes after allogeneic transplantation was observed in the canonical apoptosis signaling Up-regulated genes are highlighted in red the down-regulated genes are highlighted in green.

Apoptosis signaling
the action of the immune response. IL10 was initially characterized as a cytokine produced by Th2-cells, which inhibited the production of cytokines such as IL2, TNFα and IFNγ [42]. IL10 also modulates the expression of chemokines, chemokine receptors and MHC class II molecules [43]. Thus, the induction of IL10 after allogeneic transplantation indicates a possible feed back mechanism of the immune response (Table 1, Fig. 2).
The complement system is part of the innate immune response. The expression profile we demonstrate in allogeneic transplants for a number of complement components indicates the activation of the complement cascade by acute rejection (Fig. 7). Its major function is the elimination of pathogens and the stimulation of an inflammatory response [44,45]. Local gene expression for complement components has been shown for C3 and C1 during acute rejection in kidney and heart transplantation models, respectively [36,46]. One of the novel genes in our study with the strongest expression after allogeneic transplantation was the C4b binding protein alpha (C4bpa; 219 fold up-regulated). C4bpa is a cofactor for factor 1 in the degradation of C4b in the complement cascade [47]. The degradation of C4b inhibits the formation of the C4b2a complex, a key step in producing an inflammatory response [48]. Figure 6 Interferon signaling. Significant up-regulation of genes was also observed in the canonical interferon signaling pathway. Upregulated genes are highlighted in red.

Interferon signaling
Recently regulatory T-cells have been recognized as a promising tool for the induction of tolerance [49]. In these studies of our rat renal transplants we also find transcriptional evidence for the presence of regulatory T-cells within the graft. The function of regulatory T-cells is to regulate or suppress the activation, regulation and function of effector T-cells [50]. The development of regulatory T-cells appears to be controlled by FoxP3 [18], which we find is up-regulated after allogeneic transplantation (Fig. 2). The intragraft expression of FoxP3 has also been demonstrated in cardiac allograft patients [51]. The presence of regulatory T-cells within the grafts we studied is further supported by the concomitant increase in expression of CD25 [49] and tryptophan hydroxylase [52].

Conclusion
We have used a well-defined rat kidney transplantation model to demonstrate global changes in gene expression levels after syngeneic and allogeneic transplantation with the emphasis on early acute rejection in the allogeneic transplants. Using DNA microarrays and real time PCR analysis, we were able to identify major immune effectorrelated pathways that are activated after allogeneic transplantation and clearly linked by multiple lines of published evidence to acute rejection mechanisms of tissue injury. That we also found a list of genes that are equally highly differentially expressed in acute rejection but not presently linked in the literature represents opportunities for additional studies and discovery. In parallel, we have also identified a number of up-regulated genes that are linked to tissue protective and immune counter-regulation mechanisms. These observations fit with the evolving view of immunity as a balance between effector and regulatory mechanisms.

Kidney transplantation
Male Lewis-Brown-Norway (LBN) and Lewis (LEW) rats (250-300 g, Charles River, Sulzfeld, Germany) with free access to standard rat chow (Ssniff, Soest, Germany) and tap water were used. Experiments were approved by a governmental committee on animal welfare and were performed in accordance with national animal protection guidelines. Renal transplantation was performed as published before [17,53]. For the present study, all recipients were bilaterally nephrectomized immediately before transplantation. In brief, the left kidney including ureter, renal artery, a piece of aorta and renal vein was transplanted into the recipient. For the acute rejection-model, kidneys of LBN-rats (n = 5) were transplanted into LEWrats and, for the syngeneic transplantation-model, kidneys obtained from LBN-rats (n = 5) were transplanted into LBN-rats. The acute rejection-model leads to marked histological changes typical for acute transplant rejection Complement cascade Figure 7 Complement cascade. Significant regulated genes after allogeneic transplantation were associated with the canonical complement cascade using the Ingenuity Pathways Knowledge Base to visualize genes differently expressed in the selected pathway. Up regulated genes are highlighted in red the down regulated genes are highlighted in green.
[53] For the control group, we used native untreated kidneys from LBN-rats (n = 5).

RNA isolation, labeling, hybridization and scanning
On day four after transplantation, the transplanted kidneys were removed and total RNA was isolated using RNeasy-kit (Qiagen, Hilden, Germany) incubated with 10 U DNase I (Promega, Heidelberg, Germany) to digest genomic DNA. RNA quality was measured using the Agilent 2100 BioAnalyzer (Agilent, Palo Alto, CA, USA). Total RNA was used to prepare biotinylated target RNA. Briefly, 10 μg of total RNA was used to generate first-strand cDNA by using a T7-linked oligo (dT) primer. After second strand synthesis, in vitro transcription was performed with biotinylated UTP and CTP (Enzo Diagnostics, New York, NY, USA).
Target cDNAs generated from each sample were then processed as per manufacturer's recommendation using an Affymetrix GeneChip instrument system. Labeled samples were hybridized to the Rat Genome 230 2.0 Array. Arrays were washed and stained with streptavidin-phycoerythrin before being scanned on an Affymetrix GeneChip scanner. Data were analyzed using Affymetrix GCOS array analysis software. The data discussed in this publication have been deposited in NCBIs Gene Expression Omnibus [54] and are accessible through GEO Series accession number GSE6497.

Statistical analysis
Experiment reports were inspected to assure equal experiment quality. GeneData Refiner (GeneData, Basel, Switzerland) was used to check quality and import experiments with detection and masking of outliers (the masked area was 0.026 %, 0.01 % and 0.02 % for the control, the syngeneic and allogeneic group respectively). The array defects, assessment using 3'/5' ratios of housekeeping controls, glyseraldehyde-3-phosphate dehydrogenase (GAPDH) and β-Actin, and condensation of intensities to signal values using the "Affymetrix Statistical (MAS5)" method. In GeneData Analyst normalization onto arithmetic mean 400 was performed. Correlation coefficients of all pair wise scatter plots were calculated and summarized in a colored matrix. Clustering was done using the Euclidean distance measure and the complete linkage method. Correlation coefficients and clustering enable a visualization and assessment of the experiment structure based on global expression behavior. Identification of significantly different expressed genes between the control, allogeneic and syngeneic groups were identified using class comparison with the BRB ArrayTools developed by Dr. Richard Simon and Amy Peng [55]. The type of univariate test used was a two-sample T-test. Exact multivariate permutations test was computed based on 126 available permutations. Nominal significance level of each univariate test was set to 0.001. Confidence level of false discovery rate assessment used was 90 % and the maximum allowed numbers of false-positive genes were set to 10.

Functional annotation and pathway analysis
Up-and down-regulated genes in the allogeneic and syngeneic groups compared to control group were separately analyzed to identify enrichment of different functional pathways and gene ontology (GO) terms [56] using the Database for Annotation, Visualization, and Integrated Discovery (DAVID) [22]. The separate analysis of the down-or up-regulated genes was intended to provide information on suppressed or induced functions, respectively. The enrichment analysis was performed on the third GO level, because this gave the best compromise of specificity and coverage of the gene list. We performed the analysis with and without setting the filter on the third GO level. The Rat 230 2.0 gene chip served as the background list (Fishers exact test with corrected p-value < 0.05, Benjamini & Hochberg [57]). To illustrate the pathways with affected genes in the acute rejection group the data were analyzed through the use of Ingenuity Pathways Analysis (IPA; Ingenuity ® Systems, Redwood City, CA, USA).

Real time PCR
Expression profiles for selected genes were analyzed by real time PCR. Total RNA (10 μg) isolated from the same All instruments and reagents were purchased by Applied Biosystems (Darmstadt, Germany). Relative gene expression values were evaluated with the 2 -ΔΔCt method using GAPDH or 18s-RNA as housekeeping genes [58].