A case-control study of rheumatoid arthritis identifies an associated single nucleotide polymorphism in the NCF4 gene, supporting a role for the NADPH-oxidase complex in autoimmunity

Abstract

Rheumatoid arthritis (RA) is a chronic inflammatory disease with a heritability of 60%. Genetic contributions to RA are made by multiple genes, but only a few gene associations have yet been confirmed. By studying animal models, reduced capacity of the NADPH-oxidase (NOX) complex, caused by a single nucleotide polymorphism (SNP) in one of its components (the NCF1 gene), has been found to increase severity of arthritis. To our knowledge, however, no studies investigating the potential role played by reduced reactive oxygen species production in human RA have yet been reported. In order to examine the role played by the NOX complex in RA, we investigated the association of 51 SNPs in five genes of the NOX complex (CYBB, CYBA, NCF4, NCF2, and RAC2) in a Swedish case-control cohort consisting of 1,842 RA cases and 1,038 control individuals. Several SNPs were found to be mildly associated in men in NCF4 (rs729749, P = 0.001), NCF2 (rs789181, P = 0.02) and RAC2 (rs1476002, P = 0.05). No associations were detected in CYBA or CYBB. By stratifying for autoantibody status, we identified a strong association for rs729749 (in NCF4) in autoantibody negative disease, with the strongest association detected in rheumatoid factor negative men (CT genotype versus CC genotype: odds ratio 0.34, 95% confidence interval 0.2 to 0.6; P = 0.0001). To our knowledge, this is the first genetic association identified between RA and the NOX complex, and it supports previous findings from animal models of the importance of reactive oxygen species production capacity to the development of arthritis.

Figure 1

Haplotype blocks in NCF4. Three haplotype blocks were identified in NCF4. Block 2 contains the associated single nucleotide polymorphism (SNP) rs729749. Colour scheme of the linakge disequilibrium LD map is based on the standard D'/LOD option in the Haploview software. The LD blocks are calculated based on the CI method.