Rheumatology Advance Access published online on October 4, 2008
Rheumatology, doi:10.1093/rheumatology/ken374
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Interleukin-1 promoter region polymorphism role in rheumatoid arthritis: a meta-analysis of IL-1B-511A/G variant reveals association with rheumatoid arthritis
1Institute of Musculoskeletal Sciences and 2Cancer Research UK Epidemiology Unit, University of Oxford, Oxford, UK
Correspondence to:
Correspondence to: P. Harrison, Botnar Research Centre, Nuffield Orthopaedic Centre, Windmill Road, Oxford OX3 7LD, UK. E-mail: pille.harrison{at}ndos.ox.ac.uk
| Abstract |
|---|
|
|
|---|
Objectives. IL-1 has a central role mediating inflammation and joint destruction in RA. Single nucleotide polymorphisms (SNPs) and haplotype structure in the promoter region can modulate IL-1 function. This study examined the effects of four common promoter SNPs in the IL-1 region on susceptibility and clinical characteristics of RA in British Caucasian patients and assessed the risk of RA by meta-analysis of published studies.
Methods. Using PCR-based methods, 756 RA patients and 625 healthy controls (HCs) were genotyped for IL-1A (–889 C/A, rs17561), IL-1B (–511 A/G, rs16944), IL-1B (–1464 C/G, rs1143623) and IL-1B (–3737 G/A, rs4848306) SNPs. Further meta-analysis was performed for IL-1B (–511 A/G) incorporating 3712 RA patients and 2331 HC from six association studies.
Results. The IL-1B (–1464 C/G) G allele was found to be less common in the RA group [P = 0.01; odds ratio (OR) 1.24; 95% CI 1.04, 1.48]. There was no association between IL-1 SNPs and the presence of HLA-DRB1 shared epitope, RF or clinical characteristics. Meta-analysis revealed statistically significant association between IL-1B (–511 A/G) and RA (P = 0.02; pooled OR 1.13; 95% CI 1.02, 1.26).
Conclusions. There may be a protective effect in RA from the IL-1B (–1464 C/G) G variant. No direct association between the polymorphisms studied and clinical severity characteristics were observed. Further meta-analysis revealed IL-1B (–511 A/G) to be associated with increased susceptibility to RA.
KEY WORDS: Rheumatoid arthritis, Meta-analysis, Interleukin-1, Interleukin-1A, Interleukin-1B, Association study, Major histocompatibility complex, Radiographic severity, Genetics, Polymorphism
| Introduction |
|---|
|
|
|---|
Cytokines are small regulatory proteins that act as intercellular mediators in the generation and control of immune and inflammatory responses. In RA, two main cytokines, TNF and IL-1, are particularly involved. The most studied of the IL-1 proteins are IL-1
and IL-1β. They function as agonist molecules that activate target cells by binding to IL-1 receptor. The production and activity of IL-1
and IL-1β are regulated by transcriptional and post-translational mechanisms as well as from a natural competitive antagonist. The IL-1 gene cluster (comprising 10 related genes) maps to 2q12. Both IL-1A (MIM 147760
[OMIM]
) and IL-1B (MIM 147720
[OMIM]
) may have arisen through gene duplications, and hence encode similar protein products [1]. Strong linkage disequilibrium (LD) has been demonstrated across this gene cluster in Caucasian populations [2]. A genetic linkage study in European RA families provided suggestive evidence of linkage of RA to the IL-1 region [3]. The reported associations of certain IL-1 gene cluster polymorphisms with susceptibility and severity, especially in the more erosive form of RA, have been variable [4–8]. These studies have generally been conducted on relatively small RA patient samples, apart from the one recent exception [8]. In the latter, the authors performed a comprehensive analysis of IL-1A and IL-1B of a large number of North American and Spanish RA populations but identified no major contribution to RA from these genes.
We examined four common single nucleotide polymorphisms (SNPs) and haplotypes in the IL-1A and IL-1B promoter regions and their association with susceptibility and severity with RA in well-characterized British Caucasian patients. We performed meta-analysis of published data on IL-1B-511 to evaluate the effect of this polymorphism on susceptibility to RA among Caucasians.
| Patients and methods |
|---|
|
|
|---|
Patients and study design
A total of 756 British Caucasian RA patients (542 females, 72%) of Northern European descent attending the Nuffield Orthopaedic Centre in Oxford, UK, were recruited. All satisfied the 1987 ACR criteria for RA and mean age of disease onset was 48 yrs. The controls (n = 625) were healthy British Caucasian blood donors from the same geographical area. All subjects gave informed written consent and approval was granted by the Oxford Research Ethics Committee.
Medical records were examined for the RF status of each patient. Patients who were positive at two different time-points were considered seropositive. Where data were not available, RF was measured using nephelometric technique and a titre over 40 IU/ml was considered positive. Overall, 586 patients (78%) deemed to be RF positive.
Severity of RA was determined by the presence of rheumatoid nodules or the requirement of major joint replacement as established from medical records. The association with radiographic outcome was examined in 206 RA patients with minimum disease duration of 3 yrs. Hand radiographs were scored using Rau–Ratingen method [9].
Genotyping
Genomic DNA was extracted from whole blood. Cases and healthy controls (HCs) were genotyped for IL-1A (–889 C/A, rs17561), IL-1B (–511 A/G, rs16944), IL-1B (–1464 C/G, rs1143623) and IL-1B (–3737 G/A, rs4848306) SNPs by PCR-based methods (details available on request). Limited HLA-DRB1 and further DRB1*04 subtyping was performed as published [10]. All samples were genotyped in duplicates, and positive and negative controls were included. The HLA-DRB1 shared epitope (SE) was defined by the presence of one or two copies of the HLA DRB1*01,*0401, *0404, *0405, *0408 and *10 alleles.
Statistical analysis
Case–control association analysis
Statistical analysis for association and Hardy–Weinberg equilibrium were evaluated using
2-test. P < 0.05 was considered statistically significant. Corrections for multiple testing were applied as necessary [11]. For continuous variables, SPSS version 15.0 for Windows (SPSS Inc., IL, USA) was used. Haplotype estimations and frequency comparison were performed using PHASE v. 2 (available at http://www.stat.washington.edu/stephens/software.html) and CLUMP software (available at http://linkage.rockefeller.edu/soft/clump.html) [12, 13].
Analysis of radiological progression
To allow comparison of radiographic scores in RA patients with different disease duration, we made use of the almost linear relationship between the radiographic score and the square root of disease duration (SQRD) [14]. To ensure that the predicted severity scores on all percentiles remained within the real range of 0–140, a logit transformation of radiographic score {ln [score/(140–score)]} was entered as the dependent variable with SQRD as independent variable. The final regression equation was: logit score = (0.66 x SQRD) – 4.10. Associations between particular IL-1 gene cluster polymorphism and radiographic severity score were assessed by comparing the frequencies of IL-1 SNPs in subjects with mild (score <20) or severe disease (score >80) using Fisher's exact probability test.
Meta-analysis
Search strategy and study selection
Meta-analysis was conducted for IL-1B (-511 A/G) SNPs based on the PubMed database using the terms rheumatoid arthritis and interleukin-1. Publications on Caucasian populations in English language were included. Study selection criteria were: (i) RA patients fulfilled 1987 ACR classification criteria; (ii) controls were drawn from the same geographic area and ethnic background as patients; (iii) authors provided original genotype frequencies; and (iv) RA patient and control groups did not overlap between studies. Six studies met the selection criteria [5, 7, 8, 15–17]; all were original articles with case–control study design. The meta-analysis included 3712 RA patients and 2331 controls.
Evaluation of the statistical association
Though previous studies describe the IL-1B-511 SNP as a G
A change, in the dbSNP database the A allele is ancestral, defining the variant allele as G. However, the G allele occurs more commonly (allele frequencies among HC are 0.67). In the meta-analysis, we compared the variant G allele with ancestral A as well as recessive (G/G vs G/A + AA) and dominant (G/G + G/A vs AA) models for the G allele. Genotype-specific risks such as G/G vs A/A and G/A vs A/A were also analysed. Statistical analyses were performed using StatsDirect software. Single odd ratios (OR) and 95% CIs were estimated for each study. Two-tailed P-values were calculated using the
2-test of independence. Between-study heterogeneity was assessed with Cochran's Q-test statistics. The pooled OR and 95% CI for the variant allele in RA patients against HC were calculated from fixed-effects model, as no significant heterogeneity was observed. Forest plots were generated in StatsDirect to display the OR and 95% CI.
Genetic power calculations for the meta-analysis were done using software program Quanto v1.2.3 (available at http://hydra.usc.edu/gxe) assuming a disease prevalence of 1% and G allele frequency of 0.67. Under log-additive model, we calculated the effect of 0.05 with OR of 1.2 (average OR of all combined studies).
| Results |
|---|
|
|
|---|
Case–control association study
The genotyping results and allele counts for RA patients and HCs are presented in Table 1. IL-1B (–1464 C/G) G allele was less common among RA patients (P = 0.01; OR 1.24; 95% CI 1.04, 1.48). The minor allele frequencies for cases and controls were 24 and 29%, respectively. The association remained significant after applying corrections for false discovery rate.
|
Haplotype analysis revealed a global difference in haplotype frequencies between the RA and HC (P = 0.01). This effect was entirely due to the IL-1B (–1464 C/G) and there was no single over-represented haplotype.
Each of the polymorphisms was investigated in the RA cohort according to RF, SE status and presence of high-risk HLA-DRB1*04 (DRB1*0401, 0404, 0405 or 0408) or a low-risk allele (DRB1*01, 0406, 0407 and *10). No significant difference was observed between the groups (data not shown).
Analysis according to the severity features of RA did not reveal any independent association with the presence of nodules, requirement for joint replacement or radiographic progression (data not shown).
Meta-analysis
Genotyping data in six studies adhered to Hardy–Weinberg equilibrium. No significant heterogeneity was observed between the studies (moment-based estimate of between studies variance = 0, I2 = 0). The meta-analysis of the IL-1B (–511 A/G) showed significant association of the G/G genotype (P = 0.02; OR 1.13; 95% CI 1.02, 1.26) (Fig. 1) and variant allele G (P = 0.02; OR 1.10; 95% CI 1.02, 1.19) with RA. The genotype-specific analysis showed a significant dose effect between the G/G vs A/A (P = 0.04; OR 1.20; 95% CI 1.01, 1.44) but not between G/G vs GA or G/A vs AA genotypes (data not shown).
|
| Discussion |
|---|
|
|
|---|
We have demonstrated a novel association of the IL-1B (–1464 C/G) SNP with RA in a Caucasian population. This association should be replicated in an independent cohort in order to confirm the result. In our data set, IL-1B (–1464 C/G) is in near complete LD with IL-1B (–511 A/G) (D' > 0.90, r2 > 0.80). However, the minor allele frequencies in our cohort for IL-1B (–511 A/G) and IL-1B (–1464 C/G) are dissimilar (0.31 and 0.24 for RA and 0.34 and 0.29 for controls, respectively), indicating that independent genotyping is warranted in order to investigate the IL-1B (–1464 C/G) association. It remains unknown whether the IL-1B (–1464 C/G) SNP is the causal variant. This polymorphism has been shown to exhibit allele-specific difference in gene expression as IL-1B-1464G has shown weaker promoter activity than IL-1B-1464C on electromobility gel shift assay [18]. This promoter region harbours a core-binding motif for GATA-binding family proteins and IL-1B (–1464 C/G) could affect protein binding [18]. It is possible that the real protective variant is elsewhere and exhibits different LD with IL-1B (–511 A/G) and IL-1B (–1464 C/G). Further studies investigating genetic variations across whole gene regulatory regions and more detailed functional analysis are essential.
Although both IL-1
and IL-1β mRNA have been detected in RA synovial fluid, it is IL-1
that has been found to be strongly expressed at the cartilage–pannus junction [19]. In a study of affected sib pair RA families, Cox et al. [4] found four IL-1 gene cluster markers associated with erosive RA, particularly in sib pair families with 0 probability of sharing two HLA-DRB1 alleles identical by descent. Evidence of less severe disease progression after 20 yrs in Swiss RA patients carrying one or two copies of the IL-1B-511 minor allele has also been reported [7]. Using several different analyses, we were unable to demonstrate any significant relationship between IL-1 polymorphisms and radiographic progression. Similar conclusions were reached in the recent study of Johnsen et al. [8] where no association was observed with erosive disease in a large sample of the Brigham Rheumatoid Arthritis Sequential Study cohort.
Meta-analysis allows us to maximize the sample size, increase statistical power and provide more robust findings. The meta-analysis presented here is the largest to have investigated the effects of IL-1B (–511) in RA. We demonstrated a significant association of the GG genotype and G allele with RA. The power of this meta-analysis was over 90% to detect modest effect (OR 1.2 at a significance level of 0.05). However, the magnitude of the genetic risk associated with this polymorphism is relatively small and further studies on a much larger RA sample are desirable to confirm the true association.
In summary, the role of the genetic factors in RA is complex and the magnitude of the effects is likely to vary greatly. We have demonstrated an association of IL-1B (–511) with RA patients of Caucasian origin by conducting a meta-analysis. The potential protective influence of IL-1B (–1464) found in this study will require further confirmation in other populations.
| Acknowledgements |
|---|
|
|
|---|
We would like to thank the many patients who participated in this study and our colleagues Dr Paul Bowness, Dr Joel David and Dr Raashid Luqmani for allowing us access to their patients.
Funding: This study was funded by the University of Oxford.
Disclosure statement: The authors have declared no conflicts of interest.
| References |
|---|
|
|
|---|
- Nicklin MJ, Barton JL, Nguyen M, FitzGerald MG, Duff GW, Kornman K. A sequence-based map of the nine genes of the human interleukin-1 cluster. Genomics (2002) 79:718–25.[CrossRef][Web of Science][Medline]
- Nicklin MJ, Weith A, Duff GW. A physical map of the region encompassing the human interleukin-1 alpha, interleukin-1 beta, and interleukin-1 receptor antagonist genes. Genomics (1994) 19:382–4.[CrossRef][Web of Science][Medline]
- Barrera P, Faure S, Prudhomme JF, et al. European Consortium on Rheumatoid Arthritis Families (ECRAF). European genetic study on rheumatoid arthritis: is there a linkage of the interleukin-1 (IL-1), IL-10 or IL-4 genes to RA? Clin Exp Rheumatol (2001) 19:709–14.[Web of Science][Medline]
- Cox A, Camp NJ, Cannings C, et al. Combined sib-TDT and TDT provide evidence for linkage of the interleukin-1 gene cluster to erosive rheumatoid arthritis. Hum Mol Genet (1999) 8:1707–13.
[Abstract/Free Full Text] - Cantagrel A, Navaux F, Loubet-Lescoulie P, et al. Interleukin-1beta, interleukin-1 receptor antagonist, interleukin-4, and interleukin-10 gene polymorphisms: relationship to occurrence and severity of rheumatoid arthritis. Arthritis Rheum (1999) 42:1093–100.[CrossRef][Web of Science][Medline]
- Jouvenne P, Chaudhary A, Buchs N, Giovine FS, Duff GW, Miossec P. Possible genetic association between interleukin-1alpha gene polymorphism and the severity of chronic polyarthritis. Eur Cytokine Netw (1999) 10:33–6.[Web of Science][Medline]
- Genevay S, di Giovine FS, Perneger TV, et al. Association of interleukin-4 and interleukin-1B gene variants with Larsen score progression in rheumatoid arthritis. Arthritis Rheum (2002) 47:303–9.[CrossRef][Web of Science][Medline]
- Johnsen AK, Plenge MR, Butty V, et al. A broad analysis of IL1 polymorphism and rheumatoid arthritis. Arthritis Rheum (2008) 58:1947–57.[CrossRef][Web of Science][Medline]
- Rau R, Wassenberg S, Herborn G, Stucki G, Gebler A. A new method of scoring radiographic change in rheumatoid arthritis. J Rheumatol (1998) 25:2094–107.[Web of Science][Medline]
- Jaakkola E, Herzberg I, Laiho K, et al. Finnish HLA studies confirm the increased risk conferred by HLA-B27 homozygosity in ankylosing spondylitis. Ann Rheum Dis (2006) 65:775–80.
[Abstract/Free Full Text] - Sabatti C, Service S, Freimer N. False discovery rate in linkage and association genome screens for complex disorders. Genetics (2003) 164:829–33.
[Abstract/Free Full Text] - Stephens M, Smith NJ, Donnelly P. A new statistical method for haplotype reconstruction from population data. Am J Hum Genet (2001) 68:978–89.[CrossRef][Web of Science][Medline]
- Sham PC, Curtis D. Monte Carlo tests for associations between disease and alleles at highly polymorphic loci. Ann Hum Genet (1995) 59:97–105.[Web of Science][Medline]
- de Carvalho A, Graudal H, Jørgensen B. Evaluation of the progression of rheumatoid arthritis. Significance of age at onset and sex. Acta Radiol Diagn (1980) 21:545–50.[Medline]
- Crilly A, Maiden N, Capell HA, Madhok R. Predictive value of interleukin 1 gene polymorphism for surgery. Ann Rheum Dis (2000) 59:695–9.
[Abstract/Free Full Text] - Buchs N, di Giovine FS, Silvestri T, Vannier E, Duff GW, Miossec P. IL-1B and IL-1Ra gene polymorphisms and disease severity in rheumatoid arthritis: interaction with their plasma levels. Genes Immun (2001) 2:222–8.[CrossRef][Web of Science][Medline]
- Tolusso B, Pietrapertosa D, Morelli A, et al. IL-1B and IL-1RN gene polymorphisms in rheumatoid arthritis: relationship with protein plasma levels and response to therapy. Pharmacogenomics (2006) 7:683–95.[CrossRef][Web of Science][Medline]
- Lee KA, Ki CS, Kim HJ, et al. Novel interleukin 1beta polymorphism increased the risk of gastric cancer in a Korean population. J Gastroenterol (2004) 39:429–33.[CrossRef][Web of Science][Medline]
- Deleuran BW, Chu CQ, Field M, et al. Localization of interleukin-1 alpha, type 1 interleukin-1 receptor and interleukin-1 receptor antagonist in the synovial membrane and cartilage/pannus junction in rheumatoid arthritis. Br J Rheumatol (1992) 31:801–9.
[Abstract/Free Full Text]
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
