Open Access Highly Accessed Research article

Epistasis with HLA DR3 implicates the P2X7 receptor in the pathogenesis of primary Sjögren's syndrome

Susan Lester1, Leanne Stokes23, Kristen K Skarratt2, Ben J Gu24, Kathy L Sivils5, Christopher J Lessard5, James S Wiley24 and Maureen Rischmueller16*

Author Affiliations

1 Department of Rheumatology, The Queen Elizabeth Hospital, Woodville South, South Australia, Australia

2 Sydney Medical School Nepean, University of Sydney, Nepean Hospital, Penrith, NSW, Australia

3 Health Innovations Research Institute, School of Medical Sciences, RMIT University, Bundoora, Victoria, Australia

4 Florey Neuroscience Institutes, University of Melbourne, Parkville, Victoria, Australia

5 Arthritis and Clinical Immunology Research Program, Oklahoma Medical Research Foundation, Oklahoma City, OK, USA

6 Discipline of Medicine, University of Adelaide, South Australia, Australia

For all author emails, please log on.

Arthritis Research & Therapy 2013, 15:R71  doi:10.1186/ar4248


The electronic version of this article is the complete one and can be found online at: http://arthritis-research.com/content/15/4/R71


Received:6 November 2012
Revisions received:2 May 2013
Accepted:2 June 2013
Published:2 June 2013

© 2013 Lester et al.; licensee BioMed Central Ltd.

This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Introduction

The aim of this study was to examine the association between functional polymorphisms in the pro-inflammatory P2X7 receptor and the Ro/La autoantibody response in primary Sjögren's syndrome (pSS).

Methods

Twelve functional P2RX7 polymorphisms were genotyped in 114 pSS patients fulfilling the Revised American-European Consensus Criteria for pSS, and 136 controls. Genotyping of the A1405G (rs2230912) polymorphism was performed on a replication cohort consisting of 281 pSS patients and 534 controls. P2X7 receptor function in lymphocytes and monocytes was assessed by measurement of ATP-induced ethidium+ uptake. Serum IL-18 levels were determined by ELISA.

Results

The minor allele of P2RX7 A1405G is a tag for a common haplotype associated with gain in receptor function, as assessed by ATP-induced ethidium+ uptake. A positive association between 1405G and anti-Ro±La seropositive pSS patients was observed in Cohort 1. Although not replicated in Cohort 2, there was a consistent, significant, negative epistatic interaction effect with HLA-DR3 in seropositive pSS patients from both cohorts, thereby implicating this gain of function variant in the pathogenesis of pSS. Serum IL-18 was elevated in seropositive pSS patients, but was not influenced by P2RX7 A1405G.

Conclusions

The P2RX7 1405G gain-of-function haplotype may be a risk factor for seropositive pSS in a subset of subjects who do not carry HLA risk alleles, but has no effect in subjects who do (epistasis). Potential mechanisms relate to autoantigen exposure and inflammatory cytokine expression. The observed elevation of IL-18 levels is consistent with P2X7 receptor activation in seropositive pSS patients. Collectively these findings implicate P2X7 receptor function in the pathogenesis of pSS.

Introduction

Primary Sjögren's syndrome (pSS) is a systemic autoimmune disease characterized by lymphocytic and monocytic infiltration into the salivary and lacrimal glands with inflammation, destruction of acini and failure of exocrine secretion. A wide range of extraglandular features also manifest in a subset of patients. Autoantibodies targeting the Ro and La components of a ribonucleoprotein (RNP) complex are highly specific for pSS and constitute one of the classification criteria for this disease [1]. Once considered an epiphenomenon, these autoantibodies are now recognized to be involved in the systemic inflammation of pSS [2].

A substantial polygenic component underlies susceptibility to the Ro/La autoantibody response in pSS [3,4], however the strongest individual association is observed with HLA-DRB1*03:01 (DR3, Chr 6p21.3) [5]. Additional genes are also likely to be involved, and one candidate is the P2X7 receptor (Chr 12q24.31), whose chromosomal location lies within mapped human and murine systemic lupus erythematosus (SLE) susceptibility loci [6]. The purinergic P2X7 receptor has been implicated in several murine models of autoimmune disease, including experimental autoimmune encephalomyelitis [7], rheumatoid arthritis (RA) [8] and autoimmune hepatitis [9].

The P2X7 receptor belongs to a two-transmembrane family of receptors, which are trimeric in the membrane and are attached via non-muscle myosin to the actin cytoskeleton [10-12]. Activation of P2X7, in response to the binding of extracellular ATP, opens a non-selective cation channel allowing an immediate influx of Ca2+ and Na+, and efflux of K+ ions. Continued activation, over tens of seconds, allows a larger permeability state to develop (termed pore formation), which is accompanied by extensive membrane blebbing, and ultimately cell death, with features of both apoptosis and necrosis [13,14]. Extracellular ATP is normally maintained at low nanomolar levels by a ubiquitous ecto-ATPDase (CD39), although in an inflammatory milieu, ATP may be released from dead or injured cells, as well as by cell swelling or from autonomic nerves [15].

The P2X7 receptor is found on many cells and tissues throughout the body. The receptor is most highly expressed on cells of the macrophage lineage, including dendritic cells and microglia [13,16], and both P2X7 receptors and HLA Class II molecules are upregulated in parallel when monocytes differentiate to these mature cell types. P2X7 receptors are also abundant on both acinar and ductal cells from submandibular and other salivary glands [17,18]. A special role for this receptor in salivary glands has been highlighted by its direct involvement in the regulation of fluid secretion [19,20]. In cells of the monocyte-macrophage lineage, there is much evidence that P2X7 is pro-inflammatory [16,21], and P2X7 gene-deleted mice show major reductions in cytokine responses to various inflammatory stimuli [8]. Activation of the P2X7 receptor results in maturation and release of the pro-inflammatory cytokines, IL-1β and IL-18 [21-23], a response that has been linked to formation of the cryopyrin (NALP3) inflammasome and subsequent caspase-1 activation [21,24]. Elevated serum levels of both IL-1β and IL-18 have been reported in pSS [25,26], and these cytokines have been implicated in salivary gland pathology in both mouse models [27] and human disease [28,29].

There is considerable variability in P2X7 function between individuals, which can be largely attributed to polymorphisms of the P2RX7 gene. To date, one dozen polymorphisms have been identified in P2RX7, which alter receptor function [30,31]. We and others have shown that polymorphisms that increase or decrease P2X7 receptor function confer a corresponding increase or decrease in ATP-stimulated secretion of IL-1β or IL-18 from monocytes primed with lipopolysaccharide [31-33]. We have now examined these functional P2RX7 SNPs in patients with pSS, and have identified an epistatic interaction between the major P2RX7 gain-of-function haplotype and HLA-DR3 in patients with seropositive pSS.

Materials and methods

Study participants

Cohort 1 included 114 Caucasian pSS patients (88% female) from the South Australian Sjögren's Syndrome research registry, and 136 Caucasian population-based controls. The replication cohort (cohort 2) included 281 Caucasian pSS patients (95% female) obtained from the Oklahoma Medical Research Foundation (OMRF) plus 534 controls.

All patients met the revised 2002 American-European consensus research classification criteria for pSS [1]. Anti-Ro±La status was performed as part of routine diagnosis, and was available for all patients. Autoantibody testing using recombinant protein was performed by ELISA in cohort 1 [4], and by immunodiffusion in cohort 2 [34]. The majority of patients were autoantibody-positive, with 82% of pSS patients from the first cohort, and 71% of pSS patients from the second cohort positive for anti-Ro±La.

HLA DR3 status was available for all patients and controls for cohort 1 (medium resolution DRB1 genotyping performed by the Tissue Typing Laboratory, Adelaide Red Cross Blood Transfusion Service, South Australia), and for the majority of pSS patients from cohort 2 (a DR3 sequence-specific priming assay). A DR3 proxy SNP, rs2187668 [35], was used for analyses of the controls from cohort 2. The use of the minor allele of rs2187668 as a proxy for DR3 was validated against individuals with known DR3 status in each cohort, with observed kappa agreements of 0.95, and 0.92 respectively (data not shown).

The study was approved by the Human Ethics Committees of The Queen Elizabeth Hospital, the Sydney West Area Health Service and OMRF, in accordance with the Helsinki Declaration, and all participants gave informed, written consent.

P2RX7 SNP genotyping

Twelve SNPs (one intron splice site, and eleven non-synonymous) throughout the P2RX7 gene, most of which have been previously characterised as exerting either loss- or gain-of- function over the wild-type receptor, were examined in pSS patients and controls from cohort 1 (Table 1). Eight SNPs were analysed by high-throughput TaqManTM genotyping assays performed on an ABI Prism ABI Prism®7900HT sequence detection system (Applied Biosystems, Foster City, CA, USA) at the SUPAMAC Facility, Royal Prince Alfred Hospital, Sydney, Australia. G1068A was genotyped using the commercial TaqManTM assay, C_11704039_10 (Applied Biosystems). SNPs G151T, G474A, G946A, C1096G, A1405G, A1513C and T1729A were genotyped using primers and probes as listed in Additional file 1, Table S1. Four SNPs, T253C, C489T, G835A and G853A were genotyped by a homogeneous mass extension assay at the Australian Genome Research Facility (St Lucia, Queensland).

Table 1. P2RX7 SNPs genotyped in pSS patients and normal subjects

Additional file 1. Table S1 showing custom Taqman genotyping primers and probes. Table S2 showing P2RX7 minor allele frequency in pSS patients (cohort 1, n = 114) and normal subjects (n = 136). Table S3 showing Genotyping results for P2RX7 A1405G for seropositive pSS, seronegative pSS and control subjects in cohort 1 and cohort 2.

Format: DOCX Size: 20KB Download fileOpen Data

Genotype data for the P2X7 A1405G SNP (rs2230912) and the DR3 proxy, rs2187668 were obtained for cohort 2 using the Illumina OMNI-1 Quad array (San Diego, CA, USA) following manufacturer's protocols. These subjects were also evaluated for possible population stratification using principal component analysis implemented in Eiganstrat [36] utilizing approximately 30,000 independent SNPs (r2 <0.20) scattered thought the genome. Samples were plotted alongside the HapMap populations to determine which ones were genetic outliers.

Serum IL-18 assays

Serum IL-18 levels were assayed by sandwich ELISA (Human IL-18 ELISA Kit, Medical and Biological Laboratories Company Limited, Nagaoya, Japan, catalogue number 7620), according to the manufacturer's instructions.

P2X7 function measured by ATP-induced ethidium uptake (pore assay)

Peripheral venous blood was obtained from healthy volunteers and mononuclear cells were isolated by Ficoll-Paque gradient centrifugation. P2X7 receptor function was measured by the ATP-induced ethidium+ uptake in a mixed cell (B and T lymphocytes and monocytes) suspension using time-resolved two-colour flow cytometry [37]. Briefly, mononuclear cells were labelled with a phycoerythrin (PE)-conjugated monoclonal antibody to CD3 and an allophycocyanin (APC)-conjugated monoclonal antibody to CD14 (Dako). Ethidium+ (25 μM) was added to a cell suspension maintained at 37°C with stirring in a time-zero module, and 1 mM ATP was added 40 s later. The cells were analysed at 1000 events/s on a FACSCalibur flow cytometer and gated by forward and side scatter and by cell type-specific monoclonal antibody. The linear mean channel of fluorescence intensity for each gated subpopulation over successive 10-s intervals was analysed by WinMDI version 2.8 software [38] and plotted against time to estimate the maximum rate of ethidium+ uptake. The linear portion of the ethidium+ uptake plot was used for slope calculations over a period of 1 minute.

Statistical analysis

Differences in allele frequencies between patients and controls were analysed by logistic regression using the default additive model for allele coding and odds ratios (OR) were derived from exponentiation of the regression coefficients. Epistasis was analysed using a two-factor interaction logistic regression model, with the significance of the epistasis determined by the significance of the coefficient for the interaction term. Analyses for cohort 2 were ancestry-adjusted by the inclusion of the principal components as covariates. Serum IL-18 levels were log-normally distributed and were therefore analysed by lognormal regression and data presented as medians. These analyses were performed using Statistica v6.1 (StatSoft Inc, Tulsa, Oklahoma, USA). Haplotype inference was performed using Phase version 2 software [39].

Results

Functional P2RX7 polymorphism in pSS patients and population based controls

The minor allele frequencies for each of the 12 P2RX7 functional SNPs in pSS patients and controls from cohort 1 are presented in Additional file 1, Table S2. The minor alleles of five SNPs (G151T, G474A, G853A, G946A and T1729A) were rare (frequency <5%). There was no evidence of an association between any of the 12 functional SNPs and pSS. However, the frequency of the minor allele of the A1405G SNP (rs2230912) was quite discrepant between seronegative and seropositive pSS patients in cohort 1, and was in fact increased in seropositive pSS patients compared to controls (OR = 1.78, 95% CI 1.09, 2.92, P = 0.02) (Table 2). However this association was not confirmed in cohort 2.

Table 2. Frequency of the minor allele of P2RX7 A1405 G allele in pSS patients and controls

The minor allele of P2RX7 A1405G exhibits a negative epistatic interaction with HLA DR3 in Ro±La seropositive pSS patients

We further investigated the role of the A1405G allele in anti-Ro±La seropositive pSS by considering epistasis with HLA DR3. The HLA region confers a strong genetic risk for seropositive pSS, and DRB1 alleles are markers of this risk. In Caucasian populations, the strongest association is observed with HLA DR3, although other DRB1 alleles are also involved, and the associated alleles are known to differ across different racial groups.

The joint distribution of HLA DR3 and the P2RX7 1405G allele was examined by a two-factor interaction logistic regression model in both cohorts. Differences in the available data in relation to ascertainment of DR3 status necessitated several analyses, and the results of the regression coefficients for these analyses are depicted in Table 3.

Table 3. Epistasis between P2RX7 1405G and HLA DR3 in anti Ro±La seropositive primary Sjögren's syndrome (pSS)

These analyses confirm that HLA DR3 is a strong risk factor for seropositive, but not seronegative pSS. However this association with DR3 appears somewhat greater in cohort 1 compared to cohort 2, which used a proxy SNP for DR3 status. For example, the regression coefficients (log-OR) reported in Table 3 equate to OR of 9.3 (95% CI 4.5, 19.2) for cohort 1, compared to 2.6 (95% CI 2.6, 5.8) for cohort 2. Further, there is remarkably consistent evidence of a statistically significant negative epistatic interaction between HLA DR3 (or proxy) and the P2RX7 1405G allele in seropositive pSS patients, whether comparing these patients to controls (in both cohorts), or seropositive pSS to seronegative pSS patients (sufficient data for cohort 2 only) (Table 3). This negative epistasis, which is assessed by β3 coefficient in Table 3, implies that the proportion of seropositive pSS patients who carry both alleles is significantly less than expected. However, the association between the P2RX7 1405G allele and seropositive pSS patients who do not carry DR3 (the β2 coefficient in Table 3) observed in cohort 1 (OR 3.0, 95% CI 1.4, 6.4) was not replicated in cohort 2 (OR 1.3, 95% CI 0.8, 2.0). This epistasis effect with DR3 was specific for the P2RX7 1405G allele, and was not observed with any other P2RX7 SNPs genotyped in cohort 1 (data not shown).

The most plausible interpretation of these findings is that the P2RX7 1405G allele is a risk factor for seropositive pSS in a specific subgroup of patients who do not carry HLA risk alleles, but confers no additional risk in individuals who do. In our analysis, the use of HLA DR3 would only partially adjust for the HLA risk, as other DRB1 alleles are also involved [3]. Given that the OR for the association with DR3 are greater in cohort 1 than in cohort 2 (Table 3), this implies that other DRB1 alleles may be relatively more important in cohort 2, yet the analysis could not adjust for this. This difference in HLA risk adjustment between cohorts may be sufficient to explain the lack of replication of the P2RX7 1405G allele association in cohort 2 (Table 3).

The 1405G allele is a tag for a common P2RX7 haplotype

There was evidence of strong linkage disequilibrium between the P2RX7 SNPs. Haplotype reconstruction of pSS patients in cohort 1 and controls identified six major haplotypes (individual frequencies greater than 5%), with a combined frequency of 70% (Figure 1). These haplotypes and frequencies are comparable to those identified in a large Australian Caucasian cohort [30,31]. Importantly, the 1405G minor allele was observed on only one major haplotype (haplotype 4, frequency 13%), and this minor allele is therefore a tag for haplotype-4 which also carries minor alleles for both C489T and G1068A.

thumbnailFigure 1. Common (>5% frequency) P2RX7 haplotypes identified in the study population. Grey squares represent the major alleles and black squares the minor alleles, for each single nucleotide polymorphism. Numbering of polymorphisms are based on the original mRNA sequence [Y09561.1; GenBank]. Six major haplotypes were identified with a combined frequency of 70%.

The 1405G allele is associated with gain-of-function of the P2X7 receptor

P2X7 receptor function associated with the A1405G polymorphism was studied in monocytes and lymphocytes from peripheral blood of healthy individuals by measuring membrane pore formation (assessed by ethidium uptake), following receptor activation by ATP. Monocytes from individuals homozygous or heterozygous for the 1405G allele exhibited an increase in ATP-induced ethidium uptake compared with monocytes from homozygous 1405AA individuals (P = 0.018, n = 5 to 18 individuals) (Figure 2A, C), indicative of a gain in receptor function associated with the 1405G allele. Similar results were obtained with T-lymphocytes, where cells from 1405AG heterozygous individuals, exhibited an increased ATP-induced ethidium uptake (P = 0.015, n = 5 to 16 subjects) (Figure 2B, D). These results confirm previous reports that the 1405G allele is associated with gain of P2X7 receptor function [40].

thumbnailFigure 2. Human peripheral blood mononuclear cells were analysed for ATP-induced ethidium uptake using time-resolved flow cytometry. Monocytes were identified by an anti-CD14 APC antibody and T lymphocytes were identified by an anti-CD3-PE antibody. Ethidium uptake (25 μM) was measured in response to P2X7 receptor stimulation by 1mM ATP. (A) and (C) are representative ethidium uptake curves from genotyped individuals. (B) and (D) are scattergrams of collated data from 5 to 18 subjects from each genotype group (number of subjects shown in brackets) with ATP-induced ethidium uptake calculated from the slope of the linear portion of the uptake plot over 1 minute. *P < 0.05 from one-way analysis of variance with Dunnett's post hoc test.

Serum IL-18 levels are increased in pSS patients with Ro/La autoantibodies

IL-18 is a downstream cytokine released following P2X7 receptor activation, and serum IL-18 levels were measured in 72 pSS patients (54 seropositive for anti-Ro±La autoantibodies) and 36 age- and gender-matched controls from cohort 1 (Figure 3). When analysed by log-normal regression, mean serum IL-18 levels (in pg/mL) were 339, 318, and 498 for the controls, seronegative pSS and seropositive pSS patients, respectively. Serum IL-18 levels were almost 50% higher in Ro/La autoantibody-positive pSS patients compared to normal subjects (ratio = 1.47, 95% CI 1.22, 4.26, P < 0.001), consistent with previous reports. In contrast, the serum IL-18 levels in seronegative pSS patients were not different to that observed in normal controls (ratio = 1.12, 95% CI = 0.87, 1.45, P = 0.4). Among the seropositive pSS patients, there was no difference between the median IL-18 level in patients who carried the 1405G allele (549 pg/mL, n = 25) compared to patients with 1405A (460 pg/mL, n = 29, P = 0.13). There was no evidence of any influence of HLA DR3 status on IL-18 levels in these patients.

thumbnailFigure 3. Boxplots of serum IL-18 levels (pg/mL) in controls (n = 36), autoantibody negative pSS (n = 18), and autoantibody positive primary Sjögren's syndrome (pSS) (n = 54). Serum IL-18 levels were significantly different between the three groups (P = 0.0003, Kruskal-Wallis test), and highest in autoantibody-positive pSS.

Discussion

Autoantibodies, which target intracellular antigens, are a feature of systemic autoimmune diseases such as pSS. These autoantibodies react with autoantigen exposed on apoptotic or necrotic cells, or their debris, and resulting immune complexes promote chronic systemic inflammation via upregulation of the type I interferon system [2]. Delayed or defective clearance of apoptotic/necrotic cellular debris is understood to be an important factor in the breakdown of self-tolerance and ongoing immune complex mediated inflammation [41]. Ro±La autoantibodies are typically present in the majority of pSS patients (70 to 82% in the present study) and are associated with higher levels of glandular inflammation, and extraglandular features such as vasculitis, interstitial lung disease and lymphoma [42]. There is increasing evidence that autoantibody-positive pSS has different genetic risk factors to autoantibody-negative pSS [3,4], which possibly reflects different mechanisms of disease pathogenesis.

In this study we examined the interaction of the pro-inflammatory P2X7 receptor in the pathogenesis of Ro±La autoantibody-positive pSS (seropositive pSS) in two cohorts. The primary focus was on the P2RX7 1405G allele, which is a tag for a common haplotype conferring gain-of-function on the P2X7 receptor [30,31]. This haplotype also carries minor alleles from two other SNPS, 489T and 1068A (Table 1, Figure 1), and analysis of the functional effects of each in isolation suggests that they contribute, in an additive way, to the increased receptor function observed with this haplotype [30,43]. We have demonstrated an interaction between this gain of function SNP/haplotype and seropositive pSS, whereas no associations were observed with loss of function SNPs/haplotypes. Previous studies of P2RX7 polymorphism and systemic autoimmune disease have reported no association between the A1513C loss-of-function SNP in Caucasian patients with SLE or RA [44,45], which is consistent with the results of the present study.

We observed a remarkably consistent negative epistatic interaction between P2RX7 1405G and HLA DR3, the primary genetic risk factor for pSS, in two cohorts of seropositive pSS patients. This negative epistatic interaction was present when seropositive pSS patients were either compared to controls, or to seronegative pSS patients, so that the combined risk for individuals who carry both genetic factors is somewhat less than expected. In cohort 1 we also observed that the P2RX7 1405G allele is a risk factor for seropositive pSS in individuals who do not carry HLA DR3, but this effect was not replicated in cohort 2. The analysis in cohort 2 included principal components (derived from unrelated genetic markers) as covariates to adjust for possible population stratification, which can cause spurious associations in disease studies [36]. Therefore, the observed epistasis cannot be attributed to systematic ancestry differences between cases and controls.

The most plausible interpretation of these findings is that the P2RX7 1405G allele is a risk factor for seropositive pSS in a specific subgroup of patients who do not carry HLA risk alleles, but confers no additional risk in individuals who do, and that differences in results between the cohorts may reflect differences in adjustment for the HLA associated risk, which is only partially reflected by DR3 status. For example, it is unclear whether the use of a surrogate SNP for determining HLA-DR3 status in the cohort 2 controls may have influenced the results by underestimating the prevalence of HLA-DR3; further, HLA-DR2 is also associated with autoantibody-positive pSS [3], but this information was not available for cohort 2.

The epistatic interaction between DR3 and the P2RX7 1405G allele in seropositive pSS implies overlap in pathogenic mechanisms associated with these alleles. Unfortunately, the pathogenic mechanisms underlying HLA-associated disease susceptibility are not understood, therefore, the interpretation is highly speculative. We propose that the observed epistasis between DR3 and the P2RX7 1405G allele may reflect increased autoantigen exposure as a susceptibility mechanism for seropositive pSS.

Complement deficiency is implicated in the defective clearance of apoptotic/necrotic debris and, therefore, increased autoantigen exposure, in systemic autoimmune disease [46]. It has been proposed that this mechanism may contribute to some of the pSS disease susceptibility associated with HLA DR3, which is in linkage disequilibrium with a C4 null allele [47]. P2X7 receptor activation may also contribute to increased autoantigen exposure in systemic autoimmunity. In the presence of the ATP ligand, P2X7 receptor activation may result in the induction of cell membrane blebbing and release of microparticles. Prolonged activation of the P2X7 receptor by ATP results in cell death with both apoptotic/necrotic features [13,14], which has also been demonstrated in rat parotid acinar ParC5 cells [48], accompanied by autoantigen cleavage [49]. Importantly, the A1513C loss of function SNP is associated with a decrease in ATP-induced apoptosis in CD4(+) lymphocytes isolated from SLE patients [45], although the gain-of-function polymorphism has not yet been evaluated. Membrane microparticles released by the P2X7-induced blebbing process contain not only P2X7 receptors but also Ro-52, a component of the P2X7 membrane complex both in THP-1 macrophages and transfected HEK-293 epithelial cells [10], and also a target of autoantibodies in pSS. Further, we have recently demonstrated that P2X7, in the absence of ATP, acts as a scavenging receptor for apoptotic cellular debris [50], although the downstream effects of this have not been elucidated.

In addition to the genetic epistasis, our data suggest that there is P2X7 receptor activation in seropositive pSS. Activation of P2X7 is known to lead to processing and secretion of pro-inflammatory cytokines IL-1β and IL-18 from monocyte/macrophages via activation of the NALP3 inflammasome, which is thought to play a role in a spectrum of inflammatory diseases [51], including rheumatic inflammatory diseases [52]. We have demonstrated elevated serum IL-18 levels in Ro±La autoantibody-seropositive pSS patients, and activation of the NALP3 inflammasome, via P2X7, may be one of the mechanisms involved. There was no evidence for an effect of DR3 or 1405G status on IL-18 levels, although patient numbers were limited. Previous studies in systemic inflammatory diseases have demonstrated that the ATP-induced monocyte P2X7 receptor response (as measured by IL-1β secretion) is attenuated in SLE, but enhanced in RA [45,53], and that this response is influenced by P2RX7 polymorphism. However, the A1405G SNP has not been evaluated in this context.

Aside from its pro-inflammatory effects, ATP activation of the P2X7 receptor has also emerged as an important regulator of autonomic function, including salivary gland secretion. While P2X7 receptor activation can stimulate secretion in mouse submandibular bands, activation also results in the inhibition of muscarinic receptor-induced fluid secretion [19]. Therefore, activation of P2X7 in the milieu of salivary gland inflammation in pSS may contribute to the failure of fluid secretion irrespective of the extent of glandular destruction, through interference of autonomic signalling, a well-documented phenomenon in pSS. This dysregulation of autonomic function may also be implicated in the clustering of autonomic symptoms we have previously reported in pSS [54], and further research into the role of P2X7 and autonomic dysfunction in pSS is warranted.

Conclusions

In summary, we have identified that the P2RX7 1405G gain-of-function allele exhibits a negative epistatic interaction with HLA DR3 on the risk of autoantibody-positive pSS. This, coupled with other studies, implicates P2X7 receptor function in the pathogenesis of primary Sjogren's syndrome.

Abbreviations

APC: allophycocyanin; ELISA: enzyme-linked immunosorbent assay; IL: interleukin; NALP3: cryopyrin; OMRF: Oklahoma Medical Research Foundation; OR: odds ratio; PE: phycoerythrin; pSS: primary Sjögren's syndrome; RA: rheumatoid arthritis; RNP: ribonucleoprotein; SLE: systemic lupus erythematosus; SNP: single nucleotide polymorphism.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

SL: study design, genotyping, data analysis, interpretation, manuscript preparation; LS: ethidium uptake assays, analysis and manuscript preparation; KS: genotyping, IL-18 assays, manuscript preparation; BG: study design, manuscript preparation; KS: genotyping, interpretation and manuscript preparation; CL: genotyping, interpretation and manuscript preparation; JW: study concept and design, interpretation and manuscript preparation; MR: study concept and design, interpretation and manuscript preparation. All authors have read and approved the final manuscript.

Acknowledgements

The authors gratefully acknowledge the support of the patients who participated in this study.

This study was funded by support from Arthritis Australia, National Health and Medical Research Council (NHMRC, grant number 509184), and The Queen Elizabeth Hospital Research Foundation

References

  1. Vitali C, Bombardieri S, Jonsson R, Moutsopoulos HM, Alexander EL, Carsons SE, Daniels TE, Fox PC, Fox RI, Kassan SS, Pillemer SR, Talal N, Weisman MH: Classification criteria for Sjogren's syndrome: a revised version of the European criteria proposed by the American-European Consensus Group.

    Ann Rheum Dis 2002, 61:554-558. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  2. Nordmark G, Alm GV, Ronnblom L: Mechanisms of Disease: primary Sjogren's syndrome and the type I interferon system.

    Nat Clin Pract Rheumatol 2006, 2:262-269. PubMed Abstract | Publisher Full Text OpenURL

  3. Rischmueller M, Lester S, Chen Z, Champion G, Van Den Berg R, Beer R, Coates T, McCluskey J, Gordon T: HLA class II phenotype controls diversification of the autoantibody response in primary Sjogren's syndrome (pSS).

    Clin Exp Immunol 1998, 111:365-371. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Lester S, McLure C, Williamson J, Bardy P, Rischmueller M, Dawkins RL: Epistasis between the MHC and the RCA alpha block in primary Sjogren syndrome.

    Ann Rheum Dis 2008, 67:849-854. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. Cruz-Tapias P, Rojas-Villarraga A, Maier-Moore S, Anaya JM: HLA and Sjogren's syndrome susceptibility. A meta-analysis of worldwide studies.

    Autoimmun Rev 2012, 11:281-287. PubMed Abstract | Publisher Full Text OpenURL

  6. Elliott JI, McVey JH, Higgins CF: The P2X7 receptor is a candidate product of murine and human lupus susceptibility loci: a hypothesis and comparison of murine allelic products.

    Arthritis Res Ther 2005, 7:R468-475. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  7. Chen L, Brosnan CF: Regulation of immune response by P2X7 receptor.

    Crit Rev Immunol 2006, 26:499-513. PubMed Abstract | Publisher Full Text OpenURL

  8. Labasi JM, Petrushova N, Donovan C, McCurdy S, Lira P, Payette MM, Brissette W, Wicks JR, Audoly L, Gabel CA: Absence of the P2X7 receptor alters leukocyte function and attenuates an inflammatory response.

    J Immunol 2002, 168:6436-6445. PubMed Abstract | Publisher Full Text OpenURL

  9. Kawamura H, Aswad F, Minagawa M, Govindarajan S, Dennert G: P2X7 receptors regulate NKT cells in autoimmune hepatitis.

    J Immunol 2006, 176:2152-2160. PubMed Abstract | Publisher Full Text OpenURL

  10. Gu BJ, Rathsam C, Stokes L, McGeachie AB, Wiley JS: Extracellular ATP dissociates nonmuscle myosin from P2X(7) complex: this dissociation regulates P2X(7) pore formation.

    Am J Physiol Cell Physiol 2009, 297:C430-439. PubMed Abstract | Publisher Full Text OpenURL

  11. Kim M, Jiang LH, Wilson HL, North RA, Surprenant A: Proteomic and functional evidence for a P2X7 receptor signalling complex.

    Embo J 2001, 20:6347-6358. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Pfeiffer ZA, Aga M, Prabhu U, Watters JJ, Hall DJ, Bertics PJ: The nucleotide receptor P2X7 mediates actin reorganization and membrane blebbing in RAW 264.7 macrophages via p38 MAP kinase and Rho.

    J Leukoc Biol 2004, 75:1173-1182. PubMed Abstract | Publisher Full Text OpenURL

  13. Khakh BS, North RA: P2X receptors as cell-surface ATP sensors in health and disease.

    Nature 2006, 442:527-532. PubMed Abstract | Publisher Full Text OpenURL

  14. Taylor SR, Gonzalez-Begne M, Dewhurst S, Chimini G, Higgins CF, Melvin JE, Elliott JI: Sequential shrinkage and swelling underlie P2X7-stimulated lymphocyte phosphatidylserine exposure and death.

    J Immunol 2008, 180:300-308. PubMed Abstract | Publisher Full Text OpenURL

  15. Gourine AV, Wood JD, Burnstock G: Purinergic signalling in autonomic control.

    Trends Neurosci 2009, 32:241-248. PubMed Abstract | Publisher Full Text OpenURL

  16. Lister MF, Sharkey J, Sawatzky DA, Hodgkiss JP, Davidson DJ, Rossi AG, Finlayson K: The role of the purinergic P2X7 receptor in inflammation.

    J Inflamm (Lond) 2007, 4:5. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  17. Alzola E, Perez-Etxebarria A, Kabre E, Fogarty DJ, Metioui M, Chaib N, Macarulla JM, Matute C, Dehaye JP, Marino A: Activation by P2X7 agonists of two phospholipases A2 (PLA2) in ductal cells of rat submandibular gland. Coupling of the calcium-independent PLA2 with kallikrein secretion.

    J Biol Chem 1998, 273:30208-30217. PubMed Abstract | Publisher Full Text OpenURL

  18. Pochet S, Metioui M, Grosfils K, Gomez-Munoz A, Marino A, Dehaye JP: Regulation of phospholipase D by muscarinic receptors in rat submandibular ductal cells.

    Cell Signal 2003, 15:103-113. PubMed Abstract | Publisher Full Text OpenURL

  19. Nakamoto T, Brown DA, Catalan MA, Gonzalez-Begne M, Romanenko VG, Melvin JE: Purinergic P2X7 Receptors Mediate ATP-induced Saliva Secretion by the Mouse Submandibular Gland.

    J Biol Chem 2009, 284:4815-4822. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Novak I: Purinergic signalling in epithelial ion transport: regulation of secretion and absorption.

    Acta Physiol (Oxf) 2011, 202:501-522. PubMed Abstract | Publisher Full Text OpenURL

  21. Di Virgilio F: Liaisons dangereuses: P2X(7) and the inflammasome.

    Trends Pharmacol Sci 2007, 28:465-472. PubMed Abstract | Publisher Full Text OpenURL

  22. Ferrari D, Chiozzi P, Falzoni S, Dal Susino M, Melchiorri L, Baricordi OR, Di Virgilio F: Extracellular ATP triggers IL-1 beta release by activating the purinergic P2Z receptor of human macrophages.

    J Immunol 1997, 159:1451-1458. PubMed Abstract | Publisher Full Text OpenURL

  23. Perregaux DG, McNiff P, Laliberte R, Conklyn M, Gabel CA: ATP acts as an agonist to promote stimulus-induced secretion of IL-1 beta and IL-18 in human blood.

    J Immunol 2000, 165:4615-4623. PubMed Abstract | Publisher Full Text OpenURL

  24. Mariathasan S, Weiss DS, Newton K, McBride J, O'Rourke K, Roose-Girma M, Lee WP, Weinrauch Y, Monack DM, Dixit VM: Cryopyrin activates the inflammasome in response to toxins and ATP.

    Nature 2006, 440:228-232. PubMed Abstract | Publisher Full Text OpenURL

  25. Bombardieri M, Barone F, Pittoni V, Alessandri C, Conigliaro P, Blades MC, Priori R, McInnes IB, Valesini G, Pitzalis C: Increased circulating levels and salivary gland expression of interleukin-18 in patients with Sjogren's syndrome: relationship with autoantibody production and lymphoid organization of the periductal inflammatory infiltrate.

    Arthritis Res Ther 2004, 6:R447-456. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  26. Szodoray P, Alex P, Brun JG, Centola M, Jonsson R: Circulating cytokines in primary Sjogren's syndrome determined by a multiplex cytokine array system.

    Scand J Immunol 2004, 59:592-599. PubMed Abstract | Publisher Full Text OpenURL

  27. Bulosan M, Pauley KM, Yo K, Chan EK, Katz J, Peck AB, Cha S: Inflammatory caspases are critical for enhanced cell death in the target tissue of Sjogren's syndrome before disease onset.

    Immunol Cell Biol 2009, 87:81-90. PubMed Abstract | Publisher Full Text OpenURL

  28. Gilboa-Geffen A, Wolf Y, Hanin G, Melamed-Book N, Pick M, Bennett ER, Greenberg DS, Lester S, Rischmueller M, Soreq H: Activation of the Alternative NFkappaB Pathway Improves Disease Symptoms in a Model of Sjogren's Syndrome.

    PLoS One 2011, 6:e28727. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Sakai A, Sugawara Y, Kuroishi T, Sasano T, Sugawara S: Identification of IL-18 and Th17 cells in salivary glands of patients with Sjogren's syndrome, and amplification of IL-17-mediated secretion of inflammatory cytokines from salivary gland cells by IL-18.

    J Immunol 2008, 181:2898-2906. PubMed Abstract | Publisher Full Text OpenURL

  30. Fuller SJ, Stokes L, Skarratt KK, Gu BJ, Wiley JS: Genetics of the P2X7 receptor and human disease.

    Purinergic Signal 2009, 5:257-262. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  31. Stokes L, Fuller SJ, Sluyter R, Skarratt KK, Gu BJ, Wiley JS: Two haplotypes of the P2X(7) receptor containing the Ala-348 to Thr polymorphism exhibit a gain-of-function effect and enhanced interleukin-1beta secretion.

    Faseb J 2010, 24:2916-2927. PubMed Abstract | Publisher Full Text OpenURL

  32. Sluyter R, Dalitz JG, Wiley JS: P2X7 receptor polymorphism impairs extracellular adenosine 5'-triphosphate-induced interleukin-18 release from human monocytes.

    Genes Immun 2004, 5:588-591. PubMed Abstract | Publisher Full Text OpenURL

  33. Sluyter R, Shemon AN, Wiley JS: Glu496 to Ala polymorphism in the P2X7 receptor impairs ATP-induced IL-1 beta release from human monocytes.

    J Immunol 2004, 172:3399-3405. PubMed Abstract | Publisher Full Text OpenURL

  34. Scofield AK, Radfar L, Ice JA, Vista E, Anaya JM, Houston G, Lewis D, Stone DU, Chodosh J, Hefner K, Lessard CJ, Moser KL, Scofield RH: Relation of sensory peripheral neuropathy in Sjogren syndrome to anti-Ro/SSA.

    J Clin Rheumatol 2012, 18:290-293. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. Field J, Browning SR, Johnson LJ, Danoy P, Varney MD, Tait BD, Gandhi KS, Charlesworth JC, Heard RN, Stewart GJ, Kilpatrick TJ, Foote SJ, Bahlo M, Butzkueven H, Wiley J, Booth DR, Taylor BV, Brown MA, Rubio JP, Stankovich J: A polymorphism in the HLA-DPB1 gene is associated with susceptibility to multiple sclerosis.

    PLoS One 2010, 5:e13454. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  36. Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D: Principal components analysis corrects for stratification in genome-wide association studies.

    Nat Genet 2006, 38:904-909. PubMed Abstract | Publisher Full Text OpenURL

  37. Jursik C, Sluyter R, Georgiou JG, Fuller SJ, Wiley JS, Gu BJ: A quantitative method for routine measurement of cell surface P2X7 receptor function in leucocyte subsets by two-colour time-resolved flow cytometry.

    J Immunol Methods 2007, 325:67-77. PubMed Abstract | Publisher Full Text OpenURL

  38. The Scripps Research Institute: Windows Multiple Document Interface for Flow Cytometry (WinMDI), Version 2.8., 2000. [http://facs.scripps.edu/software.html] webcite

  39. Stephens M, Smith NJ, Donnelly P: A new statistical method for haplotype reconstruction from population data. [http://stephenslab.uchicago.edu/software.html] webcite

    Am J Hum Genet 2001, 68:978-989. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  40. Denlinger LC, Coursin DB, Schell K, Angelini G, Green DN, Guadarrama AG, Halsey J, Prabhu U, Hogan KJ, Bertics PJ: Human P2X7 pore function predicts allele linkage disequilibrium.

    Clin Chem 2006, 52:995-1004. PubMed Abstract | Publisher Full Text OpenURL

  41. Mahoney JA, Rosen A: Apoptosis and autoimmunity.

    Curr Opin Immunol 2005, 17:583-588. PubMed Abstract | Publisher Full Text OpenURL

  42. Beer RG, Rischmueller M, Coates T, Purcell AW, Keech CL, McCluskey J, Gordon TP: Nonprecipitating anti-La(SS-B) autoantibodies in primary Sjogren's syndrome.

    Clin Immunol Immunopathol 1996, 79:314-318. PubMed Abstract | Publisher Full Text OpenURL

  43. Portales-Cervantes L, Nino-Moreno P, Salgado-Bustamante M, Garcia-Hernandez MH, Baranda-Candido L, Reynaga-Hernandez E, Barajas-Lopez C, Gonzalez-Amaro R, Portales-Perez DP: The His155Tyr (489C>T) single nucleotide polymorphism of P2RX7 gene confers an enhanced function of P2X7 receptor in immune cells from patients with rheumatoid arthritis.

    Cell Immunol 2012, 276:168-175. PubMed Abstract | Publisher Full Text OpenURL

  44. Forchap SL, Anandacoomarasamy A, Wicks J, Di Virgilio F, Baricordi OR, Rubbini M, Trotta F, Wiley J, Manolios N: P2X7 gene polymorphisms do not appear to be a susceptibility gene locus in sporadic cases of systemic lupus erythematosus.

    Tissue Antigens 2008, 72:487-490. PubMed Abstract | Publisher Full Text OpenURL

  45. Portales-Cervantes L, Nino-Moreno P, Doniz-Padilla L, Baranda-Candido L, Garcia-Hernandez M, Salgado-Bustamante M, Gonzalez-Amaro R, Portales-Perez D: Expression and function of the P2X(7) purinergic receptor in patients with systemic lupus erythematosus and rheumatoid arthritis.

    Hum Immunol 2010, 71:818-825. PubMed Abstract | Publisher Full Text OpenURL

  46. Cook HT, Botto M: Mechanisms of Disease: the complement system and the pathogenesis of systemic lupus erythematosus.

    Nat Clin Pract Rheumatol 2006, 2:330-337. PubMed Abstract | Publisher Full Text OpenURL

  47. Candore G, Lio D, Colonna Romano G, Caruso C: Pathogenesis of autoimmune diseases associated with 8.1 ancestral haplotype: effect of multiple gene interactions.

    Autoimmun Rev 2002, 1:29-35. PubMed Abstract | Publisher Full Text OpenURL

  48. Hwang SM, Koo NY, Choi SY, Chun GS, Kim JS, Park K: P2X7 Receptor-mediated Membrane Blebbing in Salivary Epithelial Cells.

    Korean J Physiol Pharmacol 2009, 13:175-179. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  49. Hwang SM, Li J, Koo NY, Choi SY, Lee SJ, Oh SB, Castro R, Kim JS, Park K: Role of purinergic receptor in alpha fodrin degradation in Par C5 cells.

    J Dent Res 2009, 88:927-932. PubMed Abstract | Publisher Full Text OpenURL

  50. Gu BJ, Saunders BM, Petrou S, Wiley JS: P2X(7) is a scavenger receptor for apoptotic cells in the absence of its ligand, extracellular ATP.

    J Immunol 2011, 187:2365-2375. PubMed Abstract | Publisher Full Text OpenURL

  51. Masters SL, Simon A, Aksentijevich I, Kastner DL: Horror autoinflammaticus: the molecular pathophysiology of autoinflammatory disease (*).

    Annu Rev Immunol 2009, 2:7:621-668. OpenURL

  52. Sidiropoulos PI, Goulielmos G, Voloudakis GK, Petraki E, Boumpas DT: Inflammasomes and rheumatic diseases: evolving concepts.

    Ann Rheum Dis 2008, 67:1382-1389. PubMed Abstract | Publisher Full Text OpenURL

  53. Al-Shukaili A, Al-Kaabi J, Hassan B: A comparative study of interleukin-1beta production and p2x7 expression after ATP stimulation by peripheral blood mononuclear cells isolated from rheumatoid arthritis patients and normal healthy controls.

    Inflammation 2008, 31:84-90. PubMed Abstract | Publisher Full Text OpenURL

  54. Cai FZ, Lester S, Lu T, Keen H, Boundy K, Proudman SM, Tonkin A, Rischmueller M: Mild autonomic dysfunction in primary Sjogren's syndrome: a controlled study.

    Arthritis Res Ther 2008, 10:R31. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL