Multiplicative interaction of functional inflammasome genetic variants in determining the risk of gout

Introduction

The acute gout flare results from a localised self-limiting innate immune response to monosodium urate (MSU) crystals deposited in joints in hyperuricaemic individuals. Activation of the caspase recruitment domain-containing protein 8 (CARD8) NOD-like receptor pyrin-containing 3 (NLRP3) inflammasome by MSU crystals and production of mature interleukin-1β (IL-1β) is central to acute gouty arthritis. However very little is known about genetic control of the innate immune response involved in acute gouty arthritis. Therefore our aim was to test functional single nucleotide polymorphism (SNP) variants in the toll-like receptor (TLR)-inflammasome-IL-1β axis for association with gout.

Methods

1,494 gout cases of European and 863 gout cases of New Zealand (NZ) Polynesian (Māori and Pacific Island) ancestry were included. Gout was diagnosed by the 1977 ARA gout classification criteria. There were 1,030 Polynesian controls and 10,942 European controls including from the publicly-available Atherosclerosis Risk in Communities (ARIC) and Framingham Heart (FHS) studies. The ten SNPs were either genotyped by Sequenom MassArray or by Affymetrix SNP array or imputed in the ARIC and FHS datasets. Allelic association was done by logistic regression adjusting by age and sex with European and Polynesian data combined by meta-analysis. Sample sets were pooled for multiplicative interaction analysis, which was also adjusted by sample set.

Results

Eleven SNPs were tested in the TLR2, CD14, IL1B, CARD8, NLRP3, MYD88, P2RX7, DAPK1 and TNXIP genes. Nominally significant (P < 0.05) associations with gout were detected at CARD8 rs2043211 (OR = 1.12, P = 0.007), IL1B rs1143623 (OR = 1.10, P = 0.020) and CD14 rs2569190 (OR = 1.08; P = 0.036). There was significant multiplicative interaction between CARD8 and IL1B (P = 0.005), with the IL1B risk genotype amplifying the risk effect of CARD8.

Conclusion

There is evidence for association of gout with functional variants in CARD8, IL1B and CD14. The gout-associated allele of IL1B increases expression of IL-1β – the multiplicative interaction with CARD8 would be consistent with a synergy of greater inflammasome activity (resulting from reduced CARD8) combined with higher levels of pre-IL-1β expression leading to increased production of mature IL-1β in gout.

The immediate cause of gout is the deposition of monosodium urate (MSU) crystals in and around body tissues, particularly joints [1]. Initially, these deposits trigger a localised and self-limiting inflammatory response (acute gouty arthritis), which becomes increasingly frequent and severe, involving multiple joints and associated with fever. Monosodium urate crystals form under hyperuricaemic conditions when serum urate levels exceed the physiological saturation level (approximately 6.8 mg/dL; approximately 0.41 mM). The most significant biological cause of hyperuricaemia is relatively low renal clearance of uric acid [2, 3]. This is consistent with findings from genome-wide association studies in which 28 loci associated with serum urate levels have been identified, some of which are in genes involved in renal uric acid handling [4, 5]. Predictably most, but not all, of the 28 loci have been associated with gout [4, 6].

Although hyperuricaemia is a prerequisite for MSU formation, only a relatively small proportion of individuals with hyperuricaemia develop gout [7]. This indicates that beside genetic variants associated with urate metabolism and excretion, other factors contribute to the pathogenesis of gout. MSU crystals play an important role in activation of the innate immune system [8] and the recognition of gout as an auto-inflammatory disorder is consistent with the results of functional studies [9, 10]. Variations within genes of the innate immune system may therefore determine whether MSU crystals trigger an inflammatory reaction in susceptible individuals, leading to acute gout; while in others, no inflammation is elicited. Genetic variants that influence the activation and function of the NOD-like Receptor Pyrin containing 3 (NLRP3) inflammasome are candidate genes in this context [11]. The multi-protein inflammasome complex, comprising the NLRP3 polypeptide, ASC or PYCARD (apoptosis-associated speck-like protein containing a CARD) and caspase-1 [12] forms when monocytes and macrophages encounter damaged and pathogen-associated molecular pattern proteins (DAMPs and PAMP; e.g., bacterial lipopolysaccharide or MSU crystals) and leads to activation of caspase-1. Active caspase-1 processes the pro-interleukin (IL)-1β to the mature pro-inflammatory cytokine IL-1β that is then secreted [12].

CARD8 (also known as TUCAN or Cardinal) is a protein with a caspase-domain that interacts with caspase-1 and inhibits its activation [13] and also with a FIIND domain that binds to NLRP3 preventing its recruitment into the active inflammasome complex [14, 15]. Genetic associations between variants of CARD8 and autoimmune diseases have been previously reported (reviewed in [16]). The T allele of CARD8 rs2043211 (C10X) has been associated with increased risk of gout in Chinese [17], and rs2149356 in toll-like receptor 4 (TLR4), a receptor functionally implicated in MSU-stimulated inflammation [18], has also been associated with serum IL-1β levels and the risk of gout in Han Chinese [19]. However, using a haplotype tagging approach, there is no evidence of association between NLRP3 and gout in Chinese [20].^{Footnote 1}

Our aim was to extend the findings from the Han Chinese population [17] and to test genetic variants influencing inflammasome function for association with gout in other population groups. Eleven functional variants were tested in eight genes involved in the MSU crystal-mediated activation of the NLRP3-inflammasome and production of mature IL-1β for association with gout in people of European and New Zealand (NZ) Polynesian (Māori and Pacific Island) ancestry. The prevalence of gout in the NZ Polynesian population is 6–8 % (compared to 3 % in NZ European), exhibiting the highest prevalence worldwide [21, 22].

Participants, ethics and consent

The study was carried out on sample sets comprising: a NZ Polynesian sample set consisting of 1,893 individuals of Samoan, Tongan, Niuean, Cook Island Māori and NZ Māori descent (Table 1; 863 cases and 1,030 controls); and a sample set consisting of 1,684 cases of European ancestry (957 recruited from NZ and Australia and 727 recruited from Europe (‘Eurogout’) [23]) and 882 NZ European controls. All people with gout met the 1977 preliminary American Rheumatism Association classification criteria for gout [24]. New Zealand gout cases were recruited from the Auckland and Christchurch regions of NZ, Australian gout cases from an outpatient clinic in Adelaide, Australian private practice rheumatologists and from a previously reported pharmacogenetic study [25]. European gout cases were recruited from outpatient clinics in Edinburgh, Lausanne, Dresden, Arnhem and Nijmegen. The NZ controls did not have any self-reported history of arthritis, were >17 years of age and were convenience sampled from the Auckland, Christchurch and Otago regions of NZ. The New Zealand Multi-Region Ethics Committee (MEC/105/10/130) and these institutional committees in Europe and Australia granted ethical approval: Research Ethics Committee, University of New South Wales; Ethikkommission, Technische Universität Dresden (EK 8012012); South East Scotland Research Ethics Committee (04/S1102/41); Commission Cantonale (VD) D'éthique de la Recherche sur l'être Humain, Université de Lausanne; Commissie Mensgebonden Onderzoek regio Arnhem Nijmegen. All subjects gave written and informed consent.

Control subjects of European ancestry, who had been genotyped genome-wide, were also included from the Atherosclerosis Risk in Communities (ARIC) study (n = 6,989) and the Framingham Heart Study (FHS; n = 3,131). Participants with self-reported gout were excluded. The sample sets were also screened to remove all but one representative of any closely related family grouping (full/half siblings and parent/child duos and trios). The ARIC study and FHS analyses (project #834) were approved by the relevant Database of Genotype and Phenotype (dbGaP; [26]) Data Access Committees.

Genotyping and imputation

A Sequenom MassARRAY System was used for genotyping the 4,269 individuals without genome-wide genoytpe data available as previously described [27]. Publicly available genome-wide genotype data from ARIC and FHS were imputed as required to obtain genotypes from the eleven selected SNPs. Imputations were carried out by IMPUTE V2.2 software based on the combined 1000 Genomes V3 population reference set. Individuals and SNPs with a call rate of <0.98 were excluded from analysis, as were monomorphic markers and those with a minor allele frequency of <0.01, and a post-imputation quality threshold of 0.30 was used. Imputations were successful for all eleven SNPs and all imputed SNPs were in Hardy-Weinberg equilibrium (P >0.005) with the exception of NLRP3 rs7512998 in ARIC (P _HWE = 1.3 × 10^-27) - these data were excluded from analysis. Where necessary, imputed data were strand-adjusted to match the strand interrogated by the Sequenom data prior to analysis.

Association testing

Allelic association analysis was carried out using R software and adjusting for age and sex. For Polynesian samples, a genetically estimated proportion of Polynesian ancestry, calculated as previously described [27], was included as an additional covariate. All odds ratios are reported relative to the minor allele present in the genotyped NZ European control sample set. Meta-analyses were carried out using METAL [28] with weighting based on standard error and log (odds ratio (OR)) as the effect variable. Significant association was declared if P was <0.05 in the European and Polynesian meta-analysis. No correction for multiple testing was applied because there is prior functional evidence for each variant tested.

Power curves are shown in Additional file 1. The Polynesian sample set was adequately powered to detect common effects (minor allele frequency (MAF) >0.1) of OR ≥1.4, whereas the European sample set was adequately powered to detect common effects of MAF >0.2 at OR >1.2 and stronger effects (OR >1.4) at lower MAF (>0.05).

Interaction analysis

As reviewed by Cordell [29] there are a number of ways to test for interaction or a departure from additivity. We used Stata 13.1 software to carry out logistic regression analysis comparing the disease risk for heterozygosity and minor allele homozygosity for each locus individually and in combination. Previously we reported a large allele frequency difference at ABCG2 rs2231142 and heterogeneity in association with gout between Eastern (EP) and Western (WP) Polynesian sample sets [30]. Similarly, there are differences in allele and genotype frequencies between the two groups (Additional file 2), and heterogeneity will be magnified in combined genotypes. Therefore, in the interaction analysis adjustment was made separately for EP and WP, along with European data. An interaction term was included and a P value <0.017 (Bonferroni-adjusted for number of interaction analyses performed) was considered to indicate significant multiplicative interaction between genetic variants.

Eleven functional genetic variants in NLRP3, CARD8, IL1B, DAPK1, TXNIP, TLR2, P2XR7, MYD88 and CD14 were selected from the literature (Table 2) and genotyped; genotype distributions are presented in Additional file 2. There was nominal allelic association (P <0.05) for three variants in the combined European and Polynesian analysis (Table 3) - IL1B, CARD8 and CD14 (OR = 1.10, 1.12 and 1.08, respectively). CARD8 rs2043211 was also associated with gout in Europeans (OR = 1.11).

Because of reported interactions between NLRP3 and CARD8 in other auto-inflammatory conditions [30, 31], and the biological interaction between the inflammasome and IL-1β we tested for pairwise multiplicative interaction between NLRP3/CARD8, NLRP3/IL1B and CARD8/IL1B (Table 4), using only rs10754558 of NLRP3 owing to the low MAF of rs35829419 in both Europeans and Polynesians and the low MAF of rs7512998 in Polynesians. There was evidence for interaction between CARD8 and IL1B (P = 0.005, P _c = 0.015), driven by amplification of the risk conferred by the CARD8 rs2043211 T allele in the presence of the rs1143623 IL1B minor allele homozygous (GG) genotype. Nominally significant interaction between NLRP3/IL1B (P _Nominal  = 0.048) was also observed, although this was not significant after adjustment for multiple testing (P _c = 0.144).

TLR signalling via the NLRP3 inflammasome has been implicated in gout susceptibility and pathology in vivo and in vitro [10]; for example, MSU uptake and IL-1β production by bone marrow-derived macrophages derived from TLR2, TLR4 or Myd88 knockout mice is significantly reduced, as is neutrophil influx, in response to subcutaneous injection of MSU in whole animals [18]. To further elucidate the role of the TLR-inflammasome-IL-1β cascade in gout pathogenesis, eleven candidate genetic variants that functionally impact on this pathway (reviewed in [10]) were tested for association with gout in a sample set of 2,357 cases, adequately powered to detect association with common variants having an effect size of odds ratio 1.4 or greater (Additional file 1). As discussed below, the nominal evidence for association between gout and CD14, CARD8 and IL1B, and the multiplicative interaction between CARD8 and IL1B in determining the risk of gout, support the considerable evidence that TLR-mediated activation of the inflammasome and subsequent release of active IL-1β is a central causal pathogenic pathway of gout [10, 32].

Variants rs2043211 (CARD8), rs1143623 (IL1B) and rs2569190 (CD14), which were associated with gout, are functional variations in genes directly involved in the NLRP3 signaling pathway, and as such are likely to represent genuine disease-susceptibility loci. CARD8 is an adaptor protein that regulates IL-1β secretion by inhibiting NFKβ signaling (required for the expression of pro-IL-1β) and/or interacting with caspase 1 or NLRP3 to inhibit the generation of active IL-1β from inactive pro-IL-1β [13, 15]. The effect size of CARD8 SNP rs2043211 was consistent between the European and NZ Polynesian samples sets (OR = 1.11, P = 0.023 and OR = 1.15, P = 0.078, respectively; combined OR = 1.12, P = 0.007) and the Chinese sample set reported by Chen et al. (OR = 1.19, P = 0.08) [17]. Collectively our results and the study by Chen et al. [17] do suggest that the association of the minor allele of rs2043211 with gout is not a false positive one. SNP rs2043211 encodes a missense protein variation (C10X or F52I depending on transcript, [33]) with the minor allele increasing the risk of gout. Although the functional effect of this variation has not been specifically evaluated, the SNP is within an expression quantitative trait locus peak and carriage of the minor allele is associated with decreased CARD8 expression [34]. It has been inconsistently associated with other auto-inflammatory phenotypes, with some evidence for epistatic interaction with NLRP3 rs35829419 in determining risk of inflammatory bowel disease and abdominal aortic aneurysms [31, 35]. However, we found no evidence of interaction between NLRP3 and CARD8 in determining risk of gout (Table 4) (although we did not specifically analyse rs35829419), and it has also been suggested that under certain conditions, ASC-dependent IL1-β production in response to MSU stimulation can occur in the absence of NLRP3 [36].

IL1B SNP rs1143623 is within a promoter GATA transcription factor family binding site. Although the minor allele exhibits enhanced protein-binding [37] and decreased expression in vitro and in vivo, the effect of this variation seems to be influenced both by the wider promoter haplotype [38–40], and the identity of the stimulatory signal; the minor allele shows increased rather than decreased expression in vitro in response to TNF-α [40], and has also been associated with increased post-prandial triglyceride and IL6 (an effector of IL-1β) levels [41]. The minor allele is over-represented in gout cases compared to controls (OR = 1.10, P = 0.020). If this is replicated it would be consistent with an etiological role for increased IL1B expression in gout. There was evidence of multiplicative interaction with CARD8 rs2043211 in which the IL1B rs1143623 minor (risk) allele homozygous genotype appeared to amplify the effect of the minor allele of rs2043211. This would be consistent with a synergy of greater inflammasome activity (resulting from reduced CARD8) combined with higher levels of pre-IL-1β expression leading to increased production of mature IL-1β in gout.

The final variation nominally associated with gout was CD14 SNP rs2569190. This SNP is in the 5’UTR of one of the two CD14 splice variants, with the minor allele (that increases risk of gout) increasing expression in monocytes by decreasing affinity for the inhibitory Sp3 transcription factor [42], enhancing the loading of RNA polymerase II [43] and is associated with increased soluble CD14 levels in healthy individuals [44, 45]. Membrane-bound CD14 forms functional complexes with TLR2 or TLR4 and leukocyte β2-integrins, which could mediate TLR dimerization and optimize the innate immune response to MSU crystals [10].

In conclusion we provide evidence for association of gout with functional innate immune system variants in CARD8, IL1B and CD14, and multiplicative interaction between IL1B and CARD8. The findings involving IL1B are consistent with genetically determined levels of IL-1β being important in gout.

1. 1.

   Note that the Meng et al paper [20] reported positive association with single nucleotide polymorphism (SNP) rs7512998, however the P value for this SNP had been incorrectly calculated - the true P value for rs7512998 was 0.43 - as reported in a corrigendum [61].

ARA:

American Rheumatology Association

ARIC:

Atherosclerosis Risk in Communities

ASC:

apoptosis-associated speck-like protein containing a CARDD

CARD8:

caspase recruitment domain-containing protein 8

DAMP:

damage-associated molecular pattern

EP:

Eastern Polynesian

FHS:

Framingham Heart Study

IL:

interleukin

MAF:

minor allele frequency

MSU:

monosodium urate

NALP3:

NACHT, LRR and PYD domains-containing protein 3

NLRP3:

NOD-like Receptor Pyrin containing 3

NOD:

nucleotide-binding oligomerisation domain

NZ:

New Zealand

OR:

odds ratio

PAMP:

pathogen-associated molecular pattern

UTR:

untranslated region

RNA:

ribonucleic acid

SNP:

single nucleotide polymorphism

TLR:

toll-like receptor

TNF:

tumour necrosis factor

WP:

Western Polynesian

This work was supported by the Health Research Council of New Zealand, Arthritis New Zealand, New Zealand Lottery Health and the University of Otago. The authors would like to thank Jill Drake (Canterbury District Health Board), Roddi Laurence, Chris Franklin, Meaghan House (all University of Auckland) and Gabrielle Sexton (University of Otago) for recruitment. Matthew Brown, Linda Bradbury and The Arthritis Genomics Recruitment Initiative in Australia network are acknowledged. The European Crystal Network was formed after the first European Crystal Workshop in Paris, March 2010 (Prof Frédéric Lioté, Paris, and Prof Alexander So, Lausanne, convenors). The Atherosclerosis Risk in Communities Study is carried out as a collaborative study supported by National Heart, Lung, and Blood Institute contracts N01-HC-55015, N01-HC-55016, N01-HC-55018, N01-HC-55019, N01-HC-55020, N01-HC-55021, N01-HC-55022, R01HL087641, R01HL59367 and R01HL086694; National Human Genome Research Institute contract U01HG004402; and National Institutes of Health contract HHSN268200625226C. The authors thank the staff and participants of the ARIC study for their important contributions. Infrastructure was partly supported by Grant Number UL1RR025005, a component of the National Institutes of Health and NIH Roadmap for Medical Research. The Framingham Heart Study and the Framingham SHARe project are conducted and supported by the National Heart, Lung, and Blood Institute (NHLBI) in collaboration with Boston University. The Framingham SHARe data used for the analyses described in this manuscript were obtained through dbGaP. This manuscript was not prepared in collaboration with investigators of the Framingham Heart Study and does not necessarily reflect the opinions or views of the Framingham Heart Study, Boston University, or the NHLBI.

Affiliations

1. Department of Biochemistry, University of Otago, Box 56, Dunedin, 9054, New Zealand

   • Cushla McKinney
   • , Ruth K. Topless
   •  & Tony R. Merriman

2. Department of Medicine, University of Otago, Christchurch, PO Box 4345, Christchurch, New Zealand

   • Lisa K. Stamp

3. Department of Medicine, University of Auckland, Auckland, New Zealand

   • Nicola Dalbeth

4. School of Medical Sciences, University of New South Wales, Sydney, Australia

   • Richard O. Day
   • , Diluk RW Kannangara
   •  & Kenneth M. Williams

5. Department of Clinical Pharmacology & Toxicology, St Vincent’s Hospital, Sydney, Australia

   • Richard O. Day
   • , Diluk RW Kannangara
   •  & Kenneth M. Williams

6. Department of Rheumatology, Rijnstate Hospital, Arnhem, The Netherlands

   • Matthijs Janssen

7. Department of IQ HealthCare, VieCuri Medical Centre, Venlo, The Netherlands

   • Timothy L. Jansen

8. Scientific Institute of Quality in HealthCare, Radboud University Medical Centre, Nijmegen, The Netherlands

   • Timothy L. Jansen

9. Department of Internal Medicine and Radboud Institute of Molecular Life Science, Radboud University Medical Center, Nijmegen, The Netherlands

   • Leo A. Joosten

10. Department of Rheumatology and Clinical Immunology, Laboratory of Translational Immunology, University Medical Centre Utrecht, PO Box 85500, 3508, GA, Utrecht, The Netherlands

    • Timothy R. Radstake

11. Department of Immunology, University Medical Centre Utrecht, PO Box 85500, 3508, GA, Utrecht, The Netherlands

    • Timothy R. Radstake

12. Rheumatic Diseases Unit, Institute of Genetics and Molecular Medicine, University of Edinburgh, Edinburgh, UK

    • Philip L. Riches

13. Department of Rheumatology, University Clinic Carl-Gustav-Carus”, Dresden, Germany

    • Anne-Kathrin Tausche

14. INSERM, UMR-S 1132, Hospital Lariboisière, F-75010, Paris, France

    • Frederic Lioté

15. University Paris Diderot (UFR de Médecine), Sorbonne Paris Cité, F-75205, Paris, France

    • Frederic Lioté

16. DAL, Service of Rheumatology, Laboratory of Rheumatology, University of Lausanne, CHUV, Nestlé 05-5029, 1011, Lausanne, Switzerland

    • Alexander So

Corresponding author

Correspondence to Tony R. Merriman.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

CM, RKT and TRM helped to design the study, oversee its execution, and prepare the manuscript. LKS, ND, ROD, DRWK, KMW, MJ, TLJ, LAJ, TRR, PLR, A-KT, FL and AS helped to provide clinical recruitment and prepare the manuscript. All authors read and approved the final manuscript.

Reprints and Permissions