Association of the interleukin-1 gene cluster with osteoarthritis of the hip: comment on the article by Meulenbelt et al and the letter by Smith et al

Interleukin-17 and Toll-like Receptor 10 genetic polymorphisms and susceptibility to large joint osteoarthritis

Primary osteoarthritis (OA) is the most common type of a joint disease. It has a polygenic risk inheritance pattern and affects older people. The etiology of this disease is not fully understood. The aim of this study was to investigate the associations between polymorphisms in pro-inflammatory interleukin-17 (IL17A and IL17F) and anti-inflammatory Toll-like Receptor 10 (TLR10) genes with the risk for development of advanced stage hip and knee primary OA in the Croatian population. A total of 500 OA patients and 597 controls were genotyped for IL17A SNP (rs2275913), IL17F SNPs (rs763780 and rs1889570), and TLR10 (rs11096957) genes. The allelic and genotypic frequencies of IL17F SNP (rs763780) showed statistically significant differences in comparisons of controls with hip-but not knee-OA patients. The major allele (T) of rs763780 was associated with the lower risk for developing hip OA (p = 7.9 × 10^-4 , OR = 0.45, 95%CI = 0.27-0.74), whereas the minor allele (C) was associated with susceptibility to hip OA (p = 7.9 × 10^-4 , OR = 2.24, 95%CI = 1.35-3.72). The genotype T/T was associated with the protection to hip OA (p = 3.9 × 10^-4 , OR = 0.41, 95%CI = 0.24-0.70), and, lastly, the genotype T/C was associated with the higher risk to acquiring hip OA (p = 2.6 × 10^-4 , OR = 2.50, 95%CI = 1.47-4.25). TLR10 SNP rs11096957 was found significantly associated with predisposition to hip OA (p = 0.04, OR = 1.41, 95%CI = 1.02-1.94) but not knee OA. Our findings suggest that hip OA in Croatian population might have a different genetic risk regarding the IL17 and TLR10 gene locus than knee OA. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:1684-1693, 2018.

Relationships of common polymorphisms in IL-6, IL-1A, and IL-1B genes with susceptibility to osteoarthritis: a meta-analysis

Observational and experimental studies have arrived at inconsistent conclusions about whether common polymorphisms in IL-6, IL-1A, and IL-1B genes are associated with an increased risk of osteoarthritis (OA). Therefore, we undertook a comprehensive meta-analysis to more systematically summarize the relationships of IL-6, IL-1A, and IL-1B genetic polymorphisms with susceptibility to OA. We screened the PubMed, Embase, Web of Science, Cochrane Library, CISCOM, CINAHL, Google Scholar, China BioMedicine (CBM), and China National Knowledge Infrastructure (CNKI) databases up to 31 March 2014. We used STATA software to analyze statistical data. Odds ratios (ORs) and their corresponding 95 % confidence intervals (95 % CIs) were calculated. Seventeen independent case-control studies were included in this meta-analysis with a total number of 7,491 subjects, comprised of 3,293 OA patients and 4,729 healthy controls. Our results indicate that IL-6, IL-1A, and IL-1B genetic polymorphisms are statistically correlated with an increased risk of OA under the allele and dominant models. According to a subgroup analysis based on disease, a higher frequency of IL-6 genetic polymorphisms was observed among knee OA and hand OA patients, but not among hip OA and DIP OA patients. A higher frequency of IL-1A genetic polymorphisms were found among hip OA patients, hand OA, hip OA and DIP OA patients. Furthermore, we observed a higher IL-1B polymorphism frequency among knee OA and hip OA patients, but not among hand OA patients. Our findings provide evidence that IL-6, IL-1A, and IL-1B genetic polymorphisms may be correlated with susceptibility to OA.

Hip osteoarthritis susceptibility is associated with IL1B -511(G>A) and IL1 RN (VNTR) genotypic polymorphisms in Croatian Caucasian population

Among the predisposing factors to osteoarthritis (OA), a frequent destructive joint disease, is the complex genetic heritage including the interleukin-1 family members like the IL1β (IL1B) and the IL1 receptor antagonist (IL1RN) genes. The aim of this study was to investigate allelic and genotypic frequencies of the IL1B gene single nucleotide polymorphism (SNP) at -511(G>A) and the variable number tandem repeat (VNTR) in the IL1RN gene in a Croatian Caucasian population of hip OA (HOA) cases and healthy controls. A total of 259 HOA patients with total hip replacement (THR) and 518 healthy blood donors as controls were genotyped for IL1B gene SNP -511(G>A) and the VNTR in the IL1RN gene associated with HOA. The genotype G/A (1/2) at IL1B was significantly associated with the protection of the HOA (p < 0.036, OR = 0.72, 95% CI = 0.52-0.99). The genotype G/G (1/1) had only a trend towards the susceptibility (p = 0.053, OR = 1.35, 95% CI = 0.98-1.86) to disease. None of the haplotypes IL1B -511(G>A) and IL1RN (VNTR) were found associated with the HOA. The haplotype 1-2 at these loci had only a trend to susceptibility (p = 0.065). Haplotype 1-3 had a significant male bias in diseased. Furthermore, genotype comprising 2-1/2-2 haplotypes was found significantly associated with predisposition to HOA (p = 0.027, OR = 2.23, 95% CI = 1.03-4.88), whereas genotype 1-1/2-2 with protection to disease (p = 0.028, OR = 0.65, 95% CI = 0.43-0.97). Our findings suggest that HOA in Croatian population might have a different genetic risk regarding the IL1 locus than has been reported for other Caucasian populations previously.

Large-scale meta-analysis of interleukin-1 beta and interleukin-1 receptor antagonist polymorphisms on risk of radiographic hip and knee osteoarthritis and severity of knee osteoarthritis

OBJECTIVE

To clarify the role of common genetic variation in the Interleukin-1β (IL1B) and Interleukin-1R antagonist (IL1RN) genes on risk of knee and hip osteoarthritis (OA) and severity of knee OA by means of large-scale meta-analyses.

METHODS

We searched PubMed for articles assessing the role of IL1B and IL1RN polymorphisms/haplotypes on the risk of hip and/or knee OA. Novel data were included from eight unpublished studies. Meta-analyses were performed using fixed- and random-effects models with a total of 3595 hip OA and 5013 knee OA cases, and 6559 and 9132 controls respectively. The role of ILRN haplotypes on radiographic severity of knee OA was tested in 1918 cases with Kellgren-Lawrence (K/L) 1 or 2 compared to 199 cases with K/L 3 or 4.

RESULTS

The meta-analysis of six published studies retrieved from the literature search and eight unpublished studies showed no evidence of association between common genetic variation in the IL1B or IL1RN genes and risk of hip OA or knee OA (P>0.05 for rs16944, rs1143634, rs419598 and haplotype C-G-C (rs1143634, rs16944 and rs419598) previously implicated in risk of hip OA). The C-T-A haplotype formed by rs419598, rs315952 and rs9005, previously implicated in radiographic severity of knee OA, was associated with reduced severity of knee OA (odds ratio (OR)=0.71 95%CI 0.56-0.91; P=0.006, I(2)=74%), and achieved borderline statistical significance in a random-effects model (OR=0.61 95%CI 0.35-1.06 P=0.08).

CONCLUSION

Common genetic variation in the Interleukin-1 region is not associated with prevalence of hip or knee OA but our data suggest that IL1RN might have a role in severity of knee OA.

Interleukin-1 gene cluster variants with innate cytokine production profiles and osteoarthritis in subjects from the Genetics, Osteoarthritis and Progression Study

OBJECTIVE

To assess whether genetic variation in the interleukin-1 (IL-1) gene cluster contributes to familial osteoarthritis (OA) by influencing innate ex vivo production of IL-1beta or IL-1 receptor antagonist (IL-1Ra).

METHODS

Innate ex vivo IL-1beta and IL-1Ra production upon lipopolysaccharide (LPS) stimulation of whole blood cells was measured in subjects from the Genetics, Osteoarthritis and Progression (GARP) Study, which includes sibling pairs in which at least one sibling has symptomatic OA at multiple sites. Radiographic OA (ROA) was assessed by Kellgren/Lawrence score. Subjects from the GARP Study and controls from the Rotterdam Study were genotyped for 7 single-nucleotide polymorphisms (SNPs) encompassing the IL-1 gene cluster on chromosome 2q13. Linkage disequilibrium analysis and genotype and haplotype association analysis were performed to assess the relationship between the IL-1 gene cluster SNPs, innate ex vivo cytokine production, and OA.

RESULTS

Among subjects in the GARP Study, the haplotype variable-number tandem repeat in intron 2/T+8006C/T+11100C 2/2/1 of the IL1RN gene was significantly associated with reduced innate ex vivo bioavailability of IL-1beta upon LPS stimulation (P = 0.026) and with ROA at the highest number of joint locations.

CONCLUSION

These results show that genetic variation at the IL-1 gene cluster is associated with lower IL-1beta bioavailability and with OA at a large number of joint locations. The data further indicate that, among subjects with OA affecting the highest number of joints, the innate immune system may be activated, thereby obscuring possible underlying mechanisms.

A genome-wide linkage scan reveals CD53 as an important regulator of innate TNF-alpha levels

Cytokines are major immune system regulators. Previously, innate cytokine profiles determined by lipopolysaccharide stimulation were shown to be highly heritable. To identify regulating genes in innate immunity, we analyzed data from a genome-wide linkage scan using microsatellites in osteoarthritis (OA) patients (The GARP study) and their innate cytokine data on interleukin (IL)-1beta, IL-1Ra, IL-10 and tumor necrosis factor (TNF)alpha. A confirmation cohort consisted of the Leiden 85-Plus study. In this study, a linkage analysis was followed by manual selection of candidate genes in linkage regions showing LOD scores over 2.5. An single-nucleotide polymorphism (SNP) gene tagging method was applied to select SNPs on the basis of the highest level of gene tagging and possible functional effects. QTDT was used to identify the SNPs associated with innate cytokine production. Initial association signals were modeled by a linear mixed model. Through these analyses, we identified 10 putative genes involved in the regulation of TNFalpha. SNP rs6679497 in gene CD53 showed significant association with TNFalpha levels (P=0.001). No association of this SNP was observed with OA. A novel gene involved in the innate immune response of TNFalpha is identified. Genetic variation in this gene may have a role in diseases and disorders in which TNFalpha is closely involved.

Interleukin-1 region meta-analysis with osteoarthritis phenotypes

OBJECTIVE

Several research groups have examined osteoarthritis (OA) association with Interleukin-1 (IL-1) region markers and haplotypes. The results have been suggestive for hand OA, negative for knee OA, and conflicting for hip OA.

DESIGN

Our aim was to address conflicts employing meta-analytical methods on data from 1238 European-descent cases with various OA phenotypes and 1269 European-descent controls from four study centers. We imputed some missing genotype data and reconstructed IL-1 region extended haplotypes. A previously reported 7-marker IL1A-IL1B-IL1RN extended risk haplotype was tested for association with each specific index phenotype.

RESULTS

For hip OA, data from three centers showed heterogeneity of extended-risk-haplotype effect, two panels showing trend toward risk and another showing protection, with overall odds ratio (OR) 1.24 (95% Confidence interval (CI) 0.45-3.41, P 0.67). The heterogeneity fell partly along control ascertainment lines, chiefly between controls ascertained as spouses of arthroplasty patients and controls identified through population radiographic survey. For knee OA, the results showed no heterogeneity and no significant extended-risk-haplotype effect. For hand OA, the results showed little heterogeneity and a modest trend toward positive association (summary OR 1.34, 95% CI 0.83-2.17 P 0.23). Using a Bayesian partition modeling approach, the 7-marker extended haplotypes showed no significant effect on any OA phenotype examined. A 3-single-nucleotide polymorphism (SNP) IL1B-IL1RN haplotype rs1143627-rs16944-rs419598 showed a trend toward hand OA association (posterior probability of association 0.72) with the most prominent feature being protection from a specific haplotype representing a partial mirror image of the extended risk haplotype (OR estimated at 0.46).

CONCLUSIONS

The meta-analysis data do not confirm but only suggest that some hand and hip OA risk could be associated with the IL-1 region, particularly centered in IL1B and possibly also IL1RN.

Association Between Interleukin 1 Gene Cluster Polymorphisms and Bilateral Distal Interphalangeal Osteoarthritis

Objective.

To examine the association of the interleukin 1 gene (IL1) cluster polymorphisms and their haplotypes with bilateral distal interphalangeal joint osteoarthritis (DIP OA).

Methods.

Radiographs of both hands of 295 dentists and 248 teachers were examined and classified for the presence of OA using reference images. Bilateral DIP OA was defined by the presence of radiographic findings of grade 2 or more in at least 1 symmetrical pair of the DIP joints. We genotyped 10 single-nucleotide polymorphisms (SNP) in theIL1R1, IL1RL2, IL1A, IL1B,andIL1RNgenes using polymerase chain reaction-based methods. Haplotypes were statistically reconstructed using the PHASE program. The association between the genotypes/diplotypes and bilateral DIP OA was examined with logistic regression analysis.

Results.

TwoIL1BSNP (rs1143634 and rs1143633) were associated with bilateral DIP OA. The carriers of theIL1Brs1143634 minor allele had an increased OA risk [odds ratio (OR) 1.6; 95% confidence interval (CI) 1.08–2.26] compared to the noncarriers. The association was stronger in the dentists. The distribution of theIL1Brs1143633 genotype fit a recessive mode of inheritance (OR 3.03, 95% CI 1.35–6.83, p = 0.006). TwoIL1B-IL1RNextended haplotype alleles (211-1 and 121-1) were associated with bilateral DIP OA. An interaction between theIL1Brs1143634 and theIL1R1-IL1RL2andIL1B-IL1RNextended haplotypes and occupation (increased risk of OA among dentists only) was observed.

Conclusion.

Our results provide further evidence for the role ofIL1gene cluster polymorphisms in the etiology of OA and suggest that some of these may predispose DIP joints to the effects of mechanical overload.

Potential influence of IL1B haplotype and IL1A-IL1B-IL1RN extended haplotype on hand osteoarthritis risk

OBJECTIVE

To assess osteoarthritis (OA) association with the human interleukin-1 (IL-1) region.

DESIGN

Sixty-four European-descent cases with radiographic hand OA and 48 European-descent controls were genotyped at nine single nucleotide polymorphism (SNP), one variable-number-of-tandem-repeat (VNTR), and one microsatellite marker extending across loci for IL-1alpha (IL1A), IL-1beta (IL1B), and IL-1 receptor antagonist (IL1RN). The genotype data were used to reconstruct individual locus haplotypes, and then locus haplotypes were used as superalleles for extended haplotype reconstruction.

RESULTS

Nine different extended IL1A-IL1B-IL1RN haplotypes occurred at a frequency 0.05 or greater in either cases or controls. Only two IL1A-IL1B-IL1RN extended haplotypes were consistent with previously described extended risk haplotypes and totaled n=9 in cases and n=3 in controls [odds ratio (OR) 2.1, Haldane's chi(2) 1.67, one-sided P 0.1]. Our prior report showed hand OA association with homozygous IL1B rs1143633 minor allele genotype. All except one extended risk haplotype copy also had the IL1B rs1143633 minor allele. The rs1143633 genotype association was explained by one common six-SNP IL1B haplotype bearing rs1143633 minor allele and also risk alleles at rs1143634, rs1143627, and rs16944, component markers of the previously described extended risk haplotypes. The IL1B haplotype bearing all three risk alleles was found in 16 haplotype-homozygous hand OA cases and in four haplotype-homozygous controls and conferred OR 3.4 among homozygotes (nominal P value 0.006).

CONCLUSION

Our evidence broadly supports the genetic association of OA phenotypes with an IL-1 region extended risk haplotype and specifically IL1B genotype. The extended risk haplotype previously associated with hip OA appears to be less frequent and has weaker genetic effect in hand OA. Hand OA risk is conferred by homozygous state for the IL1B haplotype characteristic of the extended risk haplotype.