Skip to main content
  • Research article
  • Open access
  • Published:

Analysis of Fcγ receptor haplotypes in rheumatoid arthritis: FCGR3A remains a major susceptibility gene at this locus, with an additional contribution from FCGR3B

Abstract

The Fcγ receptors play important roles in the initiation and regulation of many immunological and inflammatory processes, and genetic variants (FCGR) have been associated with numerous autoimmune and infectious diseases. The data in rheumatoid arthritis (RA) are conflicting and we previously demonstrated an association between FCGR3A and RA. In view of the close molecular proximity with FCGR2A, FCGR2B and FCGR3B, additional polymorphisms within these genes and FCGR haplotypes were examined to refine the extent of association with RA. Biallelic polymorphisms in FCGR2A, FCGR2B and FCGR3B were examined for association with RA in two well characterized UK Caucasian and North Indian/Pakistani cohorts, in which FCGR3A genotyping had previously been undertaken. Haplotype frequencies and linkage disequilibrium were estimated across the FCGR locus and a model-free analysis was performed to determine association with RA. This was followed by regression analysis, allowing for phase uncertainty, to identify the particular haplotype(s) that influences disease risk. Our results reveal that FCGR2A, FCGR2B and FCGR3B were not associated with RA. The haplotype with the strongest association with RA susceptibility was the FCGR3A–FCGR3B 158V-NA2 haplotype (odds ratio 3.18, 95% confidence interval 1.13–8.92 [P = 0.03] for homozygotes compared with all genotypes). The association was stronger in the presence of nodules (odds ratio 5.03, 95% confidence interval 1.44–17.56; P = 0.01). This haplotype was also more common in North Indian/Pakistani RA patients than in control individuals, but not significantly so. Logistic regression analyses suggested that FCGR3A remained the most significant gene at this locus. The increased association with an FCGR3A–FCGR3B haplotype suggests that other polymorphic variants within FCGR3A or FCGR3B, or in linkage disequilibrium with this haplotype, may additionally contribute to disease pathogenesis.

Introduction

Rheumatoid arthritis (RA) is a heterogeneous disease characterised by a chronic, fluctuating, peripheral, symmetrical and erosive polyarthritis. It has been reported throughout the world, with a prevalence rate of approximately 1% in most populations [1]. Persistent synovial inflammation leads to progressive joint destruction, which in turn produces deformity and significant disability. In addition, RA is a systemic disease and some patients develop subcutaneous rheumatoid nodules, secondary Sjögren's syndrome, episcleritis and scleritis, interstitial lung disease, pericardial involvement, systemic vasculitis and Felty's syndrome. These extra-articular manifestations of RA appear to be rare in the absence of rheumatoid factor (RF), and IgG RF titres correlate with articular disease severity and with the extra-articular manifestations of RA [2]. In addition to RF, many other IgG autoantibodies are found in RA, most notably anti-cyclic citrullinated peptide [3] and anti-type II collagen antibodies [4].

The Fcγ receptors (FcγRs), which bind these IgG autoantibodies and IgG-containing immune complexes, have been shown to play important roles in the initiation and regulation of many immunological and inflammatory processes. In humans, there are three classes of FcγRs (I-III) encoded by eight genes, which produce at least 15 different membrane-bound and soluble isoforms that vary in their cellular distribution and affinity for different IgG isotypes. This molecular and expression diversity restricts specific biological properties to certain cell types. Activating (FcγRIIa, FcγRIIIa and FcγRIIIb) and inhibitory (FcγRIIb) FcγRs are frequently coexpressed on the same cell, thus providing a means for regulating signalling thresholds [5, 6]. Furthermore, the absolute level of receptor expression is modulated by proinflammatory and anti-inflammatory cytokines [7]. Activating functions include uptake and clearance of immune complexes (complement dependent and independent mechanisms), activation of phagocytes (trigger the oxidative burst, cytotoxic granule and cytokine release), antigen presentation and antibody-dependent cellular cytotoxicity [6]. Conversely, FcγRIIb contains an inhibitory motif in the cytoplasmic tail and abrogates cellular activation. FcγRIIb may also play a role in maintaining peripheral B cell tolerance and prevention of autoimmunity [8]. Single and multiple FcγR knockout mouse models have demonstrated that the balance between activating and inhibitory FcγRs influences the development of both immune complex-mediated and collagen-induced arthritis [9, 10].

Polymorphic variants that increase the expression or affinity of these IgG receptors or that enhance their ability to bind specific IgG isotypes may therefore play an important role in determining the severity and persistence of inflammation to IgG (auto)antibodies and immune complexes in RA. Our previous studies have supported an association between FCGR3A and RA. The higher affinity FCGR3A-158V allele was associated with an increased susceptibility to RA, with homozygotes demonstrating a 1.5-fold to twofold increased risk for RA and a twofold to fourfold increase in nodules [11, 12]. In keeping with several other RA susceptibility genes, this association has not been replicated in all populations [1318]. There was a trend toward an increased frequency of the FCGR3A-158V allele in Norwegian and Dutch RA patients [16, 18], with a skewing toward the FCGR3A-158F allele in Spanish, Japanese and Indian populations [1315, 17]. Many reasons for the lack of reproducibility of association studies have been proposed and include differences in the design, power and accuracy of genotyping strategies in the various studies. The apparently conflicting results may also be a consequence of a genuine difference in the genetic and environmental susceptibility factors according to the precise ethnic group or disease phenotype under investigation. Nonreplication may thus indicate a real biological difference between populations. Alternatively, if FCGR3A is in linkage disequilibrium with the true RA-susceptibility locus then, because the extent of linkage disequilibrium varies between different populations, the association may only exist in certain populations and contribute to nonreplication of findings [12, 19].

FCGR3A lies in a 200 kilobase FCGR gene cluster on 1q22-23 in close molecular proximity to FCGR2A, FCGR3B and FCGR2B [20]. Functional single nucleotide polymorphisms (SNPs) in the latter genes have been investigated as genetic susceptibility and severity factors in multiple infectious and autoimmune diseases [5, 6]. We therefore undertook additional genotyping in our original RA cohorts and examined SNPs in FCGR2A, FCGR2B and FCGR3B; examined the extent of linkage disequilibrium at this locus; and analyzed FCGR haplotypes for association with disease in order to investigate the possibility that there are other RA susceptibility variants at this locus. We demonstrate an increased association with a FCGR3A–FCGR3B haplotype, which suggests that other polymorphic variants within FCGR3A or FCGR3B, or in linkage disequilibrium with this haplotype may additionally contribute to disease pathogenesis.

Materials and methods

Rheumatoid arthritis patients and control individuals

This was an allelic association study conducted to examine FCGR2A, FCGR2B and FCGR3B, and FCGR haplotypes in two well characterized RA cohorts in which an association between FCGR3A and RA was previously identified [11]. The recruitment and clinical characteristics of these two RA and control populations, resident in Birmingham, UK have previously been described [11, 21]. They comprise 294 UK Caucasian individuals (150 RA patients and 144 healthy control individuals) and 256 North Indian/Pakistani individuals (126 RA patients and 130 healthy control individuals). Ethical approval was obtained from the respective local research ethics committees.

Elucidation of the FCGRgene order

Computational assemblies of 1q23 at the National Center for BioInformatics [22], Ensembl [23] and Oak Ridge National Laboratory [24] websites resulted in several variations in FCGR gene order, thus necessitating the construction of an electronic contig. Genomic exon sequences of the class II (FCGR2A, FCGR2B and FCGR2C) and class III (FCGR3A and FCGR3B) genes were aligned to enable identification of homologous regions and specific nucleotides that distinguished the FCGR genes. All available sequence data containing the class II and III FCGR genes was identified by performing BLAST (basic local alignment search tool) sequence homology searches at the National Centre for BioInformatics [22] using five 30-base-pair homologous regions from each receptor class. The bacterial artificial chromosome (BAC) sequence fragments were primarily aligned by the identification of complete gene sequences and overlapping sequence data from neighbouring BACs, and the final assembly was facilitated by utilizing published restriction enzyme maps of this locus [20, 25].

FCGR2Bsequence analysis

At the time when this study was performed, several putative polymorphic variants of FCGR2B had been identified in a cDNA library [26], but these had not been substantiated by sequencing genomic DNA. These included a 2-base-pair deletion at positions 208–210 at the start of exon 3, a G→A substitution at position 685 in exon 4, a functionally significant T→G substitution at position 855 in exon 6 [27] and a G→A substitution at position 1206 in the 3'-untranslated region (UTR). Two potential polymorphic sites (rs844 and rs1043) in STS accession G06355 (UniSTS:73835) were also identified from the reference SNP database, which corresponded to the 3'-UTR of FCGR2B.

The genomic sequence alignments of the class II FCGR genes (FCGR2A, FCGR2C and FCGR2B), generated as described above, were used to design FCGR2B-specific PCRs to facilitate direct sequencing of exons 3 and 6 and the 3'-UTR. The primer sequences and annealing temperatures of the different PCR reactions are shown in Table 1. Briefly, 20 μl PCRs were performed using 100 ng DNA, 200 nmol/l of each primer, 40 μmol/l each of 4 dNTPs, 1.5 mmol/l MgCl2 and 0.5 units of Taq DNA polymerase (Promega, Southampton, UK). The PCR reaction was performed in 30 individuals using a Techne Genius PCR machine (Techne [Cambridge] Ltd, Ducksford, Cambridge, UK) and the PCR conditions were 95°C for 5 minutes followed by 38 cycles of 95°C for 30 s, annealing temperature for 60 s and 72°C for 60 s, with a final extension step of 72°C for 10 minutes. Fluorescent automated cycle sequencing of the PCR products was performed using a dRhodamine terminator reaction kit (PE Biosystems, Warrington, UK). Electrophoresis was performed on polyacrylamide gels using the ABI PRISM® 377 DNA Sequencer (Applied Biosystems, Foster City, CA, USA) and the sequence analyzed utilizing ABI PRISM® 377 sequencing software.

Table 1 Primer pairs used to amplify specific FCGR sequencing templates

FCGRgenotyping

FCGR2A

The FCGR2A-131H/R functional polymorphism has a G→A substitution at nucleotide 519, which results in a switch from arginine (R) to histidine (H) at amino acid position 131 in the immunoglobulin-binding domain [28]. Genotyping was performed using a nested amplification refractory mutation system (ARMS) PCR approach. A 322 bp PCR product was amplified (30 cycles) using a combination of previously published FCGR2A-specific primer sequences (PCR1 and 4INM) [28]. A 1:500 dilution served as a template for two separate nested ARMS PCRs that utilized 400 nmol/l of the published IIA-R and IIA-H primers [29]. The 322 bp product served as a positive control and a 246 bp product indicated the presence of either the FCGR2A-131H or R allele, according to the ARMS primer used. Specific amplification of FCGR2A, rather than the highly homologous FCGR2B and FCGR2C, and the +519 SNP, was confirmed in 40 individuals by direct sequencing.

FCGR2B

The FCGR2B-1206G/A polymorphism was genotyped using a nested RFLP assay. The 330 bp FCGR2B-specific sequencing PCR product (see above) was diluted 1:200 and served as a template for a second PCR that used IIB-UTRSF and a mutated reverse primer (IIB-UTRGR), which introduced an allele-specific HaeIII restriction site. The resultant PCR product was incubated at 37°C for 1 hour with 6U HaeIII (Promega) and the products visualized using a 3.5% agarose gel.

FCGR3A

Additional DNA samples that had not yielded reliable results on direct sequencing [11] were genotyped using our single-stranded conformational polymorphism assay [12] and were included in the haplotype analysis.

FCGR3B

The functional neutrophil antigen (NA)1 and NA2 alleles of FCGR3B were genotyped using minor modifications to a previously published ARMS PCR assay, which included both allele-specific primers and a distinct internal control [30]. The NA1 assay included the two ARMS primers IIIB-NA1F and IIIB-NA1R, and an internal control (fragment of FCGR3A and FCGR3B) was amplified by IIIB-NA1PF and IIIB-NA1PR. The NA2 assay similarly included the IIIB-NA2F and IIIB-NA2R ARMS primers and internal control primers IIIB-NA2PR and IIIB-NA2PR from the gene MCSD1.

Statistical analyses

Statistical analyses were performed using the Stata statistical software (Stata Statistical Software, release 8.0; Stata Corporation, College Station, TX, USA) unless otherwise stated. Hardy–Weinberg equilibrium was investigated in each control population using a goodness-of-fit test to check whether the observed pattern of genotype frequencies was consistent with expectations. Allele and genotype frequencies were compared using 2 × 2 and 3 × 2 contingency tables, respectively. Nodules are only rarely present during the early stages of RA and their absence does not indicate that they may not develop in the future. The control population was therefore felt to be the most appropriate reference group for analysis of the subgroup with nodular RA [31].

Haplotype frequencies were estimated pair-wise across the FCGR locus using the Estimating Haplotypes PLUS (EHPLUS) program [32]. A pair-wise measure of linkage disequilibrium (D') was also calculated for each pair of FCGR genes. Association with disease was tested for by comparing the haplotype frequencies estimated from cases and controls separately with estimates based on the combined sample, using a likelihood ratio test. A permutation procedure implemented in the EHPLUS program was used to assess statistical significance based on 1,000 permutations [32].

If an individual is heterozygous at two loci, then the phase (for example, which variants are inherited from the same parent) is unknown. Association of FCGR3A–FCGR3B haplotypes with RA was investigated further using the haplotype trend regression (HTR) approach proposed by Zaykin and coworkers [33] for dealing with uncertain phase. In this method logistic regression can be used to predict disease status from an individual's haplotypes; where these are not known with certainty, all haplotypes consistent with the genotypes are included as predictors, weighted by their probabilities. This approach estimates the effect on risk for each haplotype, assuming that each of the individual's two haplotypes can have an independent effect. Stepwise regression analyses were also used to investigate the joint effect of the FCGR genes [34].

Results

Chromosomal order of the FCGRgenes

The FCGR genes were located on three BAC clones: RP11-474I16 (EMBL: AL359541) contained FCGR2B; RP11-25I17 (EMBL: AC021370) contained FCGR2B, FCGR2C, FCGR3A and FCGR3B; and the final clone RP11-5K23 (EMBL: AC013307) contained FCGR3A and FCGR2A. The complete genomic sequence of each FCGR (with the exception of the 5' part of FCGR2C) was identified on these BAC clones. Alignment of the sequence fragments demonstrated the FCGR gene order from centromere to telomere at chromosome 1q23 as FCGR2A, FCGR3A, FCGR2C, FCGR3B, FCGR2B. Similar work is in agreement with this gene order [35].

FCGR2Bsequencing

The rs844 G→A substitution was confirmed and was identical to that described previously at position 1206 [26], and this SNP was therefore designated FCGR2B-1206G/A. No further polymorphisms were identified in exons 3 or 6 or the 3'-UTR of FCGR2B in 30 control individuals and RA patients.

Association of FCGR2A, FCGR3B and FCGR2Bwith rheumatoid arthritis

Genotyping was complete on 274 Caucasian individuals (147 cases and 127 controls) and 249 North Indian/Pakistani individuals (122 cases and 127 controls). FCGR2A, FCGR3A and FCGR2B were in Hardy–Weinberg equilibrium in both control groups. FCGR3B was not in Hardy–Weinberg equilibrium (for the UK Caucasian group: P = 0.01; for the North Indian/Pakistani group: P = 0.002). We subsequently sequenced more than 200 individuals with 100% agreement with our genotyping assays to exclude genotyping error as an explanation.

No significant differences in the allele or genotype distributions were seen for FCGR2A, FCGR3B, or FCGR2B in either RA group compared with controls (Table 2). The results for FCGR3A in our expanded UK Caucasian cohort were consistent with our previous findings [11]. Homozygosity for the FCGR3A-158V allele demonstrated a trend toward and association with RA (odds ratio [OR] 2.1, 95% confidence interval [CI] 1.0–4.7; P = 0.06) and significant association with nodular RA (OR 4.3, 95% CI 1.5–12.3; P = 0.005).

Table 2 Genotype frequencies in UK Caucasian and North Indian/Pakistani rheumatoid arthritis patients and healthy controls

Linkage disequilibrium at the FCGRgenetic locus

There was evidence of weak linkage disequilibrium (D' = 0.30, P = 0.01) between FCGR2A and FCGR3A in the UK Caucasian but not the North Indian/Pakistani control populations. Highly significant linkage disequilibrium was seen between FCGR3B and FCGR2B in both populations (UK Caucasian: D' = -0.68, P = 0.0001; North Indian/Pakistani:D' = -0.52, P = 0.0001). The negative D' values indicate linkage disequilibrium between the common allele of one gene and the rare allele of the second gene. No significant linkage disequilibrium was detected between FCGR3A and FCGR3B in either ethnic group, although D' = 0.40 in the UK Caucasian group (Table 3).

Table 3 Pair-wise linkage disequilibrium measures (D') calculated from the control groups in the two populations

Association of FCGRhaplotypes with rheumatoid arthritis

The distributions of two locus FCGR haplotypes were compared between the RA cohorts and their control populations, with a difference approaching statistical significance for FCGR3A–FCGR3B (Table 4). Compared with the control frequency of 24%, the FCGR3A–FCGR3B 158V-NA2 haplotype was found at increased frequency in UK Caucasian RA patients (31%) and even higher frequency in those with nodular RA (37%).

Table 4 Estimated pairwise haplotype frequencies in rheumatoid arthritis patients and healthy controls

From the HTR analysis of FCGR3AFCGR3B haplotypes, the 158V-NA2 haplotype was found to have a significant effect on the risk for RA in UK Caucasians (OR 1.77, 95% CI 1.09–2.87; P = 0.02), taking the most common haplotype (158F-NA2) as baseline (Table 5). The effect was stronger in the small subgroup of UK Caucasian individuals with nodules (OR 2.51, 95% CI 1.15–5.49; P = 0.02).

Table 5 Haplotype trend regression in rheumatoid arthritis patients and healthy controls for FCGR3A-FCGR3B haplotypes

These are estimates of the effect of each haplotype, assuming that each of the individual's two haplotypes has an independent effect on risk with a combined multiplicative effect. The effect of the 158V-NA2 haplotype was found to be largely confined to those with two copies of this haplotype (data not shown). To estimate the effect under a recessive model, homozygosity for this haplotype was compared in the control population (frequency 4%) with the total RA population (11%), giving an OR of 3.18 (95% CI 1.13–8.92; P = 0.03) when comparing homozygotes with all others. Again, the effect of homozygosity was stronger in those with nodular RA (frequency 16%; OR 5.03, 95% CI 1.44–17.56; P = 0.01).

For the North Indian/Pakistani cohort the same haplotype was found to be at increased frequency in RA patients compared with controls (OR 1.52, 95% CI 0.82–2.80, from the HTR analysis) but the difference was not statistically significant (P = 0.19; Table 5). Homozygosity for this haplotype was seen in approximately 4% of the RA and 1.5% of the control population.

Stepwise logistic regression analyses of FCGR3A and FCGR3Bin the Caucasian group

Considering each locus separately, FCGR3A is associated with RA in Caucasians [11] but FCGR3B is not (Table 2). However comparing the model containing both genotypes with the model containing FCGR3A only, there is an improved fit with inclusion of FCGR3B2 = 6.27, 2 degrees of freedom, from likelihood ratio test; P = 0.04). When the cohort with nodules was examined, the inclusion of FCGR3B did not significantly improve the model (P = 0.22).

Contribution of FCGRhaplotypes and shared epitope alleles in rheumatoid arthritis susceptibility

One advantage of the HTR framework for analysis of haplotypes is that other factors can be included in the model. The analysis was repeated including the RA-associated 'shared epitope' (SE) alleles (positive or negative) in the models. As expected, the SE itself was still highly predictive of RA in these models (OR 3.16, 95% CI 1.75–5.71 in UK Caucasians; OR 3.94, 95% CI 2.17–7.18 in the North Indian/Pakistani group). There was evidence of a multiplicative joint effect between SE and the FCGR haplotypes, consistent with both of these two genetic factors contributing to the risk for disease. Thus, the risk for RA in SE-positive UK Caucasian individuals homozygous for the FCGR3AFCGR3B 158V-NA2 haplotype is increased tenfold compared with those with other FCGR3 genotypes who are SE negative.

Discussion

Haplotype analyses have started to assume an increased importance in genetic studies of human disease because they can be more informative in their ability to identify unique chromosomal segments that are likely to harbour disease predisposing genes. They may also provide additional evidence for the presence of further unidentified polymorphic variants that are the true disease-susceptibility variants [36]. We have demonstrated an increased level of association between FCGR3A–FCGR3B haplotypes and RA compared with FCGR3A alone. The effects were stronger in the subset of RA patients with nodules. The two-locus haplotype showing the strongest association with RA susceptibility in each group was the FCGR3A–FCGR3B 158V-NA2 haplotype. UK Caucasian individuals who were homozygous for this haplotype were estimated to be at threefold risk for disease compared with all others (OR 3.18, 95% CI 1.13–8.92; P = 0.03), and this analysis does not depend on inference of uncertain phase because individuals that are homozygous for a haplotype are unambiguously identified by their genotype. The relative importance of the FCGR3A-158V and FCGR3B-NA2 polymorphic variants were assessed further using stepwise regression analyses [34]. These analyses showed that, although only FCGR3A has been shown to be associated with RA [11], the model including both loci provided an improved fit for RA susceptibility, but not necessarily so for the development of nodules. In the North Indian/Pakistani population there was a 6% increase in the FCGR3AFCGR3B 158V-NA2 haplotype, but this failed to reach statistical significance.

Several methods have been proposed for the analysis of such data in the absence of family data and the consequent presence of phase uncertainty. The HTR method we have chosen to use has various advantages. The method uses all the data, including that from individuals with uncertain phase. Because the data are analyzed in a regression framework, the usual regression diagnostics are available, other factors can be included in the model and tests for interaction can be performed. A potential disadvantage of the method is that weights used in the regression analysis are based on estimated haplotype frequencies, and the uncertainty inherent in the estimates is ignored in the model. This leads to an anti-conservative test and could lead to false-positive results. However, Stram and coworkers [37] evaluated the method in comparison with other more sophisticated approaches and found it to perform well in most situations.

Haplotype frequencies were estimated using the expectation-maximization algorithm, which has been shown to perform well, even in the presence of some departure from Hardy–Weinberg equilibrium [38]. Errors due to sampling are generally of much greater concern than inaccuracies due to the estimation process.

We acknowledge that none of our results are highly statistically significant when a consideration is made for multiple tests. However, the consistent pattern of results, taken in conjunction with prior findings in relation to these genes and their known biological functions in humans and mice, gives additional credence to them. Replication of these findings in other populations will ultimately be required.

The FcγRs play important roles in the initiation and propagation of many different immunological and inflammatory processes. Consequently, they may act as susceptibility factors for RA through a variety of mechanisms. FCGR3A was the most significant gene in this study, and we have previously discussed the role that this receptor may play in RA pathogenesis [11, 12]. In humans, FcγRIIIa is expressed on natural killer cells, macrophages, γδ T cells, a subset of monocytes and cultured mast cells [5, 6]. Higher levels of FcγRII and FcγRIII expression have been demonstrated in synovial biopsy specimens from RA patients compared with control individuals [39]. Similarly, an increase in the expression level and proportion of circulating FcγRIIIa-positive monocytes has been observed in RA and may correlate with disease activity [40, 41]. In addition, in vitro derived macrophages from RA patients expressed more FcγRII and FcγRIIIa, and released higher levels of tumour necrosis factor-α and matrix degrading enzymes in response to heat-aggregated IgG [39] compared with controls. These findings are supportive of our own work whereby the higher affinity genetic variant of FCGR3A may sensitize FcγR-bearing cells to IgG-containing immune complexes. FcγRIIIa may also play an important role delivering (auto)antigens, and activation and maturation signals to dendritic cells [42]. This may provide an explanation for the over tenfold increased risk for RA in SE-positive Caucasian individuals homozygous for the FCGR3AFCGR3B 158V-NA2 haplotype, which has been a consistent finding by a number of groups [1113, 15, 18].

FcγRIIIb is selectively expressed on neutrophils and eosinophils, and has a low affinity for IgG. It is linked to the membrane by a glycosylphosphatidylinositol anchor and does not appear to associate with the known transmembrane adapter molecules [5, 6]. However, FcγRIIIb appears to interact with FcγRIIa in the phagocytosis of immune complexes and subsequent cellular activation, with signalling being mediated through the ITAM (immunoreceptor tyrosine-based activation motif) of FcγRIIa [43, 44]. FCGR3B has two common polymorphic forms, namely NA1 and NA2, which differ in five nucleotides that produce four amino acid differences. This alters the number of glycosylation sites, and neutrophils from individuals homozygous for the FCGR3B-NA2 allele have been found consistently to exhibit lower levels of phagocytosis than FCGR3B-NA1 homozygotes [45]. This polymorphism has important biological consequences, especially in the development of blood transfusion reactions, autoimmune neutropenias and the severity of renal disease in systemic vasculitis [6, 46]. Individuals with duplications and deletions of FCGR3B have been reported [30, 47], with the estimated frequency of the FCGR3B deletion being 0.001–0.08 in various Caucasian populations [48]. Standard genotyping assays, as performed in the present study, do not allow a calculation of the gene copy number. This may provide an explanation for a failure of our control populations to conform to Hardy–Weinberg equilibrium and the previously reported non-Mendelian segregation in some Caucasian families [49].

FcγRIIb plays a crucial role in the regulation of antibody production and susceptibility to several spontaneous and induced murine autoimmune diseases [5052]. We found no evidence of an association between FCGR2B- or FCGR2B-containing haplotypes and RA in our cohorts, unlike previous observations in a Japanese cohort in which an alternative SNP in FCGR2B was investigated [15].

Conclusion

There is good data that FcγRs may be critical modulators of inflammation within the synovium and that subtle changes in either expression or structure of these receptors may influence both the susceptibility to RA and the development of nodules. The analyses performed in this study have strengthened our original observation that the FCGR genetic locus is associated with RA, particularly in a UK Caucasian population with nodular disease. Our haplotype data, together with the stepwise regression analysis, suggest that additional polymorphic variants within FCGR3A or in linkage disequilibrium with the FCGR3A–FCGR3B 158V-NA2 haplotype may contribute to RA pathogenesis.

Abbreviations

ARMS:

amplification refractory mutation system

BAC:

bacterial artificial chromosome

bp:

base pairs

BLAST:

basic local alignment search tool

CI:

confidence interval

FcγR:

Fcγ receptor

HTR:

haplotype trend regression

NA:

neutrophil antigen

OR:

odds ratio

PCR:

polymerase chain reaction

RA:

rheumatoid arthritis

RF:

rheumatoid factor

SE:

shared epitope

SNP:

single nucleotide polymorphism

UTR:

untranslated region.

References

  1. Lin JP, Cash JM, Doyle SZ, Peden S, Kanik K, Amos CI, Bale SJ, Wilder RL: Familial clustering of rheumatoid arthritis with other autoimmune diseases. Hum Genet. 1998, 103: 475-482. 10.1007/s004390050853.

    Article  CAS  PubMed  Google Scholar 

  2. Williams DG: Autoimmunity in rheumatoid arthritis. Rheumatology. Edited by: Klippel JH, Dieppe PA. 1998, London: Mosby-Year Book Europe, 3.9.1-3.9.14.

    Google Scholar 

  3. Lindqvist E, Eberhardt K, Bendtzen K, Heinegard D, Saxne T: Prognostic laboratory markers of joint damage in rheumatoid arthritis. Ann Rheum Dis. 2005, 64: 196-201. 10.1136/ard.2003.019992.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Smolen JS, Steiner G: Are autoantibodies active players or epiphenomena?. Curr Opin Rheumatol. 1998, 10: 201-206. 10.1097/00002281-199805000-00008.

    Article  CAS  PubMed  Google Scholar 

  5. Salmon JE, Pricop L: Human receptors for immunoglobulin G: key elements in the pathogenesis of rheumatic disease. Arthritis Rheum. 2001, 44: 739-750. 10.1002/1529-0131(200104)44:4<739::AID-ANR129>3.0.CO;2-O.

    Article  CAS  PubMed  Google Scholar 

  6. Dijstelbloem HM, van de Winkel JG, Kallenberg CG: Inflammation in autoimmunity: receptors for IgG revisited. Trends Immunol. 2001, 22: 510-516. 10.1016/S1471-4906(01)02014-2.

    Article  CAS  PubMed  Google Scholar 

  7. Pricop L, Redecha P, Teillaud JL, Frey J, Fridman WH, Sautès-Fridman C, Salmon JE: Differential modulation of stimulatory and inhibitory Fcγ receptors on human monocytes by Th1 and Th2 cytokines. J Immunol. 2001, 166: 531-537.

    Article  CAS  PubMed  Google Scholar 

  8. Ravetch JV, Lanier LL: Immune inhibitory receptors. Science. 2000, 290: 84-89. 10.1126/science.290.5489.84.

    Article  CAS  PubMed  Google Scholar 

  9. Nabbe KC, Blom AB, Holthuysen AE, Boross P, Roth J, Verbeek S, Van Lent PL, van den Berg WB: Coordinate expression of activating Fcγ receptors I and III and inhibiting Fcγ receptor type II in the determination of joint inflammation and cartilage destruction during immune complex-mediated arthritis. Arthritis Rheum. 2003, 48: 255-265. 10.1002/art.10721.

    Article  CAS  PubMed  Google Scholar 

  10. Kleinau S, Martinsson P, Heyman B: Induction and suppression of collagen-induced arthritis is dependent on distinct Fcγ receptors. J Exp Med. 2000, 191: 1611-1616. 10.1084/jem.191.9.1611.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Morgan AW, Griffiths B, Ponchel F, Montague BM, Ali M, Gardner PP, Gooi HC, Situnayake RD, Markham AF, Emery P, Isaacs JD: Fcγ receptor type IIIA is associated with rheumatoid arthritis in two distinct ethnic groups. Arthritis Rheum. 2000, 43: 2328-2334. 10.1002/1529-0131(200010)43:10<2328::AID-ANR21>3.0.CO;2-Z.

    Article  CAS  PubMed  Google Scholar 

  12. Morgan AW, Keyte VH, Babbage SJ, Robinson JI, Ponchel F, Barrett JH, Bhakta BB, Bingham SJ, Buch MH, Conaghan PG, et al: FcγRIIIA-158V and rheumatoid arthritis: a confirmation study. Rheumatology (Oxford). 2003, 42: 528-533. 10.1093/rheumatology/keg169.

    Article  CAS  Google Scholar 

  13. Nieto A, Cáliz R, Pascual M, Matarán L, García S, Martin J: Involvement of Fcγ receptor IIIA genotypes in susceptibility to rheumatoid arthritis. Arthritis Rheum. 2000, 43: 735-739. 10.1002/1529-0131(200004)43:4<735::AID-ANR3>3.0.CO;2-Q.

    Article  CAS  PubMed  Google Scholar 

  14. Milicic A, Misra R, Agarwal S, Aggarwal A, Brown MA, Wordsworth BP: The F158V polymorphism in FcγRIIIA shows disparate associations with rheumatoid arthritis in two genetically distinct populations. Ann Rheum Dis. 2002, 61: 1021-1023. 10.1136/ard.61.11.1021.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Kyogoku C, Tsuchiya N, Matsuta K, Tokunaga K: Studies on the association of Fcγ receptor IIA, IIB, IIIA and IIIB polymorphisms with rheumatoid arthritis in the Japanese: evidence for a genetic interaction between HLA-DRB1 and FCGR3A. Genes Immun. 2002, 3: 488-493. 10.1038/sj.gene.6363921.

    Article  CAS  PubMed  Google Scholar 

  16. Brun JG, Madland TM, Vedeler CA: Immunoglobulin G Fc-receptor (FcγR) IIA, IIIA, and IIIB polymorphisms related to disease severity in rheumatoid arthritis. J Rheumatol. 2002, 29: 1135-1140.

    CAS  PubMed  Google Scholar 

  17. Martinez A, Pascual M, Pascual-Salcedo D, Balsa A, Martin J, de la Concha EG: Genetic polymorphisms in Spanish rheumatoid arthritis patients: an association and linkage study. Genes Immun. 2003, 4: 117-121. 10.1038/sj.gene.6363931.

    Article  CAS  PubMed  Google Scholar 

  18. Radstake TR, Petit E, Pierlot C, van de Putte LB, Cornelis F, Barrera P: Role of Fcγ receptors IIA, IIIA, and IIIB in susceptibility to rheumatoid arthritis. J Rheumatol. 2003, 30: 926-933.

    CAS  PubMed  Google Scholar 

  19. Cardon LR, Bell JI: Association study designs for complex diseases. Nat Rev Genet. 2001, 2: 91-99. 10.1038/35052543.

    Article  CAS  PubMed  Google Scholar 

  20. Su Y, Brooks DG, Li L, Lepercq J, Trofatter JA, Ravetch JV, Lebo RV: Myelin protein zero gene mutated in Charcot-Marie-tooth type 1B patients. Proc Natl Acad Sci USA. 1993, 90: 10856-10860. 10.1073/pnas.90.22.10856.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. Griffiths B, Situnayake RD, Clark B, Tennant A, Salmon M, Emery P: Racial origin and its effect on disease expression and HLA-DRB1 types in patients with rheumatoid arthritis: a matched cross-sectional study. Rheumatology (Oxford). 2000, 39: 857-864. 10.1093/rheumatology/39.8.857.

    Article  CAS  Google Scholar 

  22. National Center for Biotechnology Information (NCBI). [http://www.ncbi.nlm.nih.gov]

  23. Ensembl. [http://www.ensembl.org/index.html]

  24. Computational Biology at Oak Ridge National Laboratory (ORNL). [http://genome.ornl.gov]

  25. Qiu WQ, de Bruin D, Brownstein BH, Pearse R, Ravetch JV: Organisation of the human and mouse low-affinity FcγR genes: duplication and recombination. Science. 1990, 248: 732-735. 10.1126/science.2139735.

    Article  CAS  PubMed  Google Scholar 

  26. Warmerdam PAM, van den Herik-Oudijk IE, Parren PW, Westerdaal NA, van de Winkel JG, Capel PJ: Interaction of a human FcγRIIb1 (CD32) isoform with murine and human IgG subclasses. Int Immunol. 1993, 5: 239-247.

    Article  CAS  PubMed  Google Scholar 

  27. van den Herik-Oudijk IE, Westerdaal NA, Henriquez NV, Capel PJ, van de Winkel JG: Functional analysis of human FcγRII (CD32) isoforms expressed in B lymphocytes. J Immunol. 1994, 152: 574-585.

    CAS  PubMed  Google Scholar 

  28. Norris CF, Pricop L, Millard SS, Taylor SM, Surrey S, Schwartz E, Salmon JE, McKenzie SE: A naturally occurring mutation in FcγRIIA: a Q to K127 change confers unique IgG binding properties to the R131 allelic form of the receptor. Blood. 1998, 91: 656-662.

    CAS  PubMed  Google Scholar 

  29. Smyth LJ, Snowden N, Carthy D, Papasteriades C, Hajeer A, Ollier WE: FcγRIIa polymorphism in systemic lupus erythematosus. Ann Rheum Dis. 1997, 56: 744-746.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  30. de Haas M, Kleijer M, van Zwieten R, Roos D, von dem Borne AE: Neutrophil FcγRIIIb deficiency, nature, and clinical consequences: a study of 21 individuals from 14 families. Blood. 1995, 86: 2403-2413.

    CAS  PubMed  Google Scholar 

  31. Morgan AW, Griffiths B, Barrett JH, Markham AF, Emery P, Isaacs JD: Fcγ receptor type IIIA is associated with rheumatoid arthritis in two distinct ethnic groups. Reply. Arthritis Rheum. 2002, 46: 557-559. 10.1002/art.10123.

    Article  Google Scholar 

  32. Zhao JH, Curtis D, Sham PC: Model-free analysis and permutation tests for allelic associations. Hum Hered. 2000, 50: 133-139. 10.1159/000022901.

    Article  CAS  PubMed  Google Scholar 

  33. Zaykin DV, Westfall PH, Young SS, Karnoub MA, Wagner MJ, Ehm MG: Testing association of statistically inferred haplotypes with discrete and continuous traits in samples of unrelated individuals. Hum Hered. 2002, 53: 79-91. 10.1159/000057986.

    Article  PubMed  Google Scholar 

  34. Cordell HJ, Clayton DG: A unified stepwise regression procedure for evaluating the relative effects of polymorphisms within a gene using case/control or family data: application to HLA in Type 1 diabetes. Am J Hum Genet. 2002, 70: 124-141. 10.1086/338007.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Su K, Wu J, Edberg JC, McKenzie SE, Kimberly RP: Genomic organization of classical human low-affinity Fcgamma receptor genes. Genes Immun. 2002, 3 Suppl 1: S51-S56. 10.1038/sj.gene.6363879.

    Article  CAS  PubMed  Google Scholar 

  36. Fallin D, Cohen A, Essioux L, Chumakov I, Blumenfeld M, Cohen D, Schork NJ: Genetic analysis of case/control data using estimated haplotype frequencies: application to APOE locus variation and Alzheimer's disease. Genome Res. 2001, 11: 143-151. 10.1101/gr.148401.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. Stram DO, Leigh Pearce C, Bretsky P, Freedman M, Hirschhorn JN, Altshuler D, Kolonel LN, Henderson BE, Thomas DC: Modeling and E-M estimation of haplotype-specific relative risks from genotype data for a case-control study of unrelated individuals. Hum Hered. 2003, 55: 179-190. 10.1159/000073202.

    Article  PubMed  Google Scholar 

  38. Fallin D, Schork NJ: Accuracy of haplotype frequency estimation for biallelic loci, via the expectation-maximisation algorithm for unphased diploid genotype data. Am J Hum Genet. 2000, 67: 947-959. 10.1086/303069.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  39. Blom AB, Radstake TR, Holthuysen AE, Slöetjes AW, Pesman GJ, Sweep FG, van de Loo FA, Joosten LA, Barrera P, Van Lent PL, van den Berg WB: Increased expression of Fcγ receptors II and III on macrophages of rheumatoid arthritis patients results in higher production of tumor necrosis factor α and matrix metalloproteinase. Arthritis Rheum. 2003, 48: 1002-1014. 10.1002/art.10871.

    Article  CAS  PubMed  Google Scholar 

  40. Hepburn AL, Mason JC, Davies KA: Expression of Fcγ and complement receptors on peripheral blood monocytes in systemic lupus erythematosus and rheumatoid arthritis. Rheumatology (Oxford). 2004, 43: 547-554. 10.1093/rheumatology/keh112.

    Article  CAS  Google Scholar 

  41. Wijngaarden S, van Roon JA, Bijlsma JW, van de Winkel JG, Lafeber FP: Fcγ receptor expression levels on monocytes are elevated in rheumatoid arthritis patients with high erythrocyte sedimentation rates who do not use anti-rheumatic drugs. Rheumatology (Oxford). 2003, 42: 681-688. 10.1093/rheumatology/keg174.

    Article  CAS  Google Scholar 

  42. Regnault A, Lankar D, Lacabanne V, Rodriguez A, Théry C, Rescigno M, Saito T, Verbeek S, Bonnerot C, Ricciardi-Castagnoli P, Amigorena S: Fcγ receptor-mediated induction of dendritic cell maturation and major histocompatibility complex class I-restricted antigen presentation after immune complex internalisation. J Exp Med. 1999, 189: 371-380. 10.1084/jem.189.2.371.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  43. Chuang FY, Sassaroli M, Unkeless JC: Convergence of Fcγ receptor IIA and Fcγ receptor IIIB signaling pathways in human neutrophils. J Immunol. 2000, 164: 350-360.

    Article  CAS  PubMed  Google Scholar 

  44. Salmon JE, Millard SS, Brogle NL, Kimberly RP: Fcγ receptor IIIb enhances Fcγ receptor IIa function in an oxidant-dependent and allele-sensitive manner. J Clin Invest. 1995, 95: 2877-2885.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  45. Salmon JE, Edberg JC, Kimberly RP: Fc gamma receptor III on human neutrophils. Allelic variants have functionally distinct capacities. J Clin Invest. 1990, 85: 1287-1295.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  46. Kimberly RP, Salmon JE, Edberg JC: Receptors for immunoglobulin G. Molecular diversity and implications for disease. Arthritis Rheum. 1995, 38: 306-314.

    Article  CAS  PubMed  Google Scholar 

  47. Koene HR, Kleijer M, Roos D, de Haas M, Von dem Borne AE: Fc gamma RIIIB gene duplication: evidence for presence and expression of three distinct Fc gamma RIIIB genes in NA(1+,2+)SH(+) individuals. Blood. 1998, 91: 673-679.

    CAS  PubMed  Google Scholar 

  48. Muniz-Diaz E, Madoz P, de la Calle Martin O, Puig L: The polymorphonuclear neutrophil Fc gamma RIIIb deficiency is more frequent than hitherto assumed. Blood. 1995, 86: 3999-

    CAS  PubMed  Google Scholar 

  49. Kyogoku C, Tsuchiya N, Wu H, Tsao BP, Tokunaga K: Association of Fcgamma receptor IIA, but not IIB and IIIA, polymorphisms with systemic lupus erythematosus: A family-based association study in Caucasians. Arthritis Rheum. 2004, 50: 671-673. 10.1002/art.20029.

    Article  CAS  PubMed  Google Scholar 

  50. Jiang Y, Hirose S, Abe M, Sanokawa-Akakura R, Ohtsuji M, Mi X, Li N, Xiu Y, Zhang D, Shirai J, et al: Polymorphisms in IgG Fc receptor IIB regulatory regions associated with autoimmune susceptibility. Immunogenetics. 2000, 51: 429-435. 10.1007/s002510050641.

    Article  CAS  PubMed  Google Scholar 

  51. Pritchard NR, Cutler AJ, Uribe S, Chadban SJ, Morley BJ, Smith KG: Autoimmune-prone mice share a promoter haplotype associated with reduced expression and function of the Fc receptor FcγRII. Curr Biol. 2000, 10: 227-230. 10.1016/S0960-9822(00)00344-4.

    Article  CAS  PubMed  Google Scholar 

  52. Johansson AC, Sundler M, Kjellén P, Johannesson M, Cook A, Lindqvist AK, Nakken B, Bolstad AI, Jonsson R, Alarcón-Riquelme M, Holmdahl R: Genetic control of collagen-induced arthritis in a cross with NOD and C57BL/10 mice is dependent on gene regions encoding complement factor 5 and FcγRIIb and is not associated with loci controlling diabetes. Eur J Immunol. 2001, 31: 1847-1856. 10.1002/1521-4141(200106)31:6<1847::AID-IMMU1847>3.0.CO;2-F.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This work was supported by grants from the Arthritis Research Campaign and the Medical Research Council, UK. In addition, the authors would like to acknowledge Dr Philip Gardner for performing some DNA extractions and helpful discussions with Dr Ian Carr regarding some laboratory aspects of this project.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Ann W Morgan.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

AWM participated in the design of the study, undertook all database searches, oversaw all aspects of the laboratory work, analyzed the data and prepared the manuscript. JHB gave additional statistical support and performed the haplotype analysis. BG, RDS and PE participated in the collection of clinical data and the recruitment of patients into the study. DS, JR and VK undertook some of the genotyping assays on DNA prepared in the laboratory of EAJ and RWO, who participated in the original design of the study. FP and MA gave invaluable advice during the retrieval of sequence data from the public databases and during the optimization of some genotyping assays. AWB, AFM, PE and JDI participated in the design of the study, interpretation of the results and writing of the final manuscript.

Rights and permissions

Reprints and permissions

About this article

Cite this article

Morgan, A.W., Barrett, J.H., Griffiths, B. et al. Analysis of Fcγ receptor haplotypes in rheumatoid arthritis: FCGR3A remains a major susceptibility gene at this locus, with an additional contribution from FCGR3B. Arthritis Res Ther 8, R5 (2005). https://doi.org/10.1186/ar1847

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/ar1847

Keywords