Cellular adhesion gene SELP is associated with rheumatoid arthritis and displays differential allelic expression.

In rheumatoid arthritis (RA), a key event is infiltration of inflammatory immune cells into the synovial lining, possibly aggravated by dysregulation of cellular adhesion molecules. Therefore, single nucleotide polymorphisms of 14 genes involved in cellular adhesion processes (CAST, ITGA4, ITGB1, ITGB2, PECAM1, PTEN, PTPN11, PTPRC, PXN, SELE, SELP, SRC, TYK2, and VCAM1) were analyzed for association with RA. Association analysis was performed consecutively in three European RA family sample groups (Nfamilies = 407). Additionally, we investigated differential allelic expression, a possible functional consequence of genetic variants. SELP (selectin P, CD62P) SNP-allele rs6136-T was associated with risk for RA in two RA family sample groups as well as in global analysis of all three groups (ptotal = 0.003). This allele was also expressed preferentially (p<10-6) with a two- fold average increase in regulated samples. Differential expression is supported by data from Genevar MuTHER (p1 = 0.004; p2 = 0.0177). Evidence for influence of rs6136 on transcription factor binding was also found in silico and in public datasets reporting in vitro data. In summary, we found SELP rs6136-T to be associated with RA and with increased expression of SELP mRNA. SELP is located on the surface of endothelial cells and crucial for recruitment, adhesion, and migration of inflammatory cells into the joint. Genetically determined increased SELP expression levels might thus be a novel additional risk factor for RA.

The interferon regulatory factor 5 gene confers susceptibility to rheumatoid arthritis and influences its erosive phenotype.

BACKGROUND: Increased expression of type I IFN genes, also referred to as an IFN signature, has been detected in various autoimmune diseases including rheumatoid arthritis (RA). Interferon regulatory factors, such as IRF5, coordinate type I IFN expression. Multiple IRF5 variants were suggested as autoimmunity susceptibility factors. OBJECTIVE: As the linkage proof remains important to establish fully any genetic RA susceptibility factor, the authors took advantage of the largest reported European trio family resource dedicated to RA to test for linkage IRF5 and performed a genotype-phenotype analysis. METHODS: 1140 European Caucasian individuals from 380 RA trio families were genotyped for IRF5 rs3757385, rs2004640 and rs10954213 single nucleotide polymorphisms (SNP). RESULTS: Single marker analysis provided linkage evidence for each IRF5 SNP investigated. IRF5 linked to RA with two haplotypes: the CTA risk haplotype 'R' (transmission (T)=60.6%, p=23.1×10(-5)) and the AGG protective haplotype 'P' (T=39.6%, p=0.0015). Linkage was significantly stronger in non-erosive disease for both IRF5 R and P haplotypes (T=73.9%, p=4.20×10(-5) and T=19.6%, p=3.66×10(-5), respectively). Multivariate logistic regression analysis found IRF5 linked to RA independently of the rheumatoid factor status. IRF5 RR and PP haplotypic genotypes were associated with RA, restricted to the non-erosive phenotype: p=1.68×10(-4), OR 4.80, 95% CI 2.06 to 11.19; p=0.003, OR 0.17, 95% CI 0.05 to 0.57, respectively. CONCLUSION: This study provides the 'association and linkage proof' establishing IRF5 as a RA susceptibility gene and the identification of a genetic factor that seems to contribute to the modulation of the erosive phenotype. Further studies are warranted to clarify the role of IRF5 in RA and its subphenotypes.

Association of MICA with rheumatoid arthritis independent of known HLA-DRB1 risk alleles in a family-based and a case control study.

INTRODUCTION: The gene MICA encodes the protein major histocompatibility complex class I polypeptide-related sequence A. It is expressed in synovium of patients with rheumatoid arthritis (RA) and its implication in autoimmunity is discussed. We analyzed the association of genetic variants of MICA with susceptibility to RA. METHODS: Initially, 300 French Caucasian individuals belonging to 100 RA trio families were studied. An additional 100 independent RA trio families and a German Caucasian case-control cohort (90/182 individuals) were available for replication. As MICA is situated in proximity to known risk alleles of the HLA-DRB1 locus, our analysis accounted for linkage disequilibrium either by analyzing the subgroup consisting of parents not carrying HLA-DRB1 risk alleles with transmission disequilibrium test (TDT) or by implementing a regression model including all available data. Analysis included a microsatellite polymorphism (GCT)n and single-nucleotide polymorphisms (SNPs) rs3763288 and rs1051794. RESULTS: In contrast to the other investigated polymorphisms, the non-synonymously coding SNP MICA-250 (rs1051794, Lys196Glu) was strongly associated in the first family cohort (TDT: P = 0.014; regression model: odds ratio [OR] 0.46, 95% confidence interval [CI] 0.25 to 0.82, P = 0.007). Although the replication family sample showed only a trend, combined family data remained consistent with the hypothesis of MICA-250 association independent from shared epitope (SE) alleles (TDT: P = 0.027; regression model: OR 0.56, 95% CI 0.38 to 0.83, P = 0.003). We also replicated the protective association of MICA-250A within a German Caucasian cohort (OR 0.31, 95% CI 0.1 to 0.7, P = 0.005; regression model: OR 0.6, 95% CI 0.37 to 0.96, P = 0.032). We showed complete linkage disequilibrium of MICA-250 (D' = 1, r2= 1) with the functional MICA variant rs1051792 (D' = 1, r2= 1). As rs1051792 confers differential allelic affinity of MICA to the receptor NKG2D, this provides a possible functional explanation for the observed association. CONCLUSIONS: We present evidence for linkage and association of MICA-250 (rs1051794) with RA independent of known HLA-DRB1 risk alleles, suggesting MICA as an RA susceptibility gene. However, more studies within other populations are necessary to prove the general relevance of this polymorphism for RA.

Testing for the association of the KIAA1109/Tenr/IL2/IL21 gene region with rheumatoid arthritis in a European family-based study.

INTRODUCTION: A candidate gene approach, in a large case-control association study in the Dutch population, has shown that a 480 kb block on chromosome 4q27 encompassing KIAA1109/Tenr/IL2/IL21 genes is associated with rheumatoid arthritis. Compared with case-control association studies, family-based studies have the added advantage of controlling potential differences in population structure. Therefore, our aim was to test this association in populations of European origin by using a family-based approach. METHODS: A total of 1,302 West European white individuals from 434 trio families were genotyped for the rs4505848, rs11732095, rs6822844, rs4492018 and rs1398553 polymorphisms using the TaqMan Allelic discrimination assay (Applied Biosystems). The genetic association analyses for each SNP and haplotype were performed using the Transmission Disequilibrium Test and the genotype relative risk. RESULTS: We observed evidence for association of the heterozygous rs4505848-AG genotype with rheumatoid arthritis (P = 0.04); however, no significance was found after Bonferroni correction. In concordance with previous findings in the Dutch population, we observed a trend of undertransmission for the rs6822844-T allele and rs6822844-GT genotype to rheumatoid arthritis patients. We further investigated the five SNP haplotypes of the KIAA1109/Tenr/IL2/IL21 gene region. We observed, as described in the Dutch population, a nonsignificant undertransmission of the AATGG haplotype to rheumatoid arthritis patients. CONCLUSIONS: Using a family-based study, we have provided a trend for the association of the KIAA1109/Tenr/IL2/IL21 gene region with rheumatoid arthritis in populations of European descent. Nevertheless, we failed to replicate a significant association of this region in our rheumatoid arthritis family sample. Further investigation of this region, including detection and testing of all variants, is required to confirm rheumatoid arthritis association.

No evidence of major effects in several Toll-like receptor gene polymorphisms in rheumatoid arthritis.

INTRODUCTION: The objective was to study the potential genetic contribution of Toll-like receptor (TLR) genes in rheumatoid arthritis (RA). TLRs bind to pathogen-associated molecular patterns, and TLR genes influence both proinflammatory cytokine production and autoimmune responses. Host-pathogen interactions are involved in RA physiopathology. METHODS: We tested SNPs of five TLR genes (TLR9, TLR2, TLR6, TLR1, and TLR4) in a cohort of 100 French families with RA. Genotypes were analyzed using the transmission disequilibrium test. As TLR2, TLR6, and TLR1 are located on chromosome 4, we determined the haplotype relative risk. Analyses were performed in subgroups defined by status for rheumatoid factor, anti-cyclic citrullinated peptide autoantibodies, and erosions. RESULTS: We found no disequilibrium in allele transmission for any of the SNPs of the five TLR genes. In subgroup analyses, no associations were detected linking TLR9, TLR2, or TLR9/TLR2 to rheumatoid factor, anti-cyclic citrullinated peptide autoantibodies, or erosions. Haplotype analysis of the polymorphisms showed no haplotype associations in any of the subgroups. CONCLUSIONS: We found no evidence of major effects of TLR gene polymorphisms in RA, although we tested different TLR phenotypes. Moreover, no associations were noted with autoantibody production or erosions.

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.

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

The integrin alpha(v)beta3, whose alpha(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% versus 50% 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.

Linkage proof for PTPN22, a rheumatoid arthritis susceptibility gene and a human autoimmunity gene.

The tyrosine phosphatase PTPN22 allele 1858T has been associated with rheumatoid arthritis (RA) and other autoimmune diseases. RA is the most frequent of those multifactorial diseases. The RA association was usually restricted to serum rheumatoid factor positive disease (RF+). No interaction was shown with HLA-DRB1, the first RA gene. Many case-control studies replicated the RA association, showing an allele frequency increase of approximately 5% on average and large variations of population allele frequencies (2.1-15.5%). In multifactorial diseases, the final proof for a new susceptibility allele is provided by departure from Mendel's law (50% transmission from heterozygous parents). For PTPN22-1858T allele, convincing linkage proof was available only for type 1 diabetes. We aimed at providing this proof for RA. We analyzed 1,395 West European Caucasian individuals from 465 "trio" families. We replicated evidence for linkage, demonstrating departure from Mendel's law in this subset of early RA onset patients. We estimated the overtransmission of the 1858T allele in RF+ families: T = 63%, P < 0.0007. The 1858T allele frequency increased from 11.0% in controls to 17.4% in RF+ RA for the French Caucasian population and the susceptibility genotype (1858T/T or T/C) from 20.2% to 31.6% [odds ratio (OR) = 1.8 (1.2-2.8)]. In conclusion, we provided the linkage proof for the PTPN22-1858T allele and RF+ RA. With diabetes and RA, PTPN22 is therefore a "linkage-proven" autoimmunity gene. PTPN22 accounting for approximately 1% of the RA familial aggregation, many new genes could be expected that are as many leads to definitive therapy for autoimmune diseases.

Validation of the reshaped shared epitope HLA-DRB1 classification in rheumatoid arthritis.

Recently, we proposed a classification of HLA-DRB1 alleles that reshapes the shared epitope hypothesis in rheumatoid arthritis (RA); according to this model, RA is associated with the RAA shared epitope sequence (72-74 positions) and the association is modulated by the amino acids at positions 70 and 71, resulting in six genotypes with different RA risks. This was the first model to take into account the association between the HLA-DRB1 gene and RA, and linkage data for that gene. In the present study we tested this classification for validity in an independent sample. A new sample of the same size and population (100 RA French Caucasian families) was genotyped for the HLA-DRB1 gene. The alleles were grouped as proposed in the new classification: S1 alleles for the sequences A-RAA or E-RAA; S2 for Q or D-K-RAA; S3D for D-R-RAA; S3P for Q or R-R-RAA; and X alleles for no RAA sequence. Transmission of the alleles was investigated. Genotype odds ratio (OR) calculations were performed through conditional logistic regression, and we tested the homogeneity of these ORs with those of the 100 first trio families (one case and both parents) previously reported. As previously observed, the S2 and S3P alleles were significantly over-transmitted and the S1, S3D and X alleles were under-transmitted. The latter were grouped as L alleles, resulting in the same three-allele classification. The risk hierarchy of the six derived genotypes was the same: (by decreasing OR and with L/L being the reference genotype) S2/S3P, S2/S2, S3P/S3P, S2/L and S3P/L. The homogeneity test between the ORs of the initial and the replication samples revealed no significant differences. The new classification was therefore considered validated, and both samples were pooled to provide improved estimates of RA risk genotypes from the highest (S2/S3P [OR 22.2, 95% confidence interval 9.9-49.7]) to the lowest (S3P/L [OR 4.4, 95% confidence interval 2.3-8.4]).

Rheumatoid arthritis seropositive for the rheumatoid factor is linked to the protein tyrosine phosphatase nonreceptor 22-620W allele.

The protein tyrosine phosphatase nonreceptor type 22 (PTPN22) gene encodes for lymphoid tyrosine phosphatase LYP, involved in the negative regulation of early T-cell activation. An association has recently been reported between the PTPN22-620W functional allele and rheumatoid factor-positive (RF+) rheumatoid arthritis (RA), among other autoimmune diseases. Expected linkage proof for consistency cannot be definitely produced by an affected sib-pair (ASP) analysis. Our aim was therefore to search for linkage evidence with the transmission disequilibrium test. DNA from the French Caucasian population was available for two samples of 100 families with one RA patient and both parents, and for 88 RA index cases from RA ASP families. Genotyping was carried out by PCR-restriction fragment length polymorphism. The analysis was performed using the transmission disequilibrium test, genotype relative risk and ASP-based analysis. The transmission disequilibrium test of the PTPN22-620W allele revealed linkage and association for RF+ RA (61% of transmission, P = 0.037). The genotype relative risk showed the risk allele in 34% of RF+ RA patients and in 24% of controls derived from nontransmitted parental chromosomes (P = 0.047, odds ratio = 1.69, 95% confidence interval = 1.03-2.78). The ASP investigation showed no enriched risk allele in RA multiplex families, resulting in a lack of power of ASP analysis, explaining the published negative results. This study is the first to show linkage of PTPN22 to RF+ RA, consistent with PTPN22 as a new RA gene.

New classification of HLA-DRB1 alleles supports the shared epitope hypothesis of rheumatoid arthritis susceptibility.

OBJECTIVE: The shared epitope hypothesis was formulated to explain the involvement of HLA-DRB1 in rheumatoid arthritis (RA). However, several studies, which considered only the HLA-DRB1 alleles shown to be associated with RA risk, rejected this hypothesis. In this report, we propose that a different classification of HLA-DRB1 alleles be considered, based on the amino acid sequence at position 70-74. METHODS: The fit of both HLA-DRB1 classifications was tested in 2 groups of RA patients. All subjects were recruited through the European Consortium on Rheumatoid Arthritis Families, and included 100 patients with isolated RA and 132 patients with at least 1 affected sibling. RESULTS: The new classification produced risk estimates that fit all of the observed data, i.e., the distribution of the HLA-DRB1 genotype in the 2 patient groups, and the distribution of parental alleles shared by affected sibpairs. The risk of developing RA under this new classification depends on whether the RAA sequence occupies position 72-74 but is modulated by the amino acid at position 71 (K confers the highest risk, R an intermediate risk, A and E a lower risk) and by the amino acid at position 70 (Q or R confers a higher risk than D). CONCLUSION: A new classification based on amino acid sequence allows us to show that the shared epitope RAA sequence at position 72-74 explains the data, with the risk of developing RA modulated by the amino acids at positions 70 and 71.

Dense genome-wide linkage analysis of rheumatoid arthritis, including covariates.

OBJECTIVE: Rheumatoid arthritis (RA) is a heterogeneous disease that exhibits a complex genetic component. Previous RA genome scans confirmed the involvement of the HLA region and generated data on suggestive signals at non-HLA regions, albeit with few overlaps in findings between studies. The present study was undertaken to detect potential RA gene regions and to estimate the number of true RA gene regions, taking into account the heterogeneity of RA, through performance of a dense genome scan. METHODS: In a study of 88 French Caucasian families (105 RA sibpairs), 1,088 microsatellite markers were genotyped (3.3-cM genome scan), and a multipoint model-free linkage analysis was performed. The statistical assessment of the results relied on 10,000 computer simulations. A covariate-based multipoint model-free linkage analysis was performed on the locations of regions with suggestive evidence for linkage. RESULTS: Involvement of the HLA region was strongly confirmed (P = 6 x 10(-5)), and 19 non-HLA regions showed suggestive evidence for linkage (P < 0.05); 9 of these overlapped with regions suggested in other published RA genome scans. A routine 12-cM genome scan with the same families would have detected only 7 of the 19 regions, including only 4 of the 9 overlapping regions. From the 10,000 computer simulations, we estimated that 8 +/- 4 regions (mean +/- SD) were true-positives. RA covariate-based analysis provided additional linkage evidence for 3 regions, with age at disease onset, erosions, and HLA-DRB1 shared epitope as covariates. CONCLUSION: The results of this study provide evidence of 19 non-HLA RA gene regions, with an estimate of 8 +/- 4 as true-positives, and provide additional evidence for 3 regions from covariate-based analysis.

Association between tumor necrosis factor receptor II and familial, but not sporadic, rheumatoid arthritis: evidence for genetic heterogeneity.

OBJECTIVE: Tumor necrosis factor alpha (TNFalpha) binds the receptors TNFRI and TNFRII. Results of genome scans have suggested that TNFR2 is a candidate rheumatoid arthritis (RA) locus. A case-control study in a UK Caucasian population has shown an association between a TNFR2 genotype (196R/R in exon 6) and familial, but not sporadic, RA. The present study was undertaken to test this association in the French Caucasian population. METHODS: To test for an association in sporadic RA, 100 families were genotyped for the 196M/R polymorphism and analyzed using the transmission disequilibrium test and haplotype relative risk. To test for an association in familial RA, RA index cases from 100 affected sibpair (ASP) families were genotyped for 196M/R. Linkage analysis was performed with 3 TNFR2 microsatellite markers. RESULTS: The TNFR2 196R/R genotype was not associated with sporadic RA (odds ratio [OR] 0.59, P = 0.72), but was associated with familial RA (OR 4.0, P = 0.026). The association was most marked in the context of TNFR2 "twin-like" RA sibs (affected sibs sharing both TNFR2 haplotypes) (OR 9.2, P = 0.0017). Linkage analysis results were consistent with the association; most of the TNFR2 linkage evidence was found in the subgroup of families with 196R/R ASP index cases. CONCLUSION: This study is the first to replicate evidence of the involvement of TNFR2 in RA genetic heterogeneity. Our data refine the initial hypothesis, to suggest that a TNFR2 recessive factor, in linkage disequilibrium with the 196R allele, plays a major role in a subset of families with multiple cases of RA.