HLA-E gene polymorphism associates with ankylosing spondylitis in Sardinia
- Fabiana Paladini†1,
- Francesca Belfiore†1, 2,
- Elisa Cocco1,
- Carlo Carcassi3,
- Alberto Cauli4,
- Alessandra Vacca4,
- Maria Teresa Fiorillo1,
- Alessandro Mathieu4,
- Isabella Cascino2 and
- Rosa Sorrentino1, 5Email author
© Paladini et al.; licensee BioMed Central Ltd. 2009
Received: 7 August 2009
Accepted: 13 November 2009
Published: 13 November 2009
Ankylosing spondylitis (AS) is a severe, chronic inflammatory disease strongly associated with HLA-B27. The presence of additional HLA risk factors has been suggested by several studies. The aim of the current study is to assess the occurrence of an additional HLA susceptibility locus in the region between HLA-E and HLA-C in the Sardinian population.
200 random controls, 120 patients with AS and 175 HLA-B27 positive controls were genotyped for six single nucleotide polymorphisms (SNPs) spanning the HLA region between HLA-E and HLA-C loci previously shown to harbour an additional susceptibility locus for AS. Allele, genotype and haplotype frequencies were compared.
The data confirm our previous finding of a significant increase in patients with AS of allele A at SNP rs1264457 encoding for an Arg at the functional HLA-E polymorphism (Arg128/Gly128). This was due to a remarkable increase in the frequency of genotype A/A in patients vs HLA-B27-matched controls (51% vs 29%; P for trend: 5 × 10-5). Genotype distribution of three other SNPs mapping in genes (GNL1, PRR3 and ABCF-1) close to HLA-E and showing high LD with it, was also significantly skewed. Accordingly, haplotype distribution was also remarkably different. The frequency of the haplotype AAGA, is 42% in random controls, increases to 53% in the HLA-B27-positive controls, and reaches 68% in patients with AS (P values: 2 × 10-11 vs random and 3 × 10-4 vs HLA-B27 controls).
There is a strong association between the presence of a haplotype in genes mapping between HLA-E and HLA-C and AS due to an increase of homozygous markers in patients. The strongest association however, is with the HLA-E functional polymorphism rs1264457. Since HLA-E is the ligand for the NKG2A receptor, these data point to the natural killer (NK) activity as possible player in the pathogenesis of AS.
Ankylosing spondylitis (AS) is a severe, chronic inflammatory disease whose strong association with human leukocyte antigen (HLA)-B27 has been known for many years. However, the presence of HLA-B27 far from accounts for the entire genetic load in AS. Other genes within and outside the HLA region have been evoked as additional factors contributing to the disease [1–4]. Although the association of genes outside the HLA region, such as IL-1, IL23R or endoplasmic reticulum aminopeptidase 1 (ERAP1), has been demonstrated, it appears more difficult to single out independent HLA risk factors, due to the strong Linkage Disequilibrium (LD) among genes mapping within the HLA region. The Sardinian is an outlier population and, due to the relative low gene flow from outside populations , is particularly suitable for case-control association studies in common complex diseases. In previous work, by comparing HLA-B27 -positive patients with AS and controls, we have reported that the region between HLA-E and major histocompatibility class I polypeptide-related sequence A (MICA) contains genes showing different allele frequencies between the two cohorts, and that there was an excess of homozygous markers in patients . This allowed us to speculate that the co-occurring chromosome was contributing to the disease providing a double dose of the predisposing allele. Interestingly, among the markers analysed, we found that the functional single nucleotide polymorphism (SNP) in the HLA-E gene (rs1264457A/G; determining the amino acid variation Arg128Gly) showed a significantly different genotype distribution in the two cohorts. Here, we analyse a larger cohort of HLA-B27-positive patients (n = 120) and HLA-B27-positive controls (n = 175) for the distribution of rs1264457 and five other SNPs mapping in genes close to the HLA-E locus. This further analysis confirms the strong association between allele A at rs1264457 in the HLA-E gene and AS.
Materials and methods
In this study, 200 random controls, 120 HLA-B27-positive patients with AS and 175 HLA-B27-matched controls were analysed. Seventy-nine HLA-B27-positive patients and 77 controls (24 of which were typed as B*2709) have been previously described . A further 41 patients with AS and 98 HLA-B27-positive controls (16 of which typed as B*2709) were enrolled in this study. The final subtype distribution was as follow: 108 B*2705 (90%), 6 B*2702 (5%), 4 B*2707 (3%) and 2 B*2713 (2%) in the cohort of patients with AS; and 121 B*2705 (69%), 40 B*2709 (23%), 9 B*2707 (5%) and 5 B*2702 (3%) in the cohort of HLA-B27-positive controls. Population stratification had been controlled by analyzing the distribution of a SNP mapping at the 5' end (rs4988235) of the gene for lactase and no significant difference in genotype distribution between patients and controls was found . DNA samples from patients and controls were collected at the University of Cagliari, Cagliari, Sardinia. The B*2709 haplotype, not associated with AS , was not different from B*2705 for the markers analyzed here  and therefore B*2709 controls were included in the study. Genotype distribution of all six SNPs was in Hardy-Weinberg equilibrium. All patients were diagnosed according to the modified New York criteria . All subjects gave their informed consent to the study. The study was approved by the local ethics committee.
A 20 ng sample of genomic DNA in a volume of 15 ul, containing 20 mM Tris-HCl (pH 8.4), 50 mM KCl, 2.5 mM MgCl2, 0.2 mM dNTPs, 5 pmoles each of forward and reverse primer and 0.75 U of Taq Polymerase (Biotaq DNA Polymerase, Bioline (Bioline Ltd, London, UK) was used. The mixture was incubated at 94°C for three minutes followed by 25 to 30 cycles of denaturation at 94°C for 30 seconds, annealing at a range of 55 to 62°C (depending on the analysed marker) for 30 seconds, extension at 72°C for one minute, with a final step of 72°C for 10 minutes.
Five SNPs mapping in the region encompassing about 100 kb between the HLA-E and HLA-C genes have been analysed: rs2105960 in the member RAS oncogene family pseudogene 1 (RANP1) gene, rs1264457 in the HLA-E gene, rs2074505 in the guanine nucleotide binding protein-like 1 (GNL1) gene, rs2074503 in the proline-rich polypeptide 3 (PRR3) gene, and rs 2269709 and rs1264439 in the ATP-binding cassette, subfamily F (GCN20), member 1 (ABCF-1) gene. SNPs were typed by minisequencing as previously reported . Briefly, the amplified products were treated with 0.5 U of shrimp alkaline phosphatase (SAP; Roche, Gaillard, France) and 2.5 U of exonuclease I (EXO I; BioLabs, Frankfurt am Main, Germany) at 37°C for 120 minutes and at 75°C for 15 minutes for inactivation. For minisequencing reactions, the commercial fluorescent-based minisequencing kit SNaPshot multiplex (Applied Biosystems, Lincoln Centre Drive, Foster City, USA) was used, with 2 pmoles of primer. PCR for single base extension reaction was as follows: initial denaturation at 94°C for three minutes, followed by 50 cycles of primer extension at 96°C for 15 seconds, annealing at 55 to 59°C (depending on the oligonucleotide annealing temperature) for 15 seconds, extension at 60°C for one minute. After the extension and labelling reaction, the unincorporated dNTPs were removed by enzymatic treatment using 0.5 U of SAP at 37°C for 120 minutes and at 75°C for 15 minutes for inactivation. The products were analyzed on the ABI prism 310 Genetic Analysers (Applied Biosystems, Lincoln Centre Drive, Foster City, USA). Genotyping was performed by 310 ABI Prism GeneScan 2.1 software (Applied Biosystems,, Lincoln Centre Drive, Foster City, USA).
Primers were as follows: rs2105960: FOR: AGACGCTTCTGGAAGGAACA; REV: GGCTGTGTCCCATACATTGA; SEQ: ATGGTGGTACTGGAAAAACG. rs1264457: FOR: 5'-GTGGGCGGGACTGACTAAG; REV: 5'-GGTCCTCATTCAGGGTGAGA; SEQ: 5'-TGCGAGCTGGGGCCCGAC. rs2074505: FOR: CACAGTGATGGCTTCACAGGC; REV: 5'-CCTATAGGAAGCCCTCTGAGAC; SEQ: CCCCCGCACCCCACAGGA. rs2074503: FOR: GCACAAGAAACAGCACCAAA; REV: AGTGACACACCCATCCCATC; SEQ: AGGAAGCAAAGATGGGTCCC. rs2269709: FOR: TGCCGCCTCCTTAGTAACTC; REV: TTAGATGCCATTCCGAGGAG; SEQ: TGAACTGGTCAAGATGGAGG. rs1264439: FOR: GATTACAGGCGTGACCCATT; REV: CCCTTTTCCTCTTTCCCTTC; SEQ: CTCACCACCAAGTAGGGAGA.
Two-tailed Fisher's Exact test was used to assess the differences of proportion of polymorphic alleles and association between random controls and HLA-B27 patients, HLA-B*2705 and HLA-B*2709 controls. The Cochrane-Armitage trend test was performed using StatXact7 software (Massachusetts Ave, Cambridge, USA). P values less than 0.05 were considered significant. LD was calculated using Haploview software . Hardy-Weinberg equilibrium testing was performed using Pearson's goodness-of-fit chi-squared test with one degree of freedom.
Allele frequencies and Cochrane-Armitage test for trend of genotype distribution in patients with AS and HLA-B27-matched controls in Sardinia
n = 175
HLA-B27+ patients with AS
n = 120
Position on chromosome 6
1 × 10-4
1 × 10-3
4 × 10-4
5 × 10 -5
4 × 10-4
4 × 10-3
4 × 10 -4
9 × 10-3
7 × 10 -3
5 × 10-3
4 × 10 -3
Haplotype distribution in patients with AS and controls
patients with AS
n = 135
n = 40
9 × 10-4
2 × 10-11
3 × 10-4
1 × 10-4
7 × 10-4
The data reported here validate our previous observation of a strong association between a functional polymorphism of the HLA-E gene and AS in Sardinia. It has been shown that the signal peptide of the HLA-B27 molecules binds with low affinity to the HLA-E Arg128 variant  and therefore the absence of the alternative variant, Gly128 in patients, could result in a reduced expression of the HLA-E molecules on the cell surface. This, in turn, could lead to an imbalance of the inhibitory function of HLA-E on natural killer (NK) activity or to an alteration of regulatory functions as it has been shown in animal models of autoimmune diseases . More and more reports are involving the HLA-E-NKG2A pathway in autoimmune/inflammatory diseases. Recently, a polymorphism mapping in the HLA-E gene has been found to be associated with Kawasaki disease  and the HLA-E-NKG2A pathway is likely to play a role in celiac disease . From the data reported here, it seems unlikely that a HLA-E neighboring gene rather than HLA-E itself, could be the real factor contributing to the disease. The data reported here strongly suggest that NK activity could be involved in the pathogenesis of AS  although the HLA-E molecules themselves are not necessarily responsible for this in all populations. As an example, two studies performed in a Spanish and an Asian population reports an imbalance of KIR alleles between patients with AS and HLA-B27-matched controls but these data have not been replicated in a further study in an English cohort [16–18]. This might indicate that different pathways have been selected in different populations based on their genetic background and environmental pressure, leading to the impairment of a crucial function, which can be imbalanced towards an exceeding effector or a defective regulatory function. If this is the case, the more promiscuous the population is, the more difficult it can be to identify some of the minor predisposing genes. For this reason, isolated populations which have undergone strong selection by environmental factors as it occurred in Sardinia, might be particularly informative to single out genes that predispose to complex diseases particularly in those cases in which a gene with a strong effect, such as HLA-B27, plays a leading role .
There is a strong association between a gene mapping in the region between HLA-E and HLA-C and AS in Sardinia. This gene is most likely to be the HLA-E itself because the functional polymorphism rs1264457 encoding for an Arg at position 128 found to be functionally relevant shows the strongest association.
ATP-binding cassette, subfamily F (GCN20), member 1
endoplasmic reticulum aminopeptidase 1
guanine nucleotide binding protein-like 1
human leukocyte antigen
polymerase chain reaction
proline-rich polypeptide 3
member RAS oncogene family pseudogene 1
shrimp alkaline phosphatase
single nucleotide polymorphism.
The authors wish to thank the patients for their cooperation, and Federica Lucantoni and Silvana Caristi for their excellent technical assistance. This work was partially supported by Fondazione Pasteur-Cenci Bolognetti.
- Wellcome Trust Case Control Consortium; Australo-Anglo-American Spondylitis Consortium (TASC): Association scan of 14,500 nonsynonymous SNPs in four diseases identifies autoimmunity variants. Nat Genet. 2007, 39: 1329-1337. 10.1038/ng.2007.17.View ArticleGoogle Scholar
- Zhu X, Wang Y, Sun L, Song Y, Sun F, Tang L, Huo Z, Li J, Yang Z: A novel gene variation of TNFalpha associated with ankylosing spondylitis: a reconfirmed study. Ann Rheum Dis. 2007, 66: 1419-1422. 10.1136/ard.2006.068528.PubMed CentralView ArticlePubMedGoogle Scholar
- Brown MA: Genetics and the pathogenesis of ankylosing spondylitis. Curr Opin Rheumatol. 2009, 21: 318-323. 10.1097/BOR.0b013e32832b3795.View ArticlePubMedGoogle Scholar
- Sims AM, Barnardo M, Herzberg I, Bradbury L, Calin A, Wordsworth BP, Darke C, Brown MA: Non-B27 MHC associations of ankylosing spondylitis. Genes Immun. 2007, 8: 115-123. 10.1038/sj.gene.6364362.View ArticlePubMedGoogle Scholar
- Contu D, Morelli L, Santoni F, Foster JW, Francalacci P, Cucca F: Y-chromosome based evidence for pre-neolithic origin of the genetically homogeneous but diverse Sardinian population: inference for association scans. PLoS One. 2008, 3 (): e1430-10.1371/journal.pone.0001430.PubMed CentralView ArticlePubMedGoogle Scholar
- Cascino I, Paladini F, Belfiore F, Cauli A, Angelini C, Fiorillo MT, Mathieu A, Sorrentino R: Identification of previously unrecognized predisposing factors for ankylosing spondylitis from analysis of HLA-B27 extended haplotypes in Sardinia. Arthritis Rheum. 2007, 56: 2640-2651. 10.1002/art.22820.View ArticlePubMedGoogle Scholar
- Paladini F, Cocco E, Cauli A, Cascino I, Vacca A, Belfiore F, Fiorillo MT, Mathieu A, Sorrentino R: A functional polymorphism of the vasoactive intestinal peptide receptor 1 gene correlates with the presence of HLA-B*2705 in Sardinia. Genes Immun. 2008, 9: 659-667. 10.1038/gene.2008.60.View ArticlePubMedGoogle Scholar
- D'Amato M, Fiorillo MT, Carcassi C, Mathieu A, Zuccarelli A, Bitti PP, Tosi R, Sorrentino R: Relevance of residue 116 of HLA-B27 in determining susceptibility to ankylosing spondylitis. Eur J Immunol. 1995, 25: 3199-3201. 10.1002/eji.1830251133.View ArticlePubMedGoogle Scholar
- Linden Van der SM, Valkenburg HA, Cats A: Evaluation of diagnostic criteria for ankylosing spondylitis: a proposal for modification of the New York criteria. Arthritis Rheum. 1984, 27: 361-368. 10.1002/art.1780270401.View ArticlePubMedGoogle Scholar
- Haploview software. [http://www.broad.mit.edu/mpg/haploview/]
- Strong RK, Holmes MA, Li P, Braun L, Lee N, Geraghty E: HLA-E allelic variants. Correlating differential expression, peptide affinities, crystal structures, and thermal stabilities. J Biol Chem. 2003, 278: 5082-5090. 10.1074/jbc.M208268200.View ArticlePubMedGoogle Scholar
- Lu L, Kim HJ, Werneck MB, Cantor H: Regulation of CD8+ regulatory T cells: Interruption of the NKG2A-Qa-1 interaction allows robust suppressive activity and resolution of autoimmune disease. Proc Natl Acad Sci USA. 2008, 105: 19420-19425. 10.1073/pnas.0810383105.PubMed CentralView ArticlePubMedGoogle Scholar
- Lin YJ, Wan L, Wu JY, Sheu JJ, Lin CW, Lan YC, Lai CH, Hung CH, Tsai Y, Tsai CH, Lin TH, Lin JG, Hsueh KC, Huang YM, Chang JS, Tsai FJ: HLA-E gene polymorphism associated with susceptibility to Kawasaki disease and formation of coronary artery aneurysms. Arthritis Rheum. 2009, 60: 604-610. 10.1002/art.24261.View ArticlePubMedGoogle Scholar
- Bhagat G, Naiyer AJ, Shah JG, Harper J, Jabri B, Wang TC, Green PH, Manavalan JS: Small intestinal CD8+TCRgammadelta+NKG2A+ intraepithelial lymphocytes have attributes of regulatory cells in patients with celiac disease. J Clin Invest. 2008, 118: 281-293. 10.1172/JCI30989.PubMed CentralView ArticlePubMedGoogle Scholar
- Kollnberger S, Chan A, Sun MY, Chen LY, Wright C, di Gleria K, McMichael A, Bowness P: Interaction of HLA-B27 homodimers with KIR3DL1 and KIR3DL2, unlike HLA-B27 heterotrimers, is independent of the sequence of bound peptide. Eur J Immunol. 2007, 37: 1313-1322. 10.1002/eji.200635997.View ArticlePubMedGoogle Scholar
- Lopez-Larrea C, Blanco-Gelaz MA, Torre-Alonso JC, Bruges Armas J, Suarez-Alvarez B, Pruneda L, Couto AR, Gonzalez S, Lopez-Vázquez A, Martinez-Borra J: Contribution of KIR3DL1/3DS1 to ankylosing spondylitis in human leukocyte antigen-B27 Caucasian populations. Arthritis Res Ther. 2006, 8: R101-10.1186/ar1858.PubMed CentralView ArticlePubMedGoogle Scholar
- Díaz-Peña R, Blanco-Gelaz MA, Suárez-Alvarez B, Martínez-Borra J, López-Vázquez A, Alonso-Arias R, Brüges-Armas J, Vidal-Castiñeira JR, López-Larrea C: Activating KIR genes are associated with ankylosing spondylitis in Asian populations. Hum Immunol. 2008, 69: 437-442. 10.1016/j.humimm.2008.04.012.View ArticlePubMedGoogle Scholar
- Harvey D, Pointon JJ, Sleator C, Meenagh A, Farrar C, Sun JY, Senitzer D, Middleton D, Brown MA, Wordsworth BP: Analysis of killer immunoglobulin-like receptor genes in ankylosing spondylitis. Ann Rheum Dis. 2009, 68: 595-598. 10.1136/ard.2008.095927.View ArticlePubMedGoogle Scholar
- Mathieu A, Paladini F, Vacca A, Cauli A, Fiorillo MT, Sorrentino R: The interplay between the geographic distribution of HLA-B27 alleles and their role in infectious and autoimmune diseases: a unifying hypothesis. Autoimmun Rev. 2009, 8: 420-425. 10.1016/j.autrev.2009.01.003.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.