Genome-wide association study of determinants of anti-cyclic citrullinated peptide antibody titer in adults with rheumatoid arthritis

Deep immunophenotyping reveals circulating activated lymphocytes in individuals at risk for rheumatoid arthritis

Rheumatoid arthritis (RA) is a systemic autoimmune disease with currently no universally highly effective prevention strategies. Identifying pathogenic immune phenotypes in 'At-Risk' populations prior to clinical disease onset is crucial to establishing effective prevention strategies. Here, we applied mass cytometry to deeply characterize the immunophenotypes in blood from At-Risk individuals identified through the presence of serum antibodies to citrullinated protein antigens (ACPA) and/or first-degree relative (FDR) status (n=52), as compared to established RA (n=67), and healthy controls (n=48). We identified significant cell expansions in At-Risk individuals compared with controls, including CCR2+CD4+ T cells, T peripheral helper (Tph) cells, type 1 T helper cells, and CXCR5+CD8+ T cells. We also found that CD15+ classical monocytes were specifically expanded in ACPA-negative FDRs, and an activated PAX5 low naïve B cell population was expanded in ACPA-positive FDRs. Further, we developed an "RA immunophenotype score" classification method based on the degree of enrichment of cell states relevant to established RA patients. This score significantly distinguished At-Risk individuals from controls. In all, we systematically identified activated lymphocyte phenotypes in At-Risk individuals, along with immunophenotypic differences among both ACPA+ and ACPA-FDR At-Risk subpopulations. Our classification model provides a promising approach for understanding RA pathogenesis with the goal to further improve prevention strategies and identify novel therapeutic targets.

'RA and the microbiome: do host genetic factors provide the link?

Rheumatoid arthritis (RA) is a chronic autoimmune disease, characterised by painful synovium inflammation, bony erosions, immune activation and the circulation of autoantibodies. Despite recent advances in therapeutics enabling disease suppression, there is a considerable demand for alternative therapeutic strategies as well as optimising those available at present. The relatively low concordance rate between monozygotic twins, 20–30% contrasts with heritability estimates of ∼65%, indicating a substantive role of other risk factors in RA pathogenesis. There is established evidence that RA has an infective component to its aetiology. More recently, differences in the commensal microbiota in RA compared to controls have been identified. Studies have shown that the gut, oral and lung microbiota is different in new onset treatment naïve, and established RA patients, compared to controls. Key taxonomic associations are an increase in abundance of Porphyromonas gingivalis and Prevotella copri in RA patients, compared to healthy controls. Host genetics may provide the link between disease and the microbiome. Genetic influence may be mediated by the host immune system; a differential response to RA associated taxa is suggested. The gut microbiome contains elements which are as much as 30% heritable. A better understanding of the influence of host genetics will shed light onto the role of the microbiome in RA. Here we review the role of the microbiome in RA through the lens of host genetics, and consider future research areas addressing microbiome study design and bioinformatics approaches.

•

Rheumatoid arthritis (RA) affects 1% of the population and is highly debilitating.

•

RA is ~65% heritable, yet the concordance rate between monozygotic twins is just 20–30%, indicating a substantive role of other risk factors.

•

Studies have shown that the gut, oral and lung microbiome is different in treatment naïve and established RA patients, compared to controls.

•

Current findings suggest an important influence of host genetics on the microbiome, which may contribute to RA via the host immune system.

•

Associations of the microbiome with RA described thus far are confounded by host genetics, and future studies need to take account of this.

HLA-DRB1 Analysis Identified a Genetically Unique Subset within Rheumatoid Arthritis and Distinct Genetic Background of Rheumatoid Factor Levels from Anticyclic Citrullinated Peptide Antibodies

OBJECTIVE

HLA-DRB1 is the most important locus associated with rheumatoid arthritis (RA) and anticitrullinated protein antibodies (ACPA). However, fluctuations of rheumatoid factor (RF) over the disease course have made it difficult to define fine subgroups according to consistent RF positivity for the analyses of genetic background and the levels of RF.

METHODS

A total of 2873 patients with RA and 2008 healthy controls were recruited. We genotyped HLA-DRB1 alleles for the participants and collected consecutive data of RF in the case subjects. In addition to RF+ and RF- subsets, we classified the RF+ subjects into group 1 (constant RF+) and group 2 (seroconversion). We compared HLA-DRB1 alleles between the RA subsets and controls and performed linear regression analysis to identify HLA-DRB1 alleles associated with maximal RF levels. Omnibus tests were conducted to assess important amino acid positions.

RESULTS

RF positivity was 88%, and 1372 and 970 RF+ subjects were classified into groups 1 and 2, respectively. RF+ and RF- showed similar genetic associations to ACPA+ and ACPA- RA, respectively. We found that shared epitope (SE) was more enriched in group 2 than 1, p = 2.0 × 10^-5, and that amino acid position 11 showed a significant association between 1 and 2, p = 2.7 × 10^-5. These associations were independent of ACPA positivity. SE showed a tendency to be negatively correlated with RF titer (p = 0.012). HLA-DRB1*09:01, which reduces ACPA titer, was not associated with RF levels (p = 0.70).

CONCLUSION

The seroconversion group was shown to have distinct genetic characteristics. The genetic architecture of RF levels is different from that of ACPA.

Effects of Amerindian Genetic Ancestry on Clinical Variables and Therapy in Patients with Rheumatoid Arthritis

OBJECTIVE

To define whether Amerindian genetic ancestry correlates with clinical and therapeutic variables in admixed individuals with rheumatoid arthritis (RA) from Latin America.

METHODS

Patients with RA (n = 1347) and healthy controls (n = 1012) from Argentina, Mexico, Chile, and Peru were included. Samples were genotyped for the Immunochip v1 using the Illumina platform. Clinical data were obtained through interviews or the clinical history.

RESULTS

Percentage of Amerindian ancestry was comparable between cases and controls. Morning stiffness (p < 0.0001, OR 0.05), rheumatoid factor (RF; p < 0.0001, OR 0.22), radiographic changes (p < 0.0001, OR 0.05), and higher number of criteria were associated with lower Amerindian ancestry after Bonferroni correction. Higher Amerindian ancestry correlated only with weight loss (p_Bonferroni < 0.0001, OR 2.85). Increased Amerindian ancestry correlated with higher doses of azathioprine (p < 0.0001, OR 163.6) and sulfasalazine (p < 0.0001, OR 48.6), and inversely with methotrexate (p = 0.001, OR 0.35), leflunomide (p = 0.001, OR 0.16), and nonsteroidal antiinflammatory drugs (p_Bonferroni = 0.001, OR 0.37). Only the presence of RF and weight loss were modified after confounders adjustment.

CONCLUSION

Amerindian ancestry protects against most major clinical criteria of RA, but regarding the association of RF with increased European ancestry, age, sex, and smoking are modifiers. Ancestry also correlates with the therapeutic profiles.

Genetic and environmental risk factors for rheumatoid arthritis

Multiple genetic and environmental factors have been associated with an increased risk for rheumatoid arthritis (RA). Of these, the strongest associations have been seen with female sex, a family history of RA, the genetic factor the "shared epitope," and exposure to tobacco smoke. There is also renewed interest in mucosal inflammation and microbial factors as contributors to the development of RA. However, the identification of a "preclinical" period of RA that can be defined as local or systemic autoimmunity as measured by autoantibodies and other biomarkers prior to the development of clinically apparent synovitis suggests that the risk factors for RA are acting long prior to first clinical evidence of IA. As such, a major challenge to the field will be to investigate the full spectrum of the development of RA, from initiation and propagation of autoimmunity during preclinical RA and transition to clinically apparent synovitis and classifiable RA, to determine which genetic and environmental factors are important at each stage of disease development. Understanding the exact role and timing of action of risk factors for RA is especially important given the advent of prevention trials in RA, and the hope that a full understanding of genetic and environmental factors in RA could lead to effective preventive interventions.

Interactions between amino acid-defined major histocompatibility complex class II variants and smoking in seropositive rheumatoid arthritis

OBJECTIVE

To define the interaction between cigarette smoking and HLA polymorphisms in seropositive rheumatoid arthritis (RA), in the context of a recently identified amino acid-based HLA model for RA susceptibility.

METHODS

We imputed Immunochip data on HLA amino acids and classical alleles from 3 case-control studies (the Swedish Epidemiological Investigation of Rheumatoid Arthritis [EIRA] study [1,654 cases and 1,934 controls], the Nurses' Health Study [NHS] [229 cases and 360 controls], and the Korean RA Cohort Study [1,390 cases and 735 controls]). We examined the interaction effects of heavy smoking (>10 pack-years) and the genetic risk score (GRS) of multiple RA-associated amino acid positions (positions 11, 13, 71, and 74 in HLA-DRβ1, position 9 in HLA-B, and position 9 in HLA-DPβ1), as well as the interaction effects of heavy smoking and the GRS of HLA-DRβ1 4-amino acid haplotypes (assessed via attributable proportion due to interaction [AP] using the additive interaction model).

RESULTS

Heavy smoking and all investigated HLA amino acid positions and haplotypes were associated with RA susceptibility in the 3 populations. In the interaction analysis, we found a significant deviation from the expected additive joint effect between heavy smoking and the HLA-DRβ1 4-amino acid haplotype (AP 0.416, 0.467, and 0.796, in the EIRA, NHS, and Korean studies, respectively). We further identified the key interacting variants as being located at HLA-DRβ1 amino acid positions 11 and 13 but not at any of the other RA risk-associated amino acid positions. For residues in positions 11 and 13, there were similar patterns between RA risk effects and interaction effects.

CONCLUSION

Our findings of significant gene-environment interaction effects indicate that a physical interaction between citrullinated autoantigens produced by smoking and HLA-DR molecules is characterized by the HLA-DRβ1 4-amino acid haplotype, primarily by positions 11 and 13.

Genetic Susceptibility to Vitiligo: GWAS Approaches for Identifying Vitiligo Susceptibility Genes and Loci

Vitiligo is an autoimmune disease with a strong genetic component, characterized by areas of depigmented skin resulting from loss of epidermal melanocytes. Genetic factors are known to play key roles in vitiligo through discoveries in association studies and family studies. Previously, vitiligo susceptibility genes were mainly revealed through linkage analysis and candidate gene studies. Recently, our understanding of the genetic basis of vitiligo has been rapidly advancing through genome-wide association study (GWAS). More than 40 robust susceptible loci have been identified and confirmed to be associated with vitiligo by using GWAS. Most of these associated genes participate in important pathways involved in the pathogenesis of vitiligo. Many susceptible loci with unknown functions in the pathogenesis of vitiligo have also been identified, indicating that additional molecular mechanisms may contribute to the risk of developing vitiligo. In this review, we summarize the key loci that are of genome-wide significance, which have been shown to influence vitiligo risk. These genetic loci may help build the foundation for genetic diagnosis and personalize treatment for patients with vitiligo in the future. However, substantial additional studies, including gene-targeted and functional studies, are required to confirm the causality of the genetic variants and their biological relevance in the development of vitiligo.

An association between amino acid position 74 of HLA-DRB1 and anti-citrullinated protein antibody levels in Japanese patients with anti-citrullinated protein antibody-positive rheumatoid arthritis

OBJECTIVE

Anti-citrullinated protein antibodies (ACPAs) are highly specific to rheumatoid arthritis (RA), and strong associations between HLA-DRB1 alleles and ACPA levels have been detected in RA patients. We undertook this study to elucidate the associations between particular amino acid positions in HLA-DRB1 and ACPA levels in patients with RA.

METHODS

We analyzed ACPA data on a total of 4,371 Japanese ACPA-positive RA patients in whom HLA-DRB1 allele genotyping had been performed. Generalized linear regression analysis and omnibus testing were carried out to determine associations of HLA-DRB1 alleles, amino acid residues, or amino acid positions with levels of ACPA.

RESULTS

HLA-DRB1*09:01 and HLA-DR15 were confirmed to be associated with ACPA levels. HLA-DRB1*08:03 and DRB1*14:06 were associated with reduced and increased ACPA levels, respectively. We detected a strong association between ACPA levels and amino acid position 74 (P = 1.9 × 10(-51) ). The association was mainly conferred by alanine residue (P = 4.5 × 10(-51) ). After adjustment for position 74, amino acid positions 60 and 57 were found to be associated with ACPA levels. Amino acid positions 74 and 57 had previously been reported to be associated with susceptibility to ACPA-positive RA in Asians. Combinations of the amino acid residues at position 74 and position 60 or 57 could induce improvement in Akaike's information criterion comparable to that induced by the 5 significant HLA-DRB1 alleles (HLA-DRB1*08:03, DRB1*09:01, DRB1*14:06, DRB1*15:01, and DRB1*15:02).

CONCLUSION

Amino acid position 74 in HLA-DRB1 is strongly associated with ACPA levels in ACPA-positive RA, as well as with RA susceptibility. The mechanisms of ACPA production and susceptibility to ACPA-positive RA seem to partly overlap.

The relationship of antibodies to modified citrullinated vimentin and markers of bone and cartilage destruction in rheumatoid arthritis

Objective. To make individualised decisions regarding treatment is one of the most important challenges in clinical practise, and identification of sensitive and specific markers of prognosis is an important research question. The main objective of this study was to evaluate relationships between the level of autoantibodies, radiographic changes and laboratory markers of bone, and cartilage destruction. Methods. A total of 114 RA patients were examined. The serum concentration of IgM RF, antibodies to cyclic citrullinated peptide (anti-CCP), modified citrullinated vimentin (anti-MCV), matrix metalloproteinase 3 (MMP-3), and cartilage oligomeric matrix protein (COMP, ng/mL) were measured. The van der Heijde-modified Sharp Score was used to quantify the radiologic changes. Results. Among the patients who were high-positive for anti-MCV, the value of total modified Sharp score (mTSS) (96.5; 66-120) was higher as well as the joint space narrowing (82; 60.5-105.5), and a higher level of MMP-3 was recorded more frequently (56%) in comparison with negative/low-positive patients (57; 31-88, 50; 29-82, 31% resp., P < 0.05). The level of COMP was also higher among patients high-positive for anti-MCV (9.7; 8.1-13.1 and 6.8; 5.4-10.7, resp., P = 0.02). Conclusion. A high positive level of anti-MCV as contrasted with anti-CCP and IgM RF is associated with more pronounced destructive changes in the joints.

The influence of polygenic risk scores on heritability of anti-CCP level in RA

The objective of this study was to study genetic factors that influence quantitative anticyclic citrullinated peptide (anti-CCP) antibody levels in RA patients. We carried out a genome-wide association study (GWAS) meta-analysis using 1975 anti-CCP+ RA patients from three large cohorts, the Brigham Rheumatoid Arthritis Sequential Study (BRASS), North American Rheumatoid Arthritis Consortium (NARAC) and the Epidemiological Investigation of RA (EIRA). We also carried out a genome-wide complex trait analysis (GCTA) to estimate the heritability of anti-CCP levels. GWAS-meta-analysis showed that anti-CCP levels were most strongly associated with the human leukocyte antigen (HLA) region with a P-value of 2 × 10(-11) for rs1980493. There were 112 SNPs in this region that exceeded the genome-wide significance threshold of 5 × 10(-8), and all were in linkage disequilibrium (LD) with the HLA- DRB1*03 allele with LD r(2) in the range of 0.25-0.88. Suggestive novel associations outside of the HLA region were also observed for rs8063248 (near the GP2 gene) with a P-value of 3 × 10(-7). None of the known RA risk alleles (∼52 loci) were associated with anti-CCP level. Heritability analysis estimated that 44% of anti-CCP variation was attributable to genetic factors captured by GWAS variants. In summary, anti-CCP level is a heritable trait, and HLA-DR3 and GP2 are associated with lower anti-CCP levels.

Very high levels of anti-citrullinated protein antibodies are associated with HLA-DRB1*15 non-shared epitope allele in patients with rheumatoid arthritis

OBJECTIVE

Production of anti-citrullinated protein antibodies (ACPAs) is an important biomarker for rheumatoid arthritis (RA). We undertook this study to determine whether genetic factors (HLA-DRB1 alleles) are associated with extreme ACPA levels in individuals with ACPA-positive RA, and to ascertain whether there are any phenotypic characteristics associated with these subgroups of RA.

METHODS

HLA-DRB1 allelic groups were genotyped in 1,073 ACPA-positive RA patients from the Swedish Epidemiological Investigation of Rheumatoid Arthritis study. We found that 283 patients (26.4%) had high ACPA levels (defined as >1,500 units/ml using the Euro-Diagnostica anti-CCP2 test), while the rest of the patients had moderate ACPA levels and served as the comparison group. A replication group consisted of 235 RA patients.

RESULTS

No significant differences in baseline disease activity were observed between patients with high and those with moderate ACPA levels. However, the HLA-DRB1*15 allele was associated with high ACPA levels (P=0.0002). A similar trend was detected in HLA-DRB1*15-positive patients in the replication cohort, with meta-analysis of the discovery and replication cohorts demonstrating an overall effect of HLA-DRB1*15 on development of high ACPA levels in both the discovery and replication cohorts (P<0.0001 by Mantel-Haenszel test with a fixed-effects model).

CONCLUSION

Our data indicate that HLA-DRB1*15 may promote the production of high ACPA levels. Due to the high value of ACPA level scores in the 2010 American College of Rheumatology/European League Against Rheumatism classification criteria for RA, the presence of HLA-DRB1*15 may, at least in part, contribute to fulfilling the criteria for RA. This illustrates the complex nature of the genetic regulation of ACPA levels. Additional mechanistic studies of the regulation of ACPAs and ACPA-positive RA are pending.

Evidence for PTPN22 R620W polymorphism as the sole common risk variant for rheumatoid arthritis in the 1p13.2 region

OBJECTIVE

The PTPN22 rs2476601 genetic variant has been associated with rheumatoid arthritis (RA) and other autoimmune diseases. Some reports suggest that this single-nucleotide polymorphism (SNP) may not be the only causal variant in the region of PTPN22. Our aim was to identify new independent RA-associated common gene variants in the PTPN22 region.

METHODS

We analyzed Wellcome Trust Case-Control Consortium genome-wide association study data for associations in the 397.2 kb PTPN22 region and selected 9 associated SNP (with p < 5 × 10(-3)) for replication and dependence analysis. The replication cohorts comprised 2857 patients with RA and 2994 controls from Spain, Netherlands, and Norway.

RESULTS

We found that 6 of the 9 selected SNP were associated in the Spanish cohort. Of these, 4 were also associated in the Dutch and Norwegian cohorts, and all 6 were associated with RA in the combined analysis. Conditional analyses showed that none of these associations was independent of rs2476601.

CONCLUSION

The SNP rs2476601 located in the PTPN22 gene is the sole common genetic variant associated with RA in the 1p13.2 region, suggesting that neighbor genes of PTPN22 do not have a major influence in RA.

Genome-wide analysis identifies a quantitative trait locus in the MHC class II region associated with generalized vitiligo age of onset

Generalized vitiligo is a common autoimmune disease in which acquired patchy depigmentation of skin, hair, and mucous membranes results from loss of melanocytes from involved areas. Previous genetic analyses have focused on vitiligo susceptibility, and have identified a number of genes involved in disease risk. Age of onset of generalized vitiligo also involves a substantial genetic component, but has not previously been studied systematically. In this study, we report a genome-wide association study of vitiligo age of onset in 1,339 generalized vitiligo patients, with replication in an independent cohort of 677 cases. We identified a quantitative trait locus for vitiligo age of onset in the major histocompatibility complex (MHC) class II region, located near c6orf10-BTNL2 (rs7758128; P=8.14 × 10(-11)), a region that is also associated with generalized vitiligo susceptibility. In contrast, there was no association of vitiligo age of onset with any other MHC or non-MHC loci that are associated with vitiligo susceptibility. These findings highlight the differing roles played by genes involved in vitiligo susceptibility versus vitiligo age of onset, and illustrate that genome-wide analyses can be used to identify genes involved in quantitative aspects of disease natural history, as well as disease susceptibility per se.

Compilation of a comprehensive gene panel for systematic assessment of genes that govern an individual’s drug responses

Aims: Polymorphisms of genes involved in the pharmacokinetic and pharmacodynamic processes underlie the divergent drug responses among individuals. Despite some successes in identifying these polymorphisms, the candidate gene approach suffers from insufficient gene coverage whereas the genome-wide association approach is limited by less than ideal coverage of SNPs in some important genes. To expand the potential of the candidate approach, we aim to delineate a comprehensive network of drug-response genes for in-depth genetic studies. Materials & methods: Pharmacologically important genes were extracted from various sources including literatures and web resources. These genes, along with their homologs and regulatory miRNAs, were organized based on their pharmacological functions and weighted by literature evidence and confidence levels. Their coverage was evaluated by analyzing three commercial SNP chips commonly used for genome-wide association studies: Affymetrix SNP array 6.0, Illumina HumanHap1M and Illumina Omni. Results: A panel of drug-response genes was constructed, which contains 923 pharmacokinetic genes, 703 pharmacodynamic genes and 720 miRNAs. There are only 16.7% of these genes whose all known SNPs can be directly or indirectly (r2 > 0.8) captured by the SNP chips with coverage of more than 80%. This is possibly because these SNPs chips have notably poor performance over rare SNPs and miRNA genes. Conclusion: We have compiled a panel of candidate genes that may be pharmacologically important. Using this knowledgebase, we are able to systematically evaluate genes and their variants that govern an individual’s response to a given pharmaceutical therapy. This approach can serve as a necessary complement to genome-wide associations.

Using genetic and clinical data to understand response to disease-modifying anti-rheumatic drug therapy: data from the Brigham and Women's Hospital Rheumatoid Arthritis Sequential Study

The objective of this review is to report on the progress of the Brigham and Women's Hospital Rheumatoid Arthritis Sequential Study (BRASS) Registry data collection and summarize previous research in understanding therapeutic response to DMARDs using clinical and genetic data. The BRASS Registry, established in 2003, is a large, single-centre, prospective and observational cohort of 1100 RA patients. Patients with either new-onset or established RA disease are recruited from the practices of rheumatologists. Annual visits collect information on demographics, 28-joint DAS-CRP3 (DAS-28-CRP3), medication use, comorbidities and functional status (Modified Health Assessment Questionnaire, Short Form Health Survey 12). Two published studies have utilized BRASS to examine genetic predictors of treatment response. In a cross-sectional study, examining the association between candidate single nucleotide polymorphisms (SNPs) and disease activity in a subset of 120 RA patients on MTX monotherapy, the minor allele of ATIC rs4673993 was associated with low disease activity (P=0.01, DAS-28-CRP3≤3.2). In an international collaboration, 55 BRASS patients receiving anti-TNF therapy were genotyped for 31 SNPs associated with the risk of RA. With our collaborators, we discovered an SNP at the protein tyrosine phosphatase, receptor type, C (PTPRC) gene locus that was associated with EULAR 'good response'. With accurate data collection and the capacity to run genome-wide association studies and SNP analyses, the BRASS Registry has the ability to determine the contribution of genetic variants to disease onset and to assess their usefulness as biomarkers for treatment response and drug toxicity.

The rheumatoid arthritis HLA-DRB1 shared epitope

PURPOSE OF REVIEW

To update progress made between December 2008 and November 2009 on the role of the rheumatoid arthritis (RA)-shared epitope in the cause and pathogenesis of RA.

RECENT FINDINGS

New evidence has been recently presented to suggest that noninherited human leukocyte antigens (HLAs) originating through pregnancy or exposure to maternal antigens in utero could contribute to RA development in shared epitope-negative women. An interaction between smoking and shared epitope-coding non-*04 HLA-DRB1 alleles (particularly HLA-DRB1*01 and HLA-DRB1*10) was formally established for the first time. Progress has been made in determining the relative contributions and the interaction of the shared epitope, PTPN22 and smoking in conferring the risk of anticitrullinated protein antibodies-positive and negative RA. The autoantigen that anticitrullinated protein antibodies recognize in a significant number of RA patients has been identified as citrullinated alpha-enolase and the importance of genetic factors in anticitrullinated protein antibodies-negative RA has been highlighted. Additionally, associations of RA risk with several new genetic markers have been reported. Among them: two new major histocompatibility complex, non-DRB1 loci, a polymorphism marker in major histocompatibility complex class I polypeptide-related sequence A, an allele of the Fcgamma receptor, a polymorphism marker in the beta2-adrenergic receptor and a low-inducible allele of the cytochrome P450 subtype 1A2.

SUMMARY

Although the mechanistic basis of shared epitope-RA association remains an enigma, observations made during the last year shed new light on the conditions in which the shared epitope - alone or in combination with other genes or environmental factors - affects the risk of RA and the phenotype of the disease.

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.