- Research article
- Open Access
- Open Peer Review
Allelic variants of IL1R1gene associate with severe hand osteoarthritis
BMC Medical Genetics volume 11, Article number: 50 (2010)
In search for genes predisposing to osteoarthritis (OA), several genome wide scans have provided evidence for linkage on 2q. In this study we targeted a 470 kb region on 2q11.2 presenting the locus with most evidence for linkage to severe OA of distal interphalangeal joints (DIP) in our genome wide scan families.
We genotyped 32 single nucleotide polymorphisms (SNPs) in this 470 kb region comprising six genes belonging to the interleukin 1 superfamily and monitored for association with individual SNPs and SNP haplotypes among severe familial hand OA cases (material extended from our previous linkage study; n = 134), unrelated end-stage bilateral primary knee OA cases (n = 113), and population based controls (n = 436).
Four SNPs in the IL1R1 gene, mapping to a 125 kb LD block, provided evidence for association with hand OA in family-based and case-control analysis, the strongest association being with SNP rs2287047 (p-value = 0.0009).
This study demonstrates an association between severe hand OA and IL1R1 gene. This gene represents a highly relevant biological candidate since it encodes protein that is a known modulator of inflammatory processes associated with joint destruction and resides within a locus providing consistent evidence for linkage to hand OA. As the observed association did not fully explain the linkage obtained in the previous study, it is plausible that also other variants in this genome region predispose to hand OA.
Osteoarthritis (OA) is a complex late onset disease characterized by the destruction of joints and it can occur in one or multiple joint groups as hand and knee OA often do [1–3]. The most frequent clinical findings in OA are pain associated with loading and limited range of motion in the joints at later stages. Typical radiographic findings, narrowing of the intra-articular space and osteophytes, are considered superior to clinical manifestations in the diagnostic definition . There is a strong genetic effect on hand and knee OA, with heritability ranging between 39%-65% independent of known environmental or demographic confounders . Based on a Finnish study , the genetic effects seem to be more prominent in females.
Our genome wide linkage analysis in OA families identified chromosome 2q11-q21 locus as the major locus for severe hand OA in Finnish families . The same locus has been identified in screens of rare multiplex pedigrees and sibpairs from more heterogeneous populations [7, 8], and in a large population-based twin cohort . The 2q11-q21 locus harbours the interleukin 1 (IL-1) gene cluster. The marker providing strongest evidence for linkage in Finnish distal interphalangeal joint (DIP) OA families resides within the IL1R1 gene (LOD score 2.34), and some evidence for a shared haplotype among the affected individuals was observed with markers D2S2264, IL1R1, D2S373, and D2S1789 (p-value of 0.012) . Further, a handful of association studies have shown a potential role for interleukin-1 gene family in knee and hip OA [10–12].
In the human genome, IL1R2, IL1R1, IL1RL2, IL1RL1, IL18R1 and IL18RAP are located in a 470 kb cluster on 2q11. Genes coding for most other IL-1 ligands are also clustered on chromosome 2 . IL-1 is a pleiotropic cytokine involved in the initiation of inflammatory and immune responses. Using gene names also for their protein products, a complex of IL-1 (IL1A or IL1B), IL1RI and IL1RAP are needed for signal transduction. There are other interleukin ligand-receptor complexes leading to NF-κB activation too [14–16]. Therapeutic effect of IL-1 receptor antagonist (IL1RN) gene transfer has been confirmed in different experimental models of OA .
Our aim was to study more closely the 2q11-q21 locus, major locus for severe hand OA in our previous study, using SNPs and SNP haplotypes in extended hand OA families (including siblings with a severe, radiologically defined bilateral DIP OA used in our genome scan ), and in a set of severe primary bilateral knee OA patients. We targeted a 470 kb region in chromosome 2q11.2 by genotyping 32 SNP markers in total of 220 DNA samples from individuals affected with OA and in 436 population controls.
Our hand OA material was based on the set of radiologically verified severe DIP OA families from our previous study described in more detail in Leppävuori et al.  and summarised in Table 1. DNA samples of 86 severe radiologically verified bilateral DIP OA cases were available from our previous study . A radiologist had graded the osteoarthritic changes in the hand x-rays for DIP, PIP, MCP joints (according to ). Data for reported height and weight were available. Subjects with rheumatoid arthritis (RA) were excluded. Inclusion criteria were 3rd or 4th degree Kellgren and Lawrence (K/L) radiographic OA in DIP joints (bilateral) and a positive family history of OA, with at least one sibling with DIP OA. The other siblings were included as affected individuals if they had a K/L score ≥ 2 in DIP joints. Additionally, 21 unaffected family members derived from 40 families were available for this study. Study set was further extended by 17 radiologically verified severe DIP OA cases collected from the ORTON Orthopaedic Hospital (Helsinki, Finland) with the same radiographic criteria as used in selecting the index cases in initial material (Table 1) and also 4 affected siblings. The average age at grading of hand X-rays was 62 years (SD 9 y). For more about diagnostics and epidemiology of hand OA, see Haara et al..
To collect knee OA material, we reviewed case reports of 852 patients visiting ORTON Orthopaedic Hospital between 1994-2001 having knee OA and identified 220 patients with primary bilateral knee OA severe enough to fulfill the criteria for knee arthroplasty: pain, walking disability and radiologically at least stage 3/4 osteoarthritic changes according to Kellgren and Lawrence classification . Of the 220 patients, we reached 113 subjects, who had not had a major knee trauma and volunteered to participate by providing a DNA sample (Table 1). Knee OA patients' pain or other OA symptoms began at a mean age of 52 y (SD 12 y), and mean age at first arthroplasty was 67 y (SD 8). 28 patients had physician-diagnosed hand OA and were also included in the hand OA material to increase statistical power.
The control subjects (n = 436, aged 35-70) were selected from the Finnish twin cohort . Approximately half of this control material (n = 210) was selected from the Finnish twin cohort  with inclusion criteria of born in 1938-1941, responded to a questionnaire in 1996-1997, did not have diagnosis of OA and gave DNA samples for analyses. One twin from each twin pair was included in the control group if neither twin had physician diagnosed OA or RA, and neither twin reported that their mother, father, co-twin or any other sibling had OA or RA. The remaining half (n = 226) of the controls were subjects of Twin Spine study which is part of the Finnish twin cohort study. The Twin Spine study consists of 300 male twin pairs selected on the basis of putative risk factors for low back pain and disc degeneration (not particularly suffering from either) as described elsewhere [22, 23]. Subjects with a physician-diagnosed OA were excluded, and only one twin from a pair was included. The study was approved by the ethics committee of the Helsinki metropolitan hospital region and all individuals gave their informed consent.
Methods and statistical analyses
A set of 49 SNPs with a minimum heterozygosity of 0.2 in the Caucasian population were selected from SeattleSNP  and Snpper CHIP Bioinformatics Tools  databases to cover the 470 kb target region with approximately 10-kb SNP spacing (Fig. 1). These SNPs were genotyped using the MassARRAY® system (Sequenom, San Diego, CA, USA) as recommended by the manufacturer with additional quality assessment steps described in more detail in Silander et al. . The genotyping was done in multiplexes of 1-5 SNPs using hME assays in 384-well plates. Each SNP multiplex was validated by genotyping 81 trio samples (i.e. 27 families with two parents and a child) to control for Mendelian incompatibilities using Pedcheck 1.1 program . Each 384-well plate contained 8 duplicate samples and 8 non-template wells for quality assessment. If the genotypes could not be reliably called, if there were unexplainable errors in Mendelian inheritance or between duplicate samples, or genotypes in the non-template wells, SNP was excluded from the study. Further, only SNPs with call rate of over 90% were included. A total of 32 SNP passed the preset quality criteria and were included in the further analysis. All included SNPs showed HW p-values of more than 0.1.
Additionally using a microarray based allele specific primer extension method , we genotyped eight SNPs in more studied area in chromosome 2q harbouring the TRAP1, IL1A, IL1B, IL1F5 and IL1RN genes (rs1030878, rs2071375, rs1800587 (SNP genotyped approximately in half of the controls, n = 169), rs1143634, rs1530551, rs2234677, rs315952, and rs9005) three of which were previously analysed by Smith and colleagues . In this additional analysis we lacked 67 individuals (one hand OA case, four unaffected hand OA family members, 21 knee OA cases and 41 controls).
Deviation from the Hardy-Weinberg equilibrium was tested using Pearson's χ2 test in the control sample. For hand OA families, Pedcheck 1.1  was used to rule out Mendelian errors. The degree of pairwise linkage disequilibrium (LD) was calculated for each pair of SNPs in control samples using the Haploview program . Both r2 and D' measures of LD were studied. Then, evidence for association was evaluated in hand OA and knee OA subgroups using the Pseudomarker 0.9.7 beta program which can utilize the information from family and case-control settings . We monitored for LD with the disease assuming linkage. Results were confirmed by traditional χ2 analysis (including only an index case of each hand OA family) or by Fisher's test when less than five individuals were carrying the same genotype. The Pseudomarker analysis was performed additionally for females only in both disease groups. The number of males was too small in both case groups for separate analysis.
To evaluate the power in negative findings for knee OA group, we calculated the amount of cases required to have 90% power to detect similar association as with the hand OA group and the power with 113 cases using the Genetic Power Calculator . The SNPSpD method  with modifications by Li and Ji  was used to calculate the p-value threshold for 5% significance .
Next using the Haploview program , we performed haplotype association test in the area that showed the strongest evidence for association with individual SNPs located in a single LD block measured by r2, covering the IL1R1 gene (rs1465325, rs956730, rs3917225, rs2287047, rs3771200) (index case per family in hand OA group). We further estimated haplotypes for each subject individually using the Phase2 program [36, 37] for calculating odds ratios (OR) for associated haplotypes (index case per family in hand OA group). A minimum probability of ≥ 0.60 was requested for the haplotypes, and alleles with a frequency of over 0.05 in the control sample were each analyzed separately and the less common haplotype alleles were grouped together. Hardy-Weinberg equilibrium was calculated for tested haplotypes to confirm the correctness of the haplotype allele distribution. For further information on the haplotype association in the hand OA group, we conducted haplotype association test in the associated area for SNPs in LD measured by D' solid spine of LD (rs1465325, rs956730, rs3917225, rs2287047, rs3771200, rs2241132, rs870684, rs1922290 and rs1922295) using the Haploview program .
Finally, to study the linkage between hand OA and SNPs showing the strongest association with hand OA, the Pseudomarker program was used. Then, to assess if affected individuals with a certain genotype contribute more than expected to original linkage in the region, we used the genotype-IBD Sharing Test (GIST) software package [38, 39]. Analysis was performed with the families used in the original genome-wide linkage scan and also included in this study (23 families). We used genotype data for SNPs showing the most significant association and single point NPL scores for the original "IL1R1" microsatellite marker for each family obtained from the Merlin program [40, 41].
A set of 32 SNPs mapping within a 470 kb region on 2q11.2 presenting the locus with most evidence for linkage to severe hand OA in our previous genome wide scan families (Fig. 1) was genotyped in severe hand OA cases (n = 134), end-stage bilateral primary knee OA cases (n = 113) and controls (n = 436) (Table 1). According to SNPSpD, the effective number of independent SNP loci was 23.9 thus the experiment-wide significance threshold 0.0021 was required to keep the type I error rate at 5%. The p-values of the association analysis are presented in Table 2 as non-corrected. Out of 32 SNPs, four SNPs in IL1R1 gene in 125 kb area showed some evidence for association to hand OA (p-values < 0.05), in a family-based association analysis using the Pseudomarker program and/or in the case-control analysis (Table 2). Association with severe hand OA was concentrated in one LD block (r2 > 0.34, D' > 0.94) harboring the area of long promoter and the coding region of gene IL1R1 (Fig. 2). It is notable, that association of the SNP rs2287047 to severe hand OA (odds ratio (OR) for genotypes AA vs. GG = 0.16, 95% confidence interval (CI) = 0.06 to 0.45; OR for genotypes AA vs. AG = 0.26, 95% CI = 0.09 to 0.73; protective A allele being the minor allele; p = 0.0009 dominant mode of inheritance) remained statistically significant after correction for multiple testing. This SNP is located in intron 1 of the IL1R1 gene. Association analysis using only female hand OA subjects provided comparable results (best p-value 0.0045 for SNP rs2287047, dominant model).
One SNP (rs956730) in the same LD block showed nominal association in knee OA material (p = 0.045, minor allele frequency = 0.324) and further analysis using female knee OA subjects only provided some evidence for association with SNPs in the same LD block (r2 = 0.17, D' = 0.90) (Fig. 3) mapping within the IL1RL2 gene (p-values 0.0098, 0.012 and 0.012 for SNPs rs870684, rs1922290 and rs1922295, respectively; dominant model). However, none of the associations were statistically significant after correcting for multiple testing.
At least 503 knee OA cases would have been required to have 90% power to detect association at the level of p < 0.001 for SNP rs2287047 when assuming a similar effect size as observed in the family-based severe hand OA (OR 22 vs. 11 = 0.16), whereas the power to detect association with the 113 knee OA cases was 12%. Most probably the effect size of the studied SNP (rs2287047) observed here is an over-estimate and if assuming an effect size typically seen in complex diseases (expected risk ratios (RRs) for a single SNP closer to 1.2 and 1.15, for AA and Aa, respectively), approximately 4000 cases would have been needed for detecting the association with 80% power (allele frequency of 0.75 and prevalence of 0.1) .
We were not able to detect association (Pseudomarker and Pearson's χ2 tests) with any of the SNPs in the previously more intensively studied area (TRAP1, IL1A, IL1B, IL1F5, IL1RN) (p > 0.14). Also, opposite to findings by Smith et al., we did not find LD between SNPs in IL1R1 and the locus containing IL1A, IL1B, and IL1RN in our study material (r2 ≤ 0.03, D'<0.54).
To cover most of the allelic variation on the genome region showing evidence for association to severe hand OA, we monitored for haplotype association for IL1R1 (rs1465325, rs956730, rs3917225, rs2287047, rs3771200) in hand and knee OA groups using the Haploview program  (Table 3). The frequencies of allele 3 (CAAAG) differed between hand OA cases (0.104) and controls (0.188) (p-value 0.008). SNP rs1465325 (OR CT vs. TT = 0.59, CI = 0.45 to 0.78; p = 0.022) tagged the protective allele. The protective A allele of SNP rs2287047 for hand OA was observed in two allelic backgrounds, on allelic haplotypes 3 and 5. Also, frequencies of the rare haplotype 5 (TAAAG) differed between knee OA cases (0.103) and controls (0.053) (p-value 0.004), but the haplotype 5 was the predisposing haplotype.
Next, we estimated individual haplotypes with the Phase2 program [36, 37], to calculate an odds ratio for the associated allelic haplotypes on basis of one allele versus all other alleles. We were able to estimate haplotypes fulfilling the preset criteria (probability ≥ 0.60) for 96.5% of the hand OA samples, 95.6% of the knee OA samples and 100% of the controls, and the frequencies of the estimated haplotypes were comparable to ones provided by Haploview. When comparing the hand OA material and controls, the associated protective allele 3 (CAAAG) provided an odds ratio of 0.54 (OR 3X vs. XX, 95% CI = 0.41 to 0.72) and the rare allele 5 (TAAAG) an odds ratio of 1.81 in the knee OA (OR 5X vs. XX, 95% CI = 1.34 to 2.43) where the X represents all other alleles combined. We were only able to calculate an odds ratio for one associated allele versus no associated alleles (3X vs. XX and 5X vs. XX) since there were no individuals homozygote for the associated allele either in the case groups or in the control group. The effect of allele 5 to knee OA was opposite to the effect for hand OA, but it is important to note that the allele 5 was rare.
We also performed a haplotype association test by selecting SNPs in D' LD with each other (measured by D'>0.94 for the first and the last SNP in the block; SNPs rs1465325, rs956730, rs3917225, rs2287047, rs3771200, rs2241132, rs870684, rs1922290 and rs1922295) in hand OA group and in controls. Allele CAAAGCTGC corresponding to allele 3 (CAAAG) in the previous test, showed similar allele frequencies, allele CAAAGCTGC being more rare in the hand OA cases (0.101) than in the controls (0.190) (p = 0.0058).
The four SNPs showing some association with hand OA provided also nominal evidence for linkage in the families (LOD = 0.640670, p = 0.042930 for rs1465325; LOD = 1.343799, p = 0.006440 for rs956730; LOD = 0.281669, p = 0.127359 for rs3917225; LOD = 0.841737, p = 0.024494 for rs2287047). Although all except one affected individual carried the predisposing common G-allele of rs2287047, GIST analysis resulted in non-significant p-values when using combined model of dominant, recessive and additive models (p = 0.672 for rs1465325, p = 0.499 for rs956730, p = 0.563 for rs3917225, p = 0.998 for rs2287047), suggesting that the SNPs in IL1R1 do not fully explain the originally observed linkage.
This study identifies an association between severe hand OA and the variants of the IL1R1 gene. The strongest evidence for association, which remained statistically significant after correction for multiple testing, was observed with SNP (rs2287047) located in the intron 1 of the IL1R1 gene. This gene represents a highly relevant biological candidate gene since it encodes protein that is a known modulator of inflammatory processes associated with joint destruction and resides within a locus providing consistent evidence for linkage to hand OA. IL1R1 did not show significant single SNP association with knee OA. It is possible that IL1R1 does not have strong or any affect in knee OA, or the lack of association for knee OA is due to lack of power in the analysis. Also the knee OA patients did not have a similar strong family background for the disease as the hand OA patients. However, nominal association in the same D' LD block in IL1RL2 gene was observed for knee OA in females. Further, the association between radiological bilateral DIP OA (K/L at least 2) and IL1R1 and IL1RL2 genes was recently studied also by Solovieva et al.  in study material of 206 cases and 328 controls. None of the single SNPs (rs1465325, rs956730, rs2287047) showed association to bilateral DIP OA, but some occupation dependent association was observed with a haplotype covering the IL1R1 and IL1RL2 genes. However, it is worth of noting that the phenotype utilized in the study by Solovieva and collaborators is not equivalent to the phenotype used here (severe bilateral hand OA) and that approximately 4000 cases would be required to exclude SNPs with the effect size typically observed in complex diseases.
There are some limitations in the present study. Based on the GIST analysis, the observed associations for single SNPs in this study does not fully explain the linkage obtained in the previous study using microsatellite markers. It is likely that the SNPs associated with hand OA in this study are not the causative variants but are probably in LD with the true causative variant(s) located in the IL1R1 gene region. Furthermore, it is plausible that allelic variants in several genes within this locus contribute to OA susceptibility, each with a small effect. Other studies seem to support this hypothesis: The interleukin 1 gene family has recently been a widely studied target in osteoarthritis [8, 10–12, 29, 43–46], partly due to several genome wide scans demonstrating linkage to this region [6–8], but also because of the potential role of the gene products in inflammatory responses associated with the development of OA [13, 15, 16, 47, 48]. Mostly, the previous association analyses have concentrated on the genes IL1A, IL1B and IL1RN and have left the 470 kb region targeted in this study, with less attention. However, we did not see association in the previously more intensively studied area harbouring the IL1A, IL1B and IL1RN genes using the Finnish hand and knee OA materials.
The limited number of subjects with severe hand OA (Table 1) is another limitation in this study. However, it is important to emphasize the unique nature of this study sample as its strength: the hand OA study sample is based on the same material, which provided evidence for linkage to this genome region on 2q in our genome wide scan. Additionally, an association analysis with the Pseudomarker program, which can utilize information from several family members in addition to case-control data, provided the strongest evidence for association. Using the family based analysis putatively gives us the advantage to identify population specific variants that are important in understanding the aetiology of the disease. Further, most subjects in the hand OA material have a severe end stage disease and very strict criteria was used in selecting the study sample. Sequencing of the associated region in these familial hand OA cases may identify rare, potentially causative variant(s) of the gene segregating with the common allelic variant in these families. Thus, further validation in a large sample set and additional biological evidence are required to confirm the role of IL1R1 in severe hand OA.
A complex formed by IL-1, IL1R1 and IL1RAP initiates numerous signalling pathways (see review ). IL1R1 has three alternative promoters generating alternative transcripts . The promoter region of IL1R1 has proven to be highly polymorphic, possibly due to the polymorphisms affecting the degree of expression in a wide variety of tissues . IL1R1 is not abundant in cells but less than 10 ligand-occupied receptor molecules on the cell's surface is enough to induce a strong response due to parallel signalling pathways . One of these pathways leads to activation of NF-κB, a regulator of inflammatory and immune gene expression . Lately, the contribution of NF-κB signalling components to the pathogenesis of different rheumatic diseases and the pharmacologic modulation of NF-κB has been widely studied. These studies have yielded promising results such as currently available therapeutic agents, novel small molecule inhibitors and improved antisense DNA therapy and RNA interference . It has been shown, that blocking the signal transduction of IL1R1 prevents OA progression in animal models . This would support the biological significance of our finding of IL1R1 as a potential factor in OA development. All four associated SNPs were located in the non-coding region of IL1R1, but at least three of them (rs1465325, rs3917225 and rs2287047) co-localized with predicted transcription factor binding sites so that the conserved allele matched with the consensus sequence of the transcription factor binding site while the other allele did not or changed it to other transcription factor binding site (Computational Biology Research Center CBRC, http://www.cbrc.jp/index.eng.html). Polymorphisms in IL1R1 gene has been associated with other diseases as well, for instance AIDS progression , endometriosis  and Helicobacter pylori infection .
To summarize, a set of 32 SNPs mapping within a 470 kb region surrounding the strongest linkage peak of a genome wide scan performed in Finnish severe hand OA families was genotyped in 134 hand and 113 knee OA cases and 436 controls. The family based association analysis demonstrated association to four SNPs mapping to a 125 kb DNA region comprising the IL1R1 gene. The strongest evidence for association, observed with a SNP rs2287047 located in the intron 1 of the IL1RL1 gene (p = 0.0009) remained significant after correction for multiple testing. The haplotype analysis provided further evidence for the IL1R1 gene: One haplotype allele showed association with hand OA. The IL1R1 gene encodes receptor whose activation by the binding of a specific ligand leads to activation of NF-κB, which is a modulator of inflammatory and immune gene expression. However, the variants linked and associated with hand OA in this study did not fully explain the original linkage observed in our previous study suggesting that there may also be other hand OA predisposing variants within this locus. Potentially chromosome region 2q harbors several variants in different genes affecting the disease.
Kellgren JH, Moore R: Generalized osteoarthritis and Heberden's nodes. BMJ. 1952, 1 (4751): 181-187. 10.1136/bmj.1.4751.181.
Kujala UM, Leppävuori J, Kaprio J, Kinnunen J, Peltonen L, Koskenvuo M: Joint-specific twin and familial aggregation of recalled physician diagnosed osteoarthritis. Twin Res. 1999, 2 (3): 196-202. 10.1375/136905299320565861.
Hirsch R, Lethbridge-Cejku M, Scott WW, Reichle R, Plato CC, Tobin J, Hochberg MC: Association of hand and knee osteoarthritis: evidence for a polyarticular disease subset. Ann Rheum Dis. 1996, 55 (1): 25-29. 10.1136/ard.55.1.25.
Hart D, Spector T, Egger P, Coggon D, Cooper C: Defining osteoarthritis of the hand for epidemiological studies: the Chingford Study. Ann Rheum Dis. 1994, 53 (4): 220-223. 10.1136/ard.53.4.220.
Spector TD, Cicuttini F, Baker J, Loughlin J, Hart D: Genetic influences on osteoarthritis in women: a twin study. Bmj. 1996, 312 (7036): 940-943.
Leppävuori J, Kujala U, Kinnunen J, Kaprio J, Nissilä M, Heliövaara M, Klinger N, Partanen J, Terwilliger JD, Peltonen L: Genome scan for predisposing loci for distal interphalangeal joint osteoarthritis: evidence for a locus on 2q. Am J Hum Genet. 1999, 65 (4): 1060-1067. 10.1086/302569.
Wright GD, Hughes AE, Regan M, Doherty M: Association of two loci on chromosome 2q with nodal osteoarthritis. Ann Rheum Dis. 1996, 55 (5): 317-319. 10.1136/ard.55.5.317.
Loughlin J, Mustafa Z, Smith A, Irven C, Carr AJ, Clipsham K, Chitnavis J, Bloomfield VA, McCartney M, Cox O, et al: Linkage analysis of chromosome 2q in osteoarthritis. Rheumatology (Oxford). 2000, 39 (4): 377-381. 10.1093/rheumatology/39.4.377.
Livshits G, Kato BS, Zhai G, Hart DJ, Hunter D, MacGregor AJ, Williams FM, Spector TD: Genomewide linkage scan of hand osteoarthritis in female twin pairs showing replication of quantitative trait loci on chromosomes 2 and 19. Ann Rheum Dis. 2007, 66 (5): 623-627. 10.1136/ard.2006.060236.
Moos V, Rudwaleit M, Herzog V, Hohlig K, Sieper J, Muller B: Association of genotypes affecting the expression of interleukin-1beta or interleukin-1 receptor antagonist with osteoarthritis. Arthritis Rheum. 2000, 43 (11): 2417-2422. 10.1002/1529-0131(200011)43:11<2417::AID-ANR7>3.0.CO;2-R.
Loughlin J, Dowling B, Mustafa Z, Chapman K: Association of the interleukin-1 gene cluster on chromosome 2q13 with knee osteoarthritis. Arthritis Rheum. 2002, 46 (6): 1519-1527. 10.1002/art.10260.
Meulenbelt I, Seymour AB, Nieuwland M, Huizinga TW, van Duijn CM, Slagboom PE: Association of the interleukin-1 gene cluster with radiographic signs of osteoarthritis of the hip. Arthritis Rheum. 2004, 50 (4): 1179-1186. 10.1002/art.20121.
Subramaniam S, Stansberg C, Cunningham C: The interleukin 1 receptor family. Dev Comp Immunol. 2004, 28 (5): 415-428. 10.1016/j.dci.2003.09.016.
Debets R, Timans JC, Homey B, Zurawski S, Sana TR, Lo S, Wagner J, Edwards G, Clifford T, Menon S, et al: Two novel IL-1 family members, IL-1 delta and IL-1 epsilon, function as an antagonist and agonist of NF-kappa B activation through the orphan IL-1 receptor-related protein 2. J Immunol. 2001, 167 (3): 1440-1446.
O'Neill LA, Greene C: Signal transduction pathways activated by the IL-1 receptor family: ancient signaling machinery in mammals, insects, and plants. J Leukoc Biol. 1998, 63 (6): 650-657.
Towne JE, Garka KE, Renshaw BR, Virca GD, Sims JE: Interleukin (IL)-1F6, IL-1F8, and IL-1F9 signal through IL-1Rrp2 and IL-1RAcP to activate the pathway leading to NF-kappaB and MAPKs. J Biol Chem. 2004, 279 (14): 13677-13688. 10.1074/jbc.M400117200.
Evans CH, Gouze JN, Gouze E, Robbins PD, Ghivizzani SC: Osteoarthritis gene therapy. Gene Ther. 2004, 11 (4): 379-389. 10.1038/sj.gt.3302196.
Kellgren JH: Atlas of standard radiographs. The epidemiology of chronic rheumatism. Edited by: Jeffrey MR, Ball J. 1963, Oxford: Blackwell Scientific publications, 1-9.
Haara MM, Manninen P, Kroger H, Arokoski JP, Karkkainen A, Knekt P, Aromaa A, Heliovaara M: Osteoarthritis of finger joints in Finns aged 30 or over: prevalence, determinants, and association with mortality. Ann Rheum Dis. 2003, 62 (2): 151-158. 10.1136/ard.62.2.151.
Kellgren JH, Lawrence JS: Radiological assessment of osteo-arthrosis. Ann Rheum Dis. 1957, 16 (4): 494-502. 10.1136/ard.16.4.494.
Kaprio J, Koskenvuo M: Genetic and environmental factors in complex diseases: the older Finnish Twin Cohort. Twin Res. 2002, 5 (5): 358-365. 10.1375/136905202320906093.
Battié MC, Videman T, Levälahti E, Gill K, Kaprio J: Heritability of low back pain and the role of disc degeneration. Pain. 2007, 131 (3): 272-280. 10.1016/j.pain.2007.01.010.
Videman T, Leppavuori J, Kaprio J, Battie MC, Gibbons LE, Peltonen L, Koskenvuo M: Intragenic polymorphisms of the vitamin D receptor gene associated with intervertebral disc degeneration. Spine. 1998, 23 (23): 2477-2485. 10.1097/00007632-199812010-00002.
Snpper CHIP Bioinformatics Tools. [http://snpper.chip.org]
Silander K, Axelsson T, Widen E, Dahlgren A, Palotie A, Syvanen AC: Analysis of genetic variation in the GenomEUtwin project. Twin Res. 2003, 6 (5): 391-398. 10.1375/136905203770326394.
O'Connell JR, Weeks DE: PedCheck: a program for identification of genotype incompatibilities in linkage analysis. Am J Hum Genet. 1998, 63 (1): 259-266. 10.1086/301904.
Silander K, Komulainen K, Ellonen P, Jussila M, Alanne M, Levander M, Tainola P, Kuulasmaa K, Salomaa V, Perola M, et al: Evaluating whole genome amplification via multiply-primed rolling circle amplification for SNP genotyping of samples with low DNA yield. Twin Res Hum Genet. 2005, 8 (4): 368-375. 10.1375/twin.8.4.368.
Smith AJ, Keen LJ, Billingham MJ, Perry MJ, Elson CJ, Kirwan JR, Sims JE, Doherty M, Spector TD, Bidwell JL: Extended haplotypes and linkage disequilibrium in the IL1R1-IL1A-IL1B-IL1RN gene cluster: association with knee osteoarthritis. Genes Immun. 2004, 5 (6): 451-460. 10.1038/sj.gene.6364107.
Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21 (2): 263-265. 10.1093/bioinformatics/bth457.
Göring HH, Terwilliger JD: Linkage analysis in the presence of errors IV: joint pseudomarker analysis of linkage and/or linkage disequilibrium on a mixture of pedigrees and singletons when the mode of inheritance cannot be accurately specified. Am J Hum Genet. 2000, 66 (4): 1310-1327. 10.1086/302845.
Purcell S, Cherny SS, Sham PC: Genetic Power Calculator: design of linkage and association genetic mapping studies of complex traits. Bioinformatics. 2003, 19 (1): 149-150. 10.1093/bioinformatics/19.1.149.
Nyholt DR: A simple correction for multiple testing for single-nucleotide polymorphisms in linkage disequilibrium with each other. Am J Hum Genet. 2004, 74 (4): 765-769. 10.1086/383251.
Li J, Ji L: Adjusting multiple testing in multilocus analyses using the eigenvalues of a correlation matrix. Heredity. 2005, 95 (3): 221-227. 10.1038/sj.hdy.6800717.
Single Nucleotide Polymorphism Spectral Decomposition (SNPSpD). [http://gump.qimr.edu.au/general/daleN/SNPSpD/]
Stephens M, Donnelly P: A comparison of bayesian methods for haplotype reconstruction from population genotype data. Am J Hum Genet. 2003, 73 (5): 1162-1169. 10.1086/379378.
Stephens M, Smith NJ, Donnelly P: A new statistical method for haplotype reconstruction from population data. Am J Hum Genet. 2001, 68 (4): 978-989. 10.1086/319501.
Li C, Scott LJ, Boehnke M: Assessing whether an allele can account in part for a linkage signal: the Genotype-IBD Sharing Test (GIST). Am J Hum Genet. 2004, 74 (3): 418-431. 10.1086/381712.
Abecasis GR, Cherny SS, Cookson WO, Cardon LR: Merlin--rapid analysis of dense genetic maps using sparse gene flow trees. Nat Genet. 2002, 30 (1): 97-101. 10.1038/ng786.
Solovieva S, Kämäräinen O-P, Hirvonen A, Hämäläinen S, Laitala M, Vehmas T, Luoma K, Näkki A, Riihimäki H, Ala-Kokko L, et al: Association between interleukin 1 gene cluster polymorphisms and bilateral dip osteoarthritis. The Journal of Rheumatology. 2009, 36 (9): 1977-1986. 10.3899/jrheum.081238.
Kanoh T, Hasegawa Y, Masui T, Yamaguchi J, Ishiguro N, Hamajima N: Interleukin-1beta gene polymorphism associated with radiographic signs of osteoarthritis of the knee. J Orthop Sci. 2008, 13 (2): 97-100. 10.1007/s00776-007-1205-3.
Moxley G, Han J, Stern AG, Riley BP: Potential influence of IL1B haplotype and IL1A-IL1B-IL1RN extended haplotype on hand osteoarthritis risk. Osteoarthritis Cartilage. 2007, 15 (10): 1106-1112. 10.1016/j.joca.2007.03.022.
Ni H, Shi D, Dai J, Qin J, Xu Y, Zhu L, Yao C, Shao Z, Chen D, Xu Z, et al: Genetic polymorphisms of interleukin-1beta (-511C/T) and interleukin-1 receptor antagonist (86-bpVNTR) in susceptibility to knee osteoarthritis in a Chinese Han population. Rheumatol Int. 2009, 29 (11): 1301-1305. 10.1007/s00296-009-0848-0.
Sezgin M, Erdal ME, Altintas ZM, Ankarali HC, Barlas IO, Turkmen E, Sahin G: Lack of association polymorphisms of the IL1RN, IL1A, and IL1B genes with knee osteoarthritis in Turkish patients. Clin Invest Med. 2007, 30 (2): E86-92.
Bacconnier L, Jorgensen C, Fabre S: Erosive osteoarthritis of the hand: clinical experience with anakinra. Ann Rheum Dis. 2009, 68 (6): 1078-1079. 10.1136/ard.2008.094284.
Riyazi N, Slagboom E, de Craen AJ, Meulenbelt I, Houwing-Duistermaat JJ, Kroon HM, van Schaardenburg D, Rosendaal FR, Breedveld FC, Huizinga TW, et al: Association of the risk of osteoarthritis with high innate production of interleukin-1beta and low innate production of interleukin-10 ex vivo, upon lipopolysaccharide stimulation. Arthritis Rheum. 2005, 52 (5): 1443-1450. 10.1002/art.21014.
Ye K, Vannier E, Clark BD, Sims JE, Dinarello CA: Three distinct promoters direct transcription of different 5' untranslated regions of the human interleukin 1 type I receptor: a possible mechanism for control of translation. Cytokine. 1996, 8 (6): 421-429. 10.1006/cyto.1996.0058.
Dinarello CA: Biologic basis for interleukin-1 in disease. Blood. 1996, 87 (6): 2095-2147.
Roman-Blas JA, Jimenez SA: NF-kappaB as a potential therapeutic target in osteoarthritis and rheumatoid arthritis. Osteoarthritis Cartilage. 2006, 14 (9): 839-848. 10.1016/j.joca.2006.04.008.
Fernandes J, Tardif G, Martel-Pelletier J, Lascau-Coman V, Dupuis M, Moldovan F, Sheppard M, Krishnan BR, Pelletier JP: In vivo transfer of interleukin-1 receptor antagonist gene in osteoarthritic rabbit knee joints: prevention of osteoarthritis progression. Am J Pathol. 1999, 154 (4): 1159-1169.
Do H, Vasilescu A, Carpentier W, Meyer L, Diop G, Hirtzig T, Coulonges C, Labib T, Spadoni JL, Therwath A, et al: Exhaustive genotyping of the interleukin-1 family genes and associations with AIDS progression in a French cohort. J Infect Dis. 2006, 194 (11): 1492-1504. 10.1086/508545.
D'Amora P, Sato H, Girao MJ, Silva ID, Schor E: Polymorphisms in exons 1B and 1C of the type I interleukin-1 receptor gene in patients with endometriosis. Am J Reprod Immunol. 2006, 56 (3): 178-184. 10.1111/j.1600-0897.2006.00415.x.
Hartland S, Newton JL, Griffin SM, Donaldson PT: A functional polymorphism in the interleukin-1 receptor-1 gene is associated with increased risk of Helicobacter pylori infection but not with gastric cancer. Dig Dis Sci. 2004, 49 (9): 1545-1550. 10.1023/B:DDAS.0000042262.14969.2d.
The pre-publication history for this paper can be accessed here:http://0-www.biomedcentral.com.brum.beds.ac.uk/1471-2350/11/50/prepub
We wish to thank the participating OA patients and families and Drs Markku Heliövaara and Martti Nissilä for their contribution in the identification of hand OA cases. Arja Tapio is acknowledged for the expert technical assistance. The Academy of Finland, Finnish Ministry of Education, ORTON Research Institute, Invalid Foundation, TULES Graduate School, and TBGS National Graduate School of Musculoskeletal Disorders and Biomaterials have supported the study financially. Leena Peltonen and Jaakko Kaprio were supported by the Academy of Finland Centre of Excellence in Complex Disease Genetics.
The authors declare that they have no competing interests.
JS, JK, LP and UMK designed and coordinated the study and participated in writing the manuscript. SK carried out the genotyping. AN and SK participated in designing the study, carried out all the statistical analyses and drafted the manuscript. AH, KT and UMK are responsible for the phenotyping and collection of the samples, and they participated in writing the manuscript. TV, MCB and JK are responsible for the control samples and participated in writing the manuscript. All authors read and approved the final manuscript.
Annu Näkki, Sanna T Kouhia contributed equally to this work.
About this article
Cite this article
Näkki, A., Kouhia, S.T., Saarela, J. et al. Allelic variants of IL1R1gene associate with severe hand osteoarthritis. BMC Med Genet 11, 50 (2010) doi:10.1186/1471-2350-11-50
- Linkage Disequilibrium Block
- IL1R1 Gene
- Haplotype Allele
- Immune Gene Expression
- Haploview Program