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

Supervised machine learning and logistic regression identifies novel epistatic risk factors with PTPN22 for rheumatoid arthritis

Investigating genetic interactions (epistasis) has proven difficult despite the recent advances of both laboratory methods and statistical developments. With no 'best' statistical approach available, combining several analytical methods may be optimal for detecting epistatic interactions. Using a multi-stage analysis that incorporated supervised machine learning and methods of association testing, we investigated epistatic interactions with a well-established genetic factor (PTPN22 1858T) in a complex autoimmune disease (rheumatoid arthritis (RA)). Our analysis consisted of four principal stages: Stage I (data reduction)-identifying candidate chromosomal regions in 292 affected sibling pairs, by predicting PTPN22 concordance using multipoint identity-by-descent probabilities and a supervised machine learning algorithm (Random Forests); Stage II (extension analysis)-testing detailed genetic data within candidate chromosomal regions for epistasis with PTPN22 1858T in 677 cases and 750 controls using logistic regression; Stage III (replication analysis)-confirmation of epistatic interactions in 947 cases and 1756 controls; Stage IV (combined analysis)-a pooled analysis including all 1624 RA cases and 2506 control subjects for final estimates of effect size. A total of seven replicating epistatic interactions were identified. SNP variants within CDH13, MYO3A, CEP72 and near WFDC1 showed significant evidence for interaction with PTPN22, affecting susceptibility to RA.

Linkage and association studies of joint morbidity from rheumatoid arthritis

OBJECTIVE

To investigate the relationship between genetic variations of rheumatoid arthritis (RA) susceptibility in terms of joint morbidity.

METHODS

We used data from Genetic Analysis Workshop 15. The Illumina linkage panel IV included 5858 single-nucleotide polymorphisms (SNP), with 5744 SNP passing quality control filters. The phenotypic variables analyzed were the level of rheumatoid factor (RF) and score on the Joint Alignment and Motion (JAM) scale. We modified the scale, dividing by RF values relevant to disease severity. Linkage analysis for affected sibling pairs was done using the MERLIN program, and family-based association tests were carried out using PLINK and FBAT software.

RESULTS

We found a high peak (LOD = 3.29; NPL Z = 4.07) near the HLA-DRB1 region on chromosome 6. The linkage at 6p24 at rs1410766 [LOD = 2.66; nonparametric linkage (NPL) Z = 3.23] was statistically significant. Two other regions also showed possible linkage peaks: chromosome 7q30 at rs322812 (LOD = 2.47; NPL Z = 3.39) and chromosome 15p34 at rs347117 (LOD = 1.95; NPL Z = 2.80). For the family-based association study, 7 SNP related to clinical RA severity were detected.

CONCLUSION

Genetic variations may lead to an enhanced risk of joint damage and increased levels of RF. Further studies are needed to elucidate the roles of other genes involved in RA and to explore whether the clinical signs of RA are associated with particular genetic variations.

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.

Vra4 congenic rats with allelic differences in the class II transactivator gene display altered susceptibility to experimental autoimmune encephalomyelitis

Presentation of Ag bound to MHC class II (MHC II) molecules to CD4+ T cells is a key event in adaptive immune responses. Genetic differences in MHC II expression in the rat CNS were recently positioned to allelic variability in the CIITA gene (Mhc2ta), located within the Vra4 locus on rat chromosome 10. In this study, we have examined reciprocal Vra4-congenic strains on the DA and PVGav1 backgrounds, respectively. After experimental nerve injury the strain-specific MHC II expression on microglia was reversed in the congenic strains. Similar findings were obtained after intraparenchymal injection of IFN-gamma in the brain. Expression of MHC class II was also lower on B cells and dendritic cells from the DA.PVGav1-Vra4- congenic strain compared with DA rats after in vitro stimulation with IFN-gamma. We next explored whether Vra4 may affect the outcome of experimental autoimmune disease. In experimental autoimmune encephalomyelitis induced by immunization with myelin oligodendrocyte glycoprotein, DA.PVGav1-Vra4 rats displayed a lower disease incidence and milder disease course compared with DA, whereas both PVGav1 and PVGav1.DA-Vra4 rats were completely protected. These results demonstrate that naturally occurring allelic differences in Mhc2ta have profound effects on the quantity of MHC II expression in the CNS and on immune cells and that this genetic variability also modulates susceptibility to autoimmune neuroinflammation.

A two-dimensional genome scan for rheumatoid arthritis susceptibility loci

We performed a genome-wide search for pairs of susceptibility loci that jointly contribute to rheumatoid arthritis in families recruited by the North American Rheumatoid Arthritis Consortium. A complete two-dimensional (2D) non-parametric linkage scan was carried out using 380 autosomal microsatellite markers in 511 families. At each 2D peak we obtained the most likely underlying genetic model explaining the two-locus effects, defining epistasis as a departure from an additive or a multiplicative two-locus penetrance function. The highest peak in the surface identified an epistatic interaction between loci 6p21 and 16p12 (two-locus lod score = 18.02, epistasis P < 0.012). Significant and suggestive two-locus effects were also obtained for region 6p21 in combination with loci 18q21, 8p23, 1q41, and 6p22, while the highest 2D peaks excluding region 6p21 were observed at locus pairs 8p23-18q21 and 1p21-18q21. The 2D peaks were further examined using combined microsatellite and single-nucleotide polymorphism (SNP) marker genotypes in 744 families. The two-locus evidence for linkage increased for region pairs 6p21-18q12, 6p21-16p12, 6p21-8p23, 1q41-6p21, and 6p21-6p22, but decreased for pairs of regions that did not include locus 6p21. In conclusion, we obtained evidence for multi-locus interactions in rheumatoid arthritis that are mediated by the major susceptibility locus at 6p21.

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.

The interferon induced with helicase domain 1 A946T polymorphism is not associated with rheumatoid arthritis

An important feature of autoimmune diseases is the overlap of pathophysiological characteristics. Clustering of autoimmune diseases in families suggests that genetic variants may contribute to autoimmunity. The aim of the present study was to investigate the role of the interferon induced with helicase domain 1 (IFIH1) A946T (rs1990760 A>G) variant in rheumatoid arthritis (RA), as this was recently associated with susceptibility to type 1 diabetes. A total of 965 Caucasians with RA and 988 healthy controls were genotyped for IFIH1 A946T. Gene expression of IFIH1 was measured in peripheral blood leukocytes using real-time PCR. Genotypes were equally distributed in both RA cases and healthy controls (odds ratio for allele C = 0.9, 95% confidence interval = 0.8-1.0, P = 0.3). No association was detected after stratification by sex, age at onset, rheumatoid factor status, anti-cyclic citrullinated peptide status or radiological joint damage. Levels of IFIH1 mRNA were approximately twofold higher in blood leucocytes of RA cases compared with healthy controls (P < 0.0001). These results indicate that the IFIH1 is upregulated in RA but that the A946T variant does not contribute significantly to the genetic background of RA.

Contribution of genetic studies in rodent models of autoimmune arthritis to understanding and treatment of rheumatoid arthritis

Rheumatoid arthritis (RA) is a chronic and potentially debilitating autoimmune disease. While novel therapies have emerged in recent years, disease remission is rarely achieved. RA is a complex trait, and the identifying of its susceptibility and severity genes has been anticipated to generate new targets for therapeutic intervention. However, finding those genes and understanding their function has been a challenging task. Studies in rodent intercrosses and congenics generated from inbred strains have been an important complementary strategy to identify arthritis genes, and understand how they operate to regulate disease. Furthermore, these new rodent arthritis genes will be new targets for therapeutic interventions, and will identify new candidate genes or candidate pathways for association studies in RA. In this review-opinion article I discuss RA genetics, difficulties involved in gene identification, and how rodent models can facilitate (1) the discovery of both arthritis susceptibility and severity genes, (2) studies of gene-environment interactions, (3) studies of gene-gender interactions, (4) epistasis, (5) functional characterization of the specific genes, (6) development of novel therapies and (7) how the information generated from rodent studies will be useful to understanding and potentially treating RA.

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.

The Contribution of Methotrexate Exposure and Host Factors on Transcriptional Variance in Human Liver

Long-term administration of methotrexate (MTX) for management of chronic inflammatory diseases is associated with risk of liver damage. In this study, we examined the transcriptional profiles of livers from patients treated with MTX. The possibility that expression signatures correlate with grade of fibrosis or underlying rheumatic disease was evaluated. Twenty-seven patients taking MTX were accrued for this study. Ten non-MTX-exposed normal liver specimens were used as controls. Global mRNA expression was assayed using oligonucleotide arrays. A total of 205 genes were significantly altered in MTX-exposed livers. Six of these genes were validated by qPCR. Two genes, CLN8 and ANKH that map to chromosomal locations previously associated with rheumatoid arthritis, were found to be elevated in MTX-exposed samples. Subsequent pathway analysis indicates that MTX exposure is associated with the following key alterations: (1) upregulation of lipid biosynthetic genes, consistent with MTX-induced steatosis, (2) downregulation of proinflammatory chemokines, consistent with the anti-inflammatory effects of MTX, and (3) elevation of complement pathway gene expression. Complement 5, shown earlier to be correlated with liver fibrosis in mice, was found to be elevated (twofold) in MTX-exposed livers. In conclusion, we have found the expression of a number of genes associated with rheumatic disease and/or MTX exposure to be significantly different. Differences in complement expression provide the rationale for future correlative studies between MTX-induced liver fibrosis and C5 alleles in order to identify patients with increased risk for fibrosis.

Nature's choice of genes controlling chronic inflammation

Inflammation is a physiological response that may go uncontrolled and thereby develop in a chronic way. This seems to happen in many common diseases of autoimmune, degenerative, or allergic character. Rheumatoid arthritis (RA) is by definition a chronic disease with an autoimmune inflammatory attack on diarthrodial cartilaginous joints. The development of new treatment neutralizing cytokines involved in the inflammatory attack has given relief and gives the promise of more effective treatment of already established disease. It is now time to set our eyes on a new vision to develop preventive and curative treatment based on knowledge of the unique and causative pathogenic mechanisms. To do this we believe it is important to identify the natural-selected polymorphisms that are associated with disease. These have proven to be extremely difficult to identify in complex diseases such as RA, but using animal models, this work is closer to reality. Animal models have recently been developed mimicking various aspects of the human disease. We will present an example in which a genetic polymorphism associated with the development of arthritis has been identified. On the basis of this finding, a new pathway involving control of immune tolerance by reactive oxidative species has been identified and a new class of antiinflammatory agents activating the induced oxidative burst protein complex is suggested.

Recent advances in the genetics of rheumatoid arthritis

Recent progress in defining the role of genetic factors in rheumatoid arthritis (RA) has been remarkable. Anticyclic citrullinated peptide (anti-CCP) antibody-positive disease appears to be immunogenetically distinct from anti-CCP-negative disease, with the former subgroup primarily responsible for association and linkage with the HLA-DRB1 shared epitope (SE). There is preliminary evidence that non-HLA genes contribute differentially to anti-CCP-positive and negative disease. The phenotypic differences evident in anti-CCP-positive and negative disease suggest a need to reclassify RA based on the presence or absence of this autoantibody. Some recent work also suggests marked interactions between cigarette smoking, anti-CCP antibodies, and the SE, though these relationships may vary across populations. Lastly, a recent single nucleotide polymorphism-based genome-wide linkage analysis of multicase RA families revealed novel genomic regions that likely contain genes that predispose to RA or more specific phenotypes.

High-density SNP analysis of 642 Caucasian families with rheumatoid arthritis identifies two new linkage regions on 11p12 and 2q33

We have completed a genome wide linkage scan using >5700 informative single-nucleotide polymorphism (SNP) markers (Illumina IV SNP linkage panel) in 642 Caucasian families containing affected sibling pairs with rheumatoid arthritis (RA), ascertained by the North American Rheumatoid Arthritis Consortium. The results show striking new evidence of linkage at chromosomes 2q33 and 11p12 with logarithm of odds (LOD) scores of 3.52 and 3.09, respectively. In addition to a strong and broad linkage interval surrounding the major histocompatibility complex (LOD>16), regions with LOD>2.5 were observed on chromosomes 5 and 10. Additional linkage evidence (LOD scores between 1.46 and 2.35) was also observed on chromosomes 4, 7, 12, 16 and 18. This new evidence for multiple regions of genetic linkage is partly explained by the significantly increased information content of the Illumina IV SNP linkage panel (75.6%) compared with a standard microsatellite linkage panel utilized previously (mean 52.6%). Stratified analyses according to whether or not the sibling pair members showed elevated anticyclic citrullinated peptide titers indicates significant variation in evidence for linkage among strata on chromosomes 4, 5, 6 and 7. Overall, these new linkage data should reinvigorate efforts to utilize positional information to identify susceptibility genes for RA.