Certain cytokine gene polymorphisms have been shown to correlate with outcome of human leukocyte antigen (HLA) identical sibling donor stem cell transplantation (SCT), but in unrelated donor SCT such information is scarce. We have studied the association between cytokine gene polymorphism and transplant-related mortality (TRM) in 182 unrelated SCTs performed at a single center. We found association of polymorphism in the tumor necrosis factor alpha (TNFα) and interleukin-10 (IL-10) gene and TRM. Both the TNFd4 allele and the TNFα -1031C alleles are associated with high risk for TRM. Statistical analysis showed that both polymorphisms were present on a single haplotype. This haplotype was associated with high risk of TRM when present in recipient or donor, 55% (43%-67%) compared with 21% (12%-30%) when absent from both (P < .01). A further allele associated with this haplotype, TNFa5, is also associated with increased risk of TRM. For IL-10, presence of the donor R2-G-C-C haplotype was associated with decreased risk of TRM, 61% (43%-79%) versus 34% (25%-43%), P = .01. In contrast, possession of the R3-G-C-C haplotype by the donor predicted reduced risk of TRM, 30% (19%-41%, 95% CI) versus 53% (40%-66%, 95% CI), P = .01. No independent associations of cytokine polymorphisms with acute graft-versus-host disease were shown. (Blood. 2004;103:3599-3602)

Tumor necrosis factor alpha (TNFα) is a member of the TNF receptor superfamily of ligand receptors and, as a T-helper cell 1 (Th1) cytokine, is an important component of the proinflammatory immune response. Intensive chemotherapy and irradiation given prior to stem cell transplantation (SCT) is associated with increased serum levels of TNFα.1-3  High levels of serum TNFα measured before transplantation are associated with increased probability of early transplant-related mortality (TRM) and acute graft-versus-host disease (aGVHD).2,4 

Interleukin-10 (IL-10) is a member of an expanding family of cytokines known as the interferon receptor superfamily.5  It is an important anti-inflammatory cytokine in a number of immune pathways. A reduction in TRM after SCT has been associated with increased recipient levels of IL-10 immediately before and after transplantation; however, this observation remains controversial.4,6,7 

The genes encoding TNFα (TNFSF2) and IL-10 (IL10) are both highly polymorphic. The TNFSF2 gene maps to chromosome 6 (6p21.3) and contains a number of single nucleotide polymorphisms (SNPs), predominantly in the 5′ promoter region.8 TNFSF2 gene is surrounded by a number of other genes with immunologic functions. A number of well-studied microsatellites have been described within genes flanking TNFSF2.8  TNFa, TNFb, and TNFc microsatellites are found in or near the lymphotoxin alpha gene, TNFSF1. TNFd and TNFe microsatellites are located in intron 3 of the leukocyte-specific transcript gene, LST1. The IL10 gene maps to chromosome 1 (1q31-32). The 5′ promoter region of this gene contains a number of proximal9  and distal10  SNPs as well as 2 microsatellites termed IL-10R and IL-10G.11,12  In both TNFSF2 and IL10 genes investigators have sought polymorphic profiles that can be defined as haplotypes.11,13,14 

Polymorphisms in the TNFSF2 region and IL10 genes have been associated with clinical outcome in several identical sibling stem cell transplantation (ID-SCT) studies.12,15,16  However for unrelated donors (UD-SCT), very few studies have been undertaken.17  Therefore, we investigated the impact of TNFSF2 and IL10 polymorphisms on TRM after UD-SCT by analyzing individual polymorphisms and statistically inferred haplotypes. We have identified polymorphic haplotypes in the TNFSF2 region and in IL10 that may facilitate the selection of patients and donors for UD-SCT by identifying individuals at higher risk of TRM.

Patients

One hundred eighty-two patients undergoing UD-SCT at the University of Minnesota Hospital between 1990 and 2000 were studied. Forty-six percent of patients were younger than 18 years of age and 32% received transplants for nonmalignant disease. Sixty-nine percent of those with hematologic malignancy had early disease, acute leukemia in first or second complete remission, or chronic myelogenous leukemia (CML) in chronic phase at the time of transplantation. All donors were older than 18 years of age and the median donor age was 37 years. Sixty-five percent of transplantations were performed after 1995 and the median follow-up on survivors was 4.8 years. Eighty-nine percent of patients received total body irradiation-based pretransplantation protocols and 66% received T-cell-replete marrow transplants and received cyclosporine and methotrexate after transplantation. In 34% of cases T cells were removed by elutriation. Pretransplantation serologic typing indicated all donor-recipient pairs were matched at human leukocyte antigen A (HLA-A) and B loci. HLA-C locus typing was not performed. Sequence-specific oligonucleotide probes were used to genotype for HLA-DRB1, which indicated that 7% of pairs were HLA-DRB1 mismatched. Approval was obtained from the University of Minnesota institutional review board, and informed consent was provided by all individuals according to the Declaration of Helsinki.

Retrospective cytokine and HLA typing

Genotyping was performed on stored DNA or crude cell lysate. The quantity of DNA did not allow all pairs to be analyzed at all loci. Cytokine SNPs and HLA class I alleles were genotyped using induced heteroduplex generator (IHG) technology.18  Analysis of the TNFSF2 region and IL10 microsatellites was performed using standard polymerase chain reaction (PCR) methodology.

Statistical analysis

For IL10 haplotype analysis, 2 previously described methods were used to produce a physically linked 4-locus haplotype (10G, -1082, -819, -592).12,19  A PHASE-based statistical approach was used to generate haplotypes including the IL-10R microsatellite.20  For the TNFSF2 region, the PHASE-based approach was used alone because of large distances between individual polymorphic loci. Cumulative incidence rates and their 95% confidence intervals (CIs) were used to estimate treatment-related mortality, treating relapse-related deaths as a competing risk.21  Univariate comparisons of the major end points were completed by using the log-rank statistic. In preparation for inclusion into the multiple regression models, the proportional hazards assumption was tested by the Cox regression model using the Wald chi-square test and a time-dependent covariate of the logarithm of time. Cox proportional hazard models were then used to evaluate the independent effect of cytokine polymorphism.22  In addition to all significant TNFSF2 region and IL10 genotype factors, other potential confounders were included in the Cox regression models. These factors included recipient age, donor age, sex, sex mismatch, diagnosis, disease status, T-cell depletion, pretransplantation conditioning, and HLA disparity.

Relevant allele frequencies for both recipients and donors are summarized in Table 1. No significant differences were found between allele frequencies for TNFSF2 region and IL10 in the recipient and donor cohorts. In contrast to ID-SCT studies there was a considerable level of mismatch of TNFSF2 region genotypes between recipients and donors. At the TNFSF2 -238 locus, 8.8% of pairs were mismatched for one or more alleles and at the TNFSF2 -308 locus 13.3% of pairs were mismatched. As the TNFSF2 region is in linkage disequilibrium with flanking HLA genes, this level of mismatch is consistent with the levels of HLA allele disparity found using the IHG-based HLA typing technique.18 

Table 1.

Cytokine genotypes and haplotypes associated with increased or decreased risk of TRM in unrelated donor SCT at 12 months


Locus/haplotype and allele(s)

Allele frequency*

Status

Phenotype

No. died/total

Cl of TRM at 12 months, % (95% Cl, %)

P

Cl of relapse at 12 months, % (95% Cl, %)
Overall   NA   NA   NA   73/182   40 (33-47)   NA   12 (7-17)  
TNFSF2-1031C   0.2   R   C+   25/46   54 (38-70)   .04   5 (0-11)  
TNFSF2-1031C   0.2   R   C−   36/111   32 (23-41)    13 (7-19)  
TNF d4  0.22   R   d4+   35/62   56 (42-70)   < .01   3 (0-7)  
TNF d4  0.22   R   d4−   30/99   30 (21-39)    14 (7-21)  
TNF d4  0.22   D   d4+   34/65   52 (39-65)   .02   11 (4-18)  
TNF d4  0.22   D   d4−   33/99   33 (24-42)    11 (5-17)  
TNF a5  0.06   R   a5+   12/17   71 (45-97)   < .01   12 (0-26)  
TNF a5  0.06   R   a5−   57/154   37 (29-45)    12 (7-17)  
TNF a6  0.16   D   a6+   15/50   30 (17-43)   .05   20 (9-31)  
TNF a6  0.16   D   a6−   54/118   46 (37-55)    8 (3-13)  
TNF a7  0.05   D   a7+   12/18   67 (42-92)   < .01   0  
TNF a7  0.05   D   a7−   57/150   38 (30-46)    13 (8-18)  
TNF d4/TNFSF2-1031C   NA   R or D   d4+/C+   42/77   55 (43-67)   < .01   9 (3-15)  
TNF d4/TNFSF2-1031C   NA   R or D   Others   14/68   21 (12-30)    12 (8-16)  
IL 10 R2-G-C-C   NA   D   R2-G-C-C+   22/36   61 (43-79)   .01   11 (1-21)  
IL 10 R2-G-C-C   NA   D   Others   33/97   34 (25-43)    10 (4-16)  
IL 10 R3-G-C-C   NA   D   R3-G-C-C+   21/69   30 (19-41)   .01   13 (5-21)  
IL 10 R3-G-C-C
 
NA
 
D
 
Others
 
34/64
 
53 (40-66)
 

 
8 (2-14)
 

Locus/haplotype and allele(s)

Allele frequency*

Status

Phenotype

No. died/total

Cl of TRM at 12 months, % (95% Cl, %)

P

Cl of relapse at 12 months, % (95% Cl, %)
Overall   NA   NA   NA   73/182   40 (33-47)   NA   12 (7-17)  
TNFSF2-1031C   0.2   R   C+   25/46   54 (38-70)   .04   5 (0-11)  
TNFSF2-1031C   0.2   R   C−   36/111   32 (23-41)    13 (7-19)  
TNF d4  0.22   R   d4+   35/62   56 (42-70)   < .01   3 (0-7)  
TNF d4  0.22   R   d4−   30/99   30 (21-39)    14 (7-21)  
TNF d4  0.22   D   d4+   34/65   52 (39-65)   .02   11 (4-18)  
TNF d4  0.22   D   d4−   33/99   33 (24-42)    11 (5-17)  
TNF a5  0.06   R   a5+   12/17   71 (45-97)   < .01   12 (0-26)  
TNF a5  0.06   R   a5−   57/154   37 (29-45)    12 (7-17)  
TNF a6  0.16   D   a6+   15/50   30 (17-43)   .05   20 (9-31)  
TNF a6  0.16   D   a6−   54/118   46 (37-55)    8 (3-13)  
TNF a7  0.05   D   a7+   12/18   67 (42-92)   < .01   0  
TNF a7  0.05   D   a7−   57/150   38 (30-46)    13 (8-18)  
TNF d4/TNFSF2-1031C   NA   R or D   d4+/C+   42/77   55 (43-67)   < .01   9 (3-15)  
TNF d4/TNFSF2-1031C   NA   R or D   Others   14/68   21 (12-30)    12 (8-16)  
IL 10 R2-G-C-C   NA   D   R2-G-C-C+   22/36   61 (43-79)   .01   11 (1-21)  
IL 10 R2-G-C-C   NA   D   Others   33/97   34 (25-43)    10 (4-16)  
IL 10 R3-G-C-C   NA   D   R3-G-C-C+   21/69   30 (19-41)   .01   13 (5-21)  
IL 10 R3-G-C-C
 
NA
 
D
 
Others
 
34/64
 
53 (40-66)
 

 
8 (2-14)
 

Univariate comparisons were made to analyze associations between TNFSF2 region and IL 10 polymorphisms previously implicated as predictive of TRM. Only univariate statistically significant associations are shown in the Table. Factors that remained statistically significant in the Cox regression model included TNF d4/TNFSF2 -1031C and TNFa5.

NA indicates not applicable; R, recipient; and D, donor.

*

Allele frequencies were all in Hardy-Weinberg equilibrium. IL 10 haplotypes were determined by physical linkage.19 

The association between donor and recipient TNFSF2 region and IL10 genotypes and the cumulative incidence of TRM at one year was investigated (Table 1). The presence of the TNFd4 allele in conjunction with the TNFSF2 -1031C allele in recipients or donors significantly increased the cumulative incidence of TRM (P < .01; Figure 1A). Cox proportional hazards analysis showed that the TNFd4/TNFSF2 -1031 haplotype increased the risk of TRM at one year independent of all other factors. Relative risk (RR) = 2.2, CI 1.2-4.3, P = .02; and TNFa5 RR = 2.8, CI 1.2-6.2, P = .01. The other TNFSF2 region genotypes investigated were not associated with TRM.

Figure 1.

Cumulative incidence plots. (A) Cumulative incidence plot showing associations between transplant-related mortality (TRM) and presence or absence of the TNFd4+, 1031(C+) haplotype in donors and recipients of unrelated donor stem cell transplants. The cumulative incidence of TRM at one year after transplantation when the donor or recipient were TNFd4+, 1031(C+) haplotype-positive was 55% (43%-67%) compared with 21% (12%-30%) when both were negative (P < .01). (B) Cumulative incidence plot showing associations between transplant-related mortality and donor R2-G-C-C haplotype. The cumulative incidence of TRM at one year was significantly higher when the R2-G-C-C haplotype was present in the donor (61% [43%-79%]) compared with when it was absent from the donor (34% [25%-43%]), P = .01. (C) Cumulative incidence plot showing associations between transplant-related mortality and donor R3-G-C-C haplotype. The cumulative incidence of TRM at one year after transplantation was significantly lower when the R3-G-C-C haplotype was present in the donor (30% [19%-41%]) compared with when the haplotype was absent from the donor (53% [40%-66%]), P = .01.

Figure 1.

Cumulative incidence plots. (A) Cumulative incidence plot showing associations between transplant-related mortality (TRM) and presence or absence of the TNFd4+, 1031(C+) haplotype in donors and recipients of unrelated donor stem cell transplants. The cumulative incidence of TRM at one year after transplantation when the donor or recipient were TNFd4+, 1031(C+) haplotype-positive was 55% (43%-67%) compared with 21% (12%-30%) when both were negative (P < .01). (B) Cumulative incidence plot showing associations between transplant-related mortality and donor R2-G-C-C haplotype. The cumulative incidence of TRM at one year was significantly higher when the R2-G-C-C haplotype was present in the donor (61% [43%-79%]) compared with when it was absent from the donor (34% [25%-43%]), P = .01. (C) Cumulative incidence plot showing associations between transplant-related mortality and donor R3-G-C-C haplotype. The cumulative incidence of TRM at one year after transplantation was significantly lower when the R3-G-C-C haplotype was present in the donor (30% [19%-41%]) compared with when the haplotype was absent from the donor (53% [40%-66%]), P = .01.

Close modal

Individual IL10 polymorphic loci were not associated with TRM. However haplotype analysis revealed significant associations. The cumulative incidence of TRM at one year was higher in the presence of R2-G-C-C haplotype10  (P = .01; Figure 1B), whereas the R3-G-C-C haplotype was shown to be protective (P = .01; Figure 1C).

Cox proportional hazards analysis supported the association of the R3-G-C-C haplotype with TRM (RR = 0.4, CI 0.2-0.8, P = .03), although results could not be verified for the R2-G-C-C haplotype as hazards were not proportional. Cavet et al12  reported that ID-SC transplant recipients possessing longer 10G alleles were more likely to develop grades II to IV aGVHD. We have shown R2-G-C-C, a haplotype with predominantly long 10G alleles, is a risk factor for TRM after UD-SCT. The patients included in this study were a heterogeneous group. Furthermore, it must be emphasized that these analyses involved multiple comparisons and therefore P values between .05 and .01 should be interpreted with caution.

This is the first study to describe an association of cytokine polymorphism with TRM after UD-SCT, although no independent associations with aGVHD were identified. The most important finding is that recipient or donor TNFd4/TNFSF2 -1031C haplotype is a risk factor for TRM (Table 1). The same TNFd4 allele has recently been implicated as a risk factor for moderate to severe aGVHD in ID-SCT.16  It is likely that the TNFd4 allele is part of an extended chromosome 6 haplotype comprising immune HLA and a complement gene. We have shown significant associations between polymorphisms of TNFSF2 region and IL10 with TRM after UD-SCT. Confirmation of our results in a multicenter context would rule out any possibility of a center effect being responsible for our findings. Cytokine gene polymorphism analysis could become a useful pretransplantation tool to identify high-risk donors and recipients.

Prepublished online as Blood First Edition Paper, December 30, 2003; DOI 10.1182/blood-2002-11-3568.

Supported by grants from the National Health Service Research and Development Directorate, Children with Leukaemia, and the Leukaemia Research Fund.

The publication costs of this article were defrayed in part by page charge payment. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 U.S.C. section 1734.

1
Hill GR, Teshima T, Gerbitz A, et al. Differential roles of IL-1 and TNF-alpha on graft-versus-host disease and graft versus leukemia.
J Clin Invest.
1999
;
104
:
459
-467.
2
Holler E, Kolb HJ, Mittermuller J, et al. Modulation of acute graft-versus-host-disease after allogeneic bone marrow transplantation by tumor necrosis factor alpha (TNF alpha) release in the course of pretransplant conditioning: role of conditioning regimens and prophylactic application of a monoclonal antibody neutralizing human TNF alpha (MAK 195F).
Blood.
1995
;
86
:
890
-899.
3
Holler E, Kolb HJ, Moller A, et al. Increased serum levels of tumor necrosis factor alpha precede major complications of bone marrow transplantation.
Blood.
1990
;
75
:
1011
-1016.
4
Holler E, Ertl B, Hintermeier-Knabe R, et al. Inflammatory reactions induced by pretransplant conditioning—an alternative target for modulation of acute GvHD and complications following allogeneic bone marrow transplantation?
Leuk Lymphoma.
1997
;
25
:
217
-224.
5
Akdis CA, Blesken T, Akdis M, Wurthrich B, Blaser K. Role of IL-10 in specific immunotherapy.
J Clin Invest.
1998
;
102
:
98
-106.
6
Middleton PG, Taylor PR, Jackson G, Proctor SJ, Dickinson AM. Cytokine gene polymorphisms associating with severe acute graft-versus-host disease in HLA-identical sibling transplants.
Blood.
1998
;
92
:
3943
-3948.
7
Nagler A, Bishara A, Brautbar C, Barak V. Dysregulation of inflammatory cytokines in unrelated bone marrow transplantation.
Cytokines Cell Mol Ther.
1998
;
4
:
161
-167.
8
Allen RD. Polymorphism of the human TNF-alpha promoter-random variation or functional diversity.
Mol Immunol.
1999
;
36
:
1017
-1027.
9
Udalova IA, Nedospasov SA, Webb GC, Chaplin DD, Turetskaya RL. Highly informative typing of the human TNF locus using six adjacent polymorphic markers.
Genomics.
1993
;
16
:
180
-186.
10
Eskdale J, Keijsers V, Huizinga T, Gallagher G. Microsatellite alleles and single nucleotide polymorphisms (SNP) combine to form four major haplotype families at the human interleukin-10 (IL-10) locus.
Genes Immun.
1999
;
1
:
151
-155.
11
Eskdale J, Kube D, Gallagher G. A second polymorphic dinucleotide repeat in the 5′ flanking region of the human IL10 gene.
Immunogenetics.
1996
;
45
:
82
-83.
12
Cavet J, Middleton PG, Segall M, Noreen H, Davies SM, Dickinson AM. Recipient tumor necrosis factor-alpha and interleukin-10 gene polymorphisms associate with early mortality and acute graft-versus-host disease severity in HLA-matched sibling bone marrow transplants.
Blood.
1999
;
94
:
3941
-3946.
13
Gibson AW, Edberg JC, Wu J, Westerndorp RG, Huizinga TW, Kimberly RP. Novel single nucleotide polymorphisms in the distal IL-10 promoter affect IL-10 production and enhance the risk of systemic lupus erythematosus.
J Immunol.
2001
;
166
:
3915
-3922.
14
Eskdale J, Gallagher G. A polymorphic dinucleotide repeat in the human IL-10 promoter.
Immunogenetics.
1995
;
42
:
444
-445.
15
Socie G, Loiseau P, Tamousa R, et al. Both genetic and clinical factors predict the development of graft versus host disease after allogeneic hematopoietic stem cell transplantation.
Transplantation.
2001
;
72
:
699
-706.
16
Nordlander A, Uzunel M, Mattsson J, Remberger M. The TNFd4 allele is correlated to moderate-to-severe acute graft-versus-host disease after allogeneic stem cell transplantation.
Br J Haematol.
2002
;
119
:
1133
-1136.
17
Ishikawa Y, Kashiwase K, Akaza T, et al. Polymorphisms in TNFA and TNFR2 affect outcome of unrelated bone marrow transplantation.
Bone Marrow Transplant.
2002
;
29
:
569
-575.
18
Morse HR, Olomolaiye OO, Wood NA, Keen LJ, Bidwell JL. Induced heteroduplex genotyping of TNF-alpha, IL-1 beta, IL-6 and IL-10 polymorphisms associated with transcriptional regulation.
Cytokine.
1999
;
11
:
789
-795.
19
Wood NA, Keen LJ, Tilley LA, Bidwell JL. Determination of cytokine regulating haplotypes by heteroduplex analysis of DNA.
J Immunol Methods.
2001
;
2490
:
191
-198.
20
Stephens M, Smith NJ, Donnelly P. A new statistical method for haplotype reconstruction from population data.
Am J Hum Genet.
2001
;
68
:
978
-989.
21
Lin D. Non-parametric inference for cumulative incidence functions in competing risks studies.
Stat Med.
1997
;
16
:
457
-481.
22
Cox D. Regression models and life tables.
J R Stat Soc A.
1972
;
34
:
187
-220.
Sign in via your Institution