Background The present study is an analysis of the frequencies of HLA-A and -B antigens and HLA haplotypes in two groups of individuals homozygous for the two main  HFE  mutations (C282Y and H63D) and a group heterozygous for the S65C mutation.
Background
The present study is an analysis of the frequencies of HLA-A and -B antigens and HLA haplotypes in two groups of individuals homozygous for the two main  HFE  mutations (C282Y and H63D) and a group heterozygous for the S65C mutation.HFEMethods The study population includes: 1123 healthy individuals, 100 homozygous for the C282Y mutation, 138 homozygous for the H63D mutation and 17 heterozygous for the S65C mutation.  HFE  and HLA alleles were detected using DNA-based and microlymphocytotoxicity techniques respectively.
Methods
The study population includes: 1123 healthy individuals, 100 homozygous for the C282Y mutation, 138 homozygous for the H63D mutation and 17 heterozygous for the S65C mutation.  HFE  and HLA alleles were detected using DNA-based and microlymphocytotoxicity techniques respectively.HFEResults An expected significant association between C282Y and the HLA-A3/B7 haplotype was found, but other HLA haplotypes carrying the -A3 antigen were found: HLA-A3/B62 and HLA-A3/B44. Also, a significant association between H63D mutation and HLA-A29/B44 haplotype was found, and again other HLA haplotypes carrying the HLA-A29 antigen were also found: HLA-A29/B14 and HLA-A29/B62. In addition, the S65C mutation seems to be associated with a HLA haplotype carrying the HLA-A26 antigen.
Results
An expected significant association between C282Y and the HLA-A3/B7 haplotype was found, but other HLA haplotypes carrying the -A3 antigen were found: HLA-A3/B62 and HLA-A3/B44. Also, a significant association between H63D mutation and HLA-A29/B44 haplotype was found, and again other HLA haplotypes carrying the HLA-A29 antigen were also found: HLA-A29/B14 and HLA-A29/B62. In addition, the S65C mutation seems to be associated with a HLA haplotype carrying the HLA-A26 antigen.Conclusion These findings clearly suggest that HLA-A3/B7 and HLA-A29/B44 are the ancestral haplotypes from which the C282Y and H63D mutations originated, respectively. The frequencies of these mutations in different populations, their geographical distribution, and the degree of the statistical association to the ancestral haplotypes, suggest that the H63D mutation must have occurred earlier than the C282Y mutation.
Conclusion
These findings clearly suggest that HLA-A3/B7 and HLA-A29/B44 are the ancestral haplotypes from which the C282Y and H63D mutations originated, respectively. The frequencies of these mutations in different populations, their geographical distribution, and the degree of the statistical association to the ancestral haplotypes, suggest that the H63D mutation must have occurred earlier than the C282Y mutation.

Background Hereditary hemochromatosis (HH) is an autosomal recessive disease common in northern European populations. HH is characterized by increased iron absorption and deposition in the liver, pancreas, heart, joints and pituitary gland. Without treatment, death may occur from cirrhosis, primary liver cancer, diabetes or cardiomyopathy. In 1996, the HH gene ( HFE ) was cloned and located on the short arm of chromosome 6 (6p21.3) [ 1 ], 4 megabases (Mb) telomeric to the major histocompatibility complex (MHC). A single point mutation 845 G>A (exon 4), changing cysteine at position 282 to tyrosine (C282Y) is identified in 85 to 100% of patients with HH in populations of northern European descent, but it is found in only 60% of cases from Mediterranean populations [ 2 ]. Two other mutations, 187 C>G (exon 2), a histidine to aspartate change at amino acid 63 (H63D) and 193 A>T (exon 2), a serine to cysteine change at amino acid 65 (S65C) appear to be associated with milder forms of HH and may increase risk of disease in persons heterozygous for C282Y mutation [ 3 , 4 ]. C282Y lies within a Celtic ancestral haplotype which includes the human MHC (HLA) haplotype HLA-A3/B7 [ 5 ]. The HLA-A3/B7 haplotype was reported in HH patients in many European and populations of European descent [ 5 - 8 ] but HLA-A3/B14, HLA-A3/B35 and others were also reported [ 9 - 11 ]. The predominance of the HLA-A3 associated haplotypes on hemochromatosis chromosomes, and the pattern of their distribution in the world, led Simon et al [ 5 ] to propose the founder hypothesis, postulating that the hemochromatosis mutation was a rare event that occurred once on a particular chromosome which was subsequently modified by recombinations involving both HLA-B and HLA-A alleles and population migrations, producing the varied haplotype associations that were described. In contrast, the H63D substitution is not restricted to European populations: allele frequencies greater than 5% are found in countries bordering the Mediterranean, in the Middle East, and in the Indian subcontinent. H63D seems to be associated with HLA-A29 and HLA-B44 [ 3 , 12 ]. A few studies have been performed on the distribution of the S65C mutation in Europe and other countries. In 1999, Barton et al [ 13 ] identified the S65C mutation in Alabama hemochromatosis probands and found that this mutation was linked to a haplotype characterized by HLA-A32; and recently, Couto et al [ 14 ] found that this mutation is in linkage disequilibrium with the HLA-A29/B44 haplotype. The aim of this study is to find the HLA antigens and haplotypes associated with the three main mutations in the  HFE  gene in a sample of the Spanish population.
Background
Hereditary hemochromatosis (HH) is an autosomal recessive disease common in northern European populations. HH is characterized by increased iron absorption and deposition in the liver, pancreas, heart, joints and pituitary gland. Without treatment, death may occur from cirrhosis, primary liver cancer, diabetes or cardiomyopathy. In 1996, the HH gene ( HFE ) was cloned and located on the short arm of chromosome 6 (6p21.3) [ 1 ], 4 megabases (Mb) telomeric to the major histocompatibility complex (MHC). A single point mutation 845 G>A (exon 4), changing cysteine at position 282 to tyrosine (C282Y) is identified in 85 to 100% of patients with HH in populations of northern European descent, but it is found in only 60% of cases from Mediterranean populations [ 2 ]. Two other mutations, 187 C>G (exon 2), a histidine to aspartate change at amino acid 63 (H63D) and 193 A>T (exon 2), a serine to cysteine change at amino acid 65 (S65C) appear to be associated with milder forms of HH and may increase risk of disease in persons heterozygous for C282Y mutation [ 3 , 4 ].HFEC282Y lies within a Celtic ancestral haplotype which includes the human MHC (HLA) haplotype HLA-A3/B7 [ 5 ]. The HLA-A3/B7 haplotype was reported in HH patients in many European and populations of European descent [ 5 - 8 ] but HLA-A3/B14, HLA-A3/B35 and others were also reported [ 9 - 11 ]. The predominance of the HLA-A3 associated haplotypes on hemochromatosis chromosomes, and the pattern of their distribution in the world, led Simon et al [ 5 ] to propose the founder hypothesis, postulating that the hemochromatosis mutation was a rare event that occurred once on a particular chromosome which was subsequently modified by recombinations involving both HLA-B and HLA-A alleles and population migrations, producing the varied haplotype associations that were described.In contrast, the H63D substitution is not restricted to European populations: allele frequencies greater than 5% are found in countries bordering the Mediterranean, in the Middle East, and in the Indian subcontinent. H63D seems to be associated with HLA-A29 and HLA-B44 [ 3 , 12 ].A few studies have been performed on the distribution of the S65C mutation in Europe and other countries. In 1999, Barton et al [ 13 ] identified the S65C mutation in Alabama hemochromatosis probands and found that this mutation was linked to a haplotype characterized by HLA-A32; and recently, Couto et al [ 14 ] found that this mutation is in linkage disequilibrium with the HLA-A29/B44 haplotype.The aim of this study is to find the HLA antigens and haplotypes associated with the three main mutations in the  HFE  gene in a sample of the Spanish population.HFEMethods Individuals A total of 100 unrelated individuals homozygous for the C282Y mutation, 138 unrelated individuals homozygous for the H63D mutation and 17 individuals heterozygous for the S65C mutation were selected for this study. These were subjects in whom  HFE  genotyping has been previously performed on a medical care basis because of a presumptive diagnosis of hemochromatosis. In addition, 1113 unrelated, apparently healthy subjects were used as controls for the study. In addition, HLA typing was performed in 230 individuals whom  HFE  genotyping was negative for the three mutations. HLA-A and -B typing HLA-class I typing was performed on freshly collected venous blood samples by the standard complement-dependent microlymphocytotoxicity assay using commercially available alloantisera. DNA extraction and HFE amplification Genomic DNA from whole blood samples was extracted by standard protocols. Polymerase chain reaction (PCR) with the pair of primers HEMEx2-5' (5'-CTT TGG GCT ACG TGG ATG ACC) and HEMEx2-3' (5'-CTG GCT TGA AAT TCT ACT GGA AAC C) was used to amplify exon 2 of the  HFE  gene. To amplify exon 4, a second set of oligonucleotides was used: HEMEx4-5' (5'-GGT GTC GGG CCT TGA ACT ACT ACC) and HEMEx4-3' (5'-A CAT ACC CCA GAT CAC AAT GAG G). The following conditions were used for the PCR reactions: five minutes denaturation at 94°C, 40 cycles of 15 seconds denaturation at 95°C, 15 seconds annealing at 57°C and 30 seconds extension at 72°C. PCR products coming from exons 2 and 4 were 101 and 228 base pairs (bp), respectively. Digestion with mutation-specific restriction endonuclease Following the PCR amplifications, aliquots (17 μl) of the reaction mixture were digested with the restriction endonucleases  Bcl  I (exon 2),  Hinf  I (exon 2) and  Rsa  I (exon 4) for 3 hours following the protocol recommended by the manufacturer (Promega, Madison, WI). The H63D mutation destroys the  Bcl  I site in the 101 bp PCR product, so while normal DNA is cut into two fragments of 38 and 63 bp, the mutated DNA is not cut. The S65C mutation destroys the  Hinf  I site in the 101 bp PCR product, so while normal DNA is cut into two fragments of 47 and 54 bp, the mutated DNA is not cut. The C282Y mutation creates a new  Rsa  I site, the 228 bp DNA product digested with this enzyme is cut into two fragments of 145 and 83 bp in the normal allele, while in the mutated DNA three fragments of 145, 29 and 54 bp are generated after digestion. The digested products were size resolved in 10% acrylamide gel and detected by staining with ethidium bromide. Statistical analysis Allele and haplotype frequencies were estimated using Arlequin V2.0 software [ 15 ]. The haplotype frequencies were computed using the Expectation-Maximization algorithm [ 16 ]; this procedure is an interactive process aimed at obtaining maximum-likelihood estimates of haplotype frequencies from multi-locus genotype data when the gametic phase is unknown. The existence of association between  HFE  mutations and HLA-A and -B alleles and haplotypes was calculated by 2 × 2 comparison tables and p values were corrected according to the number of alleles or haplotypes compared [ 17 ] and using Yates corrected Chi 2  and Fisher's tests. Odds ratios were calculated as previously described [ 18 ].
Methods
Individuals A total of 100 unrelated individuals homozygous for the C282Y mutation, 138 unrelated individuals homozygous for the H63D mutation and 17 individuals heterozygous for the S65C mutation were selected for this study. These were subjects in whom  HFE  genotyping has been previously performed on a medical care basis because of a presumptive diagnosis of hemochromatosis. In addition, 1113 unrelated, apparently healthy subjects were used as controls for the study. In addition, HLA typing was performed in 230 individuals whom  HFE  genotyping was negative for the three mutations.
Individuals
A total of 100 unrelated individuals homozygous for the C282Y mutation, 138 unrelated individuals homozygous for the H63D mutation and 17 individuals heterozygous for the S65C mutation were selected for this study. These were subjects in whom  HFE  genotyping has been previously performed on a medical care basis because of a presumptive diagnosis of hemochromatosis. In addition, 1113 unrelated, apparently healthy subjects were used as controls for the study. In addition, HLA typing was performed in 230 individuals whom  HFE  genotyping was negative for the three mutations.HFEHFEHLA-A and -B typing HLA-class I typing was performed on freshly collected venous blood samples by the standard complement-dependent microlymphocytotoxicity assay using commercially available alloantisera.
HLA-A and -B typing
HLA-class I typing was performed on freshly collected venous blood samples by the standard complement-dependent microlymphocytotoxicity assay using commercially available alloantisera.DNA extraction and HFE amplification Genomic DNA from whole blood samples was extracted by standard protocols. Polymerase chain reaction (PCR) with the pair of primers HEMEx2-5' (5'-CTT TGG GCT ACG TGG ATG ACC) and HEMEx2-3' (5'-CTG GCT TGA AAT TCT ACT GGA AAC C) was used to amplify exon 2 of the  HFE  gene. To amplify exon 4, a second set of oligonucleotides was used: HEMEx4-5' (5'-GGT GTC GGG CCT TGA ACT ACT ACC) and HEMEx4-3' (5'-A CAT ACC CCA GAT CAC AAT GAG G). The following conditions were used for the PCR reactions: five minutes denaturation at 94°C, 40 cycles of 15 seconds denaturation at 95°C, 15 seconds annealing at 57°C and 30 seconds extension at 72°C. PCR products coming from exons 2 and 4 were 101 and 228 base pairs (bp), respectively.
DNA extraction and HFE amplification
Genomic DNA from whole blood samples was extracted by standard protocols. Polymerase chain reaction (PCR) with the pair of primers HEMEx2-5' (5'-CTT TGG GCT ACG TGG ATG ACC) and HEMEx2-3' (5'-CTG GCT TGA AAT TCT ACT GGA AAC C) was used to amplify exon 2 of the  HFE  gene. To amplify exon 4, a second set of oligonucleotides was used: HEMEx4-5' (5'-GGT GTC GGG CCT TGA ACT ACT ACC) and HEMEx4-3' (5'-A CAT ACC CCA GAT CAC AAT GAG G).HFEThe following conditions were used for the PCR reactions: five minutes denaturation at 94°C, 40 cycles of 15 seconds denaturation at 95°C, 15 seconds annealing at 57°C and 30 seconds extension at 72°C. PCR products coming from exons 2 and 4 were 101 and 228 base pairs (bp), respectively.Digestion with mutation-specific restriction endonuclease Following the PCR amplifications, aliquots (17 μl) of the reaction mixture were digested with the restriction endonucleases  Bcl  I (exon 2),  Hinf  I (exon 2) and  Rsa  I (exon 4) for 3 hours following the protocol recommended by the manufacturer (Promega, Madison, WI). The H63D mutation destroys the  Bcl  I site in the 101 bp PCR product, so while normal DNA is cut into two fragments of 38 and 63 bp, the mutated DNA is not cut. The S65C mutation destroys the  Hinf  I site in the 101 bp PCR product, so while normal DNA is cut into two fragments of 47 and 54 bp, the mutated DNA is not cut. The C282Y mutation creates a new  Rsa  I site, the 228 bp DNA product digested with this enzyme is cut into two fragments of 145 and 83 bp in the normal allele, while in the mutated DNA three fragments of 145, 29 and 54 bp are generated after digestion. The digested products were size resolved in 10% acrylamide gel and detected by staining with ethidium bromide.
Digestion with mutation-specific restriction endonuclease
Following the PCR amplifications, aliquots (17 μl) of the reaction mixture were digested with the restriction endonucleases  Bcl  I (exon 2),  Hinf  I (exon 2) and  Rsa  I (exon 4) for 3 hours following the protocol recommended by the manufacturer (Promega, Madison, WI). The H63D mutation destroys the  Bcl  I site in the 101 bp PCR product, so while normal DNA is cut into two fragments of 38 and 63 bp, the mutated DNA is not cut. The S65C mutation destroys the  Hinf  I site in the 101 bp PCR product, so while normal DNA is cut into two fragments of 47 and 54 bp, the mutated DNA is not cut. The C282Y mutation creates a new  Rsa  I site, the 228 bp DNA product digested with this enzyme is cut into two fragments of 145 and 83 bp in the normal allele, while in the mutated DNA three fragments of 145, 29 and 54 bp are generated after digestion. The digested products were size resolved in 10% acrylamide gel and detected by staining with ethidium bromide.BclHinfRsaBclHinfRsaStatistical analysis Allele and haplotype frequencies were estimated using Arlequin V2.0 software [ 15 ]. The haplotype frequencies were computed using the Expectation-Maximization algorithm [ 16 ]; this procedure is an interactive process aimed at obtaining maximum-likelihood estimates of haplotype frequencies from multi-locus genotype data when the gametic phase is unknown. The existence of association between  HFE  mutations and HLA-A and -B alleles and haplotypes was calculated by 2 × 2 comparison tables and p values were corrected according to the number of alleles or haplotypes compared [ 17 ] and using Yates corrected Chi 2  and Fisher's tests. Odds ratios were calculated as previously described [ 18 ].
Statistical analysis
Allele and haplotype frequencies were estimated using Arlequin V2.0 software [ 15 ]. The haplotype frequencies were computed using the Expectation-Maximization algorithm [ 16 ]; this procedure is an interactive process aimed at obtaining maximum-likelihood estimates of haplotype frequencies from multi-locus genotype data when the gametic phase is unknown. The existence of association between  HFE  mutations and HLA-A and -B alleles and haplotypes was calculated by 2 × 2 comparison tables and p values were corrected according to the number of alleles or haplotypes compared [ 17 ] and using Yates corrected Chi 2  and Fisher's tests. Odds ratios were calculated as previously described [ 18 ].HFEResults The allele frequencies of the HLA-A and -B antigens found in the C282Y carriers group, in comparison with the frequencies in the control population are listed in Table  1 . As expected, significant associations were found for HLA-A3 and -B7, but HLA-B62 also shows significant association. The association of HLA-A3 is stronger than that of HLA-B7 or -B62. The frequencies of the HLA-A/B haplotypes in the homozygous group in comparison with the haplotype frequencies in the control population are listed in Table  2 . Three haplotypes are significantly associated with the C282Y mutation: HLA-A3/B7, HLA-A3/B62 and HLA-A3/B44. Table 1 Allele frequencies of HLA-A and HLA-B antigens in the C282Y homozygous group in comparison with the control group. C282Y HOMOZYGOTES CONTROLS C282Y HOMOZYGOTES CONTROLS HLA Freq (n = 200) OR p Freq (n = 2226) HLA Freq (n = 200) OR p Freq (n = 2226) A1 0.070 0.69 N.S. 0.097 B8 0.035 0.81 N.S. 0.044 A3 0.425 6.23 <10 -7 0.106 B7 0.260 3.67 <10 -7 0.089 B62 0.075 2.74 0.030 0.028 A29 0.035 0.54 N.S. 0.062 B44 0.170 1.19 N.S. 0.146 A30 0.010 0.16 0.093 0.060 B18 0.015 0.27 0.120 0.088 A2 0.195 0.67 N.S. 0.264 B35 0.095 0.77 N.S. 0.119 A11 0.060 0.75 N.S. 0.078 B51 0.055 0.60 N.S. 0.088 A26 0.020 0.42 N.S. 0.046 B49 0.025 0.67 N.S. 0.036 A28 0.025 0.59 N.S. 0.041 B60 0.015 0.42 N.S. 0.035 A32 0.020 0.57 N.S. 0.034 B14 0.010 0.28 N.S. 0.035 A23 0.030 0.95 N.S. 0.031 B38 0.010 0.33 N.S. 0.029 A33 0.010 0.36 N.S. 0.027 B27 0.025 0.88 N.S. 0.028 A25 0.020 0.97 N.S. 0.020 B50 0.005 0.21 N.S. 0.023 A31 0.025 1.56 N.S. 0.016 B65 0.015 0.64 N.S. 0.023 Freq: Frequency n: Total number of alleles in each group OR: Odds ratio p: Significance level N.S.: Not significant Table 2 Frequencies of HLA-A/B haplotypes in the two groups homozygous for H63D and C282Y in comparison with the control group. HAPLOTYPE H63D HAPLOTYPES Freq (n = 276) p OR C282Y HAPLOTYPES Freq (n = 200) p OR CONTROLS Freq (n = 2226) A1/B8 0.02513 N.S. 0.01500 N.S. 0.0246 A2/B7 0.03612 N.S. 0.02859 N.S. 0.0239 A2/B44 0.08254 N.S. 0.04037 N.S. 0.0443 A2/B51 0.04839 N.S. 0.00000 0.0404 A2/B35 0.01545 N.S. 0.03543 N.S. 0.0258 A3/B14 0.00000 0.00000 0.0005 A3/B7 0.00725 N.S. 0.20275 <10 -7 7.29 0.0343 A3/B62 0.00000 0.03379 10 -5 16.11 0.0024 A3/B44 0.00409 N.S. 0.08989 <10 - 12.13 0.0079 A11/B27 0.01812 N.S. 0.00000 0.0064 A11/B35 0.00362 N.S. 0.00862 N.S. 0.0204 A24/B35 0.00794 N.S. 0.01478 N.S. 0.0233 A29/B44 0.13261 <10 -7 3.93 0.02999 N.S. 0.0378 A29/B14 0.01449 <0.01 32.7 0.00000 0.0005 A29/B62 0.01591 <0.01 32.7 0.00000 0.0005 A30/B18 0.01087 N.S. 0.00000 0.0278 A33/B14 0.00000 0.00500 N.S. 0.0121 Freq: Frequency n: Total number of haplotypes in each group OR: Odds ratio p: Significance level N.S.: Not significant The allele frequencies of the HLA-A and -B antigens in the H63D homozygous group in comparison with the frequencies in the control population are listed in Table  3 . Significant associations were found for HLA-A29 and -B44. Again, the association of the HLA-A antigen (HLA-A29) is stronger than HLA-B44. The frequencies of the HLA-A/B haplotypes are also listed in Table  2 , and three HLA-A/B haplotypes (with the same HLA-A antigen, HLA-A29) are significantly associated with the H63D mutation: HLA-A29/B44, HLA-A29/B14 and HLA-A29/B62. Table 3 Allele frequencies of HLA-A and HLA-B antigens in the H63D homozygous group in comparison with the control group. H63D HOMOZYGOTES CONTROLS H63D HOMOZYGOTES CONTROLS HLA Freq (n = 276) OR p Freq (n = 2226) HLA Freq (n = 276) OR p Freq (n = 2226) A1 0.061 0.60 N.S. 0.097 B8 0.050 1.14 N.S. 0.044 A3 0.086 0.80 N.S. 0.106 B7 0.068 0.75 N.S. 0.089 B62 0.028 1.00 N.S. 0.028 A29 0.199 3.65 <10 -7 0.062 B44 0.264 2.10 2.10 -5 0.146 A30 0.036 0.58 N.S. 0.060 B18 0.086 0.98 N.S. 0.088 A2 0.286 1.12 N.S. 0.264 B35 0.076 0.60 N.S. 0.119 A11 0.036 0.44 N.S. 0.078 B51 0.097 1.12 N.S. 0.088 A26 0.029 0.61 N.S. 0.046 B49 0.032 0.89 N.S. 0.036 A28 0.021 0.52 N.S. 0.041 B60 0.011 0.30 N.S. 0.035 A32 0.039 1.16 N.S. 0.034 B14 0.014 0.40 N.S. 0.035 A23 0.047 1.52 N.S. 0.031 B38 0.032 1.10 N.S. 0.029 A33 0.007 0.26 N.S. 0.027 B27 0.036 1.29 N.S. 0.028 A25 0.032 1.60 N.S. 0.020 B50 0.018 0.76 N.S. 0.023 A31 0.011 0.67 N.S. 0.016 B65 0.039 1.74 N.S. 0.023 Freq: Frequency n: Total number of alleles in each group OR: Odds ratio p: Significance level N.S.: Not significant From the 17 individuals heterozygous for the S65C mutation, 7 (20%) were HLA-A26 versus 4.6% found in the control population (p = 0.02, OR = 5.29). No association with HLA-B has been found. The significant associations did not change if we used a control group of 230 individuals without  HFE  mutations.
Results
The allele frequencies of the HLA-A and -B antigens found in the C282Y carriers group, in comparison with the frequencies in the control population are listed in Table  1 . As expected, significant associations were found for HLA-A3 and -B7, but HLA-B62 also shows significant association. The association of HLA-A3 is stronger than that of HLA-B7 or -B62. The frequencies of the HLA-A/B haplotypes in the homozygous group in comparison with the haplotype frequencies in the control population are listed in Table  2 . Three haplotypes are significantly associated with the C282Y mutation: HLA-A3/B7, HLA-A3/B62 and HLA-A3/B44.Allele frequencies of HLA-A and HLA-B antigens in the C282Y homozygous group in comparison with the control group.A30.4256.23<10 -70.106B70.2603.67<10 -70.089B620.0752.740.0300.028Freq: Frequency n: Total number of alleles in each group OR: Odds ratio p: Significance level N.S.: Not significantFrequencies of HLA-A/B haplotypes in the two groups homozygous for H63D and C282Y in comparison with the control group.A3/B70.20275<10 -77.290.0343A3/B620.0337910 -516.110.0024A3/B440.08989<10 -12.130.0079A29/B440.13261<10 -73.930.0378A29/B140.01449<0.0132.70.0005A29/B620.01591<0.0132.70.0005Freq: Frequency n: Total number of haplotypes in each group OR: Odds ratio p: Significance level N.S.: Not significantThe allele frequencies of the HLA-A and -B antigens in the H63D homozygous group in comparison with the frequencies in the control population are listed in Table  3 . Significant associations were found for HLA-A29 and -B44. Again, the association of the HLA-A antigen (HLA-A29) is stronger than HLA-B44. The frequencies of the HLA-A/B haplotypes are also listed in Table  2 , and three HLA-A/B haplotypes (with the same HLA-A antigen, HLA-A29) are significantly associated with the H63D mutation: HLA-A29/B44, HLA-A29/B14 and HLA-A29/B62.Allele frequencies of HLA-A and HLA-B antigens in the H63D homozygous group in comparison with the control group.A290.1993.65<10 -70.062B440.2642.102.10 -50.146Freq: Frequency n: Total number of alleles in each group OR: Odds ratio p: Significance level N.S.: Not significantFrom the 17 individuals heterozygous for the S65C mutation, 7 (20%) were HLA-A26 versus 4.6% found in the control population (p = 0.02, OR = 5.29). No association with HLA-B has been found.The significant associations did not change if we used a control group of 230 individuals without  HFE  mutations.HFEDiscussion Populations of homozygous individuals for C282Y and H63D are optimal groups to study the HLA haplotypes in which these mutations preferentially appear. To our knowledge, Barton [ 19 ] and the present work are the only studies of associations between  HFE  mutations and HLA antigens and haplotypes in homozygous probands. The paper by Barton and Acton [ 19 ] presents haplotype frequencies assessed by family studies where phase could be set; in our paper, the haplotype frequences are estimated because the probands and controls are unrelated individuals. The low frequency of the S65C mutation makes the sampling of homozygous probands difficult and imposes the use of heterozygous individuals for the analysis. C282Y and HLA The strong association between the HLA-A3/B7 haplotype and the C282Y mutation indicates that this haplotype is the main one associated with this mutation in the Spanish population. However, other haplotypes are also associated: HLA-A3/B62 and HLA-A3/B44. This finding supports the founder hypothesis of Simon et al [ 5 ]: the ancestral haplotype where the C282Y mutation occurred on the ancestral haplotype HLA-A3/B7 and subsequent recombinations involving both HLA-B and HLA-A alleles produced the varied haplotype associations that have been found. Thus, we found many HLA-A/B haplotypes in our C282Y group, but only three HLA-A3 bearing haplotypes are statistically associated with this mutation. The two less frequent haplotypes (HLA-A3/B44 and HLA-A3/B62) have been observed in other populations in association with HH [ 20 , 21 ]. In addition, the high frequency of the HLA-A3/B7 haplotype makes other HLA antigens and haplotypes have reduced frequencies in respect to the controls. It is interesting to see that haplotypes with high-frequency in the Spanish population, such as HLA-A30/B18 and HLA-A2/B51 are absent in the C282Y homozygous group (Table  2 ). These HLA haplotypes are not  contaminated  by the C282Y mutation, and up until now, these haplotypes may be considered as  protector  haplotypes. H63D and HLA Porto et al [ 3 ] and Cardoso et al [ 22 ] found individual associations between HLA-A29 and non-classical forms of iron overload in linkage disequilibrium with H63D, and a strong linkage disequilibrium between H63D and all A29 containing haplotypes assigned in a large population of normal portuguese families. In the present work we found a strong association between the HLA-A29/B44 haplotype and the H63D mutation. Our finding agrees with the association of HLA-A29/B12(44) and hemochromatosis described in the Danish population [ 21 ]. H63D and HLA-A29-bearing haplotypes follow a pattern of associations similar to that described for C282Y and HLA-A3-bearing haplotypes. This promotes speculation that HLA-A29/B44 is the ancestral haplotype from which the H63D mutation emerged, since other HLA-A29 carrying haplotypes are also statistically associated with the mutation (HLA-A29/B14 and HLA-A29/B62, see Table  2 ), confirming the results reported by Cardoso et al [ 22 ] in the normal Portuguese population. Studies in other populations might lend support to whether HLA-A29/B44 is the ancestral haplotype of the H63D mutation, and HLA-A29/B14 and HLA-A29/B62 are specific Spanish haplotypes associated with H63D mutation. C282Y and H63D mutations: which one is older? In an attempt to establish the relative age of C282Y and H63D mutations, we have analysed the geographical distribution, allele frequencies and HLA haplotype associations for each mutation, assuming that the mutations ocurred once and that its age is directly proportional to its geographical spread, its frequency in the population, and the number of HLA haplotypes to which they are linked. On the other hand, a strong association of a particular  HFE  mutation to a particular HLA haplotype could mean that the mutation arose more recently, since the lower number of ruptures and recombinations of the original haplotype would reflect that a shorter time has passed. Merryweather-Clarke et al [ 23 ] analysed 2978 samples from probands distributed world-wide and showed that the C282Y mutation was most prevalent in northern European populations and absent from samples of non-European subjects (Africans, Asians, Australasians and Americans). In contrast, the H63D mutation is not restricted to European populations, being found in countries bordering the Mediterranean, in the Middle East and in the Indian subcontinent, and its allele frequency is higher than that of the C282Y mutation [ 23 ]. Our analysis of C282Y and H63D homozygous groups yielded a higher number of HLA haplotypes in association with the H63D mutation; and the frequency of HLA-A3/B7 in the C282Y homozygous group is 20%, while the frequency of HLA-A29/B44 in the H63D homozygous group is 13% (see table  2 ), reflecting that the HLA-A29/B44/H63D haplotype has suffered more recombinations than HLA-A3/B7/C282Y, and therefore, that HLA-A29/B44/H63D is older [ 24 ]. Altogether, these features suggest that the H63D mutation may have occurred earlier than the C282Y mutation, as has been previously proposed in studies from Italian populations [ 25 , 26 ]. S65C and HLA Few studies have been performed on the distribution, frequency and HLA association of the S65C mutation in Europe and other continents. Barton et al [ 13 ] described the linkage of the S65C mutation to a HLA-A32 haplotype in hemochromatosis probands from Alabama. Surprisingly, Couto et al [ 14 ] found the linkage of S65C (and not H63D) to the HLA-A29/B44 haplotype in a population from the Azores, although only 5 H63D homozygous and 9 S65C heterozygous individuals were studied in that work. In the present work we find that the S65C mutation seems to be linked to HLA-A26 in the Spanish population. Further studies in other populations and with more S65C-bearing haplotypes are necessary to shed light on the generation of the S65C mutation.
Discussion
Populations of homozygous individuals for C282Y and H63D are optimal groups to study the HLA haplotypes in which these mutations preferentially appear. To our knowledge, Barton [ 19 ] and the present work are the only studies of associations between  HFE  mutations and HLA antigens and haplotypes in homozygous probands. The paper by Barton and Acton [ 19 ] presents haplotype frequencies assessed by family studies where phase could be set; in our paper, the haplotype frequences are estimated because the probands and controls are unrelated individuals. The low frequency of the S65C mutation makes the sampling of homozygous probands difficult and imposes the use of heterozygous individuals for the analysis.HFEC282Y and HLA The strong association between the HLA-A3/B7 haplotype and the C282Y mutation indicates that this haplotype is the main one associated with this mutation in the Spanish population. However, other haplotypes are also associated: HLA-A3/B62 and HLA-A3/B44. This finding supports the founder hypothesis of Simon et al [ 5 ]: the ancestral haplotype where the C282Y mutation occurred on the ancestral haplotype HLA-A3/B7 and subsequent recombinations involving both HLA-B and HLA-A alleles produced the varied haplotype associations that have been found. Thus, we found many HLA-A/B haplotypes in our C282Y group, but only three HLA-A3 bearing haplotypes are statistically associated with this mutation. The two less frequent haplotypes (HLA-A3/B44 and HLA-A3/B62) have been observed in other populations in association with HH [ 20 , 21 ]. In addition, the high frequency of the HLA-A3/B7 haplotype makes other HLA antigens and haplotypes have reduced frequencies in respect to the controls. It is interesting to see that haplotypes with high-frequency in the Spanish population, such as HLA-A30/B18 and HLA-A2/B51 are absent in the C282Y homozygous group (Table  2 ). These HLA haplotypes are not  contaminated  by the C282Y mutation, and up until now, these haplotypes may be considered as  protector  haplotypes.
C282Y and HLA
The strong association between the HLA-A3/B7 haplotype and the C282Y mutation indicates that this haplotype is the main one associated with this mutation in the Spanish population. However, other haplotypes are also associated: HLA-A3/B62 and HLA-A3/B44. This finding supports the founder hypothesis of Simon et al [ 5 ]: the ancestral haplotype where the C282Y mutation occurred on the ancestral haplotype HLA-A3/B7 and subsequent recombinations involving both HLA-B and HLA-A alleles produced the varied haplotype associations that have been found. Thus, we found many HLA-A/B haplotypes in our C282Y group, but only three HLA-A3 bearing haplotypes are statistically associated with this mutation. The two less frequent haplotypes (HLA-A3/B44 and HLA-A3/B62) have been observed in other populations in association with HH [ 20 , 21 ].In addition, the high frequency of the HLA-A3/B7 haplotype makes other HLA antigens and haplotypes have reduced frequencies in respect to the controls. It is interesting to see that haplotypes with high-frequency in the Spanish population, such as HLA-A30/B18 and HLA-A2/B51 are absent in the C282Y homozygous group (Table  2 ). These HLA haplotypes are not  contaminated  by the C282Y mutation, and up until now, these haplotypes may be considered as  protector  haplotypes.contaminatedprotectorH63D and HLA Porto et al [ 3 ] and Cardoso et al [ 22 ] found individual associations between HLA-A29 and non-classical forms of iron overload in linkage disequilibrium with H63D, and a strong linkage disequilibrium between H63D and all A29 containing haplotypes assigned in a large population of normal portuguese families. In the present work we found a strong association between the HLA-A29/B44 haplotype and the H63D mutation. Our finding agrees with the association of HLA-A29/B12(44) and hemochromatosis described in the Danish population [ 21 ]. H63D and HLA-A29-bearing haplotypes follow a pattern of associations similar to that described for C282Y and HLA-A3-bearing haplotypes. This promotes speculation that HLA-A29/B44 is the ancestral haplotype from which the H63D mutation emerged, since other HLA-A29 carrying haplotypes are also statistically associated with the mutation (HLA-A29/B14 and HLA-A29/B62, see Table  2 ), confirming the results reported by Cardoso et al [ 22 ] in the normal Portuguese population. Studies in other populations might lend support to whether HLA-A29/B44 is the ancestral haplotype of the H63D mutation, and HLA-A29/B14 and HLA-A29/B62 are specific Spanish haplotypes associated with H63D mutation.
H63D and HLA
Porto et al [ 3 ] and Cardoso et al [ 22 ] found individual associations between HLA-A29 and non-classical forms of iron overload in linkage disequilibrium with H63D, and a strong linkage disequilibrium between H63D and all A29 containing haplotypes assigned in a large population of normal portuguese families.In the present work we found a strong association between the HLA-A29/B44 haplotype and the H63D mutation. Our finding agrees with the association of HLA-A29/B12(44) and hemochromatosis described in the Danish population [ 21 ].H63D and HLA-A29-bearing haplotypes follow a pattern of associations similar to that described for C282Y and HLA-A3-bearing haplotypes. This promotes speculation that HLA-A29/B44 is the ancestral haplotype from which the H63D mutation emerged, since other HLA-A29 carrying haplotypes are also statistically associated with the mutation (HLA-A29/B14 and HLA-A29/B62, see Table  2 ), confirming the results reported by Cardoso et al [ 22 ] in the normal Portuguese population. Studies in other populations might lend support to whether HLA-A29/B44 is the ancestral haplotype of the H63D mutation, and HLA-A29/B14 and HLA-A29/B62 are specific Spanish haplotypes associated with H63D mutation.C282Y and H63D mutations: which one is older? In an attempt to establish the relative age of C282Y and H63D mutations, we have analysed the geographical distribution, allele frequencies and HLA haplotype associations for each mutation, assuming that the mutations ocurred once and that its age is directly proportional to its geographical spread, its frequency in the population, and the number of HLA haplotypes to which they are linked. On the other hand, a strong association of a particular  HFE  mutation to a particular HLA haplotype could mean that the mutation arose more recently, since the lower number of ruptures and recombinations of the original haplotype would reflect that a shorter time has passed. Merryweather-Clarke et al [ 23 ] analysed 2978 samples from probands distributed world-wide and showed that the C282Y mutation was most prevalent in northern European populations and absent from samples of non-European subjects (Africans, Asians, Australasians and Americans). In contrast, the H63D mutation is not restricted to European populations, being found in countries bordering the Mediterranean, in the Middle East and in the Indian subcontinent, and its allele frequency is higher than that of the C282Y mutation [ 23 ]. Our analysis of C282Y and H63D homozygous groups yielded a higher number of HLA haplotypes in association with the H63D mutation; and the frequency of HLA-A3/B7 in the C282Y homozygous group is 20%, while the frequency of HLA-A29/B44 in the H63D homozygous group is 13% (see table  2 ), reflecting that the HLA-A29/B44/H63D haplotype has suffered more recombinations than HLA-A3/B7/C282Y, and therefore, that HLA-A29/B44/H63D is older [ 24 ]. Altogether, these features suggest that the H63D mutation may have occurred earlier than the C282Y mutation, as has been previously proposed in studies from Italian populations [ 25 , 26 ].
C282Y and H63D mutations: which one is older?
In an attempt to establish the relative age of C282Y and H63D mutations, we have analysed the geographical distribution, allele frequencies and HLA haplotype associations for each mutation, assuming that the mutations ocurred once and that its age is directly proportional to its geographical spread, its frequency in the population, and the number of HLA haplotypes to which they are linked. On the other hand, a strong association of a particular  HFE  mutation to a particular HLA haplotype could mean that the mutation arose more recently, since the lower number of ruptures and recombinations of the original haplotype would reflect that a shorter time has passed.HFEMerryweather-Clarke et al [ 23 ] analysed 2978 samples from probands distributed world-wide and showed that the C282Y mutation was most prevalent in northern European populations and absent from samples of non-European subjects (Africans, Asians, Australasians and Americans). In contrast, the H63D mutation is not restricted to European populations, being found in countries bordering the Mediterranean, in the Middle East and in the Indian subcontinent, and its allele frequency is higher than that of the C282Y mutation [ 23 ]. Our analysis of C282Y and H63D homozygous groups yielded a higher number of HLA haplotypes in association with the H63D mutation; and the frequency of HLA-A3/B7 in the C282Y homozygous group is 20%, while the frequency of HLA-A29/B44 in the H63D homozygous group is 13% (see table  2 ), reflecting that the HLA-A29/B44/H63D haplotype has suffered more recombinations than HLA-A3/B7/C282Y, and therefore, that HLA-A29/B44/H63D is older [ 24 ].Altogether, these features suggest that the H63D mutation may have occurred earlier than the C282Y mutation, as has been previously proposed in studies from Italian populations [ 25 , 26 ].S65C and HLA Few studies have been performed on the distribution, frequency and HLA association of the S65C mutation in Europe and other continents. Barton et al [ 13 ] described the linkage of the S65C mutation to a HLA-A32 haplotype in hemochromatosis probands from Alabama. Surprisingly, Couto et al [ 14 ] found the linkage of S65C (and not H63D) to the HLA-A29/B44 haplotype in a population from the Azores, although only 5 H63D homozygous and 9 S65C heterozygous individuals were studied in that work. In the present work we find that the S65C mutation seems to be linked to HLA-A26 in the Spanish population. Further studies in other populations and with more S65C-bearing haplotypes are necessary to shed light on the generation of the S65C mutation.
S65C and HLA
Few studies have been performed on the distribution, frequency and HLA association of the S65C mutation in Europe and other continents. Barton et al [ 13 ] described the linkage of the S65C mutation to a HLA-A32 haplotype in hemochromatosis probands from Alabama. Surprisingly, Couto et al [ 14 ] found the linkage of S65C (and not H63D) to the HLA-A29/B44 haplotype in a population from the Azores, although only 5 H63D homozygous and 9 S65C heterozygous individuals were studied in that work.In the present work we find that the S65C mutation seems to be linked to HLA-A26 in the Spanish population. Further studies in other populations and with more S65C-bearing haplotypes are necessary to shed light on the generation of the S65C mutation.Conclusions We have found that, in the Spanish population, the three main  HFE  mutations: C282Y, H63D and S65C, are in linkage disequilibrium with HLA haplotypes carrying the HLA-A3, -A29 and -A26 alleles, respectively. In addition, the ancestral HLA haplotypes from which C282Y and H63D mutations were originated are HLA-A3/B7 and HLA-A29/B44, respectively, and H63D is older than C282Y. Further studies in other populations using homozygous individuals for  HFE  mutations will help to identify the associated ancestral and specific haplotypes.
Conclusions
We have found that, in the Spanish population, the three main  HFE  mutations: C282Y, H63D and S65C, are in linkage disequilibrium with HLA haplotypes carrying the HLA-A3, -A29 and -A26 alleles, respectively. In addition, the ancestral HLA haplotypes from which C282Y and H63D mutations were originated are HLA-A3/B7 and HLA-A29/B44, respectively, and H63D is older than C282Y. Further studies in other populations using homozygous individuals for  HFE  mutations will help to identify the associated ancestral and specific haplotypes.HFEHFECompeting interests The author(s) declare that they have no competing interests.
Competing interests
The author(s) declare that they have no competing interests.Authors' contributions Authors AP and PM conceived the study, contributed to proband characterisation, performed statistical comparisons and edited the manuscript. Authors EM, MJR, MJC, DO, MGB and LG contributed at different times to the characterisation of probands and HLA typing.
Authors' contributions
Authors AP and PM conceived the study, contributed to proband characterisation, performed statistical comparisons and edited the manuscript. Authors EM, MJR, MJC, DO, MGB and LG contributed at different times to the characterisation of probands and HLA typing.Pre-publication history The pre-publication history for this paper can be accessed here:
Pre-publication history
The pre-publication history for this paper can be accessed here: