Rheumatoid Arthritis Susceptibility Loci Research Articles

Association of rheumatoid factor and anti-cyclic citrullinated peptide positivity, but not carriage of shared epitope or PTPN22 susceptibility variants, with anti-tumour necrosis factor response in rheumatoid arthritis

Objective:To determine whether rheumatoid factor (RF), anti-cyclic citrullinated peptide (CCP) antibodies, or carriage of shared epitope (SE) and PTPN22 genetic susceptibility variants predict response to therapy in patients with rheumatoid...

14.5 Rheumatoid Arthritis susceptibility loci; STAT4, TRAF1/C5 and 6q23 region, are also associated with Juvenile Idiopathic Arthritis

14.5 Rheumatoid Arthritis susceptibility loci; STAT4, TRAF1/C5 and 6q23 region, are also associated with Juvenile Idiopathic Arthritis

Replication of the tumor necrosis factor receptor−associated factor 1/complement component 5 region as a susceptibility locus for rheumatoid arthritis in a European family‐based study

We recently showed, using a candidate gene approach in a case-control association study, that a 65-kb block encompassing tumor necrosis factor receptor-associated factor 1 (TRAF1) and C5 is strongly associated with rheumatoid arthritis (RA). Compared with case-control association studies, family-based studies have the added advantage of controlling potential differences in population structure and are not likely to be hampered by variation in population allele frequencies, as is seen for many genetic polymorphisms, including the TRAF1/C5 locus. The aim of this study was to confirm this association in populations of European origin by using a family-based approach. A total of 1,356 western European white individuals from 452 "trio" families were genotyped for the rs10818488 polymorphism, using the TaqMan allelic discrimination assay. We observed evidence for association, demonstrating departure from Mendel's law, with an overtransmission of the rs10818488 A allele (A = 55%; P = 0.036). By taking into consideration parental phenotypes, we also observed an increased A allele frequency in affected versus unaffected parents (A = 64%; combined P = 0.015). Individuals carrying the A allele had a 1.2-fold increased risk of developing RA (allelic odds ratio 1.24, 95% confidence interval 1.04-1.50). Using a family-based study that is robust against population stratification, we provide evidence for the association of the TRAF1/C5 rs10818488 A allele and RA in populations of European descent, further substantiating our previous findings. Future functional studies should yield insight into the biologic relevance of this locus to the pathways involved in RA.

Genome‐wide association study of rheumatoid arthritis in the Spanish population: KLF12 as a risk locus for rheumatoid arthritis susceptibility

To identify new genes associated with susceptibility to rheumatoid arthritis (RA), using a 2-stage genome-wide association study. Following a liability-based study design, we analyzed 317,503 single-nucleotide polymorphisms (SNPs) in 400 patients with RA and 400 control subjects. We selected a group of candidate SNPs for replication in an independent group of 410 patients with RA and 394 control subjects. Using data from the 3 previous genome-wide association studies in RA, we also looked for genomic regions showing evidence of common association signals. Finally, we analyzed the presence of genome-wide epistasis using the binary test implemented in the PLINK program. We identified several genomic regions showing evidence of genome-wide association (P < 1 x 10(-5)). In the replication analysis, we identified KLF12 SNP rs1324913 as the most strongly associated SNP (P = 0.01). In our study, we observed that this SNP showed higher significance than PTPN22 SNP rs2476601, in both the genome-wide association studies and the replication analyses. Furthermore, the integration of our data with those from previous genome-wide association studies showed that KLF12 and PTPRT are the unique loci that are commonly associated in 3 different studies (P = 0.004 and P = 0.002 for KLF12 in the Wellcome Trust Case Control Consortium study and the Brigham and Women's Rheumatoid Arthritis Sequential Study genome-wide association study, respectively). The genome-wide epistasis analysis identified several SNP pairs close to significance after multiple test correction. The present genome-wide association study identified KLF12 as a new susceptibility gene for RA. The joint analysis of our results and those from previous genome-wide association studies showed genomic regions with a higher probability of being genuine susceptibility loci for RA.

Re-evaluation of putative rheumatoid arthritis susceptibility genes in the post-genome wide association study era and hypothesis of a key pathway underlying susceptibility

Rheumatoid arthritis (RA) is an archetypal, common, complex autoimmune disease with both genetic and environmental contributions to disease aetiology. Two novel RA susceptibility loci have been reported from recent genome-wide and candidate gene association studies. We, therefore, investigated the evidence for association of the STAT4 and TRAF1/C5 loci with RA using imputed data from the Wellcome Trust Case Control Consortium (WTCCC). No evidence for association of variants mapping to the TRAF1/C5 gene was detected in the 1860 RA cases and 2930 control samples tested in that study. Variants mapping to the STAT4 gene did show evidence for association (rs7574865, P = 0.04). Given the association of the TRAF1/C5 locus in two previous large case–control series from populations of European descent and the evidence for association of the STAT4 locus in the WTCCC study, single nucleotide polymorphisms mapping to these loci were tested for association with RA in an independent UK series comprising DNA from >3000 cases with disease and >3000 controls and a combined analysis including the WTCCC data was undertaken. We confirm association of the STAT4 and the TRAF1/C5 loci with RA bringing to 5 the number of confirmed susceptibility loci. The effect sizes are less than those reported previously but are likely to be a more accurate reflection of the true effect size given the larger size of the cohort investigated in the current study.

Association mapping of susceptibility loci for rheumatoid arthritis

We analyzed a case-control data set for chromosome 18q from the Genetic Analysis Workshop 15 to detect susceptibility loci for rheumatoid arthritis (RA). A total number of 460 cases and 460 unaffected controls were genotyped on 2300 single-nucleotide polymorphisms (SNPs) by the North American Rheumatoid Arthritis Consortium. Using a multimarker approach for association mapping under the framework of the Malecot model and composite likelihood, we identified a region showing significant association with RA (p < 0.002) and the predicted disease locus was at a genomic location of 53,306 kb with a 95% confidence interval (CI) of 53,295-53,331 kb. A common haplotype in this region was protective against RA (p = 0.002). In another region showing nominal significant association (51,585 kb, 95% CI: 51,541-51,628 kb, p = 0.037), a haplotype was also protective (p = 0.002). We further demonstrated that reducing SNP density decreased power and accuracy of association mapping. SNP selection based on equal linkage disequilibrium (LD) distance generally produced higher accuracy than that based on equal kilobase distance or tagging.

One- and two-locus models for mapping rheumatoid arthritis-susceptibility genes on chromosome 6

To evaluate whether there is evidence for two rheumatoid arthritis (RA) susceptibility genes on chromosome 6, we applied new robust methods for two-locus multipoint identical-by-descent mapping to the rheumatoid arthritis data of the Genetic Analysis Workshop 15. The software GEEARP was used to estimate the locations and the corresponding genetic effects for one locus or two linked loci in a region on chromosome 6, on the basis of affected relative pairs. These methods were applied to the data sets from the North American Rheumatoid Arthritis Consortium, Canada, France, and the first screen of United Kingdom. From the resultant 95% confidence intervals given by a robust variance estimator, a linked region, other than the well-known HLA region, was at 54.7-69.6 cM, providing evidence for a second rheumatoid arthritis susceptibility locus on chromosome 6.

Large-scale linkage analysis of 1302 affected relative pairs with rheumatoid arthritis

Rheumatoid arthritis is the most common systematic autoimmune disease and its etiology is believed to have both strong genetic and environmental components. We demonstrate the utility of including genetic and clinical phenotypes as covariates within a linkage analysis framework to search for rheumatoid arthritis susceptibility loci. The raw genotypes of 1302 affected relative pairs were combined from four large family-based samples (North American Rheumatoid Arthritis Consortium, United Kingdom, European Consortium on Rheumatoid Arthritis Families, and Canada). The familiality of the clinical phenotypes was assessed. The affected relative pairs were subjected to autosomal multipoint affected relative-pair linkage analysis. Covariates were included in the linkage analysis to take account of heterogeneity within the sample. Evidence of familiality was observed with age at onset (p << 0.001) and rheumatoid factor (RF) IgM (p << 0.001), but not definite erosions (p = 0.21). Genome-wide significant evidence for linkage was observed on chromosome 6. Genome-wide suggestive evidence for linkage was observed on chromosomes 13 and 20 when conditioning on age at onset, chromosome 15 conditional on gender, and chromosome 19 conditional on RF IgM after allowing for multiple testing of covariates.

HSPD1 is not a major susceptibility gene for rheumatoid arthritis in the French Caucasian population

The heat shock 60-kDa protein 1 (HSP60) is involved in immune and inflammatory reactions, which are hallmarks of rheumatoid arthritis (RA). HSP60 is encoded by the HSPD1 gene located on 2q33, one of the suggested RA susceptibility loci in the French Caucasian population. Our aim was to test whether HSPD1 is a major susceptibility gene by studing families from the French Caucasian population. Three single nucleotide polymorphisms (SNPs) were studied in 100 RA trio families, and 100 other families were used for replication. Genetic analyses were performed by comparing allelic frequencies, by applying the transmission disequilibrium test, and by assessing the genotype relative risk. We observed a significant RA association for the C/C genotype of rs2340690 in the first sample. However, this association was not confirmed when the second sample was added. The two other SNPs and the haplotype analysis did not give any significant results. We conclude that HSPD1 is not a major RA susceptibility gene in the French Caucasian population.

Functional characterization of NF- B inhibitor-like protein 1 (NF BIL1), a candidate susceptibility gene for rheumatoid arthritis

Several studies have implicated the NF-kappaB inhibitor-like protein 1 (NFkBIL1) gene located in the class III region of the major histocompatibility complex (MHC) as a possible susceptibility locus for rheumatoid arthritis (RA). Based on limited homology, it has been suggested to be a member of the inhibitor of NF-kappaB (IkappaB) family of proteins, but a role in mRNA processing has also been proposed. We have investigated the expression of NFkBIL1 in RA synovial tissue and characterized its function. Real-time PCR showed the two NFkBIL1 mRNA splice variants are expressed in a tissue-specific manner. Dual immunofluorescent staining of human RA synovium with polyclonal anti-NFkBIL1 antibodies and anti-CD68, anti-CD3 or anti-factor VIII showed that NFkBIL1 was expressed in the rheumatoid synovial lining and sub-lining layers and co-localized in CD68+ and CD3+, but not Factor VIII+ cells. Confocal microscopy of cultured synovial fibroblasts revealed expression in speckled nuclear and homogenous cytoplasmic distributions, suggesting shuttling between the cytoplasmic and nuclear compartments. Functional tests showed that NFkBIL1 isoforms were incapable of associating with NF-kappaB and did not inhibit it, thus disproving the hypothesis that NFkBIL1 functions as an IkappaB. Affinity purification of endogenous NFkBIL1 proteins and co-immunoprecipitation experiments showed that NFkBIL1 can associate with mRNA and with three protein partners, identified by mass spectrometry as leukophysin, translation elongation factor 1 alpha and CTP synthase I. These data support a potential role for NFkBL1 in the pathogenesis of RA and indicates that it may be involved in mRNA processing or the regulation of translation.

The ITGAVrs3738919-C allele is associated with rheumatoid arthritis in the European Caucasian population: a family-based study

The integrin αvβ3, whose αv subunit is encoded by the ITGAV gene, plays a key role in angiogenesis. Hyperangiogenesis is involved in rheumatoid arthritis (RA) and the ITGAV gene is located in 2q31, one of the suggested RA susceptibility loci. Our aim was to test the ITGAV gene for association and linkage to RA in a family-based study from the European Caucasian population.Two single nucleotide polymorphisms were genotyped by PCR-restriction fragment length polymorphism in 100 French Caucasian RA trio families (one RA patient and both parents), 100 other French families and 265 European families available for replication. The genetic analyses for association and linkage were performed using the comparison of allelic frequencies (affected family-based controls), the transmission disequilibrium test, and the genotype relative risk.We observed a significant RA association for the C allele of rs3738919 in the first sample (affected family-based controls, RA index cases 66.5% versus controls 56.7%; P = 0.04). The second sample showed the same trend, and the third sample again showed a significant RA association. When all sets were combined, the association was confirmed (affected family-based controls, RA index cases 64.6% versus controls 58.1%; P = 0.005). The rs3738919-C allele was also linked to RA (transmission disequilibrium test, 56.5% versus50% of transmission; P = 0.009) and the C-allele-containing genotype was more frequent in RA index cases than in controls (RA index cases 372 versus controls 339; P = 0.002, odds ratio = 1.94, 95% confidence interval = 1.3–2.9).The rs3738919-C allele of the ITGAV gene is associated with RA in the European Caucasian population, suggesting ITGAV as a new minor RA susceptibility gene.

Molecular Profile of Peripheral Blood Mononuclear Cells from Patients with Rheumatoid Arthritis

Rheumatoid arthritis (RA) is a chronic inflammatory arthritis. Currently, diagnosis of RA may take several weeks, and factors used to predict a poor prognosis are not always reliable. Gene expression in RA may consist of a unique signature. Gene expression analysis has been applied to synovial tissue to define molecularly distinct forms of RA; however, expression analysis of tissue taken from a synovial joint is invasive and clinically impractical. Recent studies have demonstrated that unique gene expression changes can be identified in peripheral blood mononuclear cells (PBMCs) from patients with cancer, multiple sclerosis, and lupus. To identify RA disease-related genes, we performed a global gene expression analysis. RNA from PBMCs of 9 RA patients and 13 normal volunteers was analyzed on an oligonucleotide array. Compared with normal PBMCs, 330 transcripts were differentially expressed in RA. The differentially regulated genes belong to diverse functional classes and include genes involved in calcium binding, chaperones, cytokines, transcription, translation, signal transduction, extracellular matrix, integral to plasma membrane, integral to intracellular membrane, mitochondrial, ribosomal, structural, enzymes, and proteases. A k-nearest neighbor analysis identified 29 transcripts that were preferentially expressed in RA. Ten genes with increased expression in RA PBMCs compared with controls mapped to a RA susceptibility locus, 6p21.3. These results suggest that analysis of RA PBMCs at the molecular level may provide a set of candidate genes that could yield an easily accessible gene signature to aid in early diagnosis and treatment.

Linkage analysis of rheumatoid arthritis in US and UK families reveals interactions between HLA–DRB1 and loci on chromosomes 6q and 16p

HLA is the most strongly associated locus in rheumatoid arthritis (RA), accounting for up to one-third of the genetic contribution. Conditioning on the effect of true disease loci such as HLA can lead to increased power to detect effects at other loci and, in addition, allows investigation of the underlying disease models, including interactions. The aim of this study was to detect susceptibility loci for RA by conditioning on HLA in a large sample of affected sibling pairs (ASPs) and to test for evidence of interaction between novel loci and HLA. Genotype data from 3 whole-genome linkage scans for RA in a US population and a UK population were pooled, resulting in a combined data set of 886 ASPs. This pooling of data increased the power to detect loci showing low levels of heterogeneity. Nonparametric linkage analysis was performed to identify regions of interest. Joint 2-locus analysis was then performed for HLA and each of the loci that demonstrated evidence of linkage in the 886 ASPs. Evidence for linkage was most significant at HLA (P = 4 x 10(-16)), with 7 non-HLA loci showing some evidence for linkage (P = 0.05-0.003). Joint modeling of these loci with HLA provided evidence for linkage at a genome-wide significance level for loci on 6q (P = 2.7 x 10(-6)) and 16p (P = 2 x 10(-4)). These data provide the most convincing evidence to date that 6q and 16p harbor susceptibility genes. In addition, these loci may interact with HLA, facilitating the search for candidate genes within this region.

Association of PTPN22 1858 single-nucleotide polymorphism with rheumatoid arthritis in a German cohort: higher frequency of the risk allele in male compared to female patients

The functional single-nucleotide polymorphism (SNP) of the gene PTPN22 is a susceptibility locus for rheumatoid arthritis (RA). The study presented here describes the association of the PTPN22 1858T allele with RA in a German patient cohort; 390 patients with RA and 349 controls were enrolled in the study. For 123 patients, clinical and radiographic documentation over 6 years was available from the onset of disease. Genotyping of the PTPN22 1858 SNP was performed using an restriction fragment length polymorphism PCR-based genotyping assay. The odds ratio to develop RA was 2.57 for carriers of the PTPN22 1858T allele (95% confidence interval (CI) 1.85–3.58, p < 0.001), and 5.58 for homozygotes (95% CI 1.85–16.79). The PTPN22 1858T allele was significantly associated not only with rheumatoid factor (RF) and anti-cyclic citrullinated peptide (CCP) positive RA, but also with RF and anti-CCP negative disease. The frequency of the PTPN22 1858T allele was increased disproportionately in male patients (53.8% compared to 33.0% in female patients, p < 0.001), and the resulting odds ratio for male carriers was increased to 4.47 (95% CI 2.5–8.0, p < 0.001). Moreover, within the male patient population, the rare allele was significantly associated with the HLA-DRB1 shared epitope (p = 0.01). No significant differences in disease activity or Larsen scores were detected. The results provide further evidence that the PTPN22 1858T allele is associated with RA irrespective of autoantibody production. The increased frequency of the risk allele in male patients and its association with the shared epitope indicate that the genetic contribution to disease pathogenesis might be more prominent in men.

カルシウム結合によって活性化されるヒストン修飾酵素ＰＡＤ４の構造科学的基盤

Peptidylarginine deimianse 4 (PAD4) is a Ca2+-dependent enzyme that catalyzes the conversion of both arginine and mono-methyl arginine in histones into citrullines, and regulates both histone argininine methylation level and gene activity. Its gene is susceptibility locus for rheumatoid arthritis (RA) . Here we present the crystal structure of Ca2+-free wild-type PAD4, which shows that the polypeptide chain adopts an elongated fold in which the N-terminal domain forms two immunoglobulin-like subdomains, and the C-terminal domain forms an α/β propeller structure. Five Ca2+-binding sites, none of which adopts an EF-hand motif, were identified in the structure of a Ca2+-bound inactive mutant with and without bound substrate. These structural data indicate that Ca2+binding induces conformational changes that generate the active site cleft. Our findings identify a novel mechanism for enzyme activation by Ca2+ions, and are important for understanding the mechanism of protein citrullination and for developing PAD-inhibiting drugs for the treatment of RA.

Polymorphism within the telomeric MHC and expression of rheumatoid arthritis candidate genes

Several recent studies have indicated a second susceptibility locus for rheumatoid arthritis (RA) in the telomeric MHC close to the tumour necrosis factor locus, a region containing a number of genes of unknown function including HLA-B associated transcripts 1 and 2 (BAT1 and BAT2) and lymphotoxin alpha (LTA). Several single-nucleotide polymorphisms have been described in this region; however, the influence of each on susceptibility is not known. In this study we sought to examine the association of individual single-nucleotide polymorphisms with expression of each mRNA transcript in response to an inflammatory stimulus.

Evidence for a novel rheumatoid arthritis susceptibility locus on chromosome 6p

To investigate whether the large linkage peak on chromosome 6p harbors rheumatoid arthritis (RA) susceptibility loci in addition to the well-characterized HLA-DRB1 gene. DNA samples obtained from 377 UK RA affected sibling pair (ASP) families, comprising test (181 ASPs) and replication (196 ASPs) cohorts, were used for linkage analysis. Three hundred eighty-four patients with RA derived from a subset of 192 ASPs were compared with a panel of 288 unrelated healthy controls for association studies. Samples were genotyped for 35 microsatellites and 25 single-nucleotide polymorphisms (SNPs). In the test cohort, the maximum logarithm of odds (LOD) score was obtained over D6S1260 (LOD 7.1). Evidence for linkage to the telomeric portion of the peak was increased in subsets of ASPs in which both individuals had erosive disease or both carried 2 copies of the shared epitope. HLA-A, HLA-DRB1, and 8 additional markers showed evidence of linkage in the presence of association with RA (using the extended transmission disequilibrium test [ETDT]). The positive ETDT result for 2 adjacent markers (D6S1665 and 210901-4) mapping to the telomeric end of the linked region ( approximately 11 Mb from DRB1) was replicated (for D6S1665) in the second cohort of ASPs. Haplotypic overtransmission of pairwise combinations between D6S1665*7 and 210901-4*4 was identified through the TDTPhase program. Multipoint conditional analysis showed this effect to be independent of HLA-DRB1. SNP-based association studies of the region identified a 4-marker haplotype in the DEK gene that was significantly associated with RA (P = 0.009). Evidence has been presented for an RA susceptibility locus mapping under the linkage peak on 6p, 11 Mb telomeric of HLA-DRB1. Preliminary association data implicate the gene DEK.

Structural basis for Ca(2+)-induced activation of human PAD4.

Peptidylarginine deiminase 4 (PAD4) is a Ca(2+)-dependent enzyme that catalyzes the conversion of protein arginine residues to citrulline. Its gene is a susceptibility locus for rheumatoid arthritis. Here we present the crystal structure of Ca(2+)-free wild-type PAD4, which shows that the polypeptide chain adopts an elongated fold in which the N-terminal domain forms two immunoglobulin-like subdomains, and the C-terminal domain forms an alpha/beta propeller structure. Five Ca(2+)-binding sites, none of which adopt an EF-hand motif, were identified in the structure of a Ca(2+)-bound inactive mutant with and without bound substrate. These structural data indicate that Ca(2+) binding induces conformational changes that generate the active site cleft. Our findings identify a novel mechanism for enzyme activation by Ca(2+) ions, and are important for understanding the mechanism of protein citrullination and for developing PAD-inhibiting drugs for the treatment of rheumatoid arthritis.

Whole-Genome Scan, in a Complex Disease, Using 11,245 Single-Nucleotide Polymorphisms: Comparison with Microsatellites

Despite the theoretical evidence of the utility of single-nucleotide polymorphisms (SNPs) for linkage analysis, no whole-genome scans of a complex disease have yet been published to directly compare SNPs with microsatellites. Here, we describe a whole-genome screen of 157 families with multiple cases of rheumatoid arthritis (RA), performed using 11,245 genomewide SNPs. The results were compared with those from a 10-cM microsatellite scan in the same cohort. The SNP analysis detected HLA*DRB1, the major RA susceptibility locus (P=.00004), with a linkage interval of 31 cM, compared with a 50-cM linkage interval detected by the microsatellite scan. In addition, four loci were detected at a nominal significance level (P<.05) in the SNP linkage analysis; these were not observed in the microsatellite scan. We demonstrate that variation in information content was the main factor contributing to observed differences in the two scans, with the SNPs providing significantly higher information content than the microsatellites. Reducing the number of SNPs in the marker set to 3,300 (1-cM spacing) caused several loci to drop below nominal significance levels, suggesting that decreases in information content can have significant effects on linkage results. In contrast, differences in maps employed in the analysis, the low detectable rate of genotyping error, and the presence of moderate linkage disequilibrium between markers did not significantly affect the results. We have demonstrated the utility of a dense SNP map for performing linkage analysis in a late-age-at-onset disease, where DNA from parents is not always available. The high SNP density allows loci to be defined more precisely and provides a partial scaffold for association studies, substantially reducing the resource requirement for gene-mapping studies.

Additional genetic susceptibility for rheumatoid arthritis telomeric of the DRB1 locus.

Rheumatoid arthritis (RA) has an estimated genetic contribution of 30-50%, approximately one-third of which arises from the major histocompatibility complex on 6p21.3. Many studies have implicated alleles of DRB1 that encode a shared epitope. However, several recent studies have suggested that additional telomeric genetic influences may exist. In this study, we sought to investigate whether a separate non-DRB1 effect could be detected and to determine its likely location. We typed 13 single-nucleotide polymorphisms, located mainly in the telomeric class III region of the major histocompatibility complex, in 164 British Caucasian families with RA that had at least 1 affected offspring and used unconditioned and DRB1-conditioned transmission disequilibrium tests (TDTs). Unconditioned TDTs revealed overtransmission of shared epitope alleles (P = 2.12 x 10(-5)) and an allele of the HLA-B-associated transcript 1 (BAT1) gene in the telomeric class III region (P = 0.009). Using a DRB1-conditioned TDT to assess whether an independent effect existed, we detected unequal transmission of alleles of lymphocyte-specific transcript 1 (P = 0.004), BAT1 (P = 0.003), and PG8 (P = 0.003). At least 1 additional non-DRB1 susceptibility locus for RA exists in an interval that encompasses the junction of the class III and I regions. This is a genomic segment of high linkage disequilibrium containing a large number of poorly characterized immunomodulatory genes.