Type 1 Diabetes in the Spanish Population: additional factors to Class II HLA-DR3 and -DR4
© Urcelay et al; licensee BioMed Central Ltd. 2005
Received: 15 November 2004
Accepted: 20 April 2005
Published: 20 April 2005
The Major Histocompatibility Complex is the main genetic contributor to susceptibility to type 1 diabetes (T1D); genome-wide scans have consistently mapped increased predisposition to this region. The highest disease risk has been associated with HLA-DR3 and HLA-DR4. In particular, the DR3-positive ancestral haplotype 18.2 was reported as highly diabetogenic. We aimed to corroborate whether this haplotype increases the susceptibility conferred by the DQ2-DR3 alleles in a Mediterranean population. We also searched for additional susceptibility factors to the classic DQ2-DR3 and DQ8-DR4.
Genetic MHC markers were analysed in a case-control study with 302 T1D patients and 529 ethnically matched controls. DR3-TNFa1b5 carrier rate was significantly higher in DR3-positive heterozygous T1D patients than in DR3-positive heterozygous controls (p = 0.0019; odds ratio OR [95% confidence interval CI] = 2.26 [1.3–3.93]). This data was confirmed analysing the allelic frequency, which includes the information corresponding to the DR3-homozygous individuals (p = 0.001; OR = 2.09) and by using the Arlequin software to check the DR3-positive haplotypes (p = 0.004;OR = 1.93). The present results provide strong evidence of a second susceptibility region in the ancestral haplotype 18.2 in the Spanish population.
Moreover, we searched for T1D susceptibility factors in addition to the MHC classical ones, within the DR2-DQ6/DR3-DQ2/DR4-DQ8 negative population. Several genetic markers in both MHC class II (DQA1*0101-DQB1*0501 [p = 0.007;OR = 2.81], DQA1*0201-DQB1*0202 [p = 0.03; OR = 2.35]) and III (TNFa2b1 [p = 0.01 OR = 2.74], BAT-2*2 [p = 0.004; OR = 3.19]) were found. These different alleles associated with T1D were not independent and we observed linkage disequilibrium among them leading us to describe two new risk haplotypes (DQA1*0101-DQB1*0501-TNFa2b1 and DQA1*0201-DQB1*0202- BAT-2*2). Finally, we studied a T1D susceptibility/protection marker located in extended class I, D6S2223; however, no association was observed in our population.
Our results suggest that other associated MHC haplotypes might present susceptibility factors in loci different from HLA-class II and that the class II molecules are not necessarily the universal etiologic factor in every MHC haplotype.
Type 1 diabetes is a multifactorial autoimmune disease characterised by insulin deficiency, due to the T cell mediated destruction of pancreatic β-cells . Among the genetic determinants of susceptibility, with more than 18 putative loci identified to date, a region in chromosome 6p21 (IDDM1) containing the Major Histocompatibility Complex (MHC) is the only one consistently associated with T1D in genome-wide screenings. The MHC spans approximately 4 Mb and consists of over 200 genes arranged in three subregions named class II, III and I. More than 90% of Caucasoid type 1 diabetic patients have at least one copy of the class II HLA-DR3 or DR4 allele, as compared to the 45% present in the general population and the genotype frequency of the DR3/DR4 heterozygote in T1D patients is 40% vs. 3% in controls . In fact, the strongest susceptibility haplotypes described are DQB1*0201-DQA1*0501-DRB1*03 (DQ2-DR3) and DQB1*0302-DQA1*0301-DRB1*04 (DQ8-DR4), especially when both appear together in the genotype. However, not all HLA-DR3 or -DR4 positive haplotypes are equally predisposing [3, 4]. This might suggest the role of other MHC loci responsible for modifying the susceptibility to diabetes conferred by class II genes, but their search has been difficult due to the strong linkage disequilibrium (LD) present in this chromosomal region. The term ancestral, extended haplotype was coined referring to continuous sequence derived with little, if any, change from an ancestor of all those now carrying all or part of the haplotype . A published report stated that the DQ2-DR3-B18 AH 18.2 significantly increased risk compared to other haplotypes with the same class II alleles, but LD between markers hampered a more precise localisation of the presumed susceptibility gene . Taking into account that this AH 18.2 is more frequent in Southern European populations, we decided to assess the possible role of other loci in the MHC region besides DQ-DR within this haplotype. We evaluated several genetic markers along the MHC in a case-control study with 302 Spanish T1D patients and 529 healthy controls. This study reproduced the especially strong association of the AH 18.2 with T1D in the Spanish population, supporting the existence of a second susceptibility locus within this haplotype. The present work extended the search to other T1D risk factors in the MHC besides DQ2-DR3/DQ8-DR4. Additional MHC haplotypes question the paradigm of class II genes as sole responsible for the association and open up the possibility that class III susceptibility factors would be also involved.
We first aimed to corroborate in the Spanish population the published increased risk to autoimmune diabetes of the AH 18.2, as compared to other DR3-positive haplotypes . Microsatellites TNFa and TNFb, located in the untranslated region upstream of the LTA gene, exhibit extensive polymorphism and are good haplospecific markers. In order to identify the extended haplotype AH 18.2, we selected both microsatellites TNFab, which are conveniently close to HLA-B and easier to genotype than the latter. The microsatellites TNFa1b5 tag the AH18.2 defined by the characteristic alleles of five markers (D6S273*2, BAT2*2, TNFa1b5, MICA*4) in our families with IgA deficiency  and celiac disease . These AH 18.2 markers were also confirmed by the computer program Arlequin, by linkage disequilibrium studies and by analysis of homozygous individuals. As we are dealing with patients and controls with no families at our disposal, the DR3-TNFa1b5 phase was uncertain. However, due to heavy linkage disequilibrium between these alleles, most of the times they will probably be in cis. This proportion can be estimated based upon the presence of TNFa1b5 in DR3-negative individuals: 8 carriers out of 109 DR3-negative T1D patients, and 13 carriers out of 376 DR3-negative healthy controls. These numbers yield allelic frequencies of TNFa1b5 on DR3-negative haplotypes of 7 and 3%, respectively. The assumption can therefore be made that, most of the time; a DR3 heterozygous sample carrying TNFa1b5 will in fact carry these alleles in the same haplotype.
Comparison of the TNFa1b5 allelic frequency in the DR3-positive T1D patient and control cohorts.
Frequency and number of DR3-positive haplotypes estimated by the Arlequin software.
T1D (2n = 596)
Controls (2n = 1006)
In our Spanish sample this DR3-TNFa1b5 positive haplotype was present in a 33% T1D patients (98/298) vs. 7% controls (37/504). The haplotypes ascertained by the Arlequin software yielded an AH 18.2 almost five times more frequent in the Spanish patient than control cohorts (Table 2). Moreover, almost one half (46%, 66/142) of the DR3 heterozygous diabetic patients vs. one forth (27%, 33/120) of the controls displayed a DR3-TNFa1b5 positive phenotype. On the contrary, in this subpopulation the distribution of the other main DR3-positive haplotype, AH 8.1, was similar between patients (43/142, 30%) and controls (47/120, 39%).
Therefore, these results confirmed the reported specific diabetogenic role of the AH 18.2 in a Mediterranean population and are compatible with published observations that established the presence of a second susceptibility gene in the AH 18.2 [6, 9] possibly located in either class III or class I.
Comparison of the carrier rate of HLA class III genetic markers in the DR2-DQ6/DR3-DQ2/DR4-DQ8 negative diabetic and control cohorts.
n = 29
n = 251
n = 29
n = 263
n = 29
n = 234
n = 29
n = 233
DQA1-DQB1 haplotypes in DR2-DQ6/DR3-DQ2/DR4-DQ8 negative T1D patients and controls.
T1D (n = 29)
Controls (n = 282)
Finally, we tried to ascertain the role as T1D susceptibility/protection marker of the microsatellite located 4.9 Mb telomeric of DQ in extended class I, D6S2223. In Northern European countries allele D6S2223*3 was associated with a reduction in risk independent of the class II effect [9, 10] and allele D6S2223*1 accounted for susceptibility to T1D . These associations were not observed in Sardinia . In our DR3, DR4 and DR2-negative study group we did not find any association (allele D6S2223*3 for protection: 28/29 in T1D patients vs. 223/234 in healthy controls, OR = 1.38, p = 0.61; allele D6S2223*1 for susceptibility: 0/29 vs. 2/234, p = 0.79). However, the limited number of patients after stratification decreases the statistical power of detection for a positive association in our population.
Our observations support the distinctive effect on type 1 diabetes susceptibility of the AH 18.2 among all the other DR3-positive haplotypes in the Spanish population. Both, allelic and phenotypic frequencies showed the increased risk conferred by the DR3-TNFa1b5 haplotype. Moreover, the different behaviour in terms of T1D predisposition of both DR3-positive main haplotypes, AH 8.1 and 18.2, argues in favour of a second susceptibility locus besides the standard class II molecules (DR and DQ) in the latter.
Previous studies reported the contribution to T1D susceptibility of a second region in the MHC besides DQ-DR genes. Among them, one concluded that this second region extended between HLA-B and BAT-3 in the highest risk DR3/DR4 genotype , while other mapped the critical region around the microsatellite D6S273 centromeric to TNF . TNF-α has been identified as a diabetogenic effector molecule synergistic with IFN-γ in the induction of β-cell apoptosis [14–16] and therefore, it was of interest to consider the TNF region located in class III MHC as a candidate susceptibility locus. Transfected cells containing BAT1 promoter fragments from the 8.1 AH exhibited lower reporter activity compared to the 7.1 AH. Since the 7.1 and 8.1 AH are associated with resistance and susceptibility to T1D respectively, the observed effects are claimed to have significant implications for the pathogenesis of this disease . We decided to study this region in our population. Our findings of BAT-2*2, TNFa2b1 and TNFa10b4 as genetic markers of the risk conferred by the DQ-DR locus agrees with those previous reports in terms of location. Most probably the etiological variant(s) will be in linkage disequilibrium with the BAT-2*2 allele. However, given the low recombination rate of the haplotype 18.2 (shown by the TNFa1b5 present in a DR3-negative context, see Table 3), and being BAT-2*2 a common allele, this microsatellite could mark a different susceptibility gene in the AH 18.2 and in other haplotypes. Clearly, the analysis of genetic markers in class II and class III provided an emerging pattern of susceptibility/protection haplotypes in addition to the classical DR3/DR4/DR2 determinants. In this sense, one could speculate that the traditional emphasis on class II alleles alone as responsible for increasing the T1D risk could be an oversimplification of a more complex reality where class III markers remained in the background.
Population studies provide invaluable information helping to map predisposition loci. A good example is the role of D6S2223 as a marker of autoimmune diabetes in several Nordic populations  which has been reproduced neither in Sardinia  nor in Spain (present data).
A large number of studies from different populations including the Spanish one  have indicated that common allelic variants at the class II HLA-DRB1, -DQA1 and -DQB1 loci are associated with T1D. These MHC class II molecules play an important role in the presentation of peptide antigens after intracellular processing to CD4 T- lymphocytes. The correlation between the relative T1D predisposition of class II alleles and the structure of their proteins has also been described . However, the presence of HLA class II molecules alone does not by itself cause disease. The HLA transgenic mice develop diabetes when there is an islet "insult", even if the islet "insult" is, itself, not sufficient to precipitate disease .
Our data prove that there are additional susceptibility/protection factors besides DR3/DR4/DR2 in HLA in the Spanish population and particular combinations of these minor risk factors could have important effects on susceptibility. The results of the present study provide evidence of the importance of the whole HLA in T1D risk gradient after stratification for the DR3, DR4 and DR2 determinants. The accent placed on class II in the canonical DQ2-DR3/DQ8-DR4 haplotypes is not necessarily extensive to other associated haplotypes where susceptibility might map to different MHC loci.
302 T1D patients and 529 healthy unrelated controls, both groups composed of Caucasian individuals from the same Madrid area, were consecutively recruited after informed consent obtained from the indexed subjects or their parents when the patient was a child. The T1D patients were diagnosed at two hospitals in Madrid and the control subjects were collected among healthy blood donors. The age at onset for the T1D patients (150 men and 152 women) range 1–55 years old (median onset 15 years old). Diagnosis was based on patients' clinical features and laboratory data according to the criteria of the American Diabetes Association (ADA). All subjects were insulin-dependent at the time of the study. Ab-positivity was studied for GAD, IA-2, insulin and ICA. The protocol followed the principles expressed in the Declaration of Helsinki and was approved by the Hospital Ethics Committee.
Both patients and controls were genotyped for HLA-DQB1, DQA1 and DRB1 alleles  for 6 microsatellite markers (D6S273, BAT-2, TNFa, TNFb, MICA and D6S2223) as previously described [20–25]. Microsatellite alleles were ascertained using an ABI Prism™ 310 automatic sequencer (Applied Biosystems, Foster City, CA, USA). Each sample included an internal size standard (TAMRA 500, Applied Biosystems) to achieve a highly consistent measure.
Allelic distribution and frequencies of the haplotypes present in the control and diabetic cohorts were determined using a software for population genetics data analysis (Arlequin ver 2.000. Schneider S, Roessli D and Excoffier L. University of Geneva, Switzerland). Haplotypes were also corroborated by our previous findings in families affected by other immune related diseases, by LD studies in controls, and by the recent literature .
The frequencies of each marker allele in patients versus controls were compared by a standard chi-square test using a statistical software package (EPI-INFO v. 6.02; Center for Disease Control & Prevention CDC, U.S.A.) and a result was considered significant if p < 0.05. Statistics are performed upon the individuals with available data in the markers under comparison; therefore the final number could differ among different analyses. As our working hypothesis was to verify in the Spanish population the presumed special association with T1D of the characteristic alleles of the AH 18.2, there was no need for Bonferroni corrections in these initial comparisons.
type 1 diabetes
Major Histocompatibility Complex
human leukocyte antigen
tumor necrosis factor
We thank Carmen Martínez for expert technical assistance. Elena Urcelay is recipient of a "Ramón y Cajal" contract of the Spanish Science and Technology Ministry. Alfonso Martínez is a recipient of a FIS contract (CP04/00175).
- Mathis D, Vence L, Benoist C: beta-Cell death during progression to diabetes. Nature. 2001, 414 (6865): 792-798. 10.1038/414792a.PubMedView ArticleGoogle Scholar
- Field LL: Genetic linkage and association studies of Type I diabetes: challenges and rewards. Diabetologia. 2002, 45 (1): 21-35. 10.1007/s125-002-8241-7.PubMedView ArticleGoogle Scholar
- Robinson WP, Barbosa J, Rich SS, Thomson G: Homozygous parent affected sib pair method for detecting disease predisposing variants: application to insulin dependent diabetes mellitus. Genet Epidemiol. 1993, 10 (5): 273-288.PubMedView ArticleGoogle Scholar
- Lie BA, Todd JA, Pociot F, Nerup J, Akselsen HE, Joner G, Dahl-Jorgensen K, Ronningen KS, Thorsby E, Undlien DE: The predisposition to type 1 diabetes linked to the human leukocyte antigen complex includes at least one non-class II gene. Am J Hum Genet. 1999, 64 (3): 793-800. 10.1086/302283.PubMedPubMed CentralView ArticleGoogle Scholar
- Degli-Esposti MA, Leaver AL, Christiansen FT, Witt CS, Abraham LJ, Dawkins RL: Ancestral haplotypes: conserved population MHC haplotypes. Hum Immunol. 1992, 34 (4): 242-252. 10.1016/0198-8859(92)90023-G.PubMedView ArticleGoogle Scholar
- Johansson S, Lie BA, Todd JA, Pociot F, Nerup J, Cambon-Thomsen A, Kockum I, Akselsen HE, Thorsby E, Undlien DE: Evidence of at least two type 1 diabetes susceptibility genes in the HLA complex distinct from HLA-DQB1, -DQA1 and -DRB1. Genes Immun. 2003, 4 (1): 46-53. 10.1038/sj.gene.6363917.PubMedView ArticleGoogle Scholar
- De La Concha EG, Fernandez-Arquero M, Gual L, Vigil P, Martinez A, Urcelay E, Ferreira A, Garcia-Rodriguez MC, Fontan G: MHC Susceptibility Genes to IgA Deficiency Are Located in Different Regions on Different HLA Haplotypes. J Immunol. 2002, 169 (8): 4637-4643.PubMedView ArticleGoogle Scholar
- Fernandez L, Fernandez-Arquero M, Gual L, Lazaro F, Maluenda C, Polanco I, Figueredo MA, Gomez de la Concha E: Triplet repeat polymorphism in the transmembrane region of the MICA gene in celiac disease. Tissue Antigens. 2002, 59 (3): 219-222. 10.1034/j.1399-0039.2002.590307.x.PubMedView ArticleGoogle Scholar
- Lie BA, Sollid LM, Ascher H, Ek J, Akselsen HE, Ronningen KS, Thorsby E, Undlien DE: A gene telomeric of the HLA class I region is involved in predisposition to both type 1 diabetes and coeliac disease. Tissue Antigens. 1999, 54 (2): 162-168. 10.1034/j.1399-0039.1999.540207.x.PubMedView ArticleGoogle Scholar
- Koeleman BV, De Groot KN, Van Der Slik AR, Roep BO, Giphart MJ: Association between D6S2223 and type I diabetes independent of HLA class II in Dutch families. Diabetologia. 2002, 45 (4): 598-599. 10.1007/s00125-001-0725-1.PubMedView ArticleGoogle Scholar
- Zavattari P, Lampis R, Motzo C, Loddo M, Mulargia A, Whalen M, Maioli M, Angius E, Todd JA, Cucca F: Conditional linkage disequilibrium analysis of a complex disease superlocus, IDDM1 in the HLA region, reveals the presence of independent modifying gene effects influencing the type 1 diabetes risk encoded by the major HLA-DQB1, -DRB1 disease loci. Hum Mol Genet. 2001, 10 (8): 881-889. 10.1093/hmg/10.8.881.PubMedView ArticleGoogle Scholar
- Degli-Esposti MA, Abraham LJ, McCann V, Spies T, Christiansen FT, Dawkins RL: Ancestral haplotypes reveal the role of the central MHC in the immunogenetics of IDDM. Immunogenetics. 1992, 36 (6): 345-356. 10.1007/BF00218041.PubMedView ArticleGoogle Scholar
- Moghaddam PH, Zwinderman AH, de Knijff P, Roep BO, Schipper RF, Van der Auwera B, Naipal A, Gorus F, Schuit F, Giphart MJ: TNFa microsatellite polymorphism modulates the risk of IDDM in Caucasians with the high-risk genotype HLA DQA1*0501- DQB1*0201/DQA1*0301-DQB1*0302. Belgian Diabetes Registry. Diabetes. 1997, 46 (9): 1514-1515.PubMedView ArticleGoogle Scholar
- Mandrup-Poulsen T, Bendtzen K, Dinarello CA, Nerup J: Human tumor necrosis factor potentiates human interleukin 1-mediated rat pancreatic beta-cell cytotoxicity. J Immunol. 1987, 139 (12): 4077-4082.PubMedGoogle Scholar
- Campbell IL, Iscaro A, Harrison LC: IFN-gamma and tumor necrosis factor-alpha. Cytotoxicity to murine islets of Langerhans. J Immunol. 1988, 141 (7): 2325-2329.PubMedGoogle Scholar
- Suk K, Kim S, Kim YH, Kim KA, Chang I, Yagita H, Shong M, Lee MS: IFN-gamma/TNF-alpha synergism as the final effector in autoimmune diabetes: a key role for STAT1/IFN regulatory factor-1 pathway in pancreatic beta cell death. J Immunol. 2001, 166 (7): 4481-4489.PubMedView ArticleGoogle Scholar
- Wong FS, Wen L: The Study of HLA Class II and Autoimmune Diabetes. Curr Top Med Chem. 2003, 3 (1): 1-15.View ArticleGoogle Scholar
- Escribano-de-Diego J, Sanchez-Velasco P, Luzuriaga C, Ocejo-Vinyals JG, Paz-Miguel JE, Leyva-Cobian F: HLA class II immunogenetics and incidence of insulin-dependent diabetes mellitus in the population of Cantabria (Northern Spain). Hum Immunol. 1999, 60 (10): 990-1000. 10.1016/S0198-8859(99)00077-4.PubMedView ArticleGoogle Scholar
- Cucca F, Lampis R, Congia M, Angius E, Nutland S, Bain SC, Barnett AH, Todd JA: A correlation between the relative predisposition of MHC class II alleles to type 1 diabetes and the structure of their proteins. Hum Mol Genet. 2001, 10 (19): 2025-2037. 10.1093/hmg/10.19.2025.PubMedView ArticleGoogle Scholar
- Kimura A, Sasazuki T: Eleventh International Histocompability workshop reference protocol for the HLA DNA-typing technique. HLA 1991 Proceedings of the Eleventh International Histocompatibility Workshop and Conference. Edited by: Tsuji K, Aizawa M, Sasazuki T. 1991, Oxford: Oxford University Press, 1: 397-419.Google Scholar
- Nedospasov SA, Udalova IA, Kuprash DV, Turetskaya RL: DNA sequence polymorphism at the human tumor necrosis factor (TNF) locus. Numerous TNF/lymphotoxin alleles tagged by two closely linked microsatellites in the upstream region of the lymphotoxin (TNF-beta) gene. J Immunol. 1991, 147 (3): 1053-1059.PubMedGoogle Scholar
- Martin M, Mann D, Carrington M: Recombination rates across the HLA complex: use of microsatellites as a rapid screen for recombinant chromosomes. Hum Mol Genet. 1995, 4 (3): 423-428.PubMedView ArticleGoogle Scholar
- Ota M, Katsuyama Y, Mizuki N, Ando H, Furihata K, Ono S, Pivetti-Pezzi P, Tabbara KF, Palimeris GD, Nikbin B: 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, Kazakhstan, Iranian, Saudi Arabian, Greek and Italian. Tissue Antigens. 1997, 49 (5): 448-454.PubMedView ArticleGoogle Scholar
- Vorechovsky I, Cullen M, Carrington M, Hammarstrom L, Webster AD: Fine mapping of IGAD1 in IgA deficiency and common variable immunodeficiency: identification and characterization of haplotypes shared by affected members of 101 multiple-case families. J Immunol. 2000, 164 (8): 4408-4416.PubMedView ArticleGoogle Scholar
- Foissac A, Salhi M, Cambon-Thomsen A: Microsatellites in the HLA region: 1999 update. Tissue Antigens. 2000, 55 (6): 477-509. 10.1034/j.1399-0039.2000.550601.x.PubMedView ArticleGoogle Scholar