Rheumatology Advance Access originally published online on March 14, 2006
Rheumatology 2006 45(9):1096-1100; doi:10.1093/rheumatology/kel058
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
No primary association of MICA polymorphism with systemic lupus erythematosus
1Instituto de Biomedicina, CSIC, 2Servicio Inmunología, Hospital Virgen de las Nieves, 3Servicio de Medicina Interna, Hospital San Cecilio, 4Servicio de Medicina Interna, Hospital Virgen de las Nieves, Granada, 5Servicio de Inmunología, Hospital Virgen del Rocío, 6Servicio de Medicina Interna, Hospital Virgen del Rocío, Sevilla, 7Servicio de Reumatología, Hospital Xeral-Calde, Lugo and 8Servicio de Medicina Interna, Hospital Carlos-Haya, Malaga, Spain.
Correspondence to: M. F. González-Escribano, Servicio de Inmunología, Hospital Universitario Virgen del Rocío, Avda Manuel Siurot s/n, 41013 Sevilla, Spain. E-mail: mariaf.gonzalez.sspa{at}juntadeandalucia.es
| Abstract |
|---|
|
|
|---|
Objective. To replicate the described association between MHC class I chain-related A (MICA) gene polymorphism and susceptibility to systemic lupus erythematosus (SLE).
Methods. MICA transmembrane microsatellite polymorphism was genotyped using a polymerase chain reaction (PCR)-based method. Genotyping of HLA-B* and DRB1* was performed using PCR and detection with a reverse sequence-specific oligonucleotide (SSO) probe system. Combined data for these three loci (HLA-B*, DRB1* and MICA) were obtained from a total of 333 patients and 361 healthy controls.
Results. Significant association with B*08 [P<107, odds ratio (OR) 3.17, 95% confidence interval (CI) 2.025.00], DRB1*0301 (P<107, OR 2.07, 95% CI 1.592.68) and MICA5.1 (P = 0.01, OR 1.23, 95% CI 1.041.46) was observed. The combinations DRB1*0301-MICA5.1-B8 and HLA-DRB1*0301-B*08-positive and MICA5-1-negative were more frequent among SLE patients (11.4 vs 3.3% in healthy controls, P = 3.9 x 105, OR 3.76, 95% CI 1.857.73, and 6.9 vs 1.7%, P = 0.0007, OR 4.32, 95% CI 1.6813.10, respectively). Additionally, individuals who were HLA-DRB1*0301-B*08-negative and MICA5-1-positive were less frequent among patients (22.2 vs 31.3% in healthy controls, P = 0.007, OR 0.63, 95% CI 0.440.89) and the magnitude of the OR was similar to that obtained in individuals negative for all the three factors (OR 0.69, 95% CI 0500.94). Further analysis performed to detect independent association strongly suggested that the association between MICA5.1 and SLE is secondary to the linkage disequilibrium of this allele with B*08.
Conclusions. Our results do not support an independent association of MICA gene polymorphism with susceptibility to SLE.
KEY WORDS: MICA, HLA-DRB1, HLA-B, SLE, Susceptibility.
| Introduction |
|---|
|
|
|---|
Systemic lupus erythematosus (SLE) is a systemic autoimmune disease with a complex pathogenesis involving multiple genetic and environmental factors. Production of multiple autoantibodies and immune complex deposits causing tissue damage is a characteristic of this disease. Human major histocompatibility (HLA) class I and class II molecules have been described as genetic risk factors for many autoimmune diseases [reviewed in 13]. Two HLA class II haplotypes, DRB1*0301-DQA1*0501-DQB1*0201 and DRB1*1501-DQA1*0102-DQB1*0602, have been consistently associated with susceptibility to SLE in different Caucasian populations [46]. Additionally, several studies have described an association of HLA-B8 with susceptibility to SLE [7], and other genes located in the HLA region, such as TNFA, C2 and C4 have also been associated with susceptibility to the disease [reviewed in 8, 9].
A family of non-classical HLA genes known as major histocompatibility complex class I chain-related gene A (MICA) has been located 46 kb centromeric to the HLA-B gene. The MICA gene has a full length of 11.7 kb and encodes a stress-inducible molecule with three extracellular domains (
1,
2 and
3), a transmembrane (TM) region and a cytoplasmic tail (CY) [10, 11], and it shows preferential expression in fibroblasts and epithelial cells. The MICA chain folds similarly to typical class I chains and may have the capacity to bind peptides or other short ligands. A microsatellite polymorphism consisting of a variable number of GCT repeats encoding 4, 5, 6, 7, 9 or 10 alanine residues and located on the MICA exon 5 (TM segment) has been identified. Additionally, the MICA A5.1 allele presents a nucleotide insertion (G) between the second and third triplet repeats, resulting in a premature stop codon [1214].
The location of the MICA gene in the HLA region and its association with other autoimmune diseases suggests the MICA gene as a possible candidate gene for the development of SLE. A recent study reported an independent contribution of the MICA TM polymorphism to susceptibility to SLE in an Italian population [15]. The aim of the present study was to investigate whether the MICA TM polymorphism, as well as HLA class I (locus B) and class II (locus DRB1*), influences genetic predisposition to SLE in a large Spanish cohort of patients.
| Materials and methods |
|---|
|
|
|---|
Subjects
A total of 574 SLE patients (13% men and 87% women), meeting the American College of Rheumatology (ACR) criteria for SLE [16, 17], were recruited from five Spanish hospitals: Hospital Virgen de las Nieves and Hospital Clínico San Cecilio (Granada), Hospital Virgen del Rocio (Seville), Hospital Xeral-Calde (Lugo) and Hospital Carlos-Haya (Malaga). The mean (S.D.) age of the SLE patients at diagnosis was 43±13.3 yr and the mean age at onset was 32±15 yr. A total of 934 blood bank and bone marrow donors from the corresponding cities were included as healthy controls. The patient and control groups were both of Spanish Caucasian origin and were matched for age and sex. Samples were obtained from subjects after they had given written informed consent. The study was approved by all local ethical committees of the corresponding hospitals.
MICA genotyping
Genomic DNA was isolated from anticoagulant-treated peripheral blood mononuclear cells using standard methods. Genotyping of the MICA TM microsatellite polymorphism was performed using a polymerase chain reaction (PCR)-based method as previously described [18]. Briefly, primer sequences were MICA 5F 5'-CCT TTT TTT CAG GGA AAG TGC-3' and MICA 5R 5'-CCT TAC CAT CTC CAG AAA CTG C-3'. The reverse primer was labelled at the 5' end with 6-FAM dye. PCR products were electrophoresed in an ABI Prism 3100 Genetic Analyzer (Applied Biosystems, Foster City, CA, USA) and their sizes were determined using the Genescan 672 software (Applied Biosystems).
HLA genotyping
Genotyping for HLA-B* and DRB1* was carried out using a reverse dot-blot kit with sequence-specific oligonucleotide (SSO) probes (Dynal RELITM SSO HLA-B and DRB1* typing kits; Dynal Biotech, Bromborough, UK). When necessary, high resolution typing of HLA-DRB1*03 samples was performed using Dynal AllSetTM SSP DRB1*03.
Genotyping of the three loci HLA-B*, MICA and DRB1* was performed in 333 SLE patients and 361 healthy controls.
Statistical analyses
Allelic and phenotypic frequencies of the HLA markers studied were obtained by direct counting. Statistical analysis to compare distributions was performed with the
2 test. Odds ratios (OR) and 95% confidence intervals (CIs) were calculated according to Wolf 's method. The software used was Statcalc program (Epi Info 2002; Centers for Disease Control and Prevention, Atlanta, GA, USA). P-values below 0.05 were considered statistically significant. The P-values were corrected (Pc) according to the Bonferroni method for the MICA TM alleles (n = 5). To detect the strongest association for data in Table 2, the test described by Svejgaard and Ryder was used and the P-values were corrected as recommended by these authors [19]. The correction factor used in each case is shown in parentheses.
|
| Results |
|---|
|
|
|---|
No significant differences in the frequency of the different alleles of the three genes studied were observed among patient groups or among control groups from the five cohorts. Hence we combined the five cohorts to form a SLE casecontrol group, which was used in further analyses. The control study population was found to be in HardyWeinberg equilibrium for all the loci studied.
Table 1 shows the distribution of the MICA alleles in SLE patients and controls. As has been found in other Caucasian populations [15, 18], in our population we found five alleles: A4, A5, A5-1, A6 and A9; the MICA allele displaying the highest frequency in both patient and control groups was A6 (35.8 and 38.9% respectively). There was a significant deviation in the distribution of the alleles of MICA that distinguished between patient and control groups (P = 0.007,
2 test with a 5 x 2 contingency table). This significant deviation corresponded to a higher frequency of MICA5.1 [28.1 vs 24.0% in the control group; P = 0.01, P(5) = 0.05, OR 1.23, 95% CI 1.041.46] and a lower frequency of MICA5 among SLE patients [10.9 vs 14.1% in the control group; P = 0.009, P(5) = 0.045, OR 0.74, 95% CI 0.590.94]. No significant differences in the distribution of individuals homozygous for MICA5.1 were observed when patient and control groups were compared (12.0 vs 12.3%, P>0.05).
|
Regarding HLA typing, the HLA-B* gene was genotyped in 399 SLE patients and 385 control subjects. The allelic frequency of HLA-B*08 was significantly higher in patients (11.0 vs 3.8% in the control group, P<107, OR 3.17, 95% CI 2.025.00). Additionally, the HLA-DRB1* was genotyped in 459 SLE patients and 522 healthy controls and, as in other studies, the allelic frequency of HLA-DRB1*0301 was higher in patients (20.4 vs 11.0% in controls, P<107, OR 2.07, 95% CI 1.592.68). The allelic frequency of HLA-DRB1*15 was higher in SLE patients (12.0 vs 9.4% in the control group) but the difference did not reach statistical significance (P = 0.06). The distribution of HLA-DRB1*08 was similar in patients and controls (3.4 vs 3.1%, P = 0.7).
Strong linkage disequilibrium between B*08 and both MICA5.1 and DRB1*0301 was observed in our control population [P(9) = 0.00018 and P(9)<106, respectively]. Nevertheless, no significant linkage disequilibrium between MICA5.1 and DRB1*0301 was observed in our control group (P = 0.4). No significant linkage disequilibrium between MICA5.1 and DRB1*15 was observed among B*08- or DRB1*0301-negative controls [P(9)>0.05 in all cases].
We performed further analysis to better define the contribution of each marker to susceptibility to SLE. Table 2 shows data for the four phenotypic combinations of the three factors studied two by two. The highest OR, taking into account two positive risk factors, was for the combination B*08-positiveDRB1*0301-positive (P<107, OR 4.18, 95% CI 2.407.36) and the lowest OR was for the combination B*08-negativeDRB1*0301-negative (P<107, OR 0.43, 95% CI 0.310.59). Additionally, frequencies of individuals who were B*08 positiveMICA5.1-negative and individuals who were DRB1*0301-positiveMICA5.1-negative were higher in patients (7.8 vs 1.7% in controls, P = 7.7 x 105, OR 5.07 95% CI 2.0515.01, and 19.5% vs 13.0% in controls, P = 0.02, OR 1.62, 95% CI 1.062.49). Nevertheless, frequencies of individuals who were B*08-negative-MICA5.1-positive and individuals who were DRB1*0301-negative-MICA5.1-positive were lower among patients (29.3 vs 34.1% in the control group, P>0.05, and 24.9 vs 33.0 in the control group, respectively, P = 0.02, OR 0.68 95% CI 0.480.95).
Table 3 shows data for the eight possible phenotypic combinations of the three risk factors studied. The highest OR (4.32, 95% CI 1.6813.10) found was for the combination B*08 positiveMICA5.1-negativeDRB1 0301-positive, whereas the lowest OR found was for the combination B*08-negativeMICA5.1-positiveDRB1*0301-negative (0.63, 0.440.89). Next, we investigated the association of DRB1*0301-MICA5.1-B*08 to evaluate which of these factors are likely to be the responsible of the association, using the test recommended by Sverjgaard and Ryde [19]. In this analysis, MICA5.1 did not show any association with SLE in B*08-positive (P = 0.2) or in B*08-negative (P = 0.8) individuals or in DR*0301-positive (P = 0.2) or DRB1*0301-negative (P = 0.6) individuals, suggesting that MICA5.1 is not an independent risk factor for SLE. On the contrary, B*08 was associated with SLE in MICA5.1-positive [P(6) = 0.0012, OR 2.57, 95% CI 1.575.20] and negative [P(6) = 0.00012, OR 5.70, 95% CI 2.0717.03] individuals and DRB1*0301 was also associated in both MICA5.1-positive [P(6) = 0.0042, OR 2.92, 95% CI 1.655.20] and MICA5.1-negative [P(6) = 0.05, OR 1.79, 1.132.84] individuals. Regarding the independent contribution of B*08 and DRB1*0301, the frequency of B*08 was not increased in DRB1*0301-positive individuals (P = 0.1) and was probably only increased in DRB1*0301-negative [P = 0.02, P(6) = 0.12, OR 1.56 95% CI 1.042.34] individuals. Conversely, the frequency of DRB1*0301 was increased in B*08-positive individuals [P(6) = 0.0018, OR 2.98, 95% CI 1.555.77] and probably only increased in B*08 negative individuals [P = 0.03, P(6) = 0.18, OR 2.33, 95% CI 1.005.54]. Using the test recommended by Sverjgaard and Ryde [19], associations with SLE of B*08 and DRB1*0301 were stronger than the association of MICA5.1 [P(6) = 0.00036 and P(6) = 0.024 respectively]. Nevertheless, no significant differences in the strength of the association of B*08 and DRB1*0301 was observed (P = 0.4). Finally, strong linkage disequilibrium between B*08 and both MICA5.1 and DRB1*0301 was observed also in the patient group [P(9) = 0.000072 and P(9)<106, respectively]. Nevertheless, no significant linkage disequilibrium between MICA5.1 and DRB1*0301 was observed among patients (P = 0.2) or among controls (P = 0.4), even after stratification for B*08 [P(9)>0.05 in all cases].
|
Additionally, there was no evidence for an association of MICA5.1 with SLE among B*08 positiveDRB1*0301-positive individuals ORs were very similar for MICA5.1-negative and MICA5.1-positive individuals (4.32 vs 3.76).
| Discussion |
|---|
|
|
|---|
The MHC region has been associated with susceptibility to SLE in many studies. However, the literature concerning association with HLA-SLE is often in disagreement regarding the genes involved in the susceptibility [20, 21]. This could be due to confounding factors caused by strong linkage disequilibrium in the region.
Recently, a study suggesting an independent contribution of the MICA TM polymorphism to susceptibility to SLE in an Italian population was reported [15]. This casecontrol study was performed in a relatively small group of individuals (48 patients and 158 controls) and the authors reported a positive association of susceptibility to SLE with both MICA5.1 and MICA5 alleles and a negative association with MICA9. Gambelunghe et al. [15] reported an association between HLA-B8 and SLE in the Italian population. They also found strong linkage disequilibrium between HLA-B8 and MICA5.1, indicating that HLA-B8 was not significantly and independently associated with SLE. However, they did not perform any analysis to detect which of these genes showed the strongest association. We think this may explain the discrepant results obtained in the present work. Our study failed to confirm an independent association between MICA and SLE. The MICA5.1 allele showed an association with SLE susceptibility in the one-marker analysis, although the OR (1.23) found was lower than those obtained for B*08 (3.17) and DRB1*0301 (2.07). Tables 2 and 3 show that the presence of MICA5.1 is neutral if B*08 is absent, or even protective if DRB1*0301 or both B*08 and DRB1*0301 are absent. The magnitude of the OR (0.63), as given by the simultaneous presence of MICA5.1 and absence of both B*08 and DRB1*0301, was similar to that given by the absence of these three factors (0.69). Additionally, analysis performed to detect independent association strongly suggests that the association described between MICA5.1 and SLE is secondary to the linkage disequilibrium of this allele with B*08. Regarding MICA5, we observed a lower the frequency of this allele among patients. This contradiction with the previous study could be due to a difference in linkage disequilibrium between HLA-B* and MICA in Italian and Spanish populations and, in any case, the contrary effect found for MICA5 in both populations does not support any association of this allele. No differences in the distribution of the allele MICA9 were observed to allow comparison between patient and control groups.
Our results confirm previous findings regarding the association of HLA-B*08 and DRB1*0301 with susceptibility to SLE [reviewed in 15]. Data obtained in the present study are not conclusive enough to clarify which of these two markers is more strongly associated with the disease. Tables 2 and 3 show a strong positive association of susceptibility to the disease with the presence of both factors and a strong negative association with their absence. Nevertheless, individuals positive only for one of these factors are in the minority, although they are not significantly increased among patients. Additional analysis performed to detect independent association did not provide strong evidence supporting any hypothesis because several critical P-values are not significant after correction, both markers show a similar strength of association and they are in strong linkage disequilibrium. Regarding other class II specificities associated with SLE in different populations, as already described in other Caucasian populations, only a weak increase in the frequency of HLA-DRB1*1501 was observed in our series of patients and the distribution of DRB1*08 was very similar in both patients and controls.
In conclusion, our results do not support an independent association between MICA and SLE and are in agreement with the previous finding reporting an association with B*08 and DRB1*0301.
| Acknowledgments |
|---|
This work was supported by Plan Nacional de I + D (SAF03-3460), Fondo de Investigaciones Sanitarias (PI 04 0067) and in part by Junta de Andalucía, grupos CTS-197 and CTS-180. B.T. is the recipient of a fellowship from Fundación Reina Mercedes.
| Notes |
|---|
* These authors contributed equally to this work
| References |
|---|
|
|
|---|
- Tsao BP. (2004) Update on human systemic lupus erythematosus genetics. Curr Opin Rheumatol 16:51321.[CrossRef][Web of Science][Medline]
- Yu CY and Whitacre CC. (2004) Sex, MHC and complement C4 in autoimmune diseases. Trends Immunol 25:6949.[CrossRef][Web of Science][Medline]
- Nath SK, Kilpatrick J, Harley JB. (2004) Genetics of human systemic lupus erythematosus: the emerging picture. Curr Opin Immunol 16:794800.[CrossRef][Web of Science][Medline]
- Graham RR, Ortmann WA, Langefeld CD, et al. (2002) Visualizing human leukocyte antigen class II risk haplotypes in human systemic lupus erythematosus. Am J Hum Genet 71:54353.[CrossRef][Web of Science][Medline]
- Mok CC and Lau CS. (2003) Pathogenesis of systemic lupus erythematosus. J Clin Pathol 56:48190.
[Abstract/Free Full Text] - Ramal LM, Lopez-Nevot MA, Sabio JM, et al. (2004) Systemic lupus erythematosus in southern Spain: a comparative clinical and genetic study between Caucasian and Gypsy patients. Lupus 13:93440.
[Abstract/Free Full Text] - Price P, Witt C, Allcock R, Sayer D, et al. (1999) The genetic basis for the association of the 8.1 ancestral haplotype (A1, B8, DR3) with multiple immunopathological diseases. Immunol Rev 167:25774.[CrossRef][Web of Science][Medline]
- Maddison PJ. (1999) Nature and nurture in systemic lupus erythematosus. Adv Exp Med Biol 455:713.[Web of Science][Medline]
- Reveille JD. (2005) Genetic studies in the rheumatic diseases: present status and implications for the future. J Rheumatol Suppl. 72, 103.
- Bahram S, Bresnahan M, Geraghty DE, Spies T. (1994) A second lineage of mammalian major histocompatibility complex class I genes. Proc Natl Acad Sci USA 91:625963.
[Abstract/Free Full Text] - Bahram S, Mizuki N, Inoko H, Spies T. (1996) Nucleotide sequence of the human MHC class I MICA gene. Immunogenetics 44:801.[CrossRef][Web of Science][Medline]
- Mizuki N, Ota M, Kimura M, et al. (1997) Triplet repeat polymorphism in the transmembrane region of the MICA gene: a strong association of six GCT repetitions with Behcet disease. Proc Natl Acad Sci USA 94:1298303.
[Abstract/Free Full Text] - Perez-Rodriguez M, Corell A, Arguello JR, et al. (2000) A new MICA allele with ten alanine residues in the exon 5 microsatellite. Tissue Antigens 55:1625.[CrossRef][Web of Science][Medline]
- Rueda B, Pascual M, Lopez-Nevot MA, Gonzalez E, Martin J. (2002) A new allele within the transmembrane region of the human MICA gene with seven GCT repeats. Tissue Antigens 60:5268.[CrossRef][Web of Science][Medline]
- Gambelunghe G, Gerli R, Bocci EB, et al. (2005) Contribution of MHC class I chain-related gene A (MICA) polymorphism to genetic susceptibility for systemic lupus erythematosus. Rheumatology 44:28792.
[Abstract/Free Full Text] - Tan EM, Cohen AS, Fries JF, et al. (1982) The 1982 revised criteria for the classification of systemic lupus erythematosus. Arthritis Rheum 25:12717.[Web of Science][Medline]
- Hochberg MC. (1997) Updating the American College of Rheumatology revised criteria for the classification of systemic lupus erythematosus. Arthritis Rheum 40:1725.[Web of Science][Medline]
- Ota M, Katsuyama Y, Mizuki N, et al. (1997) Trinucleotide repeat polymorphism within exon 5 of the MICA gene (MHC class I chain-related gene A): allele frequency data in the nine population groups Japanese, Northern Han, Hui, Uygur, Kazakhastan, Iranian, Saudi Arabian, Greek and Italian. Tissue Antigens 49:44854.[Web of Science][Medline]
- Svergaard A and Ryder LP. (1994) HLA and disease associations: detecting the strongest association. Tissue Antigens 43:1827.[Web of Science][Medline]
- Goldstein R and Sengar DP. (1993) Comparative study of the major histocompatibility complex in French Canadian and non-French Canadian Caucasian with systemic lupus erythematosus. Arthritis Rheum 36:11217.[Web of Science][Medline]
- Bekker-Mendez C and Yamamoto-Furusho JK, et al. (1998) Haplotype distribution of class II MHC genes in Mexican patients with systemic lupus erythematosus. Scand J Rheumatol 27:3736.[CrossRef][Web of Science][Medline]
This article has been cited by other articles:
![]() |
C. Pan, C. Wu, H. Chen, C. Dang, F. Chang, H. Liu, C. Chu, M Lin, and Y. Lee Molecular analysis of HLA-DRB1 allelic associations with systemic lupus erythematous and lupus nephritis in Taiwan Lupus, July 1, 2009; 18(8): 698 - 704. [Abstract] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
