Dense Genotyping of Immune-Related Regions Identifies Loci for Rheumatoid Arthritis Risk and Damage in African Americans.

Over 100 risk loci for rheumatoid arthritis (RA) have been identified in individuals of European and Asian descent, but the genetic basis for RA in African Americans is less well understood. We genotyped 610 African Americans with autoantibody positive RA and 933 African American controls on the ImmunoChip (iChip) array. Using multivariable regression we evaluated the association between iChip markers and the risk of RA and radiographic severity. The single nucleotide polymorphism (SNP) rs1964995 (OR = 1.97, p = 1.28 × 10-15) near HLA-DRB1 was the most strongly associated risk SNP for RA susceptibility; SNPs in AFF3, TNFSF11, and TNFSF18 loci were suggestively associated (10-4 < p < 3.1 × 10-6). Trans-ethnic fine mapping of AFF3 identified a 90% credible set containing previously studied variants including rs9653442, rs7608424, and rs6712515 as well as the novel candidate variant rs11681966; several of these likely influence AFF3 gene expression level. Variants in TNFRSF9, CTLA4, IL2RA, C5/TRAF1, and ETS1 - but no variants within the major histocompatibility complex - were associated with RA radiographic severity. Conditional regression and pairwise linkage disequilibrium (LD) analyses suggest that additional pathogenic variants may be found in ETS1 and IL2RA beyond those found in other ethnicities. In summary, we use the dense genotyping of the iChip array and unique LD structure of African Americans to validate known risk loci for RA susceptibility and radiographic severity, and to better characterize the associations of AFF3, ETS1, and IL2RA.

Increased pretreatment serum IFN-ß/α ratio predicts non-response to tumour necrosis factor α inhibition in rheumatoid arthritis.

OBJECTIVE: Studies suggest that circulating type I interferon (IFN) may predict response to biological agents in rheumatoid arthritis (RA). Prediction of response prior to initiating therapy would represent a major advancement. METHODS: We studied sera from a test set of 32 patients with RA from the Auto-immune Biomarkers Collaborative Network Consortium and a validation set of 92 patients with RA from the Treatment Efficacy and Toxicity in Rheumatoid Arthritis Database and Repository registry. The test set included those with good response or no response to tumour necrosis factor (TNF) inhibitors at 14 weeks by European League Against Rheumatism criteria. The validation set included subjects with good, moderate or no response at 12 weeks. Total serum type I IFN activity, IFN-α and IFN-ß activity were measured using a functional reporter cell assay. RESULTS: In the test set, an increased ratio of IFN-ß to IFN-α (IFN-ß/α activity ratio) in pretreatment serum associated with lack of response to TNF inhibition (p=0.013). Anti-cyclic citrullinated peptide antibody titre and class of TNF inhibitor did not influence this relationship. A receiver-operator curve supported a ratio of 1.3 as the optimal cut-off. In the validation set, subjects with an IFN-ß/α activity ratio >1.3 were significantly more likely to have non-response than good response (OR=6.67, p=0.018). The test had 77% specificity and 45% sensitivity for prediction of non-response compared with moderate or good response. Meta-analysis of test and validation sets confirmed strong predictive capacity of IFN-ß/α activity ratio (p=0.005). CONCLUSIONS: Increased pretreatment serum IFN-ß/α ratio strongly associated with non-response to TNF inhibition. This study supports further investigation of serum type I IFN in predicting outcome of TNF inhibition in RA.

Identification of a genetic variant for joint damage progression in autoantibody-positive rheumatoid arthritis.

BACKGROUND: Joint destruction is a hallmark of autoantibody-positive rheumatoid arthritis (RA), though the severity is highly variable between patients. The processes underlying these interindividual differences are incompletely understood. METHODS: We performed a genome-wide association study on the radiological progression rate in 384 autoantibody-positive patients with RA. In stage-II 1557 X-rays of 301 Dutch autoantibody-positive patients with RA were studied and in stage-III 861 X-rays of 742 North American autoantibody-positive patients with RA. Sperm-Associated Antigen 16 (SPAG16) expression in RA synovium and fibroblast-like synoviocytes (FLS) was examined using Real-Time Quantitative Polymerase Chain Reaction (RT-qPCR) and immunohistochemistry. FLS secrete metalloproteinases that degrade cartilage and bone. SPAG16 genotypes were related to matrix metalloproteinase (MMP)-3 and MMP-1 expression by FLS in vitro and MMP-3 production ex vivo. RESULTS: A cluster of single nucleotide polymorphisms (SNPs) at 2q34, located at SPAG16, associated with the radiological progression rate; rs7607479 reached genome-wide significance. A protective role of rs7607479 was replicated in European and North American patients with RA. Per minor allele, patients had a 0.78-fold (95% CI 0.67 to 0.91) progression rate over 7 years. mRNA and protein expression of SPAG16 in RA synovium and FLS was verified. FLS carrying the minor allele secreted less MMP-3 (p=1.60×10(-2)). Furthermore, patients with RA carrying the minor allele had lower serum levels of MMP-3 (p=4.28×10(-2)). In a multivariate analysis on rs7607479 and MMP-3, only MMP-3 associated with progression (p=2.77×10(-4)), suggesting that the association between SPAG16-rs7607479 and joint damage is mediated via an effect on MMP-3 secretion. CONCLUSIONS: Genetic and functional analyses indicate that SPAG16 influences MMP-3 regulation and protects against joint destruction in autoantibody-positive RA. These findings could enhance risk stratification in autoantibody-positive RA.

Gene expression patterns in peripheral blood cells associated with radiographic severity in African Americans with early rheumatoid arthritis.

Gene expression profiling may be used to stratify patients by disease severity to test the hypothesis that variable disease outcome has a genetic component. In order to define unique expression signatures in African American rheumatoid arthritis (RA) patients with severe erosive disease, we undertook a gene expression study using samples of RNA from peripheral blood mononuclear cells (PBMCs). RNA from baseline PBMC samples of 96 African American RA patients with early RA (<2 years disease duration) was hybridized to cDNA probes of the Illumina Human HT-V3 expression array. Expression analyses were performed using the ca. 25,000 cDNA probes, and then expression levels were compared to the total number of erosions in radiographs of the hands and feet at baseline and 36 months. Using a false discovery rate cutoff of Q = 0.30, 1,138 genes at baseline and 680 genes at 36 months significantly correlated with total erosions. No evidence of a signal differentiating disease progression, or change in erosion scores between baseline and 36 months, was found. Further analyses demonstrated that the differential gene expression signature was localized to the patients with the most erosive disease (>10 erosions). Ingenuity Pathway Analysis demonstrated that genes with fold change greater than 1.5 implicated immune pathways such as CTLA signaling in cytotoxic T lymphocytes. These results demonstrate that CLEAR patients with early RA having the most severe erosive disease, as compared to more mild cases (<10 erosions), may be characterized by a set of differentially expressed genes that represent biological pathways with relevance to autoimmune disease.

Beliefs and attitudes towards participating in genetic research - a population based cross-sectional study.

BACKGROUND: Biobanks have the potential to offer a venue for chronic disease biomarker discovery, which would allow for disease early detection and for identification of carriers of a certain predictor biomarker. To assess the general attitudes towards genetic research and participation in biobanks in the Long Island/Queens area of New York, and what factors would predict a positive view of such research, participants from the NSLIJ hospital system were surveyed. METHODS: Participants were recruited at six hospital centers in the NSLIJ system during the summers of 2009 and again in 2011 (n = 1,041). Those who opted to participate were given a questionnaire containing 22 questions assessing demographics, lifestyle and attitudes towards genetic research. These questions addressed individual participant's beliefs about the importance of genetic research, willingness to participate in genetic research themselves, and their views on informed consent issues. RESULTS: Respondents took a generally positive view of genetic research in general, as well as their own participation in such research. Those with reservations were most likely to cite concerns over the privacy of their medical and genetic information. Those who were married tended to view genetic research as important, while those in the younger age group viewed it as less important. Prior blood donation of respondents was found to be a predictor of their approval for genetic research. Demographic factors were not found to be predictive of personal willingness to participate in genetic research, or of approval for the opt-out approach to consent. CONCLUSIONS: While respondents were generally inclined to approve of genetic research, and those who disapproved did not do so based on an underlying moral objection to such research, there is a disconnect between the belief in the importance of genetic research and the willingness of individuals to participate themselves. This indicates a continued concern for the ways in which genetic materials are safeguarded once they are collected. Also indicated was a general lack of understanding about the various consent processes that go along with genetic research, which should be addressed further to ensure the successful continuation of biobanks.

Rheumatoid arthritis risk allele PTPRC is also associated with response to anti-tumor necrosis factor alpha therapy.

OBJECTIVE: Anti-tumor necrosis factor alpha (anti-TNF) therapy is a mainstay of treatment in rheumatoid arthritis (RA). The aim of the present study was to test established RA genetic risk factors to determine whether the same alleles also influence the response to anti-TNF therapy. METHODS: A total of 1,283 RA patients receiving etanercept, infliximab, or adalimumab therapy were studied from among an international collaborative consortium of 9 different RA cohorts. The primary end point compared RA patients with a good treatment response according to the European League Against Rheumatism (EULAR) response criteria (n = 505) with RA patients considered to be nonresponders (n = 316). The secondary end point was the change from baseline in the level of disease activity according to the Disease Activity Score in 28 joints (triangle upDAS28). Clinical factors such as age, sex, and concomitant medications were tested as possible correlates of treatment response. Thirty-one single-nucleotide polymorphisms (SNPs) associated with the risk of RA were genotyped and tested for any association with treatment response, using univariate and multivariate logistic regression models. RESULTS: Of the 31 RA-associated risk alleles, a SNP at the PTPRC (also known as CD45) gene locus (rs10919563) was associated with the primary end point, a EULAR good response versus no response (odds ratio [OR] 0.55, P = 0.0001 in the multivariate model). Similar results were obtained using the secondary end point, the triangle upDAS28 (P = 0.0002). There was suggestive evidence of a stronger association in autoantibody-positive patients with RA (OR 0.55, 95% confidence interval [95% CI] 0.39-0.76) as compared with autoantibody-negative patients (OR 0.90, 95% CI 0.41-1.99). CONCLUSION: Statistically significant associations were observed between the response to anti-TNF therapy and an RA risk allele at the PTPRC gene locus. Additional studies will be required to replicate this finding in additional patient collections.

Genome-wide association scan identifies candidate polymorphisms associated with differential response to anti-TNF treatment in rheumatoid arthritis.

The prediction of response (or non-response) to anti-TNF treatment for rheumatoid arthritis (RA) is a pressing clinical problem. We conducted a genome-wide association study using the Illumina HapMap300 SNP chip on 89 RA patients prospectively followed after beginning anti-TNF therapy as part of Autoimmune Biomarkers Collaborative Network (ABCoN [Autoimmune Bio-markers Collaborative Network]) patient cohort. Response to therapy was determined by the change in Disease Activity Score (DAS28) observed after 14 wks. We used a two-part analysis that treated the change in DAS28 as a continuous trait and then incorporated it into a dichotomous trait of "good responder" and "nonresponder" by European League Against Rheumatism (EULAR) criteria. We corrected for multiple tests by permutation, and adjusted for potential population stratification using EIGENSTRAT. Multiple single nucleotide polymorphism (SNP) markers showed significant associations near or within loci including: the v-maf musculoaponeurotic fibrosarcoma oncogene homolog B (MAFB) gene on chromosome 20; the type I interferon gene IFNk on chromosome 9; and in a locus on chromosome 7 that includes the paraoxonase I (PON1) gene. An SNP in the IL10 promoter (rs1800896) that was previously reported as associated with anti-TNF response was weakly associated with response in this cohort. Replications of these results in independent and larger data sets clearly are required. We provide a reference list of candidate SNPs (P < 0.01) that can be investigated in future pharmacogenomic studies.

The Genotype and Phenotype (GaP) registry: a living biobank for the analysis of quantitative traits.

We describe the development of the Genotype and Phenotype (GaP) Registry, a living biobank of normal volunteers who are genotyped for genetic markers related to human disease. Participants in the GaP can be recalled for hypothesis driven study of disease associated genetic variants. The GaP has facilitated functional studies of several autoimmune disease associated loci including Csk, Blk, PDRM1 (Blimp-1) and PTPN22. It is likely that expansion of such living biobank registries will play an important role in studying and understanding the function of disease associated alleles in complex disease.

Modular analysis of peripheral blood gene expression in rheumatoid arthritis captures reproducible gene expression changes in tumor necrosis factor responders.

OBJECTIVE: To establish whether the analysis of whole-blood gene expression is useful in predicting or monitoring response to anti-tumor necrosis factor (anti-TNF) therapy in patients with rheumatoid arthritis (RA). METHODS: Whole-blood RNA (using a PAXgene system to stabilize whole-blood RNA in the collection tube) was obtained at baseline and at 14 weeks from 3 independent cohorts, consisting of a combined total of 240 RA patients who were beginning therapy with anti-TNF. We used an approach to gene expression analysis that is based on modular patterns of gene expression, or modules. RESULTS: Good and moderate responders according to the European League Against Rheumatism criteria exhibited highly significant and consistent changes in multiple gene expression modules after 14 weeks of therapy, as demonstrated by hypergeometric analysis. Strikingly, nonresponders exhibited very little change in any modules, despite exposure to TNF blockade. These patterns of change were highly consistent across all 3 cohorts, indicating that immunologic changes after TNF treatment are specific to the combination of both drug exposure and responder status. In contrast, modular patterns of gene expression did not exhibit consistent differences between responders and nonresponders at baseline in the 3 study cohorts. CONCLUSION: These data provide evidence that using gene expression modules related to inflammatory disease may provide a valuable method for objective monitoring of the response of RA patients who are treated with TNF inhibitors.

Association of a single-nucleotide polymorphism in CD40 with the rate of joint destruction in rheumatoid arthritis.

OBJECTIVE: The severity of joint destruction in rheumatoid arthritis (RA) is highly variable from patient to patient and is influenced by genetic factors. Genome-wide association studies have enormously boosted the field of the genetics of RA susceptibility, but risk loci for RA severity remain poorly defined. A recent meta-analysis of genome-wide association studies identified 6 genetic regions for susceptibility to autoantibody-positive RA: CD40, KIF5A/PIP4K2C, CDK6, CCL21, PRKCQ, and MMEL1/TNFRSF14. The purpose of this study was to investigate whether these newly described genetic regions are associated with the rate of joint destruction. METHODS: RA patients enrolled in the Leiden Early Arthritis Clinic were studied (n=563). Yearly radiographs were scored using the Sharp/van der Heijde method (median followup 5 years; maximum followup 9 years). The rate of joint destruction between genotype groups was compared using a linear mixed model, correcting for age, sex, and treatment strategies. A total of 393 anti-citrullinated protein antibody (ACPA)-positive RA patients from the North American Rheumatoid Arthritis Consortium (NARAC) who had radiographic data available were used for the replication study. RESULTS: The TT and CC/CG genotypes of 2 single-nucleotide polymorphisms, rs4810485 (CD40) and rs42041 (CDK6), respectively, were associated with a higher rate of joint destruction in ACPA-positive RA patients (P=0.003 and P=0.012, respectively), with rs4810485 being significant after Bonferroni correction for multiple testing. The association of the CD40 minor allele with the rate of radiographic progression was replicated in the NARAC cohort (P=0.021). CONCLUSION: A polymorphism in the CD40 locus is associated with the rate of joint destruction in patients with ACPA-positive RA. Our findings provide one of the first non-HLA-related genetic severity factors that has been replicated.

Interaction between smoking, the shared epitope, and anti-cyclic citrullinated peptide: a mixed picture in three large North American rheumatoid arthritis cohorts.

OBJECTIVE: Recently, Swedish members of the Epidemiological Investigation of Rheumatoid Arthritis (EIRA) provided evidence that smoking may trigger RA-specific immune reactions to citrullinated protein in carriers of HLA-DR shared epitope alleles. In an effort to confirm this interaction between smoking and shared epitope alleles, we performed a case-only analysis of 3 North American RA cohorts. METHODS: A total of 2,476 white patients with RA were studied, 1,105 from the North American Rheumatoid Arthritis Consortium (NARAC) family collection, 753 from the National Inception Cohort of Rheumatoid Arthritis Patients (Inception Cohort), and 618 from the Study of New Onset Rheumatoid Arthritis (SONORA). All patients were HLA-DRB1 typed, and tested for anti-cyclic citrullinated peptide (anti-CCP) and rheumatoid factor. Information about smoking history was obtained by questionnaire. RESULTS: A significant association was found between smoking and the presence of anti-CCP in the NARAC and the Inception Cohort, but not in the SONORA. The shared epitope alleles consistently correlated with anti-CCP in all 3 populations. Using multiple logistic regression analyses, shared epitope alleles were still the most significant risk factor for anti-CCP positivity. Weak evidence of gene-environment interaction between smoking and shared epitope alleles for anti-CCP formation was found only in the NARAC. CONCLUSION: Unlike the EIRA data, we could not confirm a major gene-environment interaction for anti-CCP formation between shared epitope alleles and smoking in 3 North American RA cohorts. Our data indicate a need for further studies to address the full range of environmental factors other than smoking that may be associated with citrullination and RA.

Microarray analyses of peripheral blood cells identifies unique gene expression signature in psoriatic arthritis.

Psoriatic arthritis (PsA) is a chronic and erosive form of arthritis of unknown cause. We aimed to characterize the PsA phenotype using gene expression profiling and comparing it with healthy control subjects and patients rheumatoid arthritis (RA). Peripheral blood cells (PBCs) of 19 patients with active PsA and 19 age- and sex-matched control subjects were used in the analyses of PsA, with blood samples collected in PaxGene tubes. A significant alteration in the pattern of expression of 313 genes was noted in the PBCs of PsA patients on Affymetrix U133A arrays: 257 genes were expressed at reduced levels in PsA, and 56 genes were expressed at increased levels, compared with controls. Downregulated genes tended to cluster to certain chromosomal regions, including those containing the psoriasis susceptibility loci PSORS1 and PSORS2. Among the genes with the most significantly reduced expression were those involved in downregulation or suppression of innate and acquired immune responses, such as SIGIRR, STAT3, SHP1, IKBKB, IL-11RA, and TCF7, suggesting inappropriate control that favors proin-flammatory responses. Several members of the MAPK signaling pathway and tumor suppressor genes showed reduced expression. Three proinflammatory genes--S100A8, S100A12, and thioredoxin--showed increased expression. Logistic regression and recursive partitioning analysis determined that one gene, nucleoporin 62 kDa, could correctly classify all controls and 94.7% of the PsA patients. Using a dataset of 48 RA samples for comparison, the combination of two genes, MAP3K3 followed by CACNA1S, was enough to correctly classify all RA and PsA patients. Thus, PBC gene expression profiling identified a gene expression signature that differentiated PsA from RA, and PsA from controls. Several novel genes were differentially expressed in PsA and may prove to be diagnostic biomarkers or serve as new targets for the development of therapies.

Regulation of anti-cyclic citrullinated peptide antibodies in rheumatoid arthritis: contrasting effects of HLA-DR3 and the shared epitope alleles.

OBJECTIVE: To examine the association between HLA-DRB1 alleles and the production of anti-cyclic citrullinated peptide (anti-CCP) and rheumatoid factor (RF) autoantibodies in patients with rheumatoid arthritis (RA). METHODS: We studied 1,723 Caucasian RA patients enrolled in the North American Rheumatoid Arthritis Consortium (NARAC) family cohort and the Study of New Onset Rheumatoid Arthritis (SONORA) cohort. All patients were tested for anti-CCP antibodies (by enzyme-linked immunosorbent assay), RF (by nephelometry), and HLA-DR genotype (by polymerase chain reaction and sequence-specific oligonucleotide hybridization). RESULTS: When controlled for the presence of RF, anti-CCP positivity was strongly associated with the HLA-DRB1 shared epitope (SE). In RF+ patients, the presence of the SE was very significantly associated with anti-CCP positivity, with an odds ratio (OR) of 5.8 and a 95% confidence interval (95% CI) of 4.1-8.3. This relationship was also seen in RF- patients (OR 3.1 [95% CI 1.8-5.3]). In contrast, RF positivity was not significantly associated with presence of the SE independently of anti-CCP antibodies. Strikingly, HLA-DRB1*03 was strongly associated with reduced anti-CCP titers, even after controlling for the presence of the SE and restricting the analysis to anti-CCP+ patients. HLA-DR3 was also associated with anti-CCP- RA in our population. CONCLUSION: The HLA-DRB1 SE is strongly associated with the production of anti-CCP antibodies, but not RF. In contrast, HLA-DR3 alleles are associated with anti-CCP- disease and with lower levels of anti-CCP antibodies, even when controlling for the SE. These data emphasize the complexity of the genetic effects of the major histocompatibility complex on the RA phenotype.

Analysis of families in the multiple autoimmune disease genetics consortium (MADGC) collection: the PTPN22 620W allele associates with multiple autoimmune phenotypes.

Autoimmune disorders constitute a diverse group of phenotypes with overlapping features and a tendency toward familial aggregation. It is likely that common underlying genes are involved in these disorders. Until very recently, no specific alleles--aside from a few common human leukocyte antigen class II genes--had been identified that clearly associate with multiple different autoimmune diseases. In this study, we describe a unique collection of 265 multiplex families assembled by the Multiple Autoimmune Disease Genetics Consortium (MADGC). At least two of nine "core" autoimmune diseases are present in each of these families. These core diseases include rheumatoid arthritis (RA), systemic lupus erythematosus (SLE), type 1 diabetes (T1D), multiple sclerosis (MS), autoimmune thyroid disease (Hashimoto thyroiditis or Graves disease), juvenile RA, inflammatory bowel disease (Crohn disease or ulcerative colitis), psoriasis, and primary Sjogren syndrome. We report that a recently described functional single-nucleotide polymorphism (rs2476601, encoding R620W) in the intracellular tyrosine phosphatase (PTPN22) confers risk of four separate autoimmune phenotypes in these families: T1D, RA, SLE, and Hashimoto thyroiditis. MS did not show association with the PTPN22 risk allele. These findings suggest a common underlying etiologic pathway for some, but not all, autoimmune disorders, and they suggest that MS may have a pathogenesis that is distinct from RA, SLE, and T1D. DNA and clinical data for the MADGC families are available to the scientific community; these data will provide a valuable resource for the dissection of the complex genetic factors that underlie the various autoimmune phenotypes.