Genome-wide meta-analysis for rheumatoid arthritis

NLRP1, PTPN22 and PADI4 gene polymorphisms and rheumatoid arthritis in ACPA-positive Singaporean Chinese

Studies have shown that the genetic risk factors for rheumatoid arthritis (RA) differ substantially between Asian and Caucasian populations. Even among Asian populations, the genetic contributions of NLRP1, PTPN22 and PADI4 have been controversial. Consequently, we sought to address these separate findings and determine whether any of these proposed risk variants are associated with RA susceptibility, onset, DAS activity and erosion in a Singaporean Chinese cohort. We genotyped five SNPs within NLRP1 (rs878329 and rs6502867), PTPN22 (rs2488457 and rs6665194), and PADI4 (rs2240340) in 500 anti-cyclic citrullinated peptide antibody-positive (ACPA) patients with RA and 500 healthy controls using TaqMan assays. The CC genotype of NLRP1 rs878329 and TT genotype of PADI4 rs2240340 were associated with RA susceptibility. The risk association of the T allele of PADI4 rs2240340 with RA was confirmed through a meta-analysis based on previous reports in Asian populations. The GG genotype of PTPN22 rs6665194 (-3508A>G) was associated with significantly reduced risk of RA. No significant association was found for NLRP1 rs6502867 T/C and PTPN22 rs2488457 G/C polymorphisms. None of the five SNPs was associated with RA's clinical features. This work supports the association of the T allele of PADI4 rs2240340 with RA in Asians. The roles of NLRP1 rs878329 G/C and PTPN22 rs6665194 A/G polymorphisms were demonstrated for the first time. We also propose rs6665194 to be a promising candidate for RA risk evaluation between ethnicities.

Integrative analyses for functional mechanisms underlying associations for rheumatoid arthritis

OBJECTIVE

Extensive association analyses including genome-wide association studies (GWAS) and powerful metaanalysis studies have identified a long list of loci associated with rheumatoid arthritis (RA) in very large populations, but most of them established statistical associations of genetic markers and RA only at the DNA level, without supporting evidence of functional relevance. Our study serves as a trial to detect the functional mechanisms underlying associations for RA by searching publicly available datasets and results.

METHODS

Based on publicly available datasets and results, we performed integrative analyses (gene relationships across implicated loci analysis, differential gene expression analysis, and functional annotation clustering analysis) and combined them with the expression quantitative trait locus (eQTL) results to dissect functional mechanisms underlying the associations for RA.

RESULTS

By searching 2 GWAS, Integrator and PheGenI, we selected 98 RA association results (p < 10(-5)). Among these associations, we found that 8 single-nucleotide polymorphisms (SNP; rs1600249, rs2736340, rs3093023, rs3093024, rs4810485, rs615672, rs660895, and rs9272219) serve as cis-effect regulators of the corresponding eQTL genes (BLK and CD4 in non-HLA region; CCR6, HLA-DQA1, and HLA-DQB1 in HLA region) that also were differentially expressed in RA-related cell groups. These 5 genes are closely related with immune response in function.

CONCLUSION

Our results showed the functional mechanisms underlying the associations of 8 SNP and the corresponding genes. This study is an example of mining publicly available datasets and results in validation of significant disease-association results. Using public data resources for integrative analyses may provide insights into the molecular genetic mechanisms underlying human diseases.

Non-MHC risk alleles in rheumatoid arthritis and in the syntenic chromosome regions of corresponding animal models

Rheumatoid arthritis (RA) is a polygenic autoimmune disease primarily affecting the synovial joints. Numerous animal models show similarities to RA in humans; some of them not only mimic the clinical phenotypes but also demonstrate the involvement of homologous genomic regions in RA. This paper compares corresponding non-MHC genomic regions identified in rodent and human genome-wide association studies (GWAS). To date, over 30 non-MHC RA-associated loci have been identified in humans, and over 100 arthritis-associated loci have been identified in rodent models of RA. The genomic regions associated with the disease are designated by the name(s) of the gene having the most frequent and consistent RA-associated SNPs or a function suggesting their involvement in inflammatory or autoimmune processes. Animal studies on rats and mice preferentially have used single sequence length polymorphism (SSLP) markers to identify disease-associated qualitative and quantitative trait loci (QTLs) in the genome of F2 hybrids of arthritis-susceptible and arthritis-resistant rodent strains. Mouse GWAS appear to be far ahead of rat studies, and significantly more mouse QTLs correspond to human RA risk alleles.

Disease-promoting and -protective genomic loci on mouse chromosomes 3 and 19 control the incidence and severity of autoimmune arthritis

Proteoglycan (PG)-induced arthritis (PGIA) is a murine model of rheumatoid arthritis. Arthritis-prone BALB/c mice are 100% susceptible, whereas the major histocompatibility complex-matched DBA/2 strain is completely resistant to PGIA. To reduce the size of the disease-suppressive loci for sequencing and to find causative genes of arthritis, we created a set of BALB/c.DBA/2-congenic/subcongenic strains carrying DBA/2 genomic intervals overlapping the entire Pgia26 locus on chromosome 3 (chr3) and Pgia23/Pgia12 loci on chr19 in the arthritis-susceptible BALB/c background. Upon immunization of these subcongenic strains and their wild-type (BALB/c) littermates, we identified a major Pgia26a sublocus on chr3 that suppressed disease onset, incidence and severity via controlling the complex trait of T-cell responses. The region was reduced to 3 Mbp (11.8 Mbp with flanking regions) in size and contained gene(s) influencing the production of a number of proinflammatory cytokines. Additionally, two independent loci (Pgia26b and Pgia26c) suppressed the clinical scores of arthritis. The Pgia23 locus (∼3 Mbp in size) on chr19 reduced arthritis susceptibility and onset, and the Pgia12 locus (6 Mbp) associated with low arthritis severity. Thus, we have reached the critical sizes of arthritis-associated genomic loci on mouse chr3 and chr19, which are ready for high-throughput sequencing of genomic DNA.

Antibodies against cyclic citrullinated peptides of IgG, IgA and IgM isotype and rheumatoid factor of IgM and IgA isotype are increased in unaffected members of multicase rheumatoid arthritis families from northern Sweden

BACKGROUND

Rheumatoid factors (RFs) and antibodies against cyclic citrullinated peptides (CCPs) of IgG, IgA and IgM isotype have been shown to precede disease onset by years.

OBJECTIVE

To evaluate serological risk markers in first-degree relatives from multicase families in relation to genetic and environmental risk factors.

METHODS

51 multicase families consisting of 163 individuals with rheumatoid arthritis (RA) (mean±SD age, 60±14 years; disease duration 21 years; 71.8% female) and with 157 first-degree relatives unaffected by RA (54±17 years; 59.9% female) were recruited. Isotypes of antibodies against CCPs (IgG, IgA and IgM) and RFs (IgM and IgA) were determined using automated enzyme immunoassays. Cut-off levels were established using receiver operating characteristic curves based on values for 100 unrelated healthy controls.

RESULTS

The concentrations and frequencies of all anti-CCP and RF isotypes were significantly increased in first-degree relatives and patients with RA compared with unrelated healthy controls. The relative distribution of IgA and IgM isotypes was higher than IgG in the relatives, whereas the IgG isotype dominated in patients with RA. The patients carried human leucocyte antigen-shared epitope (HLA-SE) significantly more often than the relatives (71.4% vs 53.9%, p=0.01), while the frequency of the PTPN22 T variant was similar. HLA-SE, combined with smoking, was significantly related to all combinations of anti-CCP and RF isotypes in patients with RA. No such relationships were found for the first-degree relatives.

CONCLUSIONS

All anti-CCP and RF isotypes analysed occurred more commonly in unaffected first-degree relatives from multicase families than in controls, but with different isotype distribution from patients with RA.

Wnt signaling genes of murine chromosome 15 are involved in sex-affected pathways of inflammatory arthritis

OBJECTIVE

Sex disparities in rheumatoid arthritis (RA) are well documented despite the lack of any known major RA susceptibility genes mapped to sex chromosomes. Murine chromosome 15 carries the sex-affected Pgia8 locus that mediates proteoglycan-induced arthritis, and homologous human loci are associated with RA. This study was undertaken to identify genes/mechanisms implicated in sex disparities in arthritis.

METHODS

Gene expression analysis was performed using RNA isolated from the paws of male and female Pgia8-congenic mice with collagen antibody-induced arthritis. Results were corroborated by reverse transcription-polymerase chain reaction, and mice were also studied prior to disease onset. Ingenuity Pathways Analysis of the expression patterns and gene functions was used to discover locus-specific and sex-affected signature transcripts.

RESULTS

We found that the Pgia8 locus regulates antibody-mediated inflammatory arthritis differently in males and females. In Pgia8-congenic males, arthritis severity was 30% less (P < 0.005) than in wild-type males, but the antiinflammatory effect was similar in wild-type and congenic females. Transcriptome analysis indicated that 12 genes within the locus were significantly dysregulated in arthritic joints of congenic mice; expression of these genes was also sex specific. The genes that correlated most highly with arthritis severity included those for collagen triple-helix repeat-containing 1 (Cthrc1), metalloproteinase (Adamts12), R-spondin (Rspo2), and syndecan (Sdc2) (r = 0.87-0.91). The level of Cthrc1 message also correlated with that of the genes for the proinflammatory cytokines interleukin-1β and interleukin-6.

CONCLUSION

These results indicate that sex-specific disparities in RA are linked to transcriptional regulation of genes involved in cartilage degradation (Adamts12) and canonical and noncanonical Wnt signaling (Cthrc1, Rspo2, Sdc2).

Proteoglycan-induced arthritis and recombinant human proteoglycan aggrecan G1 domain-induced arthritis in BALB/c mice resembling two subtypes of rheumatoid arthritis

OBJECTIVE

To develop a simplified and relatively inexpensive version of cartilage proteoglycan-induced arthritis (PGIA), an autoimmunity model of rheumatoid arthritis (RA), and to evaluate the extent to which this new model replicates the disease parameters of PGIA and RA.

METHODS

Recombinant human G1 domain of human cartilage PG containing "arthritogenic" T cell epitopes was generated in a mammalian expression system and used for immunization of BALB/c mice. The development and progression of arthritis in recombinant human PG G1-immunized mice (designated recombinant human PG G1-induced arthritis [GIA]) was monitored, and disease parameters were compared with those in the parent PGIA model.

RESULTS

GIA strongly resembled PGIA, although the clinical symptoms and immune responses in mice with GIA were more uniform than in those with PGIA. Mice with GIA showed evidence of stronger Th1 and Th17 polarization than those with PGIA, and anti-mouse PG autoantibodies were produced in different isotype ratios in the 2 models. Rheumatoid factor (RF) and anti-cyclic citrullinated peptide (anti-CCP) antibodies were detected in both models; however, serum levels of IgG-RF and anti-CCP antibodies were different in GIA and PGIA, and both parameters correlated better with disease severity in GIA than in PGIA.

CONCLUSION

GIA is a novel model of seropositive RA that exhibits all of the characteristics of PGIA. Although the clinical phenotypes are similar, GIA and PGIA are characterized by different autoantibody profiles, and the 2 models may represent 2 subtypes of seropositive RA, in which more than 1 type of autoantibody can be used to monitor disease severity and response to treatment.

CIITA is not associated with risk of developing rheumatoid arthritis

The major histocompatibility complex (MHC) class II transactivator gene (CIITA) encodes an important transcription factor regulating genes required for human leukocyte antigen (HLA) class II MHC-restricted antigen presentation. MHC genes, particularly HLA class II, are strongly associated with risk of developing rheumatoid arthritis (RA). Given the strong biological relationship between CIITA and HLA class II genes, a comprehensive investigation of CIITA variation in RA was conducted. This study tested 31 CIITA single-nucleotide polymorphisms in 2542 RA cases and 3690 controls (N=6232). All individuals were of European ancestry, as determined by ancestry informative genetic markers. No evidence for association between CIITA variation and RA was observed after a correction for multiple testing was applied. This is the largest study to fully characterize common genetic variation in CIITA, including an assessment of haplotypes. Results exclude even a modest role for common CIITA polymorphisms in susceptibility to RA.

Genetic variations in genes encoding RANK, RANKL, and OPG in rheumatoid arthritis: a case-control study

OBJECTIVE

Rheumatoid arthritis (RA) is a systemic autoimmune disease characterized by chronic inflammation of the joints, which may lead to structural damage of the cartilage and bone. The receptor activator of nuclear factor-kappaB (RANK) and the osteoprotegerin (OPG) cascade system have been reported to be essential in osteoclastogenesis. Genetic variations in the genes coding for RANK, RANK ligand (RANKL), and OPG are thought to play roles in the susceptibility to RA.

METHODS

In our case-control study, genomic DNA was obtained from 534 patients with RA who fulfilled the American College of Rheumatology 1987 criteria and 516 healthy control blood donors (HC). We studied 7 single-nucleotide polymorphisms (SNP) in the genes of RANK (2 SNP: rs1805034, rs35211496), OPG (2 SNP: rs3102735, rs2073618), and RANKL (3 SNP: rs9533156, rs2277438, rs1054016) using TaqMan assay-guided polymerase chain reaction. Genotype and allelic frequencies comparing RA patients with HC were analyzed by chi-square test for 2 x 3 and 2 x 2 tables, respectively.

RESULTS

Genotype distributions of the SNP rs35211496 in the RANK gene as well as the SNP rs2277438 in the RANKL gene differed significantly between patients with RA and HC. The frequency of the minor allele of rs9533156 of RANKL was significantly higher in patients with RA than in HC (OR 0.84, 95% CI 0.71-0.99, p = 0.047). Multivariate analysis adjusted to sex and investigating SNP demonstrated a significantly elevated risk for RA associated with the major allele in the RANK SNP rs35211496 (p = 0.0231) and with the minor allele in the RANKL SNP rs2277438 (p = 0.0092). No significantly increased risk was detected in the other SNP.

CONCLUSION

The minor allele of the RANK SNP rs35211496 may be protective against RA, while the minor alleles of the RANKL SNP rs2277438 may increase susceptibility to RA.

A candidate gene study of CLEC16A does not provide evidence of association with risk for anti-CCP-positive rheumatoid arthritis

CLEC16A, a putative immunoreceptor, was recently established as a susceptibility locus for type I diabetes and multiple sclerosis. Subsequently, associations between CLEC16A and rheumatoid arthritis (RA), Addison's disease and Crohn's disease have been reported. A large comprehensive and independent investigation of CLEC16A variation in RA was pursued. This study tested 251 CLEC16A single-nucleotide polymorphisms in 2542 RA cases (85% anti-cyclic citrullinated peptide (anti-CCP) positive) and 2210 controls (N=4752). All individuals were of European ancestry, as determined by ancestry informative genetic markers. No evidence for significant association between CLEC16A variation and RA was observed. This is the first study to fully characterize common genetic variation in CLEC16A including assessment of haplotypes and gender-specific effects. The previously reported association between RA and rs6498169 was not replicated. Results show that CLEC16A does not have a prominent function in susceptibility to anti-CCP-positive RA.

Use of supplementary phenotype to identify additional rheumatoid arthritis loci in a linkage analysis of 342 UK affected sibling pair families

Background

Although rheumatoid arthritis has been shown to have moderately strong genetic component, both linked loci identified in linkage analyses and susceptibility variants from association studies are short of adequately accounting for a comprehensive catalogue of the molecular factors underlying this complex disease. The objective of this study was to use supplementary phenotype based on cumulative hazard of rheumatoid arthritis to identify linkage evidence for new and additional rheumatoid arthritis loci in a genome-wide linkage analysis of 342 affected sibling pair families from the United Kingdom.

Methods

Using proportional hazards model, we estimated cumulative hazard of rheumatoid arthritis and then used it as a quantitative trait in a non-parametric multipoint variance component linkage analysis with 353 microsatellite markers distributed across the 22 autosomal chromosomes.

Results

We identified 3 new loci with genome-wide suggestive linkage evidence for rheumatoid arthritis on 9q21.13, 15p11.1 and 20q13.33. Our results also confirmed previously reported linkage evidence in the HLA-DRB1 region on chromosome 6 and on locus 1q32.1.

Conclusion

This study demonstrates the potential for information gain through the use of supplementary phenotypes in genetic study of complex diseases to identify new and additional potential linked loci that are not detected by linkage analysis of traditional phenotypes; and our results provide further evidence of the involvement of multiple loci in the genetic aetiology of rheumatoid arthritis.

The effect of parental consanguinity on the clinical and laboratory findings of rheumatoid arthritis

AIMS

We aimed to evaluate the frequency of consanguinity among the parents of patients with rheumatoid arthritis (RA) and the influence of parental consanguinity on several clinical and laboratory parameters which reflect the severity of the disease.

METHODS AND PATIENTS

The study population consisted of 265 patients with RA which were divided into two groups with respect to the presence or absence of consanguinity between their parents. The frequency of parental consanguinity was compared with the general population. The two groups were compared with respect to family history of RA, the age of onset, the age at which RA was diagnosed, duration of the disease, the presence of rheumatoid nodules, vasculitis, serositis and the need for orthopaedic surgery, amyloidosis, the presence and level of rheumatoid factor and anti-cyclic citrullinated peptide antibodies, erosive changes on radiographs, and the need for anti-tumour necrosis factor therapy.

RESULTS

Twenty-one patients (8%) had parents who were consanguineous, which was not more frequent compared with the general population (14%). The mean age of disease onset and the mean age at which RA was diagnosed were lower in patients with parental consanguinity, although the difference was not statistically significant. The other clinical and laboratory parameters were also not different between the two groups.

CONCLUSION

The present data suggests that parental consanguinity has no effect on disease severity, and the frequency of consanguinity is not increased among the parents of patients with RA. A possible exception is the earlier disease onset and age at diagnosis which needs to be confirmed by larger studies.

Lack of association with rheumatoid arthritis of selected polymorphisms in 4 candidate genes: CFH, CD209, eotaxin-3, and MHC2TA

OBJECTIVE

To investigate associations with rheumatoid arthritis (RA) of single-nucleotide polymorphisms (SNP) in 4 candidate genes, complement factor H (CFH), CD209 or DC-SIGN, eotaxin-3, and the MHC class II Transactivator (MHC2TA) genes. These SNP have been reported as important for RA (eotaxin-3 and MHC2TA) or for other immune-mediated diseases (CFH and CD209).

METHODS

Genotypes for the 7 selected SNP were obtained from 1587 patients with RA and 1570 controls of Spanish ancestry. Analyses were carried out after stratification for sex, erosions, rheumatoid factor, shared epitope, anti-cyclic citrullinated peptide antibodies, and the R620W PTPN22 SNP.

RESULTS

None of the comparisons between patients with RA and controls or between the different strata of patients according to disease features was significant.

CONCLUSION

None of the SNP in CFH and CD209 showed evidence of association with RA. We did not replicate the association of eotaxin-3 with RA described in Koreans, or that of the MHC2T SNP.

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.

Combining multiple family-based association studies

While high-throughput genotyping technologies are becoming readily available, the merit of using these technologies to perform genome-wide association studies has not been established. One major concern is that for studies of complex diseases and traits, the whole-genome approach requires such large sample sizes that both recruitment and genotyping pose considerable challenge. Here we propose a novel statistical method that boosts the effective sample size by combining data obtained from several studies. Specifically, we consider a situation in which various studies have genotyped non-overlapping subjects at largely non-overlapping sets of markers. Our approach, which exploits the local linkage disequilibrium structure without assuming an explicit population model, opens up the possibility of improving statistical power by incorporating existing data into future association studies.

Combining linkage data sets for meta-analysis and mega-analysis: the GAW15 rheumatoid arthritis data set

We have used the genome-wide marker genotypes from Genetic Analysis Workshop 15 Problem 2 to explore joint evidence for genetic linkage to rheumatoid arthritis across several samples. The data consisted of four high-density genome scans on samples selected for rheumatoid arthritis. We cleaned the data, removed intermarker linkage disequilibrium, and assembled the samples onto a common genetic map using genome sequence positions as a reference for map interpolation. The individual studies were combined first at the genotype level (mega-analysis) prior to a multipoint linkage analysis on the combined sample, and second using the genome scan meta-analysis method after linkage analysis of each sample. The two approaches were compared, and give strong support to the HLA locus on chromosome 6 as a susceptibility locus. Other regions of interest include loci on chromosomes 11, 2, and 12.

Modeling of PTPN22 and HLA-DRB1 susceptibility to rheumatoid arthritis

In the present paper, we used the North American Rheumatoid Arthritis Consortium data provided for Genetic Analysis Workshop 15 Problem 2 to: 1) estimate the penetrances of PTPN22 and HLA-DRB1 and, 2) test the selected model of PTPN22 conditional on the rheumatoid factor status. To achieve these aims, we used the marker association segregation chi-square method, fitting simultaneously both genotype frequency and identical by descent distributions in a sample of 3690 White individuals from 604 nuclear families. A co-dominant model fitted the rs2476601 (R620W) single-nucleotide polymorphism (SNP) of the PTPN22 gene well, whereas a lack of fit for all models was observed for the HLA-DRB1 locus. Testing genetic models of rheumatoid arthritis that include the PTPN22 SNP in addition to the HLA-DRB1 locus did not affect the results, nor did subgroup analysis of PTPN22 conditional on the rheumatoid factor status. In conclusion, PTPN22 R620W SNP is a risk factor for rheumatoid arthritis. The genetic architecture of the HLA-DRB1 locus is highly complex, and more elaborate modeling of this locus is required.

Novel approach for genome scan meta-analysis of rheumatoid arthritis: a kernel-based estimation procedure

Genome scan meta-analysis (GSMA) can prove very useful in detecting genetic effects too small to be detected in an individual linkage study and can also lead to more consistent results. In this paper, we propose a new kernel-based estimation procedure for GSMA. Instead of estimating identity by descent between markers, as performed in interval mapping approaches, we estimated directly the nonparametric linkage score between markers using a kernel procedure. The GSMA is then extended to take into account the kernel estimate of the nonparametric linkage score and its variance at a given chromosomal position. The method is applied to the rheumatoid arthritis genome scan data (Genetic Analysis Workshop 15 Problem 2).

Identifying single-nucleotide polymorphisms responsible for the linkage signal of rheumatoid arthritis on chromosome 6 by joint modeling of linkage and association

This study evaluated the utility of unrelated controls and flanking markers when performing joint modeling of linkage and association by the LAMP software (version 0.0.6) [Am J Hum Genet 2005, 76:934–949; Am J Hum Genet 2006, 78:778–792]. Analyses were conducted on the simulated rheumatoid arthritis (RA) data in Genetic Analysis Workshop 15 (GAW15), using single-nucleotide polymorphisms (SNPs) on chromosome 6 over the 100 simulated replicates. We found that the LOD score for testing association in the presence of linkage dramatically increased when unrelated controls were added to affected sib pairs (ASPs), and that choosing a sufficient number of flanking markers is critical in order to distinguish between perfect linkage disequilibrium (which leads to the conclusion of a measured SNP explaining a linkage signal) and incomplete linkage disequilibrium (which leads to the conclusion of other undetected causal variants in a linkage region).

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.

Linkage analyses of rheumatoid arthritis and related quantitative phenotypes: the GAW15 experience

The group that formed on the theme of linkage analyses of rheumatoid arthritis RA and related phenotypes (Group 10) in the Genetic Analysis Workshop 15 comprised 18 sets of investigators. Two data sets were available: one was a real set provided by the North American Rheumatoid Arthritis Consortium and collaborators in Canada, France (European Consortium Of Rheumatoid Arthritis Families) and the UK; the other was a simulated data set modelled after the real data set. Whereas a majority of the investigators analyzed the RA affection status as a binary phenotype, a few contributions considered data on correlated quantitative traits such as anti-cyclic citrullinated peptide and rheumatoid factor-immunoglobulin M. The different investigators applied a wide spectrum of linkage methods. As expected, most methods could identify the human leukocyfeantigen region on chromosome 6 as a major genetic factor for RA. In addition, some novel chromosomal regions provided significant evidence of linkage in multiple contributions in the group. In this report, we discuss the different strategies explored by the different investigators with the common goal of improving the power to detect linkage.

Genetic factors in systemic sclerosis

A number of genetic loci have been identified that appear to be associated with systemic sclerosis (SSc; scleroderma). There is mounting evidence suggesting that these genetic associations may in fact be associated with distinct phenotypes in SSc based on autoantibody pattern rather than with SSc as a single disease entity. This may ultimately have implications for approaches to therapy as well as responses to therapy. The most promising candidate genes are those involved in pathways that lead to the vascular damage and fibrosis that are the hallmarks of this disease. There is uncertainty, however, regarding the nature of the key pathological mechanisms that link these two disease processes. Recent studies have focused on Fli1 (friend leukaemia integration 1), a transcription factor that is found in immune cells, fibroblasts, and endothelial cells that regulates collagen gene function and angiogenesis. Fli1 is dysregulated in SSc skin and dermal blood vessels, and appears to play a pathological role in SSc skin fibrosis and vessel degeneration. Whether this dysregulation is due to genetic polymorphisms in the Fli1 pathway or to epigenetic mechanisms is not clear.

TRAF1-C5 as a risk locus for rheumatoid arthritis--a genomewide study

BACKGROUND

Rheumatoid arthritis has a complex mode of inheritance. Although HLA-DRB1 and PTPN22 are well-established susceptibility loci, other genes that confer a modest level of risk have been identified recently. We carried out a genomewide association analysis to identify additional genetic loci associated with an increased risk of rheumatoid arthritis.

METHODS

We genotyped 317,503 single-nucleotide polymorphisms (SNPs) in a combined case-control study of 1522 case subjects with rheumatoid arthritis and 1850 matched control subjects. The patients were seropositive for autoantibodies against cyclic citrullinated peptide (CCP). We obtained samples from two data sets, the North American Rheumatoid Arthritis Consortium (NARAC) and the Swedish Epidemiological Investigation of Rheumatoid Arthritis (EIRA). Results from NARAC and EIRA for 297,086 SNPs that passed quality-control filters were combined with the use of Cochran-Mantel-Haenszel stratified analysis. SNPs showing a significant association with disease (P<1x10(-8)) were genotyped in an independent set of case subjects with anti-CCP-positive rheumatoid arthritis (485 from NARAC and 512 from EIRA) and in control subjects (1282 from NARAC and 495 from EIRA).

RESULTS

We observed associations between disease and variants in the major-histocompatibility-complex locus, in PTPN22, and in a SNP (rs3761847) on chromosome 9 for all samples tested, the latter with an odds ratio of 1.32 (95% confidence interval, 1.23 to 1.42; P=4x10(-14)). The SNP is in linkage disequilibrium with two genes relevant to chronic inflammation: TRAF1 (encoding tumor necrosis factor receptor-associated factor 1) and C5 (encoding complement component 5).

CONCLUSIONS

A common genetic variant at the TRAF1-C5 locus on chromosome 9 is associated with an increased risk of anti-CCP-positive rheumatoid arthritis.

STAT4 and the risk of rheumatoid arthritis and systemic lupus erythematosus

BACKGROUND

Rheumatoid arthritis is a chronic inflammatory disease with a substantial genetic component. Susceptibility to disease has been linked with a region on chromosome 2q.

METHODS

We tested single-nucleotide polymorphisms (SNPs) in and around 13 candidate genes within the previously linked chromosome 2q region for association with rheumatoid arthritis. We then performed fine mapping of the STAT1-STAT4 region in a total of 1620 case patients with established rheumatoid arthritis and 2635 controls, all from North America. Implicated SNPs were further tested in an independent case-control series of 1529 patients with early rheumatoid arthritis and 881 controls, all from Sweden, and in a total of 1039 case patients and 1248 controls from three series of patients with systemic lupus erythematosus.

RESULTS

A SNP haplotype in the third intron of STAT4 was associated with susceptibility to both rheumatoid arthritis and systemic lupus erythematosus. The minor alleles of the haplotype-defining SNPs were present in 27% of chromosomes of patients with established rheumatoid arthritis, as compared with 22% of those of controls (for the SNP rs7574865, P=2.81x10(-7); odds ratio for having the risk allele in chromosomes of patients vs. those of controls, 1.32). The association was replicated in Swedish patients with recent-onset rheumatoid arthritis (P=0.02) and matched controls. The haplotype marked by rs7574865 was strongly associated with lupus, being present on 31% of chromosomes of case patients and 22% of those of controls (P=1.87x10(-9); odds ratio for having the risk allele in chromosomes of patients vs. those of controls, 1.55). Homozygosity of the risk allele, as compared with absence of the allele, was associated with a more than doubled risk for lupus and a 60% increased risk for rheumatoid arthritis.

CONCLUSIONS

A haplotype of STAT4 is associated with increased risk for both rheumatoid arthritis and systemic lupus erythematosus, suggesting a shared pathway for these illnesses.

Meta-analysis of HLA-DRB1 polymorphism in Latin American patients with rheumatoid arthritis

OBJECTIVES

To estimate the common effect size of HLA-DRB1 alleles on rheumatoid arthritis (RA) susceptibility across Latin America populations through a meta-analysis combining the results of published data.

METHODS

Case-control studies on HLA-DRB1 association with RA in Latin America were searched up to October 2006. Genotype frequencies were extracted according to both shared epitope (SE) and HLA-DR4 positive or negative alleles. The effect summary odds ratio (OR) and 95% confidence intervals was obtained. Heterogeneity and publication bias were assessed.

RESULTS

Eight studies containing 684 cases and 1015 controls were included. Under the random effects model, the common OR was 3.28 (1.93, 5.60) (p<0.0001) and 3.54 (2.47, 5.05) (p=4.22 x 10(-12)) for HLA-DR4 and SE, respectively. There was no evidence of publication bias according to Funnel plot and Egger's regression test (p=0,445 for DR4 and p=0,464 for SE meta-analysis). Significant heterogeneity was observed for HLA-DR4 (I2=81.06%, Q=36.96, p=0.000005) but not for the SE meta-analysis.

CONCLUSIONS

HLA-DR4 and SE positive HLA-DRB1 alleles (mainly HLA-DRB10404) are associated with RA in Latin Americans. Heterogeneity is expected owing to the diverse degree of admixture between the examined populations. Our findings support the HLA as a major susceptibility locus for RA and validate the SE hypothesis in Latin America.

Cia25 on rat chromosome 12 regulates severity of autoimmune arthritis induced with pristane and with collagen

BACKGROUND

A genomewide scan in a DA x ACI F2 intercross studied for collagen-induced arthritis (CIA) identified the severity quantitative trait locus Cia25 on rat chromosome 12. Cia25 co-localises with loci regulating several forms of autoimmune diseases in rats, mice and humans, suggesting a common gene.

OBJECTIVE

To characterise the effects of Cia25 on severity of arthritis in congenic rats.

METHODS

DA.ACI(Cia25) congenic rats were constructed according to a genotype-guided strategy, and tested for pristane-induced arthritis (PIA) and CIA, induced with rat type II collagen (CII). A well-established scoring system previously shown to correlate with histological damage, including cartilage and bone erosions, synovial hyperplasia and synovial inflammation, was used.

RESULTS

The introgression of ACI alleles at Cia25 into DA background, as in DA.ACI(Cia25) rats, was enough to significantly reduce arthritis severity by 60% in PIA and by 40% in CIA, both in males and females compared with DA rats of the same sex. Levels of IgG anti-CII in male DA.ACI(Cia25) rats were 83% lower than in male DA. Levels of anti-CII in females were not affected by the congenic interval.

CONCLUSIONS

Cia25 contains a gene that regulates disease severity in two distinct models of autoimmune arthritis. Although both genders were protected in arthritis studies, only male congenic rats had a dramatic reduction in levels of anti-CII, suggesting the possibility of a second arthritis gene in this interval that operates via the regulation of autoantibodies in a sex-specific manner. The identification of the gene(s) accounting for Cia25 is expected to generate novel prognostic biomarkers and targets for therapy.

Genetic association between the PRKCH gene encoding protein kinase Ceta isozyme and rheumatoid arthritis in the Japanese population

OBJECTIVE

Analyses of families with rheumatoid arthritis (RA) have suggested the presence of a putative susceptibility locus on chromosome 14q21-23. This large population-based genetic association study was undertaken to examine this region.

METHODS

A 2-stage case-control association study of 950 unrelated Japanese patients with RA and 950 healthy controls was performed using >400 gene-based common single-nucleotide polymorphisms (SNPs).

RESULTS

Multiple SNPs in the PRKCH gene encoding the eta isozyme of protein kinase C (PKCeta) showed significant single-locus disease associations, the most significant being SNP c.427+8134C>T (odds ratio 0.72, 95% confidence interval 0.62-0.83, P = 5.9 x 10(-5)). Each RA-associated SNP was consistently mapped to 3 distinct regions of strong linkage disequilibrium (i.e., linkage disequilibrium or haplotype blocks) in the PRKCH gene locus, suggesting that multiple causal variants influence disease susceptibility. Significant SNPs included a novel common missense polymorphism of the PRKCH gene, V374I (rs2230500), which lies within the ATP-binding site that is highly conserved among PKC superfamily members. In circulating lymphocytes, PRKCH messenger RNA was expressed at higher levels in resting T cells (CD4(+) or CD8(+)) than in B cells (CD19(+)) or monocytes (CD14(+)) and was significantly down-regulated through immune responses.

CONCLUSION

Our results provide evidence of the involvement of PRKCH as a susceptibility gene for RA in the Japanese population. Dysregulation of PKCeta signal transduction pathway(s) may confer increased risk of RA through aberrant T cell-mediated autoimmune responses.

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.

Association of the microsatellite in the 3' untranslated region of the CD154 gene with rheumatoid arthritis in females from a Spanish cohort: a case-control study

CD40-CD154 interaction is an important mediator of inflammation and has been implicated in T helper type 1-mediated autoimmune diseases including rheumatoid arthritis (RA). Linkage studies have shown association of markers in the proximity of the CD154 gene. In the present work we investigated whether specific allele variants of the microsatellite in the 3' UTR of the CD154 gene might modulate the risk of RA. The study, in a case-control setting, included 189 patients and 150 healthy controls from the Canary Islands, Spain. The 24CAs allele was less represented in female patients than in controls (0.444 in controls versus 0.307 in patients, P = 0.006, odds ratio (OR) 0.556, 95% confidence interval (CI) 0.372 to 0.831) but not in males (0.414 versus 0.408), and only when homozygous (P = 0.012; OR 0.35, 95% CI 0.16 to 0.77). We also verified that CD154 association with RA was independent of human leukocyte antigen (HLA) phenotype. A further functional study showed that after stimulation anti-CD3, CD154 mRNA was more stable in CD4+ T lymphocytes from patients with RA bearing the 24CAs allele (mRNA half-life 208 minutes) than in patients without the 24CAs allele (109 minutes, P = 0.009). However, a lower percentage of CD154+CD4+ T lymphocytes was seen in freshly isolated peripheral blood mononuclear cells from patients carrying 24CAs alleles (mean 4.28 versus 8.12; P = 0.033), and also in CD4+ T lymphocytes stimulated with anti-CD3 (median 29.40 versus 47.60; P = 0.025). These results were concordant with the smaller amounts of CD154 mRNA isolated from stimulated T lymphocytes with 24CAs alleles. The CD154 microsatellite therefore seems to affect the expression of the gene in a complex manner that implies not only mRNA stability. These data suggest that the CD154 microsatellite contributes to the regulation of mRNA and protein expression, although further studies will be necessary to elucidate its role in disease predisposition.