|
|
||||||||

* Laboratory of Molecular Immunology, Department of Medicine, and
Laboratory of Immunology, Department of Surgery, Brigham and Womens Hospital and Harvard Medical School, Boston, Massachusetts
Correspondence: Dr. David L Perkins, Laboratory of Molecular Immunology, Department of Medicine, Brigham and Womens Hospital and Harvard Medical School, 75 Francis Street, Boston, MA 02115. E-mail: dperkins{at}rics.bwh.harvard.edu
| ABSTRACT |
|---|
|
|
|---|
Key Words: molecular biology gene regulation rodent
| INTRODUCTION |
|---|
|
|
|---|
To study these complex biologic phenomena, we used DNA microarrays with results confirmed by real-time polymerase chain reaction (PCR) to compare patterns and kinetics of gene-expression profiles of ischemia, stress, and innate and adaptive immune responses following transplantation. Considering complexity, pleiotropy, and redundancy in the immune system, our objective was not to filter for a single gene or gene family but to analyze global patterns of gene expression. We used hierarchical clustering algorithms and self-organizing maps (SOMs) to identify patterns of gene expression. In addition, the Gene Ontology (GO) database was used to assign core biological processes associated with the genes differentially regulated following transplantation [7 , 8 ].
It is interesting that our analyses demonstrate disproportionate expression of genes from categories of biological processes modulated by ischemia, systemic stress, and innate and adaptive immunity.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Vascularized heterotopic cardiac transplantation
Murine hearts were transplanted as described previously [9
]. Briefly, hearts were harvested from freshly killed donors and immediately transplanted into 8- to 12-week-old recipients anaesthetized via intraperitoneal injection with 60 mg/kg pentobarbital sodium. The donor aorta was attached to the recipient abdominal aorta by end-to-side anastomosis, and the donor pulmonary artery was attached to the recipient vena cava by end-to-side anastomosis. All surgical procedures were completed in less than 60 min from the time that the donor heart was harvested. Donor hearts that did not beat immediately after reperfusion or stopped within 1 day following transplantation were excluded (>95% of all grafts functioned at 1 day following transplantation). The native heart of the recipient was not surgically manipulated and remained functional. Donor grafts were harvested at 24 h or 7 days following transplantation and were divided into equal sections for preparation of RNA and tissue sections for histology. Altogether, organs from 37 animals were harvested.
DNA microarrays
Experimental procedures for the microarrays (GeneChip) were performed according to the Affymetrix (Santa Clara, CA) GeneChip Expression Analysis technical manual. Briefly, double-stranded cDNA was synthesized by means of the SuperScript Choice system (Gibco-BRL, Life Technologies, Rockville, MD) and a T7(dT) 24 primer (Gensetoligos, La Jolla, CA). The cDNA was purified using phenol/chloroform extraction with Phase Lock Gel (Eppendorf, Germany) and was concentrated by ethanol precipitation. In vitro transcription was performed to produce biotin-labeled cRNA using a BioArray HighYield RNA transcript labeling kit (Enzo Diagnostics, Farmingdale, NY) according to the manufacturers instructions. cRNA was linearly amplified
40-fold with T7 polymerase using double-stranded cDNA, which was synthesized. The biotinylated RNA was cleaned with an RNeasy mini kit (Qiagen, Valencia, CA), fragmented to 50200 nt, and then hybridized to an Affymetrix murine array (Mu11kB), which contains probe sets for 6500 genes and expressed sequence tags (ests). After being washed, the array was stained with streptavidin-phycoerythrin (Molecular Probes, Eugene, OR), amplified by biotinylated antistreptavidin (Vector Laboratories, Burlingame, CA), and then scanned on a HP Genearray scanner. The intensity for each feature of the array was captured with Affymetrix GeneChip software, according to standard Affymetrix procedures. Correlation coeffecient of duplicate control heart microarray results was r = 0.98.
Statistics and data analysis
Array data were analyzed with Microarray Suite 4.0.1. A single expression level for each gene was derived from the 20 probe pairs representing each gene, 20 perfectly matched (PM) and mismatched (MM) control probes. The MM probes act as specificity control, which allows the direct subtraction of background and cross-hybridization signals. Each array was normalized to a standard of 2500 units per probe set. To determine the quantitative RNA level, the average of the differences representing PM MM for each gene-specific probe set was calculated. In addition the expression of each probe set was qualitatively categorized as present, marginal, or absent. Calculations of means and variances were performed with JMP statistical software (SAS Institute, Cary, NC).
Cluster analysis
Hierarchical clustering, using Cluster and TreeView software [10
] (courtesy of Michael B. Eisen, Lawrence Livermore Radiation Laboratory, Berkeley, CA), and SOMs, using GeneCluster software [11
] (courtesy of Whitehead Institute for Biomedical Institute, Cambridge, MA), were applied to analyze gene expression.
Average difference values were analyzed by Cluster using the hierarchical clustering algorithm with average linkage clustering. Briefly, dissimilarity of gene expression between each experimental group is calculated by Pearson correlation relationships and visualized by a tree diagram (dendrogram) [10 ]. Branch lengths represent the degree of similarity between the groups.
A second algorithm, SOMs, was used to cluster genes with similar expression patterns across all experimental groups. Genes showing a less than twofold change in relative expression across the experimental groups were eliminated, and the expression level of each gene was normalized across the experiments. The remaining genes were analyzed by GeneCluster using a 2 x 3 geometry, and six clusters were generated, eliminating clusters with few genes or large standard deviations [11 ]. Each cluster contained genes with similar patterns of gene expression across the experimental groups. Based on multiple heuristic observations, increased numbers of nodes produced clusters with low numbers of genes, whereas decreased numbers of nodes produced larger SD. Increasing the number of epochs (=500) did not produce detectable changes in the clusters or SD.
Real-time RT-PCR
Primer pairs were designed using Primer Express software (Applied Biosystems, Foster City, CA). Forward and reverse primers were chosen to have a length between 18 and 22 base pairs and were designed to amplify an amplicon length of 51 base pairs. The housekeeping gene glyceraldehyde 3-phosphate dehydrogenase (GAPDH; M32599) was selected as endogeneous control. All primer pairs were tested in immune-rich tissue samples and nontemplate controls for specificity, primerdimer conformation, and reproducibility. For the real-time RT-PCR analyses, a coefficient of variance between 0.2 and 1.2% was obtained measuring GAPDH in different runs of identical replicates in the various groups. The sequences of the forward (FW) and reverse (RE) primer pairs used in the experiments are as follows: major histocompatibility complex (MHC) class II invariant (Ii; FW) GAAGCTTCCGAAATCTGCCA, (RE) GGAGTAGCCATCCGCATCTG; MHC class II antigen (I-A-b; FW) TTCATGGGCGAGTGCTACTTC, (RE) CACATATCGTATGCGCTGCG; C1q complement (C1q)-A (FW) CTGGCATCCGGACTGGTATC, (RE) CAGATTCCCCTGGGTCTCCT; C1q-B (FW) AGGACCATCAACAGCCCCTT, (RE) CTTTTCGAAGCGAATGACCTG; C1q-B (FW) TGTGGAGGGCCGATACAAAC, (RE) CGGGTGACTGTGAATACCGA; acute-phase proteins (APP; FW) GGGCTGACAAACATCAAGACG, (RE) TGCATCCATCTTCACTTCCG; lymphocyte function-associated antigen-1 (LFA-1; FW) CAGCCATCTGCCTATGACCA, (RE) CCCAACCAAGGCCTCTAGTGT; and GAPDH (FW) TGTGGAAGGGCTCATGACC, (RE) TCTTCTGGGTGGCAGTGATG.
RNA for the real-time RT-PCR studies was obtained from separate tissue samples to control for technical and biological variability. Total murine RNA was isolated from three hearts per experimental group comprising untransplanted BALB/c hearts (control group), day 1 B6 graft hearts (syngeneic group), day 1 allogeneic BALB/c-RAG-/- graft hearts (RAG-deficient group), and days 1 and 7 BALB/c graft hearts (allogeneic group). TRI reagent (Sigma-Aldrich, St. Louis, MO) was used for the RNA isolation, and samples were treated with DNase to eliminate DNA (DNase I, Amplification Grade, Invitrogen Life Technologies, Carlsbad, CA). The RNA was reverse-transcribed using SuperScript II RNase RT (Gibco, Carlsbad, CA). The GeneAmp 5700 sequence detection system (Applied Biosystems) was used to perform real-time PCR using 250 ng template cDNA, 5 µM FW and RE primer, and 10 µL 10x SYBR Green PCR Master Mix (Applied Biosystems) per well in a MicroAmp Opitcal 96-well reaction plate (Applied Biosystems). The gene-specific PCR products are continuously measured by the increase in fluorescence as a result of the binding of SYBR Green to double-stranded DNA during 40 cycles. Dye ROX, included in the SYBR Green PCR Master Mix, was used as a passive reference to normalize for non-PCR-related fluctuations in fluorescence signal. The cycle thresholds (CT) of each target gene are relative to the CT value of the endogeneous reference, GAPDH. The differences in the CT values between GAPDH and each single target gene, respectively, are calculated as fold-change in expression relative to GAPDH.
Gene classification
All genes selected by the SOMs were classified according to the biological processes to which they contribute using the GO (<www.geneontology.org>) annotation system [7
]. The GO identifiers were obtained from the mouse genome informatics database [12
]. The biological process is defined as the biological objective to which the gene contributes. Most genes could be assigned to the following five broad processes: cell communication (communication; GO:0007154), cell organization and biogenesis (organization; GO:0007010), metabolism (GO:0008152), defense response (defense; GO:0006952), and stress response (stress; GO:0006950). The remaining genes, in particular those with unknown biological process, were classified under "other". If a gene contributed to more than one distinct biological process, the most appropriate to the underlying experimental settings was chosen. A detailed description of each gene, including identifier, standard name, characterization, and relevant reference as well as related biological process according to the GO classification, is provided on our website (<perkinslab.bwh.harvard.edu>).
| RESULTS |
|---|
|
|
|---|
Cluster analysis of experimental groups at days 1 and 7
To determine the global relationships among the experimental groups at day 1, we analyzed expression profiles from DNA microarrays of allogeneic graft heart, alymphoid graft heart, native heart, ischemic heart, and untransplanted control heart using agglomerative, hierarchical clustering algorithms that depict the degree of dissimilarity in dendrograms [10
] (Fig. 1a
). Also, we analyzed untransplanted control LN to include an immune gene-rich sample. As shown in the dendrogram, the three untransplanted groups (control, native, and ischemic hearts) were the most similar; also, the two transplanted groups (allogeneic day 1 and alymphoid day 1) were only modestly dissimilar from each other. As expected, the alymphoid and syngeneic day 1 groups were highly similar (not shown). Gene expression in the LN was the most dissimilar compared with all five heart groups. These observations show that the changes in gene expression at day 1 following transplantation are not highly dissimilar in the allogeneic versus alymphoid group, suggesting that gene regulation is predominantly a result of antigen-independent mechanisms.
|
Gene modulation associated with antigen-independent mechanisms
Genes (181) were identified using a selection threshold of a greater than twofold change in relative expression across the experimental groups. These genes are listed in rank order by the absolute magnitude of increased or decreased expression, respectively (Table 1
). Supplementary information to each of the 181 modulated genes is available at our website (<perkinslab.bwh.harvard.edu>).
|
|
|
Genes modulated by innate immunity
Our experimental design includes two groups that lack an alloantigenic stimulus: the syngeneic group (B6 into B6) and the alymphoid group (BALB/c-RAG-/- into B6-RAG-/-). Cluster 4 of the SOM shows a subset of genes up-regulated in the syngeneic and alymphoid groups. In addition, cluster 5 includes genes up-regulated in the ischemic and native hearts, as well as the syngeneic and alymphoid graft hearts (Fig. 2)
. A compilation of genes based on the same selection threshold up-regulated in the syngeneic and alymphoid groups but not in the ischemic or native hearts generates 38 genes (Table 1)
. Based on the GO characterization, these genes predominantly belong to defense responses (see Fig. 4C ). They can be subclassified into immune response (for example, cytokines) or acute phase and inflammatory responses (for example, serum amyloid A 3). In addition, there is a large increase in the percentage of genes associated with communication, which includes predominantly adhesion molecules. Strikingly, 16% of all genes up-regulated in the innate experimental model were mediating cell-adhesion processes.
Genes modulated in the alloantigen-dependent, adaptive model
Based on our analysis with the clustering algorithms, the proinflammatory responses during the early phase following transplantation (day 1) were similar in syngeneic, alymphoid, and allogeneic recipients. In contrast, during the late phase (day 7), the allogeneic grafts produced a distinct profile (Fig. 1)
. To investigate the antigen-dependent, adaptive immune response during rejection, we analyzed allogeneic groups from days 1 and 7 following transplantation using SOMs (Fig. 3
). In this analysis, we also included untransplanted control heart, alymphoid graft heart, and LN. The analysis with SOMs shows the absence of a cluster-containing genes up-regulated exclusively at day 1 in the allogeneic group, which is consistent with the dendrograms produced by hierarchical clustering showing similarity among all of the groups of transplanted hearts at day 1. In contrast, a unique subset of genes was identified with modulated expression at day 7 in the allogeneic group. Specifically, cluster 2 contains genes up-regulated at day 7 in the allogeneic group and also highly expressed in LN. In addition, cluster 4 contains genes up-regulated at day 7 but not in LN or other groups. Using the same threshold criteria, compilation of genes significantly modulated in the allogeneic group at day 7 but not in the ischemic, native heart, or alymphoid groups produced 49 up-regulated and 11 down-regulated genes (Table 1)
. Based on the GO classification, the majority of genes in this subset is associated with the defense response (Fig. 4D
). The 20 genes showing the highest up-regulation include mostly genes linked to the MHC or TCR complexes or to the complement cascade (Table 1)
. Importantly, the percentage of the defense response is greatest in the adaptive group. It is interesting that in the innate and adaptive model adhesion molecules, classified in the segment communication, constitute a large group. Conversely, the metabolic and stress response components are reduced in the adaptive group.
|
|
| DISCUSSION |
|---|
|
|
|---|
To identify specific subsets of modulated genes with similar expression patterns during the late or adaptive phase of rejection, we used SOMs (Figs. 2 and 3) . This clustering technique has been successfully applied to high-dimensional and nonlinear problems in many areas of science. Supplementary information on all selected genes of the experimental groups including accession number, characterization, and relevant references is provided in our website (<perkinslab.bwh.harvard.edu>). Specific genes identified with the SOMs included five MHC genes, four antigen receptor genes, four complement genes, three antigen-processing genes, multiple cell-surface markers, and adhesion molecules. In addition, a T cell surface marker (Thy-1) and an IFNIP-30 were detected [14 ]. This subset indicates that a major component of the late response belongs to adaptive immunity.
The method of real-time RT-PCR was used as a validation method, and specifically, genes up-regulated during this late adaptive response were selected for corroboration. Despite the methodological differences between PCR amplification and microarray hybridization and the use of different RNA samples, the seven target genes show comparable expression profiles for both methods. In particular, genes associated with the adaptive immune response display the highest degree of up-regulation during the allogeneic response. In addition, confirming the findings by cluster analysis and GO attribution, this allogeneic response is associated with increased expression of markers associated with innate and inflammatory processes such as complement factors.
Our analysis of the early response following transplantation focused on three components including ischemia, systemic stress, and innate immunity. The ischemic response, which was induced in vitro to eliminate confounding effects of systemic mediators or infiltrating cells, contained 63 genes that were identified using SOMs. It is interesting that five genes (elongation factor, gelsolin, ATP synthase, crystallin, and ubiquitin) up-regulated by ischemia were also up-regulated in a previous report analyzing congestive heart failure [15 ]. Importantly, most of the additional genes identified in our ischemia group have previously reported functions consistent with a functional role during ischemia or injury induced by ischemia. For example, a recent study demonstrates that hypoxia up-regulates expression of ferritin by decreasing promoter binding of the inhibitory iron-regulatory protein [16 ]. MnSOD, which protects against reactive oxygen species, has been shown to be transcriptionally up-regulated by hypoxia and lipopolysaccharides (LPS) [17 , 18 ]. C1 esterase inhibitor is protective against myocardial ischemia and cold ischemia during organ preservation [19 , 20 ]. Crystallin, related to the small hsp, including hsp27, is induced by stress and is highly expressed in encysted artemia franciscana (brine shrimp), which can survive 3 years of anoxia [21 , 22 ]. In addition, nine genes highly expressed in cardiac tissue were down-regulated in response to ischemia.
Our experimental model included a group to analyze gene expression by systemic stress. The graft heart is transplanted heterotopically into the abdomen, whereas the recipients native heart is not surgically manipulated but is exposed to the systemic stresses of the transplant procedure. Our analysis of the native heart identified 21 modulated genes using SOMs. Supporting the conclusion that these genes were, at least in part, regulated by stress is the observation that multiple genes in this subset are regulated by other types of stress. For example, MT-I and SelP are up-regulated by exposure to heavy metals. It is interesting that MT has been shown to provide protection against electrophiles and oxidant stress and to be induced during acute liver injury and regeneration [23 , 24 ]. The ryanodine receptor, a Na/Ca exchanger, binds calmodulin and FKBP12 and is sensitive to redox status and mechanical stress [25 , 26 ]. 24p3, a gelatinase-associated lipocalin, is an acute-phase response protein that binds formylpeptides and LPS and has been suggested to function as a regulator of inflammation [27 ]. EF-2 is up-regulated by angiotensin II or ß-adrenergic agonists, both of which could be released by the transplant surgery [28 , 29 ]. 11 ß-Hydroxysteroid dehydrogenase is an important deactivator of glucocorticoids [30 ], and 1-8U is an IFN-inducible protein [31 ]. CD36, a lipoprotein receptor expressed by multiple cell types including endothelial and red blood cells, can mediate cell adherence via thrombospondin [32 , 33 ]. We detected only a single est showing down-regulation in response to stress. Together, these observations demonstrate a major stress response following transplantation that is independent of ischemia and immunity.
In addition to ischemia and systemic stress, our experimental design analyzed antigen-independent modulation of gene expression in two groups without antigen-dependent responses, the syngeneic and alymphoid (RAG-/-) groups, which were operationally defined as innate groups. Our analysis identified 38 genes up-regulated by antigen-independent signals but independent of ischemia or systemic stress. In this subset, IFN-ß, a cytokine shown to be important in innate-immune responses, was highly up-regulated [34
]. Serglycin, a widely distributed proteoglycan highly produced by endothelial cells, binds the leukocyte receptor CD44 [35
, 36
]. Osteopontin, a cytokine produced by multiple cell types including chrondocytes, has been shown to promote T helper cell type 1 responses [37
]. Serum amyloid A 3, up-regulated by glucocorticoids, is a mediator of the acute-phase response [38
]. Prothymosin
, which has receptors on monocytes, promotes cellular proliferation and enhances protective immune responses induced by vaccination [39
, 40
]. Apo E is a glycosylated protein with multiple functions. It is interesting that mice deficient in Apo E have increased susceptibility to LPS endotoxemic shock [41
]. Ly-6C, which is expressed on multiple cell types including macrophages, endothelial, natural killer (NK), and T cells, has been associated with NK cell cytotoxicity [42
43
44
45
]. Cofilin is an actin-binding protein important in actin-filament depolymerization crucial in motility and capping [46
]. The integrin ß subunit plays a fundamental role in cell division, differentiation, and cell migration by mediating cellmatrix and cellcell interactions [47
]. Calpactin I (annexin II) is induced by transforming growth factor-ß (TFG-ß) and is differentially expressed during wound repair [48
]. hsp84 and hsp47 are stress-induced proteins. It is interesting that hsp47 binds collagen and fetuin and is associated with fibrosis and atheromatous plaques [49
, 50
]. Endoglin, a TGF-ß-binding protein, is up-regulated by TGF-ß and overexpressed after arterial injury [51
, 52
]. Activated macrophages can up-regulate the production of fibronectin, which has been shown to modulate TGF-ß induction of IL-8 [53
, 54
]. TGF-ß stimulates TSC-22, a DNA-binding protein of the kruppel family [55
]. We detected 10 ests but no known genes down-regulated in the innate-immune group. These results do not identify the specific triggers of the innate response, which include multiple, potential candidates, for example, complement activation, necrotic or apoptotic cells, and activation of thrombosis among others. However, our results demonstrate that transplantation injury, even in the absence of infection, can stimulate a robust innate-immune response.
Our study analyzed the expression of
6500 genes and ests following transplantation and identified 181 genes with significantly modulated expression in at least one of the experimental groups of ischemia, stress, and innate or adaptive immunity. To understand the biological importance of changes in gene expression, we used the GO database to categorize the biological process of all modulated genes [7
, 8
, 10
]. Remarkably, greater than 80% of all modulated genes were contained within only five categories of biological processes, including defense, stress, metabolism, communication, and cell organization. Supporting the validity of this approach, the relative contribution of each biological process changed in the experimental groups consistent with our hypotheses based on previous studies. For example, greater than 70% of genes modulated during ischemia were categorized in metabolism and cell organization, whereas only 3% were attributed to defense. In the stress group, there was a major increase in the stress (15%) and defense (15%) categories, and metabolism (40%) also remained a major category. In the innate response, the defense category (24%) increased further, whereas metabolism (18%) decreased. It is interesting that the communication category (18%), which includes genes of adhesion molecules and signal transducers, markedly increased during the innate response. This response is consistent with a critical role by innate immunity in the activation and recruitment of an inflammatory response. In the adaptive response,
50% of the modulated genes are within the defense category. The categories of metabolism (18%), cell organization (8%), and communication (13%) remain substantial components, suggesting a significant contribution to an antigen-specific response.
The profiles of gene expression and associated biological processes provide a basis to differentiate major biological responses in transplantation. Despite the remarkable power of profiling gene expression with microarrays, some limitations are noteworthy. The DNA microarrays cover only a subset of the whole murine genome, the molecular and biological function of the ests is not known, and post-transcriptional changes are not detected. In particular, the genes and biological processes represented on the individual chip do not per se reflect their true distribution in the genome. In this respect, information provided by the gene chip producers on the relative distribution of genes representing the different biological processes would be most helpful. Using the GO annotations, approximately 6% of all genes on the Affymetrix array (Mu11kB) used in our study are related to the defense response, comprising immune and inflammatory genes. During the late allogeneic response, approximately 50% of all genes up-regulated were associated with the defense response, indicating a marked increase compared with the baseline of 6% on the whole chip. Thus, the underlying analysis shows characteristic shifts in the biological responses, but the inherent bias of gene selection by the individual chips always has to be taken into account.
We propose that the presented approach provides a framework for future studies investigating the complex responses induced by organ transplantation.
Understanding the interactions between the defense category, which is commonly the focus of transplant biologists, and other biological processes of metabolism, cell organization, and communication may increase our understanding of the rejection response in vivo and provide new targets for novel, therapeutic strategies.
| ACKNOWLEDGEMENTS |
|---|
Received January 30, 2002; revised September 9, 2002; accepted September 24, 2002.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
V. J. Gimino, J. D. Lande, T. R. Berryman, R. A. King, and M. I. Hertz Gene Expression Profiling of Bronchoalveolar Lavage Cells in Acute Lung Rejection Am. J. Respir. Crit. Care Med., November 15, 2003; 168(10): 1237 - 1242. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |