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–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–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 infiltration 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.
Results
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.
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.
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 control). 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 down-regulated 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.
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 TaqMan 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.
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].
Table 1
Selection of significant up-regulated expression of cytokines, cytokine receptors, chemokines, chemokine receptors, regulatory or protective genes after allogeneic transplantation compared to control group.
Fold change
Probe set
Description
Gene symbol
Cytokines
121.2
1387711_at
interleukin 10
Il10
84.1
1370790_at
interferon gamma
Ifng
41.6
1391384_at
tumor necrosis factor (TNF superfamily, member 2)
Tnf
15.4
1398256_at
interleukin 1 beta
Il1b
13.1
1368592_at
interleukin 1 alpha
Il1a
8.6
1369665_a_at
interleukin 18
Il18
8.4
1368763_at
interleukin 3
Il3
7.4
1369191_at
interleukin 6
Il6
6.1
1369208_at
interleukin 7
Il7
Cytokine receptors
17
1389092_at
interleukin 2 receptor, gamma
Il2rg
13.6
1387591_at
interleukin 2 receptor, alpha chain
Il2ra
11.8
1369697_at
interleukin 8 receptor, beta
Il8rb
10.1
1387394_at
interleukin 2 receptor, beta chain
Il2rb
3.5
1370728_at
interleukin 13 receptor, alpha 1
Il13ra1
3
1386987_at
interleukin 6 receptor
Il6r
Chemokines
139.7
1382454_at
chemokine (C-X-C motif) ligand 9
Cxcl9
122.8
1387831_at
small inducible cytokine subfamily C, member 1 (lymphotactin)
T-cells and macrophages release factors like cytokines for differentiation and proliferation and chemokines for recruitment of other immunologically active cells. Thus, 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).
Table 2
Overrepresented GO categories compared to control within the syngeneic transplantation group. The top 10 significantly overrepresented GO categories on the third level (p < 0.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 syngeneic transplantation (sTX) compared to control are shown in Table 2 The numbers indicate the quantity of genes corresponding with the GO-terms.
sTX_down
Category
Term
Count
p-value
BP
organic acid metabolism
57
1,1E-10
BP
carboxylic acid metabolism
56
1,5E-10
BP
amino acid and derivative metabolism
41
7,1E-9
BP
amino acid metabolism
31
2,2E-6
BP
nitrogen compound metabolism
40
8,3E-6
BP
amine metabolism
38
9,7E-6
BP
sulfur amino acid metabolism
10
1,8E-5
BP
amino acid derivative metabolism
18
6,0E-5
BP
amino acid catabolism
13
1,1E-4
BP
nitrogen compound catabolism
15
1,2E-4
CC
mitochondrion
59
1,9E-7
CC
microbody
14
4,0E-4
CC
peroxisome
14
4,0E-4
MF
oxidoreductase activity
59
3,2E-9
MF
catalytic activity
195
1,5E-6
MF
lyase activity
20
1,7E-3
sTX_up
BP
mitotic cell cycle
44
5,3E-7
BP
cell division
30
6,1E-7
BP
cell cycle
81
6,8E-7
BP
cell adhesion
70
2,5E-6
BP
mitosis
29
2,6E-6
BP
M phase of mitotic cell cycle
29
2,7E-6
BP
cell proliferation
83
1,1E-5
BP
M phase
33
3,9E-5
BP
regulation of cell cycle
54
7,9E-5
BP
regulation of progression through cell cycle
54
7,9E-5
CC
extracellular matrix
56
3,5E-16
CC
extracellular matrix (sensu Metazoa)
54
2,7E-15
CC
intracellular non-membrane-bound organelle
144
6,0E-13
CC
non-membrane-bound organelle
144
6,0E-13
CC
collagen
16
2,0E-7
CC
ribosome
36
2,4E-5
CC
cytoskeleton
76
4,1E-5
CC
spindle
17
4,1E-5
CC
microtubule cytoskeleton
40
4,2E-5
CC
chromosome
36
4,4E-5
MF
structural molecule activity
75
3,7E-8
MF
extracellular matrix structural constituent
23
9,6E-8
MF
structural constituent of ribosome
32
9,3E-5
MF
collagen binding
8
1,7E-2
MF
pattern binding
17
4,3E-2
MF
enzyme inhibitor activity
29
4,7E-2
Decreased expression of metabolism-related genes after syngeneic transplantation suggests the suppression of BP-related 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 down-regulated 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.
Table 3
Overrepresented GO categories compared to control within the allogeneic transplantation group. The top 10 significantly overrepresented GO categories on the third level (p < 0.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 Table 3. The numbers indicate the quantity of genes corresponding with the GO-terms.
aTX_down
Category
Term
Count
p-value
BP
nitrogen compound metabolism
119
8,1E-11
BP
catabolism
130
3,7E-6
BP
cellular metabolism
864
3,9E-4
BP
transport
408
5,1E-4
BP
establishment of localization
459
1,4E-3
BP
excretion
18
5,1E-3
CC
organelle inner membrane
110
1,7E-25
CC
mitochondrial envelope
116
1,1E-23
CC
organelle envelope
128
1,9E-15
CC
cytoplasm
695
5,8E-15
CC
organelle membrane
165
9,3E-15
CC
mitochondrial lumen
41
4,6E-7
CC
membrane fraction
171
8,8E-7
MF
electron carrier activity
45
7,6E-17
MF
oxidoreductase activity, acting on NADH or NADPH
37
1,8E-12
MF
primary active transporter activity
55
2,5E-8
MF
metal ion transporter activity
39
4,3E-7
MF
oxidoreductase activity, acting on the CH-CH group of donors
22
8,8E-6
MF
oxidoreductase activity, acting on CH-OH group of donors
38
1,1E-4
MF
coenzyme binding
21
2,0E-4
aTX_up
BP
immune response
283
1,7E-44
BP
defense response
292
8,1E-43
BP
response to pest, pathogen or parasite
180
6,4E-26
BP
response to other organism
180
2,3E-24
BP
macromolecule metabolism
797
1,4E-23
BP
cell activation
78
2,7E-14
BP
cell death
198
1,9E-13
BP
response to wounding
141
1,8E-11
BP
primary metabolism
1079
8,0E-9
BP
cell cycle
171
2,2E-8
CC
ribonucleoprotein complex
138
1,7E-22
CC
intracellular non-membrane-bound organelle
318
3,9E-18
CC
nucleus
581
3,0E-12
CC
ribosome
81
1,2E-11
CC
cytosolic small ribosomal subunit (sensu Eukaryota)
23
2,2E-8
CC
T cell receptor complex
14
1,5E-6
CC
nuclear envelope
44
2,0E-5
CC
intracellular organelle
1000
2,3E-5
CC
spliceosome complex
25
9,6E-5
CC
proteasome complex (sensu Eukaryota)
25
1,5E-4
MF
RNA binding
178
3,5E-31
MF
cytokine binding
30
8,5E-4
MF
sugar binding
36
6,8E-3
MF
translation factor activity, nucleic acid binding
33
8,6E-3
MF
translation initiation factor activity
22
1,3E-2
MF
protease activator activity
10
4,3E-2
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.
Table 4
Pathways with up- or down-regulated genes after syngeneic transplantation or allogeneic transplantation. Pathway charts were obtained using DAVID functional annotation tools (p < 0.05, Fisher's exact test, p-value with Benjamini correction). Pathways with a significantly elevated number of genes down- or up-regulated after syngeneic transplantation (sTX) or allogeneic transplantation (aTX) compared to control are shown Table 4. The numbers indicate the genes affected within the pathways.
sTX-down
sTX-up
Pathway
Count
Pathway
Count
tryptophan metabolism
13
ribosome
22
methionine metabolism
7
cell cycle
28
valine, leucine and isoleucine degradation
11
ribosome
13
lysine degradation
9
focal adhesion
31
fatty acid metabolism
11
ecm-receptor interaction
16
pyruvate metabolism
9
cell communication
16
beta-alanine metabolism
7
propanoate metabolism
8
butanoate metabolism
9
selenoamino acid metabolism
6
caprolactam degradation
5
alkaloid biosynthesis II
5
limonene and pinene degradation
5
glycine, serine and threonine metabolism
7
aTX-down
aTX-up
Pathway
Count
Pathway
Count
oxidative phosphorylation
48
antigen processing and presentation
46
fatty acid metabolism
29
ribosome
43
valine, leucine and isoleucine degradation
25
cell cycle
50
glutathione metabolism
21
proteasome
22
citrate cycle (tca cycle)
13
t cell receptor signaling pathway
45
pyruvate metabolism
17
type I diabetes mellitus
28
butanoate metabolism
20
hematopoietic cell lineage
33
glycine, serine and threonine metabolism
17
cell adhesion molecules (cams)
50
metabolism of xenobiotics by cytochrome p450
23
cytokine-cytokine receptor interaction
48
tryptophan metabolism
21
toll-like receptor signaling pathway
26
propanoate metabolism
14
reductive carboxylate cycle (CO2 fixation)
6
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 enrichment in up-regulated genes can be directly or indirectly linked to the immune response.
To illustrate the pathways activated by acute rejection related to the immune response, we used IPA to map networks 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 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 (24-fold) and TAP2 (15-fold; Fig. 2C).
Figure 3
Antigen presentation pathway. Ingenuity Pathways Knowledge Base was used to display pathways containing multiple differently expressed genes. Significant up-regulation of genes compared to control was observed in the canonical antigen presentation pathway after allogeneic transplantation. Up-regulated genes are highlighted in red the down-regulated genes are highlighted in green. The legend in Fig 3 applies also to Figs. 4-7.
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.
Figure 6
Interferon signaling. Significant up-regulation of genes was also observed in the canonical interferon signaling pathway. Up-regulated genes are highlighted in red.
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.
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.
Table 5
Potential novel transcriptional markers for acute rejection
phospholipase A2, group IIA (platelets, synovial fluid)
Pla2g2a
46.437
1369766_at
prostaglandin E receptor 2, subtype EP2
Ptger2
39.419
1392171_at
chitinase 3-like 1
Chi3l1
30.24
1370869_at
branched chain aminotransferase 1, cytosolic
Bcat1
26.936
1368413_at
amiloride binding protein 1
Abp1
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 up-regulation 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].
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 B-cell 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 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].
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 effector-related 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.
Methods
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 LEW-rats 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 [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 samples was used for cDNA synthesis with the High Capacity cDNA Reverse Transcription Kit. Real time PCR was performed using SYBR Green PCR Master Mix or TaqMan Universal PCR Master Mix on an ABI PRISM 7700 Sequence Detection System. Specific primer pairs or TaqMan Gene expression assays were used. The PCR efficiency of the primers pairs was tested in dilution series with a cDNA from the allogeneic transplanted group. A list with the corresponding primer sequences or TaqMan assay IDs is provided with this manuscript (additional file 3).
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].
Declarations
Acknowledgements
The authors thank Ute Neugebauer and Rita Schröter for excellent technical assistance.
This study was supported by grants from the Innovative Medical Research of the Medical Faculty of the University of Muenster (ED210404) to B.E. and from the Else-Kröner-Fresenius-Foundation P22/05//A43/05//F00 to E.S.
Daniel R Salomon and Sunil Kurian are supported by NIH grant U19 AI063603 and the Molly Baber Research Fund.
Authors’ Affiliations
(1)
Experimentelle Nephrologie, Medizinische Klinik und Poliklinik D, Universitätsklinikum Münster
(2)
Department of Molecular and Experimental Medicine, The Scripps Research Institute
(3)
Integrierte Funktionelle Genomik, Interdisziplinäres Zentrum für Klinische Forschung, Universitätsklinikum Münster
(4)
Hämatologie und Onkologie, Medizinische Klinik und Poliklinik A, Universitätsklinikum Münster
Kurian SM, Flechner SM, Kaouk J, Modlin C, Goldfarb D, Cook DJ, Head S, Salomon DR: Laparoscopic donor nephrectomy gene expression profiling reveals upregulation of stress and ischemia associated genes compared to control kidneys. Transplantation. 2005, 80: 1067-1071. 10.1097/01.tp.0000176485.85088.f7.PubMedView ArticleGoogle Scholar
Flechner SM, Kurian SM, Solez K, Cook DJ, Burke JT, Rollin H, Hammond JA, Whisenant T, Lanigan CM, Head SR, Salomon DR: De novo kidney transplantation without use of calcineurin inhibitors preserves renal structure and function at two years. Am J Transplant. 2004, 4: 1776-1785. 10.1111/j.1600-6143.2004.00627.x.PubMedView ArticleGoogle Scholar
Flechner SM, Kurian SM, Head SR, Sharp SM, Whisenant TC, Zhang J, Chismar JD, Horvath S, Mondala T, Gilmartin T, Cook DJ, Kay SA, Walker JR, Salomon DR: Kidney transplant rejection and tissue injury by gene profiling of biopsies and peripheral blood lymphocytes. Am J Transplant. 2004, 4: 1475-1489. 10.1111/j.1600-6143.2004.00526.x.PubMedPubMed CentralView ArticleGoogle Scholar
Kurian S, Grigoryev Y, Head S, Campbell D, Mondala T, Salomon DR: Applying genomics to organ transplantation medicine in both discovery and validation of biomarkers. Int Immunopharmacol. 2007, 7: 1948-1960. 10.1016/j.intimp.2007.07.017.PubMedPubMed CentralView ArticleGoogle Scholar
Sarwal M, Chua MS, Kambham N, Hsieh SC, Satterwhite T, Masek M, Salvatierra O: Molecular heterogeneity in acute renal allograft rejection identified by DNA microarray profiling. N Engl J Med. 2003, 349: 125-138. 10.1056/NEJMoa035588.PubMedView ArticleGoogle Scholar
Akalin E, Hendrix RC, Polavarapu RG, Pearson TC, Neylan JF, Larsen CP, Lakkis FG: Gene expression analysis in human renal allograft biopsy samples using high-density oligoarray technology. Transplantation. 2001, 72: 948-953. 10.1097/00007890-200109150-00034.PubMedView ArticleGoogle Scholar
Saiura A, Mataki C, Murakami T, Umetani M, Wada Y, Kohro T, Aburatani H, Harihara Y, Hamakubo T, Yamaguchi T, Hasegawa G, Naito M, Makuuchi M, Kodama T: A comparison of gene expression in murine cardiac allografts and isografts by means DNA microarray analysis. Transplantation. 2001, 72: 320-329. 10.1097/00007890-200107270-00027.PubMedView ArticleGoogle Scholar
Lande JD, Dalheimer SL, Mueller DL, Hertz MI, King RA: Gene expression profiling in murine obliterative airway disease. Am J Transplant. 2005, 5: 2170-2184. 10.1111/j.1600-6143.2005.01026.x.PubMedView ArticleGoogle Scholar
Einecke G, Melk A, Ramassar V, Zhu LF, Bleackley RC, Famulski KS, Halloran PF: Expression of CTL associated transcripts precedes the development of tubulitis in T-cell mediated kidney graft rejection. Am J Transplant. 2005, 5: 1827-1836. 10.1111/j.1600-6143.2005.00974.x.PubMedView ArticleGoogle Scholar
Famulski KS, Einecke G, Reeve J, Ramassar V, Allanach K, Mueller T, Hidalgo LG, Zhu LF, Halloran PF: Changes in the transcriptome in allograft rejection: IFN-gamma-induced transcripts in mouse kidney allografts. Am J Transplant. 2006, 6: 1342-1354. 10.1111/j.1600-6143.2006.01337.x.PubMedView ArticleGoogle Scholar
Stegall M, Park W, Kim D, Kremers W: Gene expression during acute allograft rejection: novel statistical analysis of microarray data. Am J Transplant. 2002, 2: 913-925. 10.1034/j.1600-6143.2002.21007.x.PubMedView ArticleGoogle Scholar
Erickson LM, Pan F, Ebbs A, Kobayashi M, Jiang H: Microarray-based gene expression profiles of allograft rejection and immunosuppression in the rat heart transplantation model. Transplantation. 2003, 76: 582-588. 10.1097/01.TP.0000081398.65568.1B.PubMedView ArticleGoogle Scholar
Velic A, Hirsch JR, Bartel J, Thomas R, Schroter R, Stegemann H, Edemir B, August C, Schlatter E, Gabriels G: Renal transplantation modulates expression and function of receptors and transporters of rat proximal tubules. J Am Soc Nephrol. 2004, 15: 967-977. 10.1097/01.ASN.0000117287.74203.89.PubMedView ArticleGoogle Scholar
Velic A, Gabriels G, Hirsch JR, Schroter R, Edemir B, Paasche S, Schlatter E: Acute rejection after rat renal transplantation leads to downregulation of NA+ and water channels in the collecting duct. Am J Transplant. 2005, 5: 1276-1285. 10.1111/j.1600-6143.2005.00890.x.PubMedView ArticleGoogle Scholar
Hori S, Nomura T, Sakaguchi S: Control of Regulatory T Cell Development by the Transcription Factor Foxp3. Science. 2003, 299: 1057-1061. 10.1126/science.1079490.PubMedView ArticleGoogle Scholar
Kristiansen M, Graversen JH, Jacobsen C, Sonne O, Hoffman HJ, Law SK, Moestrup SK: Identification of the haemoglobin scavenger receptor. Nature. 2001, 409: 198-201. 10.1038/35051594.PubMedView ArticleGoogle Scholar
Runza VL, Hehlgans T, Echtenacher B, Zahringer U, Schwaeble WJ, Mannel DN: Localization of the mouse defense lectin ficolin B in lysosomes of activated macrophages. J Endotoxin Res. 2006, 12: 120-126. 10.1179/096805106X102147.PubMedView ArticleGoogle Scholar
Olson TS, Ley K: Chemokines and chemokine receptors in leukocyte trafficking. Am J Physiol Regul Integr Comp Physiol. 2002, 283: R7-28.PubMedView ArticleGoogle Scholar
Dennis G, Sherman B, Hosack D, Yang J, Gao W, Lane HC, Lempicki R: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biology. 2003, 4: 3-10.1186/gb-2003-4-5-p3.View ArticleGoogle Scholar
Kloetzel PM: Antigen processing by the proteasome. Nat Rev Mol Cell Biol. 2001, 2: 179-187. 10.1038/35056572.PubMedView ArticleGoogle Scholar
Momburg F, Hammerling GJ: Generation and TAP-mediated transport of peptides for major histocompatibility complex class I molecules. Adv Immunol. 1998, 68: 191-256.PubMedView ArticleGoogle Scholar
Carlsson R, Persson C, Leanderson T: SPI-C, a PU-box binding ETS protein expressed temporarily during B-cell development and in macrophages, contains an acidic transactivation domain located to the N-terminus. Molecular Immunology. 2003, 39: 1035-1043. 10.1016/S0161-5890(03)00032-4.PubMedView ArticleGoogle Scholar
Carlsson R, Thorell K, Liberg D, Leanderson T: SPI-C and STAT6 can cooperate to stimulate IgE germline transcription. Biochemical and Biophysical Research Communications. 2006, 344: 1155-1160. 10.1016/j.bbrc.2006.04.026.PubMedView ArticleGoogle Scholar
Lindbom J, Ljungman AG, Lindahl M, Tagesson C: Increased gene expression of novel cytosolic and secretory phospholipase A(2) types in human airway epithelial cells induced by tumor necrosis factor-alpha and IFN-gamma. J Interferon Cytokine Res. 2002, 22: 947-955. 10.1089/10799900260286650.PubMedView ArticleGoogle Scholar
Sherman LA, Chattopadhyay S: The molecular basis of allorecognition. Annu Rev Immunol. 1993, 11: 385-402. 10.1146/annurev.iy.11.040193.002125.PubMedView ArticleGoogle Scholar
Gould DS, Auchincloss H: Direct and indirect recognition: the role of MHC antigens in graft rejection. Immunol Today. 1999, 20: 77-82. 10.1016/S0167-5699(98)01394-2.PubMedView ArticleGoogle Scholar
Graubert TA, Ley TJ: How do lymphocytes kill tumor cells?. Clin Cancer Res. 1996, 2: 785-789.PubMedGoogle Scholar
York IA, Goldberg AL, Mo XY, Rock KL: Proteolysis and class I major histocompatibility complex antigen presentation. Immunol Rev. 1999, 172: 49-66. 10.1111/j.1600-065X.1999.tb01355.x.PubMedView ArticleGoogle Scholar
Ma W, Lehner PJ, Cresswell P, Pober JS, Johnson DR: Interferon-gamma Rapidly Increases Peptide Transporter (TAP) Subunit Expression and Peptide Transport Capacity in Endothelial Cells. J Biol Chem. 1997, 272: 16585-16590. 10.1074/jbc.272.26.16585.PubMedView ArticleGoogle Scholar
Hancock WW, Gao W, Faia KL, Csizmadia V: Chemokines and their receptors in allograft rejection. Curr Opin Immunol. 2000, 12: 511-516. 10.1016/S0952-7915(00)00130-8.PubMedView ArticleGoogle Scholar
Christopher K, Mueller TF, DeFina R, Liang Y, Zhang J, Gentleman R, Perkins DL: The graft response to transplantation: a gene expression profile analysis. Physiol Genomics. 2003, 15: 52-64.PubMedView ArticleGoogle Scholar
Mueller TF, Ma C, Lederer JA, Perkins DL: Differentiation of stress, metabolism, communication, and defense responses following transplantation. J Leukoc Biol. 2003, 73: 379-390. 10.1189/jlb.0102052.PubMedView ArticleGoogle Scholar
Uren AG, Pakusch M, Hawkins CJ, Puls KL, Vaux DL: Cloning and expression of apoptosis inhibitory protein homologs that function to inhibit apoptosis and/or bind tumor necrosis factor receptor-associated factors. PNAS. 1996, 93: 4974-4978. 10.1073/pnas.93.10.4974.PubMedPubMed CentralView ArticleGoogle Scholar
Duckett CS, Li F, Wang Y, Tomaselli KJ, Thompson CB, Armstrong RC: Human IAP-Like Protein Regulates Programmed Cell Death Downstream of Bcl-xL and Cytochrome c. Mol Cell Biol. 1998, 18: 608-615.PubMedPubMed CentralView ArticleGoogle Scholar
Soares MP, Brouard S, Smith RN, Bach FH: Heme oxygenase-1, a protective gene that prevents the rejection of transplanted organs. Immunol Rev. 2001, 184: 275-285. 10.1034/j.1600-065x.2001.1840124.x.PubMedView ArticleGoogle Scholar
Walunas TL, Lenschow DJ, Bakker CY, Linsley PS, Freeman GJ, Green JM, Thompson CB, Bluestone JA: CTLA-4 can function as a negative regulator of T cell activation. Immunity. 1994, 1: 405-413. 10.1016/1074-7613(94)90071-X.PubMedView ArticleGoogle Scholar
Fiorentino DF, Bond MW, Mosmann TR: Two types of mouse T helper cell. IV. Th2 clones secrete a factor that inhibits cytokine production by Th1 clones. J Exp Med. 1989, 170: 2081-2095. 10.1084/jem.170.6.2081.PubMedView ArticleGoogle Scholar
Pestka S, Krause CD, Sarkar D, Walter MR, Shi Y, Fisher PB: Interleukin-10 and related cytokines and receptors. Annu Rev Immunol. 2004, 22: 929-979. 10.1146/annurev.immunol.22.012703.104622.PubMedView ArticleGoogle Scholar
Walport MJ: Complement. Second of two parts. N Engl J Med. 2001, 344: 1140-1144. 10.1056/NEJM200104123441506.PubMedView ArticleGoogle Scholar
Walport MJ: Complement. First of two parts. N Engl J Med. 2001, 344: 1058-1066. 10.1056/NEJM200104053441406.PubMedView ArticleGoogle Scholar
Sacks S, Zhou W: The effect of locally synthesised complement on acute renal allograft rejection. J Mol Med. 2003, 81: 404-410. 10.1007/s00109-003-0454-7.PubMedView ArticleGoogle Scholar
Fujita T, Gigli I, Nussenzweig V: Human C4-binding protein. II. Role in proteolysis of C4b by C3b-inactivator. J Exp Med. 1978, 148: 1044-1051. 10.1084/jem.148.4.1044.PubMedView ArticleGoogle Scholar
Liszewski MK, Farries TC, Lublin DM, Rooney IA, Atkinson JP: Control of the complement system. Adv Immunol. 1996, 61: 201-283.PubMedView ArticleGoogle Scholar
Zelenika D, Adams E, Humm S, Lin CY, Waldmann H, Cobbold SP: The role of CD4+ T-cell subsets in determining transplantation rejection or tolerance. Immunological Reviews. 2001, 182: 164-179. 10.1034/j.1600-065X.2001.1820113.x.PubMedView ArticleGoogle Scholar
Gabriels G, August C, Grisk O, Steinmetz M, Kosch M, Rahn KH, Schlatter E: Impact of renal transplantation on small vessel reactivity. Transplantation. 2003, 75: 689-697. 10.1097/01.TP.0000044111.12370.ED.PubMedView ArticleGoogle Scholar
Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucl Acids Res. 2002, 30: 207-210. 10.1093/nar/30.1.207.PubMedPubMed CentralView ArticleGoogle Scholar
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene Ontology: tool for the unification of biology. Nat Genet. 2000, 25: 25-29. 10.1038/75556.PubMedPubMed CentralView ArticleGoogle Scholar
Benjamini Y, Y. HOCHBERG: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B. 1995, 85: 289-300.Google Scholar
Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
