To the editor:

We read the manuscript by Parmar with interest. In an exploratory analysis, these investigators found that mismatching of major histocompatibility complex class I–related chain A (MICA) between donor and recipient was associated with gastrointestinal (GI) acute graft-versus-host disease (aGVHD).1  Based on the polymorphism of this protein, the constitutive expression on the GI epithelium and the T- and natural killer (NK)–cell immune activating function,2,3  we too have hypothesized a role for MICA in GI aGVHD. To address this we developed a high-resolution MICA typing method (see supplemental Table 1 and supplemental Figure 1, available on the Blood website; see the Supplemental Materials link at the top of the online article). MICA alleles were determined in 38 donor-recipient pairs after human leukocyte antigen (HLA) 12/12 matched unrelated donor transplantations. One-half of these recipients developed severe (grade II-IV) GI aGVHD, whereas the other half did not. All but one donor-recipient pair were matched at the MICA locus. Thus in this setting, MICA mismatching was rare. This is likely due to tight linkage disequilibrium with HLA B (D′0.99484, P < .001). Given that the single MICA mismatch occurred in a patient without aGVHD, this also suggests that MICA mismatching per se did not uniformly lead to severe GI aGVHD.

In reconsidering our hypothesis, we reasoned that the important interaction might be between the immune receptor NKG2D (on donor T and NK cells) and recipient MICA (on the GI tract). MICA amino acid substitutions at position 129 (methionine or valine) are associated with “weak” or “strong” NKG2D binding.4  We thus hypothesized that strong NKG2D:MICA binding might be associated with GI aGVHD. However no association between recipient MICA allele binding affinity and GI aGVHD was observed (not shown). An additional exploratory analysis examined whether certain recipient MICA alleles were associated with GI aGVHD protection or risk. There was a trend toward less aGVHD for recipients with MICA*008 (P = .07). Similarly, a possible increase in severe GI aGHVD was noted in recipients with an amino acid motif encoded by MICA*004, *006, *009, *044, or *049 (P = .075).

Given that MICA is tightly linked to HLA-B, we examined whether certain recipient HLA-B alleles (based on MICA linkage) are associated with GI aGVHD. Using a second cohort of 1676 recipients of myeloablative, HLA 8/8 matched unrelated donor transplantation, we tested whether recipient HLA-B alleles correlated with transplantation-associated outcomes. HLA-B alleles were divided into low and high risk based on linkage to MICA*008 or *004, *006, *009, *044, and *049, respectively (described in Table 1 footnotes). All other HLA-B alleles were considered to have intermediate aGVHD risk. As shown in Table 1, in multivariate analysis there were no differences in transplantation outcomes between the high-, intermediate-, and low-risk groups.

Table 1

Impact of HLA-B alleles on transplantation-related outcomes based on linkage to high- and low-risk HLA alleles

Clinical endpointOdds ratio (95% CI)P
GI aGVHD (grade II-IV)*  .059 
    2 low-risk alleles  
    Intermediate 1.04 (0.75-1.46)  
    ≥ 1 high-risk allele 0.63 (0.38-1.04)  
GI aGVHD (grade III-IV)*  .88 
    2 low-risk alleles  
    Intermediate 0.90 (0.58-1.40)  
    ≥ 1 high-risk allele 0.80 (0.44-1.44)  
Overall aGVHD (grade II-IV)*  .41 
    2 low-risk alleles  
    Intermediate 0.99 (0.79-1.24)  
    ≥ 1 high-risk allele 1.18 (0.89-1.57)  
Overall aGVHD (grade III-IV)*  .60 
    2 low-risk alleles  
    Intermediate 0.90 (0.68-1.19)  
    ≥ 1 high-risk allele 1.05 (0.73-1.51)  
Chronic aGVHD*  .15 
    2 low-risk alleles  
    Intermediate 1.15 (0.89-1.48)  
    ≥ 1 high-risk allele 1.38 (1.0-1.90)  
Relapse*  .17 
    2 low-risk alleles  
    Intermediate 1.35 (0.94-1.93)  
    ≥ 1 high-risk allele 1.01 (0.63-1.60  
Treatment-related mortality*  .11 
    2 low-risk alleles  
    Intermediate 0.76 (0.59-0.98)  
    ≥ 1 high-risk allele 0.80 (0.58-1.12)  
Disease-free survival  .64 
    2 low-risk alleles  
    Intermediate 0.96 (0.84-1.11)  
    ≥ 1 high-risk allele 0.91 (0.76-1.10)  
Overall survival  .45 
    2 low-risk alleles  
    Intermediate 0.95 (0.82-1.10)  
    ≥ 1 high-risk allele 0.89 (0.73-1.07)  
Clinical endpointOdds ratio (95% CI)P
GI aGVHD (grade II-IV)*  .059 
    2 low-risk alleles  
    Intermediate 1.04 (0.75-1.46)  
    ≥ 1 high-risk allele 0.63 (0.38-1.04)  
GI aGVHD (grade III-IV)*  .88 
    2 low-risk alleles  
    Intermediate 0.90 (0.58-1.40)  
    ≥ 1 high-risk allele 0.80 (0.44-1.44)  
Overall aGVHD (grade II-IV)*  .41 
    2 low-risk alleles  
    Intermediate 0.99 (0.79-1.24)  
    ≥ 1 high-risk allele 1.18 (0.89-1.57)  
Overall aGVHD (grade III-IV)*  .60 
    2 low-risk alleles  
    Intermediate 0.90 (0.68-1.19)  
    ≥ 1 high-risk allele 1.05 (0.73-1.51)  
Chronic aGVHD*  .15 
    2 low-risk alleles  
    Intermediate 1.15 (0.89-1.48)  
    ≥ 1 high-risk allele 1.38 (1.0-1.90)  
Relapse*  .17 
    2 low-risk alleles  
    Intermediate 1.35 (0.94-1.93)  
    ≥ 1 high-risk allele 1.01 (0.63-1.60  
Treatment-related mortality*  .11 
    2 low-risk alleles  
    Intermediate 0.76 (0.59-0.98)  
    ≥ 1 high-risk allele 0.80 (0.58-1.12)  
Disease-free survival  .64 
    2 low-risk alleles  
    Intermediate 0.96 (0.84-1.11)  
    ≥ 1 high-risk allele 0.91 (0.76-1.10)  
Overall survival  .45 
    2 low-risk alleles  
    Intermediate 0.95 (0.82-1.10)  
    ≥ 1 high-risk allele 0.89 (0.73-1.07)  

Allogeneic transplant recipients (n = 1676) were assigned to high-, intermediate-, or low-risk groups based on known linkage between HLA-B and MICA.5-8  MICA alleles (and the associated HLA-B alleles) that trended toward GI aGVHD protection (*008) were considered low-risk. MICA alleles associated with increased GI aGVHD risk (*004, *006, *009, *044, and *049) were considered high risk. All others were coded as intermediate risk. To investigate risk factors for MICA-associated HLA-B alleles and aGVHD in the cohort analysis, cumulative incidence rates of aGVHD (the chance a patient will have experienced aGVHD event before time t, and where death without the aGVHD event is the competing risk) were calculated using methods previously described.9  Multivariate analyses were applied to adjust for the effects of other covariates on the cumulative incidence of aGVHD using the pseudovalue approach of Klein10,11  with a logistic link function. A forward stepwise regression model using a generalized linear model for the pseudovalues was used. Shown are the P values comparing the clinical endpoint for recipients with 2 low-risk alleles to either intermediate or ≥ 1 high-risk allele.

CI indicates confidence interval.

*

Pseudovalue technique.

Cox regression.

In contrast with the data presented by Parmar et al, our studies do not support the concept that donor-recipient MICA mismatching plays a role in GI aGVHD. In addition, these data also do not suggest that certain recipient MICA alleles are associated with GI aGVHD. While the reasons for the differences between our study and that of Parmar et al are not entirely clear, potential explanations may include registry versus single-center reporting, high- versus low-resolution MICA typing, HLA disparity in the patient cohorts (matched vs potentially mismatched), and differing experimental approaches. As called for by Parmar et al, confirmation of the role of MICA in GI aGVHD in a larger cohort of patients is necessary, particularly using high-resolution typing methods.

The online version of this correspondence contains a data supplement.

Conflict-of-interest disclosure: The authors declare no competing financial interests.

Correspondence: Michael R. Verneris, Department of Pediatrics, University of Minnesota, Minneapolis, MN; e-mail: Verneris@umn.edu.

1
Parmar
 
S
de Lima
 
M
Zou
 
Y
, et al. 
Donor-recipient mismatches in MHC class I chain-related gene A in unrelated donor transplantation lead to increased incidence of acute graft-versus-host disease.
Blood
2009
, vol. 
114
 
14
(pg. 
2884
-
2887
)
2
Groh
 
V
Bahram
 
S
Bauer
 
S
Herman
 
A
Beauchamp
 
M
Spies
 
T
Cell stress-regulated human major histocompatibility complex class I gene expressed in gastrointestinal epithelium.
Proc Natl Acad Sci U S A
1996
, vol. 
93
 
22
(pg. 
12445
-
12450
)
3
Zwirner
 
NW
Dole
 
K
Stastny
 
P
Differential surface expression of MICA by endothelial cells, fibroblasts, keratinocytes, and monocytes.
Hum Immunol
1999
, vol. 
60
 
4
(pg. 
323
-
330
)
4
Steinle
 
A
Li
 
P
Morris
 
DL
, et al. 
Interactions of human NKG2D with its ligands MICA, MICB, and homologs of the mouse RAE-1 protein family.
Immunogenetics
2001
, vol. 
53
 
4
(pg. 
279
-
287
)
5
Zhang
 
Y
Lazaro
 
AM
Lavingia
 
B
Stastny
 
P
Typing for all known MICA alleles by group-specific PCR and SSOP.
Hum Immunol
2001
, vol. 
62
 
6
(pg. 
620
-
631
)
6
Gao
 
X
Single
 
RM
Karacki
 
P
Marti
 
D
O'Brien
 
SJ
Carrington
 
M
Diversity of MICA and linkage disequilibrium with HLA-B in two North American populations.
Hum Immunol
2006
, vol. 
67
 
3
(pg. 
152
-
158
)
7
Dorak
 
MT
Shao
 
W
Machulla
 
HK
, et al. 
Conserved extended haplotypes of the major histocompatibility complex: further characterization.
Genes Immun
2006
, vol. 
7
 
6
(pg. 
450
-
467
)
8
Ahmad
 
T
Marshall
 
SE
Mulcahy-Hawes
 
K
, et al. 
High resolution MIC genotyping: design and application to the investigation of inflammatory bowel disease susceptibility.
Tissue Antigens
2002
, vol. 
60
 
2
(pg. 
164
-
179
)
9
Klein
 
JP
Moeschberger
 
ML
Survival Analysis: Techniques for Censored and Truncated Data
2003
2nd ed.
New York
Springer-Verlag
10
Klein
 
JP
Anderson
 
PK
Regression modeling of competing risks data based on pseudovalues of the cumulative incidence function.
Biometrics
2005
, vol. 
61
 
1
(pg. 
223
-
229
)
11
Klein
 
JP
Modeling competing risks in cancer studies.
Statist Med
2006
, vol. 
25
 
6
(pg. 
1015
-
1034
)

Supplemental data

Sign in via your Institution